Intronic hammerhead ribozymes are ultraconserved in the human genome

Minimal twister sister-like self-cleaving ribozymes in the human genome revealed by deep mutational scanning

Despite their importance in a wide range of living organisms, self-cleaving ribozymes in the human genome are few and poorly studied. Here, we performed deep mutational scanning and covariance analysis of two previously proposed self-cleaving ribozymes (LINE-1 and OR4K15). We found that the regions essential for ribozyme activities are made of two short segments, with a total of 35 and 31 nucleotides only. The discovery makes them the simplest known self-cleaving ribozymes. Moreover, the essential regions are circular permutated with two nearly identical catalytic internal loops, supported by two stems of different lengths. These two self-cleaving ribozymes, which are shaped like lanterns, are similar to the catalytic regions of the twister sister ribozymes in terms of sequence and secondary structure. However, the nucleotides at the cleavage site have shown that mutational effects on two twister sister-like (TS-like) ribozymes are different from the twister sister ribozyme. The discovery of TS-like ribozymes reveals a ribozyme class with the simplest and, perhaps, the most primitive structure needed for self-cleavage.

Evaluating Performance of Different RNA Secondary Structure Prediction Programs Using Self-cleaving Ribozymes

Accurate identification of the correct, biologically relevant RNA structures is critical to understanding various aspects of RNA biology since proper folding represents the key to the functionality of all types of RNA molecules and plays pivotal roles in many essential biological processes. Thus, a plethora of approaches have been developed to predict, identify, or solve RNA structures based on various computational, molecular, genetic, chemical, or physicochemical strategies. Purely computational approaches hold distinct advantages over all other strategies in terms of the ease of implementation, time, speed, cost, and throughput, but they strongly underperform in terms of accuracy that significantly limits their broader application. Nonetheless, the advantages of these methods led to a steady development of multiple in silico RNA secondary structure prediction approaches including recent deep learning-based programs. Here, we compared the accuracy of predictions of biologically relevant secondary structures of dozens of self-cleaving ribozyme sequences using seven in silico RNA folding prediction tools with tasks of varying complexity. We found that while many programs performed well in relatively simple tasks, their performance varied significantly in more complex RNA folding problems. However, in general, a modern deep learning method outperformed the other programs in the complex tasks in predicting the RNA secondary structures, at least based on the specific class of sequences tested, suggesting that it may represent the future of RNA structure prediction algorithms.

The structure and catalytic mechanism of a pseudoknot-containing hammerhead ribozyme

We have determined the crystal structure of a pseudoknot (PK)-containing hammerhead ribozyme that closely resembles the pistol ribozyme, with essentially identical secondary structure and connectivity. The activity is more sensitive to deletion of the G8 2'OH than to the absence of magnesium ions, indicating that the catalytic mechanism is the same as the extended hammerhead, and distinct from the pistol ribozyme. Here we show that nucleophilic attack is almost perfectly in-line, and the G8 2'OH is well positioned to act as general acid, being directed towards the O5' leaving group, and 2.9 Å away from it. Despite the similarity in overall structure to the pistol ribozyme, the local structure close to the cleavage site differs, and the PK hammerhead retains its unique mechanistic identity and demonstrates enhanced activity over other hammerhead ribozymes under standard conditions.

