Eukaryotic penelope-like retroelements encode hammerhead ribozyme motifs

Viroids and Retrozymes: Plant Circular RNAs Capable of Autonomous Replication

Among the long non-coding RNAs that are currently recognized as important regulatory molecules influencing a plethora of processes in eukaryotic cells, circular RNAs (circRNAs) represent a distinct class of RNAs that are predominantly produced by back-splicing of pre-mRNA. The most studied regulatory mechanisms involving circRNAs are acting as miRNA sponges, forming R-loops with genomic DNA, and encoding functional proteins. In addition to circRNAs generated by back-splicing, two types of circRNAs capable of autonomous RNA-RNA replication and systemic transport have been described in plants: viroids, which are infectious RNAs that cause a number of plant diseases, and retrozymes, which are transcripts of retrotransposon genomic loci that are capable of circularization due to ribozymes. Based on a number of common features, viroids and retrozymes are considered to be evolutionarily related. Here, we provide an overview of the biogenesis mechanisms and regulatory functions of non-replicating circRNAs produced by back-splicing and further discuss in detail the currently available data on viroids and retrozymes, focusing on their structural features, replication mechanisms, interaction with cellular components, and transport in plants. In addition, biotechnological approaches involving replication-capable plant circRNAs are discussed, as well as their potential applications in research and agriculture.

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.

Internal initiation of reverse transcription in a Penelope-like retrotransposon

Eukaryotic retroelements are generally divided into two classes: long terminal repeat (LTR) retrotransposons and non-LTR retrotransposons. A third class of eukaryotic retroelement, the Penelope-like elements (PLEs), has been well-characterized bioinformatically, but relatively little is known about the transposition mechanism of these elements. PLEs share some features with the R2 retrotransposon from Bombyx mori, which uses a target-primed reverse transcription (TPRT) mechanism, but their distinct phylogeny suggests PLEs may utilize a novel mechanism of mobilization. Using protein purified from E. coli, we report unique in vitro properties of a PLE from the green anole (Anolis carolinensis), revealing mechanistic aspects not shared by other retrotransposons. We found that reverse transcription is initiated at two adjacent sites within the transposon RNA that is not homologous to the cleaved DNA, a feature that is reflected in the genomic "tail" signature shared between and unique to PLEs. Our results for the first active PLE in vitro provide a starting point for understanding PLE mobilization and biology.

In-Plant Persistence and Systemic Transport of Nicotiana benthamiana Retrozyme RNA

Retrozymes are nonautonomous retrotransposons with hammerhead ribozymes in their long terminal repeats (LTRs). Retrozyme transcripts can be self-cleaved by the LTR ribozyme, circularized, and can undergo RNA-to-RNA replication. Here, we demonstrate that the Nicotiana benthamiana genome contains hundreds of retrozyme loci, of which nine represent full-length retrozymes. The LTR contains a promoter directing retrozyme transcription. Although retrozyme RNA is easily detected in plants, the LTR region is heavily methylated, pointing to its transcriptional silencing, which can be mediated by 24 nucleotide-long retrozyme-specific RNAs identified in N. benthamiana. A transcriptome analysis revealed that half of the retrozyme-specific RNAs in plant leaves have no exact matches to genomic retrozyme loci, containing up to 13% mismatches with the closest genomic sequences, and could arise as a result of many rounds of RNA-to-RNA replication leading to error accumulation. Using a cloned retrozyme copy, we show that retrozyme RNA is capable of replication and systemic transport in plants. The presented data suggest that retrozyme loci in the N. benthamiana genome are transcriptionally inactive, and that circular retrozyme RNA can persist in cells due to its RNA-to-RNA replication and be transported systemically, emphasizing functional and, possibly, evolutionary links of retrozymes to viroids-noncoding circular RNAs that infect plants.

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.