Inhibition of Cpeb3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element-binding protein 3 (Cpeb3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the Cpeb3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Conserved untranslated regions of multipartite viruses: Natural markers of novel viral genomic components and tags of viral evolution

Viruses with split genomes are classified as being either segmented or multipartite based on whether their genomic segments occur within a single virion or across different virions. Despite variations in number and sequence during evolution, the genomic segments of many viruses are conserved within the untranslated regions (UTRs). In this study, we present a methodology that combines RNA sequencing with iterative BLASTn of UTRs (https://github.com/qq371260/Iterative-blast-v.1.0) to identify new viral genomic segments. Some novel multipartite-like viruses related to the phylum Kitrinoviricota were annotated using sequencing data from field plant samples and public databases. We identified potentially plant-infecting jingmen-related viruses (order Amarillovirales) and jivi-related viruses (order Martellivirales) with at least six genomic components. The number of RNA molecules associated with a genome of the novel viruses in the families Closteroviridae, Kitaviridae, and Virgaviridae within the order Martellivirales reached five. Several of these viruses seem to represent new taxa at the subgenus, genus, and family levels. The diversity of novel genomic components and the multiple duplication of proteins or protein domains within single or multiple genomic components emphasize the evolutionary roles of genetic recombination (horizontal gene transfer), reassortment, and deletion. The relatively conserved UTRs at the genome level might explain the relationships between monopartite and multipartite viruses, as well as how subviral agents such as defective RNAs and satellite viruses can coexist with their helper viruses.

A comparative survey of the influence of small self-cleaving ribozymes on gene expression in human cell culture

Self-cleaving ribozymes are versatile tools for synthetic biologists when it comes to controlling gene expression. Up to date, 12 different classes are known, and over the past decades more and more details about their structure, cleavage mechanisms and natural environments have been uncovered. However, when these motifs are applied to mammalian gene expression constructs, the outcome can often be unexpected. A variety of factors, such as surrounding sequences and positioning of the ribozyme influences the activity and hence performance of catalytic RNAs. While some information about the efficiency of individual ribozymes (each tested in specific contexts) is known, general trends obtained from standardized, comparable experiments are lacking, complicating decisions such as which ribozyme to choose and where to insert it into the target mRNA. In many cases, application-specific optimization is required, which can be very laborious. Here, we systematically compared different classes of ribozymes within the 3'-UTR of a given reporter gene. We then examined position-dependent effects of the best-performing ribozymes. Moreover, we tested additional variants of already widely used hammerhead ribozymes originating from various organisms. We were able to identify functional structures suited for aptazyme design and generated highly efficient hammerhead ribozyme variants originating from the human genome. The present dataset will aide decisions about how to apply ribozymes for affecting gene expression as well as for developing ribozyme-based switches for controlling gene expression in human cells.

Inhibition of CPEB3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs, and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element binding protein 3 (CPEB3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the CPEB3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Viroids: Non-Coding Circular RNAs Able to Autonomously Replicate and Infect Higher Plants

Viroids are a unique type of infectious agent, exclusively composed of a relatively small (246-430 nt), highly base-paired, circular, non-coding RNA. Despite the small size and non-coding nature, the more-than-thirty currently known viroid species infectious of higher plants are able to autonomously replicate and move systemically through the host, thereby inducing disease in some plants. After recalling viroid discovery back in the late 60s and early 70s of last century and discussing current hypotheses about their evolutionary origin, this article reviews our current knowledge about these peculiar infectious agents. We describe the highly base-paired viroid molecules that fold in rod-like or branched structures and viroid taxonomic classification in two families, Pospiviroidae and Avsunviroidae, likely gathering nuclear and chloroplastic viroids, respectively. We review current knowledge about viroid replication through RNA-to-RNA rolling-circle mechanisms in which host factors, notably RNA transporters, RNA polymerases, RNases, and RNA ligases, are involved. Systemic movement through the infected plant, plant-to-plant transmission and host range are also discussed. Finally, we focus on the mechanisms of viroid pathogenesis, in which RNA silencing has acquired remarkable importance, and also for the initiation of potential biotechnological applications of viroid molecules.

A life of research on circular RNAs and ribozymes: towards the origin of viroids, deltaviruses and life

The first examples of circular RNAs (circRNAs) were reported in the '70s as a family of minimal infectious agents of flowering plants; the viroids and viral satellites of circRNA. In some cases, these small circular genomes encode self-cleaving RNA motifs or ribozymes, including an exceptional circRNA infecting not plants but humans: the Hepatitis Delta Virus. Autocatalytic ribozymes not only allowed to propose a common rolling-circle replication mechanism for all these subviral agents, but also a tentative link with the origin of life as molecular fossils of the so-called RNA world. Despite the weak biologic connection between angiosperm plants and the human liver, diverse scientists, and most notably Ricardo Flores, firmly supported an evolutionary relationship between plant viroids and human deltavirus agents. The tireless and inspiring work done by Ricardo's lab in the field of infectious circRNAs fuelled multiple hypotheses for the origin of these entities, allowing advances in other fields, from eukaryotic circRNAs to small ribozymes in genomes from all life kingdoms. The recent discovery of a plethora of viral-like circRNAs with ribozymes in disparate biological samples may finally allow us to connect plant and animal subviral agents, confirming again that Ricardo's eye for science was always a keen eye.

The potential versatility of RNA catalysis

It is commonly thought that in the early development of life on this planet RNA would have acted both as a store of genetic information and as a catalyst. While a number of RNA enzymes are known in contemporary cells, they are largely confined to phosphoryl transfer reactions, whereas an RNA based metabolism would have required a much greater chemical diversity of catalysis. Here we discuss how RNA might catalyze a wider variety of chemistries, and particularly how information gleaned from riboswitches could suggest how ribozymes might recruit coenzymes to expand their chemical range. We ask how we might seek such activities in modern biology. This article is categorized under: RNA-Based Catalysis > Miscellaneous RNA-Catalyzed Reactions Regulatory RNAs/RNAi/Riboswitches > Riboswitches RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Kinetic Study of the Avocado Sunblotch Viroid Self-Cleavage Reaction Reveals Compensatory Effects between High-Pressure and High-Temperature: Implications for Origins of Life on Earth

A high pressure apparatus allowing one to study enzyme kinetics under pressure was used to study the self-cleavage activity of the avocado sunblotch viroid. The kinetics of this reaction were determined under pressure over a range up to 300 MPa (1-3000 bar). It appears that the initial rate of this reaction decreases when pressure increases, revealing a positive ΔV≠ of activation, which correlates with the domain closure accompanying the reaction and the decrease of the surface of the viroid exposed to the solvent. Although, as expected, temperature increases the rate of the reaction whose energy of activation was determined, it appeared that it does not significantly influence the ΔV≠ of activation and that pressure does not influence the energy of activation. These results provide information about the structural aspects or this self-cleavage reaction, which is involved in the process of maturation of this viroid. The behavior of ASBVd results from the involvement of the hammerhead ribozyme present at its catalytic domain, indeed a structural motif is very widespread in the ancient and current RNA world.

Enhancement of the Catalytic Activity of Hammerhead Ribozymes by Organic Cations

Catalytic turnover is important for the application of ribozymes to biotechnology. However, the turnover is often impaired because of the intrinsic high stability of base pairs with cleaved RNA products. Here, organic cations were used as additives to improve the catalytic performance of hammerhead ribozyme constructs that exhibit different kinetic behaviors. Kinetic analysis of substrate cleavage demonstrated that bulky cations, specifically tetra-substituted ammonium ions containing pentyl groups or a benzyl group, have the ability to greatly increase the turnover rate of the ribozymes. Thermal stability analysis of RNA structures revealed that the bulky cations promote the dissociation of cleaved products and refolding of incorrectly folded structures with small disruption of the catalytic structure. The use of bulky cations is a convenient method for enhancing the catalytic activity of hammerhead ribozymes, and the approach may be useful for advancing ribozyme technologies.

Hovlinc is a recently evolved class of ribozyme found in human lncRNA

Although naturally occurring catalytic RNA molecules-ribozymes-have attracted a great deal of research interest, very few have been identified in humans. Here, we developed a genome-wide approach to discovering self-cleaving ribozymes and identified a naturally occurring ribozyme in humans. The secondary structure and biochemical properties of this ribozyme indicate that it belongs to an unidentified class of small, self-cleaving ribozymes. The sequence of the ribozyme exhibits a clear evolutionary path, from its appearance between ~130 and ~65 million years ago (Ma), to acquiring self-cleavage activity very recently, ~13-10 Ma, in the common ancestors of humans, chimpanzees and gorillas. The ribozyme appears to be functional in vivo and is embedded within a long noncoding RNA belonging to a class of very long intergenic noncoding RNAs. The presence of a catalytic RNA enzyme in lncRNA creates the possibility that these transcripts could function by carrying catalytic RNA domains.

A cell-based ribozyme reporter system employing a chromosomally-integrated 5' exonuclease gene

Background

Bioinformatic genome surveys indicate that self-cleaving ribonucleic acids (ribozymes) appear to be widespread among all domains of life, although the functions of only a small number have been validated by biochemical methods. Alternatively, cell-based reporter gene assays can be used to validate ribozyme function. However, reporter activity can be confounded by phenomena unrelated to ribozyme-mediated cleavage of RNA.

Results

We established a ribozyme reporter system in Escherichia coli in which a significant reduction of reporter activity is manifest when an active ribozyme sequence is fused to the reporter gene and the expression of a foreign Bacillus subtilis RNaseJ1 5′ exonuclease is induced from a chromosomally-integrated gene in the same cell.

Conclusions

The reporter system could be useful for validating ribozyme function in candidate sequences identified from bioinformatics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00357-7.

Cloning and Detection of Genomic Retrozymes and Their circRNA Intermediates

Retrozymes are a novel family of non-autonomous retrotransposable elements that contain hammerhead ribozyme motifs. These retroelements are found widespread in eukaryotic genomes, with active copies present in many species, which rely on other autonomous transposons for mobilization. Contrary to other retrotransposons, transcription of retrozymes in vivo leads to the formation and accumulation of circular RNAs, which can be readily detected by RNA blotting. In this chapter, we describe the procedures needed to carry out the cloning of genomic retrozymes, and to detect by northern blot their circular RNA retrotransposition intermediates.

Co-transcriptional Analysis of Self-Cleaving Ribozymes and Their Ligand Dependence

Self-cleaving ribozymes are RNA molecules that catalyze a site-specific self-scission reaction. Analysis of self-cleavage is a crucial aspect of the biochemical study and understanding of these molecules. Here we describe a co-transcriptional assay that allows the analysis of self-cleaving ribozymes in different reaction conditions and in the presence of desired ligands and/or cofactors. Utilizing a standard T7 RNA polymerase in vitro transcription system under limiting Mg²⁺ concentration, followed by a 25-fold dilution of the reaction in desired conditions of self-cleavage (buffer, ions, ligands, pH, temperature, etc.) to halt the synthesis of new RNA molecules, allows the study of self-scission of these molecules without the need for purification or additional preparation steps, such as refolding procedures. Furthermore, because the transcripts are not denatured, this assay likely yields RNAs in conformations relevant to co-transcriptionally folded species in vivo.

A Singular and Widespread Group of Mobile Genetic Elements: RNA Circles with Autocatalytic Ribozymes

Circular DNAs, such as most prokaryotic and phage genomes, are a frequent form of nucleic acids, whereas circular RNAs had been regarded as unusual macromolecules until very recently. The first reported RNA circles were the family of small infectious genomes of viroids and circular RNA (circRNA) satellites of plant viruses, some of which contain small self-cleaving RNA motifs, such as the hammerhead (HHR) and hairpin ribozymes. A similar infectious circRNA, the unique human hepatitis delta virus (HDV), is another viral satellite that also encodes self-cleaving motifs called HDV ribozymes. Very recently, different animals have been reported to contain HDV-like circRNAs with typical HDV ribozymes, but also conserved HHR motifs, as we describe here. On the other hand, eukaryotic and prokaryotic genomes encode sequences able to self-excise as circRNAs, like the autocatalytic Group I and II introns, which are widespread genomic mobile elements. In the 1990s, the first circRNAs encoded in a mammalian genome were anecdotally reported, but their abundance and importance have not been unveiled until recently. These gene-encoded circRNAs are produced by events of alternative splicing in a process generally known as backsplicing. However, we have found a second natural pathway of circRNA expression conserved in numerous plant and animal genomes, which efficiently promotes the accumulation of small non-coding RNA circles through the participation of HHRs. Most of these genome-encoded circRNAs with HHRs are the transposition intermediates of a novel family of non-autonomous retrotransposons called retrozymes, with intriguing potential as new forms of gene regulation.

Small circRNAs with self-cleaving ribozymes are highly expressed in diverse metazoan transcriptomes

Ribozymes are catalytic RNAs present in modern genomes but regarded as remnants of a prebiotic RNA world. The paradigmatic hammerhead ribozyme (HHR) is a small self-cleaving motif widespread from bacterial to human genomes. Here, we report that most of the classical type I HHRs frequently found in the genomes of animals are contained within a novel family of non-autonomous non-LTR retrotransposons of the retrozyme class. These retroelements are expressed as abundant linear and circular RNAs of ∼170-400 nt in different animal tissues. Bioinformatic and in vitro analyses indicate an efficient self-cleavage of the HHRs harboured in most invertebrate retrozymes, whereas HHRs in retrozymes of vertebrates, such as the axolotl and other amphibians, require to act as dimeric motifs to reach higher self-cleavage rates. Ligation assays of retrozyme RNAs with a protein ligase versus HHR self-ligation indicate that, most likely, tRNA ligases and not the ribozymes are involved in the step of RNA circularization. Altogether, these results confirm the existence of a new and conserved pathway in animals and, likely, eukaryotes in general, for the efficient biosynthesis of RNA circles through small ribozymes, which opens the door for the development of new tools in the emerging field of study of circRNAs.

Accurate inference of the full base-pairing structure of RNA by deep mutational scanning and covariation-induced deviation of activity

Despite the large number of noncoding RNAs in human genome and their roles in many diseases include cancer, we know very little about them due to lack of structural clues. The centerpiece of the structural clues is the full RNA base-pairing structure of secondary and tertiary contacts that can be precisely obtained only from costly and time-consuming 3D structure determination. Here, we performed deep mutational scanning of self-cleaving CPEB3 ribozyme by error-prone PCR and showed that a library of <5 × 104 single-to-triple mutants is sufficient to infer 25 of 26 base pairs including non-nested, nonhelical, and noncanonical base pairs with both sensitivity and precision at 96%. Such accurate inference was further confirmed by a twister ribozyme at 100% precision with only noncanonical base pairs as false negatives. The performance was resulted from analyzing covariation-induced deviation of activity by utilizing both functional and nonfunctional variants for unsupervised classification, followed by Monte Carlo (MC) simulated annealing with mutation-derived scores. Highly accurate inference can also be obtained by combining MC with evolution/direct coupling analysis, R-scape or epistasis analysis. The results highlight the usefulness of deep mutational scanning for high-accuracy structural inference of self-cleaving ribozymes with implications for other structured RNAs that permit high-throughput functional selections.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

RNA Back and Forth: Looking through Ribozyme and Viroid Motifs

Current cellular facts allow us to follow the link from chemical to biochemical metabolites, from the ancient to the modern world. In this context, the "RNA world" hypothesis proposes that early in the evolution of life, the ribozyme was responsible for the storage and transfer of genetic information and for the catalysis of biochemical reactions. Accordingly, the hammerhead ribozyme (HHR) and the hairpin ribozyme belong to a family of endonucleolytic RNAs performing self-cleavage that might occur during replication. Furthermore, regarding the widespread occurrence of HHRs in several genomes of modern organisms (from mammals to small parasites and elsewhere), these small ribozymes have been regarded as living fossils of a primitive RNA world. They fold into 3D structures that generally require long-range intramolecular interactions to adopt the catalytically active conformation under specific physicochemical conditions. By studying viroids as plausible remains of ancient RNA, we recently demonstrated that they replicate in non-specific hosts, emphasizing their adaptability to different environments, which enhanced their survival probability over the ages. All these results exemplify ubiquitous features of life. Those are the structural and functional versatility of small RNAs, ribozymes, and viroids, as well as their diversity and adaptability to various extreme conditions. All these traits must have originated in early life to generate novel RNA populations.

Bulky cations greatly increase the turnover of a native hammerhead ribozyme

Methods to facilitate the catalytic turnover of ribozymes are required for advancing oligonucleotide-based technologies. This study examined tetraalkylammonium ions for their ability to increase the efficiency of catalytic turnover of a native hammerhead ribozyme. Kinetic analysis showed that large tetraalkylammonium ions significantly increased the turnover rate of the ribozyme and was much more effective than poly(ethylene glycol) (PEG) and urea. The magnitude of the rate increase depended on the concentrations of Mg²⁺ and tetrapentylammonium ions, and the rate was enhanced by more than 180-fold at the optimal concentrations of these salts. The results provide physical insights into interactions of ribozymes with large cationic molecules through electrostatic forces and steric hindrance.

Large tetraalkylammonium ions increase the turnover rate of the ribozyme derived from an intronic ribozyme in the human genome. The rate can be enhanced by more than a hundred-fold at the optimal concentrations of Mg²⁺ and TPeA ions.

Kinetic Parameters of trans Scission by Extended HDV-like Ribozymes and the Prospect for the Discovery of Genomic trans-Cleaving RNAs

Hepatitis delta virus (HDV)-like ribozymes are self-cleaving catalytic RNAs with a widespread distribution in nature and biological roles ranging from self-scission during rolling-circle replication in viroids to co-transcriptional processing of eukaryotic retrotransposons, among others. The ribozymes fold into a double pseudoknot with a common catalytic core motif and highly variable peripheral domains. Like other self-cleaving ribozymes, HDV-like ribozymes can be converted into trans-acting catalytic RNAs by bisecting the self-cleaving variants at non-essential loops. Here we explore the trans-cleaving activity of ribozymes derived from the largest examples of the ribozymes (drz-Agam-2 family), which contain an extended domain between the substrate strand and the rest of the RNA. When this peripheral domain is bisected at its distal end, the substrate strand is recognized through two helices, rather than just one 7 bp helix common among the HDV ribozymes, resulting in stronger binding and increased sequence specificity. Kinetic characterization of the extended trans-cleaving ribozyme revealed an efficient trans-cleaving system with a surprisingly high K_M', supporting a model that includes a recently proposed activation barrier related to the assembly of the catalytically competent ribozyme. The ribozymes also exhibit a very long k_off for the products (∼2 weeks), resulting in a trade-off between sequence specificity and turnover. Finally, structure-based searches for the catalytic cores of these ribozymes in the genome of the mosquito Anopheles gambiae, combined with sequence searches for their putative substrates, revealed two potential ribozyme-substrate pairs that may represent examples of natural trans-cleaving ribozymes.

Circular RNAs Biogenesis in Eukaryotes Through Self-Cleaving Hammerhead Ribozymes

Circular DNAs are frequent genomic molecules, especially among the simplest life beings, whereas circular RNAs have been regarded as weird nucleic acids in biology. Now we know that eukaryotes are able to express circRNAs, mostly derived from backsplicing mechanisms, and playing different biological roles such as regulation of RNA splicing and transcription, among others. However, a second natural and highly efficient pathway for the expression in vivo of circRNAs has been recently reported, which allows the accumulation of abundant small (100-1000 nt) non-coding RNA circles through the participation of small self-cleaving RNAs or ribozymes called hammerhead ribozymes. These genome-encoded circRNAs with ribozymes seem to be a new family of small and nonautonomous retrotransposable elements of plants and animals (so-called retrozymes), which will offer functional clues to the biology and evolution of circular RNA molecules as well as new biotechnological tools in this emerging field.

The Importance of ncRNAs as Epigenetic Mechanisms in Phenotypic Variation and Organic Evolution

Neo-Darwinian explanations of organic evolution have settled on mutation as the principal factor in producing evolutionary novelty. Mechanistic characterizations have been also biased by the classic dogma of molecular biology, where only proteins regulate gene expression. This together with the rearrangement of genetic information, in terms of genes and chromosomes, was considered the cornerstone of evolution at the level of natural populations. This predominant view excluded both alternative explanations and phenomenologies that did not fit its paradigm. With the discovery of non-coding RNAs (ncRNAs) and their role in the control of genetic expression, new mechanisms arose providing heuristic power to complementary explanations to evolutionary processes overwhelmed by mainstream genocentric views. Viruses, epimutation, paramutation, splicing, and RNA editing have been revealed as paramount functions in genetic variations, phenotypic plasticity, and diversity. This article discusses how current epigenetic advances on ncRNAs have changed the vision of the mechanisms that generate variation, how organism-environment interaction can no longer be underestimated as a driver of organic evolution, and how it is now part of the transgenerational inheritance and evolution of species.

Highly motif- and organism-dependent effects of naturally occurring hammerhead ribozyme sequences on gene expression

Recent bioinformatics studies have demonstrated a wide-spread occurrence of the hammerhead ribozyme (HHR) and similar small endonucleolytic RNA motifs in all domains of life. It is becoming increasingly evident that such ribozyme motifs participate in important genetic processes in diverse organisms. Although the HHR motif has been studied for more than three decades, only little is known about the consequences of ribozyme activity on gene expression. In the present study we analysed eight different naturally occurring HHR sequences in diverse genetic and organismal contexts. We investigated the influence of active ribozymes incorporated into mRNAs in mammalian, yeast and bacterial expression systems. The experiments show an unexpectedly high degree of organism-specific variability of ribozyme-mediated effects on gene expression. The presented findings demonstrate that ribozyme cleavage profoundly affect gene expression. However, the extent of this effect varies and depends strongly on the respective genetic context. The fast-cleaving type 3 HHRs [CChMVd(-) and sLTSV(-)] generally tended to cause the strongest effects on intracellular gene expression. The presented results are important in order to address potential functions of naturally occurring ribozymes in RNA processing and post-transcriptional regulation of gene expression. Additionally, our results are of interest for biotechnology and synthetic biology approaches that aim at the utilisation of self-cleaving ribozymes as widely applicable tools for controlling genetic processes.

Circular RNAs with hammerhead ribozymes encoded in eukaryotic genomes: The enemy at home

A new family of non-autonomous retrotransposons with self-cleaving hammerhead ribozymes, the so called retrozymes, has recently been found encoded in diverse plant genomes. These retroelements can be actively transcribed, and their RNAs accumulate in the cells as abundant non-coding circular RNAs (circRNAs) of small size (600–1000 nt). Related circRNAs with self-cleaving ribozymes had already been described in plants, and belong to a group of infectious RNA agents with an uncertain origin: the viroids and viroid-like satellites of plant RNA viruses. These pathogenic circRNAs show many structural similarities with retrozyme circRNAs, and both have been found to occur in flowering plants as heterogeneous RNA molecules of positive and negative polarities. Taking all these data together, we hypothesize that circRNAs encoded by genomic retrozymes could have given origin to infectious circRNAs with self-cleaving ribozymes. Moreover, we propose that retrozymes in time could have evolved from the ancient family of Penelope-like retroelements, which also harbour hammerhead ribozymes. Putative retrozyme sequences with hammerhead ribozymes have been detected as well in metazoan genomes, opening the door to a common occurrence of circRNAs with self-cleaving motifs among eukaryotes.

The Hammerhead Ribozyme: A Long History for a Short RNA

Small nucleolytic ribozymes are a family of naturally occurring RNA motifs that catalyse a self-transesterification reaction in a highly sequence-specific manner. The hammerhead ribozyme was the first reported and the most extensively studied member of this family. However, and despite intense biochemical and structural research for three decades since its discovery, the history of this model ribozyme seems to be far from finished. The hammerhead ribozyme has been regarded as a biological oddity typical of small circular RNA pathogens of plants. More recently, numerous and new variations of this ribozyme have been found to inhabit the genomes of organisms from all life kingdoms, although their precise biological functions are not yet well understood.

Thirty-five years of research into ribozymes and nucleic acid catalysis: where do we stand today?

Since the discovery of the first catalytic RNA in 1981, the field of ribozyme research has developed from the discovery of catalytic RNA motifs in nature and the elucidation of their structures and catalytic mechanisms, into a field of engineering and design towards application in diagnostics, molecular biology and medicine. Owing to the development of powerful protocols for selection of nucleic acid catalysts with a desired functionality from random libraries, the spectrum of nucleic acid supported reactions has greatly enlarged, and importantly, ribozymes have been accompanied by DNAzymes. Current areas of research are the engineering of allosteric ribozymes for artificial regulation of gene expression, the design of ribozymes and DNAzymes for medicinal and environmental diagnostics, and the demonstration of RNA world relevant ribozyme activities. In addition, new catalytic motifs or novel genomic locations of known motifs continue to be discovered in all branches of life by the help of high-throughput bioinformatic approaches. Understanding the biological role of the catalytic RNA motifs widely distributed in diverse genetic contexts belongs to the big challenges of future RNA research.

Retrozymes are a unique family of non-autonomous retrotransposons with hammerhead ribozymes that propagate in plants through circular RNAs

Background

Catalytic RNAs, or ribozymes, are regarded as fossils of a prebiotic RNA world that have remained in the genomes of modern organisms. The simplest ribozymes are the small self-cleaving RNAs, like the hammerhead ribozyme, which have been historically considered biological oddities restricted to some RNA pathogens. Recent data, however, indicate that small self-cleaving ribozymes are widespread in genomes, although their functions are still unknown.

Results

We reveal that hammerhead ribozyme sequences in plant genomes form part of a new family of small non-autonomous retrotransposons with hammerhead ribozymes, referred to as retrozymes. These elements contain two long terminal repeats of approximately 350 bp, each harbouring a hammerhead ribozyme that delimitates a variable region of 600–1000 bp with no coding capacity. Retrozymes are actively transcribed, which gives rise to heterogeneous linear and circular RNAs that accumulate differentially depending on the tissue or developmental stage of the plant. Genomic and transcriptomic retrozyme sequences are highly heterogeneous and share almost no sequence homology among species except the hammerhead ribozyme motif and two small conserved domains typical of Ty3-gypsy long terminal repeat retrotransposons. Moreover, we detected the presence of RNAs of both retrozyme polarities, which suggests events of independent RNA-RNA rolling-circle replication and evolution, similarly to that of infectious circular RNAs like viroids and viral satellite RNAs.

Conclusions

Our work reveals that circular RNAs with hammerhead ribozymes are frequently occurring molecules in plant and, most likely, metazoan transcriptomes, which explains the ubiquity of these genomic ribozymes and suggests a feasible source for the emergence of circular RNA plant pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1002-4) contains supplementary material, which is available to authorized users.

Secondary structure confirmation and localization of Mg2+ ions in the mammalian CPEB3 ribozyme

Most of today's knowledge of the CPEB3 ribozyme, one of the few small self-cleaving ribozymes known to occur in humans, is based on comparative studies with the hepatitis delta virus (HDV) ribozyme, which is highly similar in cleavage mechanism and probably also in structure. Here we present detailed NMR studies of the CPEB3 ribozyme in order to verify the formation of the predicted nested double pseudoknot in solution. In particular, the influence of Mg(2+), the ribozyme's crucial cofactor, on the CPEB3 structure is investigated. NMR titrations, Tb(3+)-induced cleavage, as well as stoichiometry determination by hydroxyquinoline sulfonic acid fluorescence and equilibrium dialysis, are used to evaluate the number, location, and binding mode of Mg(2+)ions. Up to eight Mg(2+)ions interact site-specifically with the ribozyme, four of which are bound with high affinity. The global fold of the CPEB3 ribozyme, encompassing 80%-90% of the predicted base pairs, is formed in the presence of monovalent ions alone. Low millimolar concentrations of Mg(2+)promote a more compact fold and lead to the formation of additional structures in the core of the ribozyme, which contains the inner small pseudoknot and the active site. Several Mg(2+)binding sites, which are important for the functional fold, appear to be located in corresponding locations in the HDV and CPEB3 ribozyme, demonstrating the particular relevance of Mg(2+)for the nested double pseudoknot structure.

The structural stability and catalytic activity of DNA and RNA oligonucleotides in the presence of organic solvents

Organic solvents and apolar media are used in the studies of nucleic acids to modify the conformation and function of nucleic acids, to improve solubility of hydrophobic ligands, to construct molecular scaffolds for organic synthesis, and to study molecular crowding effects. Understanding how organic solvents affect nucleic acid interactions and identifying the factors that dominate solvent effects are important for the creation of oligonucleotide-based technologies. This review describes the structural and catalytic properties of DNA and RNA oligonucleotides in organic solutions and in aqueous solutions with organic cosolvents. There are several possible mechanisms underlying the effects of organic solvents on nucleic acid interactions. The reported results emphasize the significance of the osmotic pressure effect and the dielectric constant effect in addition to specific interactions with nucleic acid strands. This review will serve as a guide for the selection of solvent systems based on the purpose of the nucleic acid-based experiments.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Eukaryotic penelope-like retroelements encode hammerhead ribozyme motifs

Small self-cleaving RNAs, such as the paradigmatic Hammerhead ribozyme (HHR), have been recently found widespread in DNA genomes across all kingdoms of life. In this work, we found that new HHR variants are preserved in the ancient family of Penelope-like elements (PLEs), a group of eukaryotic retrotransposons regarded as exceptional for encoding telomerase-like retrotranscriptases and spliceosomal introns. Our bioinformatic analysis revealed not only the presence of minimalist HHRs in the two flanking repeats of PLEs but also their massive and widespread occurrence in metazoan genomes. The architecture of these ribozymes indicates that they may work as dimers, although their low self-cleavage activity in vitro suggests the requirement of other factors in vivo. In plants, however, PLEs show canonical HHRs, whereas fungi and protist PLEs encode ribozyme variants with a stable active conformation as monomers. Overall, our data confirm the connection of self-cleaving RNAs with eukaryotic retroelements and unveil these motifs as a significant fraction of the encoded information in eukaryotic genomes.

The rise of regulatory RNA

Discoveries over the past decade portend a paradigm shift in molecular biology. Evidence suggests that RNA is not only functional as a messenger between DNA and protein but also involved in the regulation of genome organization and gene expression, which is increasingly elaborate in complex organisms. Regulatory RNA seems to operate at many levels; in particular, it plays an important part in the epigenetic processes that control differentiation and development. These discoveries suggest a central role for RNA in human evolution and ontogeny. Here, we review the emergence of the previously unsuspected world of regulatory RNA from a historical perspective.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

Solution structure and metal ion binding sites of the human CPEB3 ribozyme's P4 domain

Three ribozymes are known to occur in humans, the CPEB3 ribozyme, the CoTC ribozyme, and the hammerhead ribozyme. Here, we present the NMR solution structure of a well-conserved motif within the CPEB3 ribozyme, the P4 domain. In addition, we discuss the binding sites and impact of Mg(2+) and [Co(NH3)6](3+), a spectroscopic probe for [Mg(H2O)6](2+), on the structure. The well-defined P4 region is a hairpin closed with a UGGU tetraloop that shows a distinct electrostatic surface potential and a characteristic, strongly curved backbone trajectory. The P4 hairpin contains two specific Mg(2+) binding sites: one outer-sphere binding site close to the proposed CPEB3 ribozyme active site with potential relevance for maintaining a compact fold of the ribozyme core, and one inner-sphere binding site, probably stabilizing the tetraloop structure. The structure of the tetraloop resembles an RNase III recognition structure, as previously described for an AGUU tetraloop. The detailed knowledge of the P4 domain and its metal ion binding preferences thus brings us closer to understanding the importance of Mg(2+) binding for the CPEB3 ribozyme's fold and function in the cell.

Finding instances of riboswitches and ribozymes by homology search of structured RNA with Infernal

In the genomics era, computational tools are essential to extract information from sequences and annotate them to allow easy access to genes. Fortunately, many of these tools are now part of standard pipelines. As a consequence, a cornucopia of genomic features is available in multiple databases. Nevertheless, as novel genomes are sequenced and new structured RNAs are discovered, homology searches and additional analyses need to be performed. In this chapter, we propose simple ways of finding instances of riboswitches and ribozymes in databases or in unannotated genomes, as well as ways of finding variants that deviate from the typical consensus.

Circularity and self-cleavage as a strategy for the emergence of a chromosome in the RNA-based protocell

Background

It is now popularly accepted that an “RNA world” existed in early evolution. During division of RNA-based protocells, random distribution of individual genes (simultaneously as ribozymes) between offspring might have resulted in gene loss, especially when the number of gene types increased. Therefore, the emergence of a chromosome carrying linked genes was critical for the prosperity of the RNA world. However, there were quite a few immediate difficulties for this event to occur. For example, a chromosome would be much longer than individual genes, and thus more likely to degrade and less likely to replicate completely; the copying of the chromosome might start at middle sites and be only partial; and, without a complex transcription mechanism, the synthesis of distinct ribozymes would become problematic.

Results

Inspired by features of viroids, which have been suggested as “living fossils” of the RNA world, we supposed that these difficulties could have been overcome if the chromosome adopted a circular form and small, self-cleaving ribozymes (e.g. the hammer head ribozymes) resided at the sites between genes. Computer simulation using a Monte-Carlo method was conducted to investigate this hypothesis. The simulation shows that an RNA chromosome can spread (increase in quantity and be sustained) in the system if it is a circular one and its linear “transcripts” are readily broken at the sites between genes; the chromosome works as genetic material and ribozymes “coded” by it serve as functional molecules; and both circularity and self-cleavage are important for the spread of the chromosome.

Conclusions

In the RNA world, circularity and self-cleavage may have been adopted as a strategy to overcome the immediate difficulties for the emergence of a chromosome (with linked genes). The strategy suggested here is very simple and likely to have been used in this early stage of evolution. By demonstrating the possibility of the emergence of an RNA chromosome, this study opens on the prospect of a prosperous RNA world, populated by RNA-based protocells with a number of genes, showing complicated functions.

Reviewers

This article was reviewed by Sergei Kazakov (nominated by Laura Landweber), Nobuto Takeuchi (nominated by Anthony Poole), and Eugene Koonin.

The ubiquitous hammerhead ribozyme

The hammerhead ribozyme is a small catalytic RNA motif capable of endonucleolytic (self-) cleavage. It is composed of a catalytic core of conserved nucleotides flanked by three helices, two of which form essential tertiary interactions for fast self-scission under physiological conditions. Originally discovered in subviral plant pathogens, its presence in several eukaryotic genomes has been reported since. More recently, this catalytic RNA motif has been shown to reside in a large number of genomes. We review the different approaches in discovering these new hammerhead ribozyme sequences and discuss possible biological functions of the genomic motifs.

Site-specific platinum(II) cross-linking in a ribozyme active site

The function of RNA depends on its ability to adopt complex and dynamic structures, and the incorporation of site-specific cross-linking probes is a powerful method for providing distance constraints that are valuable in RNA structural biology. Here we describe a new RNA-RNA cross-linking strategy based on Pt(II) targeting of specific phosphorothioate substitutions. In this strategy cis-diammine Pt(II) complexes are kinetically recruited and anchored to a phosphorothioate substitution embedded within a structured RNA. Substitution of the remaining exchangeable Pt(II) ligand with a nucleophile supplied by a nearby RNA nucleobase results in metal-mediated cross-links that are stable during isolation. This type of cross-linking strategy was explored within the catalytic core of the Hammerhead ribozyme (HHRz). When a phosphorothioate substitution is installed at the scissile bond normally cleaved by the HHRz, Pt(II) cross-linking takes place to nucleotides G8 and G10 in the ribozyme active site. Both of these positions are predicted to be within ~8 Å of a phosphorothioate-bound Pt(II) metal center. Cross-linking depends on Mg(2+) ion concentration, reaching yields as high as 30%, with rates that indicate cation competition within the RNA three-helix junction. Cross-linking efficiency depends on accurate formation of the HHRz tertiary structure, and cross-links are not observed for RNA helices. Combined, these results show promise for using kinetically inert Pt(II) complexes as new site-specific cross-linking tools for exploring RNA structure and dynamics.

Non coding RNAs and viruses in the framework of the phylogeny of the genes, epigenesis and heredity

The origin of genes is one of the most enigmatic events in the origin of life. It has been suggested that noncoding (nc) RNA was probably a precursor in the formation of the first polypeptide, and also at the origin of the first manifestation of life and genes. ncRNAs are also becoming central for understanding gene expression and silencing. Indeed, before the discovery of ncRNAs, proteins were viewed as the major molecules in the regulation of gene expression and gene silencing; however, recent findings suggest that ncRNA also plays an important role in gene expression. Reverse transcription of RNA viruses and their integration into the genome of eukaryotes and also their relationship with the ncRNA suggest that their origin is basal in genome evolution, and also probably constitute the first mechanism of gene regulation. I am to review the different roles of ncRNAs in the framework of gene evolution, as well as the importance of ncRNAs and viruses in the epigenesis and in the non-Mendelian model of heredity and evolution.

Structure-based search and in vitro analysis of self-cleaving ribozymes

Detecting functional RNAs is increasingly accomplished through structure-based searches for patterns of conserved secondary structure. With large amounts of new sequencing data becoming available, there is a greater demand for efficient methods of identifying new RNAs. Here we present a method of identifying self-cleaving ribozymes and characterizing the in vitro activity.

Characterization of hammerhead ribozyme reactions

Hammerhead ribozymes are small catalytic RNA motifs ubiquitously present in a large variety of genomes. The reactions catalyzed by these motifs are both their self-scission and the reverse ligation reaction. Here, we describe methods for the generation of DNA templates for the subsequent in vitro transcription of hammerhead ribozymes. This is followed by a description of the preparation of suitable RNA molecules for both reaction types, and their kinetic analysis.

Processing and translation initiation of non-long terminal repeat retrotransposons by hepatitis delta virus (HDV)-like self-cleaving ribozymes

Many non-long terminal repeat (non-LTR) retrotransposons lack internal promoters and are co-transcribed with their host genes. These transcripts need to be liberated before inserting into new loci. Using structure-based bioinformatics, we show that several classes of retrotransposons in phyla-spanning arthropods, nematodes, and chordates utilize self-cleaving ribozymes of the hepatitis delta virus (HDV) family for processing their 5' termini. Ribozyme-terminated retrotransposons include rDNA-specific R2, R4, and R6, telomere-specific SART, and Baggins and RTE. The self-scission of the R2 ribozyme is strongly modulated by the insertion site sequence in the rDNA, with the most common insertion sequences promoting faster processing. The ribozymes also promote translation initiation of downstream open reading frames in vitro and in vivo. In some organisms HDV-like and hammerhead ribozymes appear to be dedicated to processing long and short interspersed elements, respectively. HDV-like ribozymes serve several distinct functions in non-LTR retrotransposition, including 5' processing, translation initiation, and potentially trans-templating.