cyPhyRNA-seq: a genome-scale RNA-seq method to detect active self-cleaving ribozymes by capturing RNAs with 2',3' cyclic phosphates and 5' hydroxyl ends

Self-cleaving ribozymes are catalytically active RNAs that cleave themselves into a 5′-fragment with a 2′,3′-cyclic phosphate and a 3′-fragment with a 5′-hydroxyl. They are widely applied for the construction of synthetic RNA devices and RNA-based therapeutics. However, the targeted discovery of self-cleaving ribozymes remains a major challenge. We developed a transcriptome-wide method, called cyPhyRNA-seq, to screen for ribozyme cleavage fragments in total RNA extract. This approach employs the specific ligation-based capture of ribozyme 5′-fragments using a variant of the Arabidopsis thaliana tRNA ligase we engineered. To capture ribozyme 3′-fragments, they are enriched from total RNA by enzymatic treatments. We optimized and enhanced the individual steps of cyPhyRNA-seq in vitro and in spike-in experiments. Then, we applied cyPhyRNA-seq to total RNA isolated from the bacterium Desulfovibrio vulgaris and detected self-cleavage of the three predicted type II hammerhead ribozymes, whose activity had not been examined to date. cyPhyRNA-seq can be used for the global analysis of active self-cleaving ribozymes with the advantage to capture both ribozyme cleavage fragments from total RNA. Especially in organisms harbouring many self-cleaving RNAs, cyPhyRNA-seq facilitates the investigation of cleavage activity. Moreover, this method has the potential to be used to discover novel self-cleaving ribozymes in different organisms.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

An ancient clade of Penelope-like retroelements with permuted domains is present in the green lineage and protists, and dominates many invertebrate genomes

Penelope-like elements (PLEs) are an enigmatic clade of retrotransposons whose reverse transcriptases (RTs) share a most recent common ancestor with telomerase RTs. The single ORF of canonical endonuclease (EN)+ PLEs encodes RT and a C-terminal GIY–YIG EN that enables intrachromosomal integration, whereas EN− PLEs lack EN and are generally restricted to chromosome termini. EN+ PLEs have only been found in animals, except for one case of horizontal transfer to conifers, whereas EN− PLEs occur in several kingdoms. Here, we report a new, deep-branching PLE clade with a permuted domain order, whereby an N-terminal GIY–YIG EN is linked to a C-terminal RT by a short domain with a characteristic CxC motif. These N-terminal EN+ PLEs share a structural organization, including pseudo-LTRs and complex tandem/inverted insertions, with canonical EN+ PLEs from Penelope/Poseidon, Neptune, and Nematis clades, and show insertion bias for microsatellites, but lack canonical hammerhead ribozyme motifs. However, their phylogenetic distribution is much broader. The Naiads, found in numerous invertebrate phyla, can reach tens of thousands of copies per genome. In spiders and clams, Naiads independently evolved to encode selenoproteins containing multiple selenocysteines. Chlamys, which lack the CCHH motif universal to PLE ENs, occur in green algae, spike mosses (targeting ribosomal DNA), and slime molds. Unlike canonical PLEs, RTs of N-terminal EN+ PLEs contain the insertion-in-fingers domain (IFD), strengthening the link between PLEs and telomerases. Additionally, we describe Hydra, a novel metazoan C-terminal EN+ clade. Overall, we conclude that PLE diversity, taxonomic distribution, and abundance are comparable with non-LTR and LTR-retrotransposons.

Biochemical analysis of cleavage and ligation activities of the pistol ribozyme from Paenibacillus polymyxa

Nine distinct classes of self-cleaving ribozymes are known to date, of which the pistol ribozyme class was discovered only 5 years ago. Self-cleaving ribozymes are able to cleave their own phosphodiester backbone at a specific site with rates much higher than those of spontaneous RNA degradation. Our study focuses on a bioinformatically predicted pistol ribozyme from the bacterium Paenibacillus polymyxa. We provide a biochemical characterization of this ribozyme, which includes an investigation of the effect of various metal ions on ribozyme cleavage and a kinetic analysis of ribozyme activity under increasing Mg²⁺ concentrations and pH. Based on the obtained results, we discuss a possible catalytic role of divalent metal ions. Moreover, we investigated the ligation activity of the P. polymyxa pistol ribozyme - an aspect that has not been previously analysed for this ribozyme class. We determined that the P. polymyxa pistol ribozyme is almost fully cleaved at equilibrium with the ligation rate constant being nearly 30-fold lower than the cleavage rate constant. In summary, we have characterized an additional representative of this recently discovered ribozyme class isolated from P. polymyxa. We expect that our biochemical characterization of a pistol representative in a cultivatable, genetically tractable organism will support our future investigation of the biological roles of this ribozyme class in bacteria.

Hepatitis delta virus-like circular RNAs from diverse metazoans encode conserved hammerhead ribozymes

Human hepatitis delta virus (HDV) is a unique infectious agent whose genome is composed of a small circular RNA. Recent data, however, have reported the existence of highly divergent HDV-like circRNAs in the transcriptomes of diverse vertebrate and invertebrate species. The HDV-like genomes described in amniotes such as birds and reptiles encode self-cleaving RNA motifs or ribozymes similar to the ones present in the human HDV, whereas no catalytic RNA domains have been reported for the HDV-like genomes detected in metagenomic data from some amphibians, fish, and invertebrates. Herein, we describe the self-cleaving motifs of the HDV-like genomes reported in newts and fish, which belong to the characteristic class of HDV ribozymes. Surprisingly, HDV-like genomes from a toad and a termite show conserved type III hammerhead ribozymes, which belong to an unrelated class of catalytic RNAs characteristic of plant genomes and plant subviral circRNAs, such as some viral satellites and viroids. Sequence analyses revealed the presence of similar HDV-like hammerhead ribozymes encoded in two termite genomes, but also in the genomes of several dipteran species. In vitro transcriptions confirmed the cleaving activity for these motifs, with moderate rates of self-cleavage. These data indicate that all described HDV-like agents contain self-cleaving motifs from either the HDV or the hammerhead class. Autocatalytic ribozymes in HDV-like genomes could be regarded as interchangeable domains and may have arisen from cellular transcriptomes, although we still cannot rule out some other evolutionary explanations.

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.

Telomere-Specialized Retroelements in Drosophila: Adaptive Symbionts of the Genome, Neutral, or in Conflict?

Linear chromosomes shorten in every round of replication. In Drosophila, telomere-specialized long interspersed retrotransposable elements (LINEs) belonging to the jockey clade offset this shortening by forming head-to-tail arrays at Drosophila telomere ends. As such, these telomeric LINEs have been considered adaptive symbionts of the genome, protecting it from premature decay, particularly as Drosophila lacks a conventional telomerase holoenzyme. However, as reviewed here, recent work reveals a high degree of variation and turnover in the telomere-specialized LINE lineages across Drosophila. There appears to be no absolute requirement for LINE activity to maintain telomeres in flies, hence the suggestion that the telomere-specialized LINEs may instead be neutral or in conflict with the host, rather than adaptive.

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.

Viroid research and its significance for RNA technology and basic biochemistry

Viroids were described 47 years ago as the smallest RNA molecules capable of infecting plants and autonomously self-replicating without an encoded protein. Work on viroids initiated the development of a number of innovative methods. Novel chromatographic and gelelectrophoretic methods were developed for the purification and characterization of viroids; these methods were later used in molecular biology, gene technology and in prion research. Theoretical and experimental studies of RNA folding demonstrated the general biological importance of metastable structures, and nuclear magnetic resonance spectroscopy of viroid RNA showed the partially covalent nature of hydrogen bonds in biological macromolecules. RNA biochemistry and molecular biology profited from viroid research, such as in the detection of RNA as template of DNA-dependent polymerases and in mechanisms of gene silencing. Viroids, the first circular RNA detected in nature, are important for studies on the much wider spectrum of circular RNAs and other non-coding RNAs.

MGERT: a pipeline to retrieve coding sequences of mobile genetic elements from genome assemblies

Background

Genomes of eukaryotes are inhabited by myriads of mobile genetic elements (MGEs) – transposons and retrotransposons - which play a great role in genome plasticity and evolution. A lot of computational tools were developed to annotate them either in genomic assemblies or raw reads using de novo or homology-based approaches. But there has been no pipeline enabling users to get coding and flanking sequences of MGEs suitable for a downstream analysis from genome assemblies.

Results

We developed a new pipeline, MGERT (Mobile Genetic Elements Retrieving Tool), that automates all the steps necessary to obtain protein-coding sequences of mobile genetic elements from genomic assemblies even if no previous knowledge on MGE content of a particular genome is available.

Conclusions

Using MGERT, researchers can easily find MGEs, their coding and flanking sequences in the genome of interest. Thus, this pipeline helps researchers to focus on the biological analysis of MGEs rather than excessive scripting and pipelining.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0163-6) contains supplementary material, which is available to authorized users.

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.

Giant Transposons in Eukaryotes: Is Bigger Better?

Transposable elements (TEs) are ubiquitous in both prokaryotes and eukaryotes, and the dynamic character of their interaction with host genomes brings about numerous evolutionary innovations and shapes genome structure and function in a multitude of ways. In traditional classification systems, TEs are often being depicted in simplistic ways, based primarily on the key enzymes required for transposition, such as transposases/recombinases and reverse transcriptases. Recent progress in whole-genome sequencing and long-read assembly, combined with expansion of the familiar range of model organisms, resulted in identification of unprecedentedly long transposable units spanning dozens or even hundreds of kilobases, initially in prokaryotic and more recently in eukaryotic systems. Here, we focus on such oversized eukaryotic TEs, including retrotransposons and DNA transposons, outline their complex and often combinatorial nature and closely intertwined relationship with viruses, and discuss their potential for participating in transfer of long stretches of DNA in eukaryotes.

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.

Giant Reverse Transcriptase-Encoding Transposable Elements at Telomeres

Transposable elements are omnipresent in eukaryotic genomes and have a profound impact on chromosome structure, function and evolution. Their structural and functional diversity is thought to be reasonably well-understood, especially in retroelements, which transpose via an RNA intermediate copied into cDNA by the element-encoded reverse transcriptase, and are characterized by a compact structure. Here, we report a novel type of expandable eukaryotic retroelements, which we call Terminons. These elements can attach to G-rich telomeric repeat overhangs at the chromosome ends, in a process apparently facilitated by complementary C-rich repeats at the 3′-end of the RNA template immediately adjacent to a hammerhead ribozyme motif. Terminon units, which can exceed 40 kb in length, display an unusually complex and diverse structure, and can form very long chains, with host genes often captured between units. As the principal polymerizing component, Terminons contain Athena reverse transcriptases previously described in bdelloid rotifers and belonging to the enigmatic group of Penelope-like elements, but can additionally accumulate multiple cooriented ORFs, including DEDDy 3′-exonucleases, GDSL esterases/lipases, GIY-YIG-like endonucleases, rolling-circle replication initiator (Rep) proteins, and putatively structural ORFs with coiled-coil motifs and transmembrane domains. The extraordinary length and complexity of Terminons and the high degree of interfamily variability in their ORF content challenge the current views on the structural organization of eukaryotic retroelements, and highlight their possible connections with the viral world and the implications for the elevated frequency of gene transfer.

Numerous small hammerhead ribozyme variants associated with Penelope-like retrotransposons cleave RNA as dimers

Hammerhead ribozymes represent the most common of the 9 natural classes of self-cleaving RNAs. The hammerhead catalytic core includes 11 highly-conserved nucleotides located largely within the unpaired regions of a junction formed by stems I, II and III. The vast majority of previously reported examples carry an additional pseudoknot or other tertiary interactions between nucleotides that precede stem I and nucleotides in the loop of stem II. These extra contacts are critical for high-speed RNA catalysis. Herein, we report the discovery of ∼150,000 additional variant hammerhead representatives that exhibit diminished stem III substructures. These variants are frequently associated with Penelope-like retrotransposons, which are a type of mobile genetic element. Kinetic analyses indicate that these RNAs form dimers to cleave RNA.

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.

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.

Hammerhead ribozymes going viral

An association between hammerhead ribozymes and non-autonomous, long terminal repeat retrotransposons is uncovered in plants, shedding light on the biological function of genomically encoded ribozymes.

An Ancient Transkingdom Horizontal Transfer of Penelope-Like Retroelements from Arthropods to Conifers

Comparative genomics analyses empowered by the wealth of sequenced genomes have revealed numerous instances of horizontal DNA transfers between distantly related species. In eukaryotes, repetitive DNA sequences known as transposable elements (TEs) are especially prone to move across species boundaries. Such horizontal transposon transfers, or HTTs, are relatively common within major eukaryotic kingdoms, including animals, plants, and fungi, while rarely occurring across these kingdoms. Here, we describe the first case of HTT from animals to plants, involving TEs known as Penelope-like elements, or PLEs, a group of retrotransposons closely related to eukaryotic telomerases. Using a combination of in situ hybridization on chromosomes, polymerase chain reaction experiments, and computational analyses we show that the predominant PLE lineage, EN(+)PLEs, is highly diversified in loblolly pine and other conifers, but appears to be absent in other gymnosperms. Phylogenetic analyses of both protein and DNA sequences reveal that conifers EN(+)PLEs, or Dryads, form a monophyletic group clustering within a clade of primarily arthropod elements. Additionally, no EN(+)PLEs were detected in 1,928 genome assemblies from 1,029 nonmetazoan and nonconifer genomes from 14 major eukaryotic lineages. These findings indicate that Dryads emerged following an ancient horizontal transfer of EN(+)PLEs from arthropods to a common ancestor of conifers approximately 340 Ma. This represents one of the oldest known interspecific transmissions of TEs, and the most conspicuous case of DNA transfer between animals and plants.

Two Divalent Metal Ions and Conformational Changes Play Roles in the Hammerhead Ribozyme Cleavage Reaction

The hammerhead ribozyme is a self-cleaving RNA broadly dispersed across all kingdoms of life. Although it was the first of the small, nucleolytic ribozymes discovered, the mechanism by which it catalyzes its reaction remains elusive. The nucleobase of G12 is well positioned to be a general base, but it is unclear if or how this guanine base becomes activated for proton transfer. Metal ions have been implicated in the chemical mechanism, but no interactions between divalent metal ions and the cleavage site have been observed crystallographically. To better understand how this ribozyme functions, we have solved crystal structures of wild-type and G12A mutant ribozymes. We observe a pH-dependent conformational change centered around G12, consistent with this nucleotide becoming deprotonated. Crystallographic and kinetic analysis of the G12A mutant reveals a Zn(2+) specificity switch suggesting a direct interaction between a divalent metal ion and the purine at position 12. The metal ion specificity switch and the pH-rate profile of the G12A mutant suggest that the minor imino tautomer of A12 serves as the general base in the mutant ribozyme. We propose a model in which the hammerhead ribozyme rearranges prior to the cleavage reaction, positioning two divalent metal ions in the process. The first metal ion, positioned near G12, becomes directly coordinated to the O6 keto oxygen, to lower the pKa of the general base and organize the active site. The second metal ion, positioned near G10.1, bridges the N7 of G10.1 and the scissile phosphate and may participate directly in the cleavage reaction.

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.