R2 retrotransposons encode a self-cleaving ribozyme for processing from an rRNA cotranscript

Structure and sequence at an RNA template 5' end influence insertion of transgenes by an R2 retrotransposon protein

R2 non-long terminal repeat retrotransposons insert site-specifically into ribosomal RNA genes (rDNA) in a broad range of multicellular eukaryotes. R2-encoded proteins can be leveraged to mediate transgene insertion at 28S rDNA loci in cultured human cells. This strategy, precise RNA-mediated insertion of transgenes (PRINT), relies on the codelivery of an mRNA encoding R2 protein and an RNA template encoding a transgene cassette of choice. Here, we demonstrate that the PRINT RNA template 5' module, which as a complementary DNA 3' end will generate the transgene 5' junction with rDNA, influences the efficiency and mechanism of gene insertion. Iterative design and testing identified optimal 5' modules consisting of a hepatitis delta virus-like ribozyme fold with high thermodynamic stability, suggesting that RNA template degradation from its 5' end may limit transgene insertion efficiency. We also demonstrate that transgene 5' junction formation can be either precise, formed by annealing the 3' end of first-strand complementary DNA with the upstream target site, or imprecise, by end-joining, but this difference in junction formation mechanism is not a major determinant of insertion efficiency. Sequence characterization of imprecise end-joining events indicates surprisingly minimal reliance on microhomology. Our findings expand the current understanding of the role of R2 retrotransposon transcript sequence and structure, and especially the 5' ribozyme fold, for retrotransposon mobility and RNA-templated gene synthesis in cells.

All-RNA-mediated targeted gene integration in mammalian cells with rationally engineered R2 retrotransposons

All-RNA-mediated targeted gene integration methods, rendering reduced immunogenicity, effective deliverability with non-viral vehicles, and a low risk of random mutagenesis, are urgently needed for next-generation gene addition technologies. Naturally occurring R2 retrotransposons hold promise in this context due to their site-specific integration profile. Here, we systematically analyzed the biodiversity of R2 elements and screened several R2 orthologs capable of full-length gene insertion in mammalian cells. Robust R2 system gene integration efficiency was attained using combined donor RNA and protein engineering. Importantly, the all-RNA-delivered engineered R2 system showed effective integration activity, with efficiency over 60% in mouse embryos. Unbiased high-throughput sequencing demonstrated that the engineered R2 system exhibited high on-target integration specificity (99%). In conclusion, our study provides engineered R2 tools for applications based on hit-and-run targeted DNA integration and insights for further optimization of retrotransposon systems.

Distinct and overlapping RNA determinants for binding and target-primed reverse transcription by Bombyx mori R2 retrotransposon protein

Eukaryotic retrotransposons encode a reverse transcriptase that binds RNA to template DNA synthesis. The ancestral non-long terminal repeat (non-LTR) retrotransposons encode a protein that performs target-primed reverse transcription (TPRT), in which the nicked genomic target site initiates complementary DNA (cDNA) synthesis directly into the genome. The best understood model system for biochemical studies of TPRT is the R2 protein from the silk moth Bombyx mori. The R2 protein selectively binds the 3' untranslated region of its encoding RNA as template for DNA insertion to its target site in 28S ribosomal DNA. Here, binding and TPRT assays define RNA contributions to RNA-protein interaction, template use for TPRT and the fidelity of template positioning for TPRT cDNA synthesis. We quantify both sequence and structure contributions to protein-RNA interaction. RNA determinants of binding affinity overlap but are not equivalent to RNA features required for TPRT and its fidelity of template positioning for full-length TPRT cDNA synthesis. Additionally, we show that a previously implicated RNA-binding protein surface of R2 protein makes RNA binding affinity dependent on the presence of two stem-loops. Our findings inform evolutionary relationships across R2 retrotransposon RNAs and are a step toward understanding the mechanism and template specificity of non-LTR retrotransposon mobility.

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

Harnessing eukaryotic retroelement proteins for transgene insertion into human safe-harbor loci

Current approaches for inserting autonomous transgenes into the genome, such as CRISPR-Cas9 or virus-based strategies, have limitations including low efficiency and high risk of untargeted genome mutagenesis. Here, we describe precise RNA-mediated insertion of transgenes (PRINT), an approach for site-specifically primed reverse transcription that directs transgene synthesis directly into the genome at a multicopy safe-harbor locus. PRINT uses delivery of two in vitro transcribed RNAs: messenger RNA encoding avian R2 retroelement-protein and template RNA encoding a transgene of length validated up to 4 kb. The R2 protein coordinately recognizes the target site, nicks one strand at a precise location and primes complementary DNA synthesis for stable transgene insertion. With a cultured human primary cell line, over 50% of cells can gain several 2 kb transgenes, of which more than 50% are full-length. PRINT advantages include no extragenomic DNA, limiting risk of deleterious mutagenesis and innate immune responses, and the relatively low cost, rapid production and scalability of RNA-only delivery.

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.

Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon

Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.

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.

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.

Impaired function of rDNA transcription initiation machinery leads to derepression of ribosomal genes with insertions of R2 retrotransposon

Eukaryotic genomes harbor hundreds of rRNA genes, many of which are transcriptionally silent. However, little is known about selective regulation of individual rDNA units. In Drosophila melanogaster, some rDNA repeats contain insertions of the R2 retrotransposon, which is capable to be transcribed only as part of pre-rRNA molecules. rDNA units with R2 insertions are usually inactivated, although R2 expression may be beneficial in cells with decreased rDNA copy number. Here we found that R2-inserted rDNA units are enriched with HP1a and H3K9me3 repressive mark, whereas disruption of the heterochromatin components slightly affects their silencing in ovarian germ cells. Surprisingly, we observed a dramatic upregulation of R2-inserted rRNA genes in ovaries lacking Udd (Under-developed) or other subunits (TAF1b and TAF1c-like) of the SL1-like complex, which is homologues to mammalian Selective factor 1 (SL1) involved in rDNA transcription initiation. Derepression of rRNA genes with R2 insertions was accompanied by a reduction of H3K9me3 and HP1a enrichment. We suggest that the impairment of the SL1-like complex affects a mechanism of selective activation of intact rDNA units which competes with heterochromatin formation. We also propose that R2 derepression may serve as an adaptive response to compromised rRNA synthesis.

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.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

Hepatitis Delta Virus (HDV) and Delta-Like Agents: Insights Into Their Origin

Hepatitis delta virus (HDV) is a human pathogen, and the only known species in the genus Deltavirus. HDV is a satellite virus and depends on the hepatitis B virus (HBV) for packaging, release, and transmission. Extracellular HDV virions contain the genomic HDV RNA, a single-stranded negative-sense and covalently closed circular RNA molecule, which is associated with the HDV-encoded delta antigen forming a ribonucleoprotein complex, and enveloped by the HBV surface antigens. Replication occurs in the nucleus and is mediated by host enzymes and assisted by cis-acting ribozymes allowing the formation of monomer length molecules which are ligated by host ligases to form unbranched rod-like circles. Recently, meta-transcriptomic studies investigating various vertebrate and invertebrate samples identified RNA species with similarities to HDV RNA. The delta-like agents may be representatives of novel subviral agents or satellite viruses which share with HDV, the self-complementarity of the circular RNA genome, the ability to encode a protein, and the presence of ribozyme sequences. The widespread distribution of delta-like agents across different taxa with considerable phylogenetic distances may be instrumental in comprehending their evolutionary history by elucidating the transition from transcriptome to cellular circular RNAs to infectious subviral agents.

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.

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.

Repression of interrupted and intact rDNA by the SUMO pathway in Drosophila melanogaster

Ribosomal RNAs (rRNAs) are essential components of the ribosome and are among the most abundant macromolecules in the cell. To ensure high rRNA level, eukaryotic genomes contain dozens to hundreds of rDNA genes, however, only a fraction of the rRNA genes seems to be active, while others are transcriptionally silent. We found that individual rDNA genes have high level of cell-to-cell heterogeneity in their expression in Drosophila melanogaster. Insertion of heterologous sequences into rDNA leads to repression associated with reduced expression in individual cells and decreased number of cells expressing rDNA with insertions. We found that SUMO (Small Ubiquitin-like Modifier) and SUMO ligase Ubc9 are required for efficient repression of interrupted rDNA units and variable expression of intact rDNA. Disruption of the SUMO pathway abolishes discrimination of interrupted and intact rDNAs and removes cell-to-cell heterogeneity leading to uniformly high expression of individual rDNA in single cells. Our results suggest that the SUMO pathway is responsible for both repression of interrupted units and control of intact rDNA expression.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.

Completion of LINE integration involves an open '4-way' branched DNA intermediate

Long Interspersed Elements (LINEs), also known as non-LTR retrotransposons, encode a multifunctional protein that reverse transcribes its mRNA into DNA at the site of insertion by target primed reverse transcription. The second half of the integration reaction remains very poorly understood. Second-strand DNA cleavage and second-strand DNA synthesis were investigated in vitro using purified components from a site-specific restriction-like endonuclease (RLE) bearing LINE. DNA structure was shown to be a critical component of second-strand DNA cleavage. A hitherto unknown and unexplored integration intermediate, an open ‘4-way’ DNA junction, was recognized by the element protein and cleaved in a Holliday junction resolvase-like reaction. Cleavage of the 4-way junction resulted in a natural primer-template pairing used for second-strand DNA synthesis. A new model for RLE LINE integration is presented.

Sequence-specific retrotransposition of 28S rDNA-specific LINE R2Ol in human cells

R2 is a long interspersed element (LINE) found in a specific sequence of the 28S rDNA among a wide variety of animals. Recently, we observed that R2Ol isolated from medaka fish, Oryzias latipes, retrotransposes sequence specifically into the target sequence of zebrafish. Because the 28S target and flanking regions are widely conserved among vertebrates, we examined whether R2Ol can also integrate in a sequence-specific manner in human cells. Using adenovirus-mediated expression of R2Ol constructs, we confirmed an accurate insertion of R2Ol into the 28S target of human 293T cells. However, the R2Ol mutant devoid of endonuclease (EN) activity showed no retrotransposition ability, suggesting that the sequence-specific integration of R2Ol into 28S rDNA occurs via the cleavage activity of EN. By introducing both R2Ol helper virus and donor plasmid in human cells, we succeeded in retrotransposing an exogenous EGFP gene into the 28S target site by the trans-complementation system, which enabled simplification of specific gene knock-in in a time-efficient manner. We believe that R2Ol may provide an alternative targeted gene knock-in method for practical applications such as gene therapy in future.

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.

Mutations in Cytosine-5 tRNA Methyltransferases Impact Mobile Element Expression and Genome Stability at Specific DNA Repeats

The maintenance of eukaryotic genome stability is ensured by the interplay of transcriptional as well as post-transcriptional mechanisms that control recombination of repeat regions and the expression and mobility of transposable elements. We report here that mutations in two (cytosine-5) RNA methyltransferases, Dnmt2 and NSun2, impact the accumulation of mobile element-derived sequences and DNA repeat integrity in Drosophila. Loss of Dnmt2 function caused moderate effects under standard conditions, while heat shock exacerbated these effects. In contrast, NSun2 function affected mobile element expression and genome integrity in a heat shock-independent fashion. Reduced tRNA stability in both RCMT mutants indicated that tRNA-dependent processes affected mobile element expression and DNA repeat stability. Importantly, further experiments indicated that complex formation with RNA could also contribute to the impact of RCMT function on gene expression control. These results thus uncover a link between tRNA modification enzymes, the expression of repeat DNA, and genomic integrity.

The Role of Sequence Duplication in Transcriptional Regulation and Genome Evolution

Sequence duplication is nowadays recognized as an important mechanism that underlies the evolution of eukaryote genomes, being indeed one of the most powerful strategies for the generation of adaptive diversity by modulating transcriptional activity. The evolutionary novelties simultaneously associated with sequence duplication and differential gene expression can be collectively referred to as duplication-mediated transcriptional regulation. In the last years, evidence has emerged supporting the idea that sequence duplication and functionalization represent important evolutionary strategies acting at the genome level, and both coding and non-coding sequences have been found to be targets of such events. Moreover, it has been proposed that deleterious effects of sequence duplication might be potentially silenced by endogenous cell machinery (i.e., RNA interference, epigenetic repressive marks, etc). Along these lines, our aim is to highlight the role of sequence duplication on transcriptional activity and the importance of both in genome evolution.

LINEs Contribute to the Origins of Middle Bodies of SINEs besides 3' Tails

Short interspersed elements (SINEs), which are nonautonomous transposable elements, require the transposition machinery of long interspersed elements (LINEs) to mobilize. SINEs are composed of two or more independently originating parts. The 5′ region is called the “head” and is derived mainly from small RNAs, and the 3′ region (“tail”) originates from the 3′ region of LINEs and is responsible for being recognized by counterpart LINE proteins. The origin of the middle “body” of SINEs is enigmatic, although significant sequence similarities among SINEs from very diverse species have been observed. Here, a systematic analysis of the similarities among SINEs and LINEs deposited on Repbase, a comprehensive database of eukaryotic repeat sequences was performed. Three primary findings are described: 1) The 5′ regions of only two clades of LINEs, RTE and Vingi, were revealed to have contributed to the middle parts of SINEs; 2) The linkage of the 5′ and 3′ parts of LINEs can be lost due to occasional tail exchange of SINEs; and 3) The previously proposed Ceph-domain was revealed to be a fusion of a CORE-domain and a 5′ part of RTE clade of LINE. Based on these findings, a hypothesis that the 5′ parts of bipartite nonautonomous LINEs, which possess only the 5′ and 3′ regions of the original LINEs, can contribute to the undefined middle part of SINEs is proposed.

R1 retrotransposons in the nucleolar organizers of Drosophila melanogaster are transcribed by RNA polymerase I upon heat shock

The ribosomal RNA genes (rDNA) of Drosophila melanogaster reside within centromere-proximal nucleolar organizers on both the X and Y chromosomes. Each locus contains between 200-300 tandem repeat rDNA units that encode 18S, 5.8S, 2S, and 28S ribosomal RNAs (rRNAs) necessary for ribosome biogenesis. In arthropods like Drosophila, about 60% of the rDNA genes have R1 and/or R2 retrotransposons inserted at specific sites within their 28S regions; these units likely fail to produce functional 28S rRNA. We showed earlier that R2 expression increases upon nucleolar stress caused by the loss of the ribosome assembly factor, Nucleolar Phosphoprotein of 140 kDa (Nopp140). Here we show that R1 expression is selectively induced by heat shock. Actinomycin D, but not α-amanitin, blocked R1 expression in S2 cells upon heat shock, indicating that R1 elements are transcribed by Pol I. A series of RT-PCRs established read-through transcription by Pol I from the 28S gene region into R1. Sequencing the RT-PCR products confirmed the 28S-R1 RNA junction and the expression of R1 elements within nucleolar rDNA rather than R1 elements known to reside in centromeric heterochromatin. Using a genome-wide precision run-on sequencing (PRO-seq) data set available at NCBI-GEO, we show that Pol I activity on R1 elements is negligible under normal non-heat shock conditions but increases upon heat shock. We propose that prior to heat shock Pol I pauses within the 5' end of R1 where we find a consensus "pause button", and that heat shock releases Pol I for read-through transcription farther into R1.

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.

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.

Allosteric Modulation of the Faecalibacterium prausnitzii Hepatitis Delta Virus-like Ribozyme by Glucosamine 6-Phosphate: The Substrate of the Adjacent Gene Product

Self-cleaving ribozymes were discovered 30 years ago and have been found throughout nature, from bacteria to animals, but little is known about their biological functions and regulation, particularly how cofactors and metabolites alter their activity. A hepatitis delta virus-like self-cleaving ribozyme maps upstream of a phosphoglucosamine mutase (glmM) open reading frame in the genome of the human gut bacterium Faecalibacterium prausnitzii. The presence of a ribozyme in the untranslated region of glmM suggests a regulation mechanism of gene expression. In the bacterial hexosamine biosynthesis pathway, the enzyme glmM catalyzes the isomerization of glucosamine 6-phosphate into glucosamine 1-phosphate. In this study, we investigated the effect of these metabolites on the co-transcriptional self-cleavage rate of the ribozyme. Our results suggest that glucosamine 6-phosphate, but not glucosamine 1-phosphate, is an allosteric ligand that increases the self-cleavage rate of drz-Fpra-1, providing the first known example of allosteric modulation of a self-cleaving ribozyme by the substrate of the adjacent gene product. Given that the ribozyme is activated by the glmM substrate, but not the product, this allosteric modulation may represent a potential feed-forward mechanism of gene expression regulation in bacteria.

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.

Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes

Fourteen well-defined ribozyme classes have been identified to date, among which nine are site-specific self-cleaving ribozymes. Very recently, small self-cleaving ribozymes have attracted renewed interest in their structure, biochemistry, and biological function since the discovery, during the last three years, of four novel ribozymes, termed twister, twister sister, pistol, and hatchet. In this review, we mainly address the structure, biochemistry, and catalytic mechanism of the novel ribozymes. They are characterized by distinct active site architectures and divergent, but similar, biochemical properties. The cleavage activities of the ribozymes are highly dependent upon divalent cations, pH, and base-specific mutations, which can cause changes in the nucleotide arrangement and/or electrostatic potential around the cleavage site. It is most likely that a guanine and adenine in close proximity of the cleavage site are involved in general acid-base catalysis. In addition, metal ions appear to play a structural rather than catalytic role although some of their crystal structures have shown a direct metal ion coordination to a non-bridging phosphate oxygen at the cleavage site. Collectively, the structural and biochemical data of the four newest ribozymes could contribute to advance our mechanistic understanding of how self-cleaving ribozymes accomplish their efficient site-specific RNA cleavages.

Integration site selection by retroviruses and transposable elements in eukaryotes

Transposable elements and retroviruses are found in most genomes, can be pathogenic and are widely used as gene-delivery and functional genomics tools. Exploring whether these genetic elements target specific genomic sites for integration and how this preference is achieved is crucial to our understanding of genome evolution, somatic genome plasticity in cancer and ageing, host-parasite interactions and genome engineering applications. High-throughput profiling of integration sites by next-generation sequencing, combined with large-scale genomic data mining and cellular or biochemical approaches, has revealed that the insertions are usually non-random. The DNA sequence, chromatin and nuclear context, and cellular proteins cooperate in guiding integration in eukaryotic genomes, leading to a remarkable diversity of insertion site distribution and evolutionary strategies.

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.

The Wide Distribution and Change of Target Specificity of R2 Non-LTR Retrotransposons in Animals

Transposons, or transposable elements, are the major components of genomes in most eukaryotes. Some groups of transposons have developed target specificity that limits the integration sites to a specific nonessential sequence or a genomic region to avoid gene disruption caused by insertion into an essential gene. R2 is one of the most intensively investigated groups of sequence-specific non-LTR retrotransposons and is inserted at a specific site inside of 28S ribosomal RNA (rRNA) genes. R2 is known to be distributed among at least six animal phyla even though its occurrence is reported to be patchy. Here, in order to obtain a more detailed picture of the distribution of R2, we surveyed R2 using both in silico screening and degenerate PCR, particularly focusing on actinopterygian fish. We found two families of the R2C lineage from vertebrates, although it has previously only been found in platyhelminthes. We also revealed the apparent movement of insertion sites of a lineage of actinopterygian R2, which was likely concurrent with the acquisition of a 28S rRNA-derived sequence in their 3' UTR. Outside of actinopterygian fish, we revealed the maintenance of a single R2 lineage in birds; the co-existence of four lineages of R2 in the leafcutter bee Megachile rotundata; the first examples of R2 in Ctenophora, Mollusca, and Hemichordata; and two families of R2 showing no target specificity. These findings indicate that R2 is relatively stable and universal, while differences in the distribution and maintenance of R2 lineages probably reflect characteristics of some combination of both R2 lineages and host organisms.

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.

Topological constraints of structural elements in regulation of catalytic activity in HDV-like self-cleaving ribozymes

Self-cleaving ribozymes fold into intricate structures, which orient active site groups into catalytically competent conformations. Most ribozyme families have distinct catalytic cores stabilized by tertiary interactions between domains peripheral to those cores. We show that large hepatitis delta virus (HDV)-like ribozymes are activated by peripheral domains that bring two helical segments, P1 and P2, into proximity – a “pinch” that results in rate acceleration by almost three orders of magnitude. Kinetic analysis of ribozymes with systematically altered length and stability of the peripheral domain revealed that about one third of its free energy of formation is used to lower an activation energy barrier, likely related to a rate-limiting conformational change leading to the pre-catalytic state. These findings provide a quantitative view of enzyme regulation by peripheral domains and may shed light on the energetics of allosteric regulation.

Site-specific non-LTR retrotransposons

Although most of non-long terminal repeat (non-LTR) retrotransposons are incorporated in the host genome almost randomly, some non-LTR retrotransposons are incorporated into specific sequences within a target site. On the basis of structural and phylogenetic features, non-LTR retrotransposons are classified into two large groups, restriction enzyme-like endonuclease (RLE)-encoding elements and apurinic/apyrimidinic endonuclease (APE)-encoding elements. All clades of RLE-encoding non-LTR retrotransposons include site-specific elements. However, only two of more than 20 APE-encoding clades, Tx1 and R1, contain site-specific non-LTR elements. Site-specific non-LTR retrotransposons usually target within multi-copy RNA genes, such as rRNA gene (rDNA) clusters, or repetitive genomic sequences, such as telomeric repeats; this behavior may be a symbiotic strategy to reduce the damage to the host genome. Site- and sequence-specificity are variable even among closely related non-LTR elements and appeared to have changed during evolution. In the APE-encoding elements, the primary determinant of the sequence- specific integration is APE itself, which nicks one strand of the target DNA during the initiation of target primed reverse transcription (TPRT). However, other factors, such as interaction between mRNA and the target DNA, and access to the target region in the nuclei also affect the sequence-specificity. In contrast, in the RLE-encoding elements, DNA-binding motifs appear to affect their sequence-specificity, rather than the RLE domain itself. Highly specific integration properties of these site-specific non-LTR elements make them ideal alternative tools for sequence-specific gene delivery, particularly for therapeutic purposes in human diseases.

Integration, Regulation, and Long-Term Stability of R2 Retrotransposons

R2 elements are sequence specific non-LTR retrotransposons that exclusively insert in the 28S rRNA genes of animals. R2s encode an endonuclease that cleaves the insertion site and a reverse transcriptase that uses the cleaved DNA to prime reverse transcription of the R2 transcript, a process termed target primed reverse transcription. Additional unusual properties of the reverse transcriptase as well as DNA and RNA binding domains of the R2 encoded protein have been characterized. R2 expression is through co-transcription with the 28S gene and self-cleavage by a ribozyme encoded at the R2 5' end. Studies in laboratory stocks and natural populations of Drosophila suggest that R2 expression is tied to the distribution of R2-inserted units within the rDNA locus. Most individuals have no R2 expression because only a small fraction of their rRNA genes need to be active, and a contiguous region of the locus free of R2 insertions can be selected for activation. However, if the R2-free region is not large enough to produce sufficient rRNA, flanking units - including those inserted with R2 - must be activated. Finally, R2 copies rapidly turnover within the rDNA locus, yet R2 has been vertically maintained in animal lineages for hundreds of millions of years. The key to this stability is R2's ability to remain dormant in rDNA units outside the transcribed regions for generations until the stochastic nature of the crossovers that drive the concerted evolution of the rDNA locus inevitably reshuffle the inserted and uninserted units, resulting in transcription of the R2-inserted units.

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.

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations

Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.

Plant subviral RNAs as a long noncoding RNA (lncRNA): Analogy with animal lncRNAs in host-virus interactions

Satellite RNAs (satRNAs) and viroids belong to the group called subviral agents and are the smallest pathogens of plants. In general, small satRNAs and viroids are 300-400 nt in size and do not encode any functional proteins; they are thus regarded as so-called long noncoding RNAs (lncRNAs). These lncRNAs are receiving great attention as a new RNA class involved in gene regulation to control important biological processes such as gene transcription and epigenetic regulation. A substantial number of lncRNAs in animal cells have been found to play important roles in the interactions between a virus and its host. We here discuss the pathogenicity of subviral RNAs (especially satRNAs) in plant cells and their functions as lncRNAs associated with viral diseases, using animal lncRNAs as an analogy. Because, unlike animal lncRNAs, plant subviral RNAs can replicate and accumulate at very high levels in infected cells, we here considered the unique possibility that the RNA silencing machinery of plants, an important defense mechanism against virus infection, may have brought about the replication ability of subviral molecules. In addition, we also discuss the possibility that satRNAs may have arisen from plant-virus interactions in virus-infected cells. Understanding the molecular functions of these unique lncRNAs in plants will enable us to reveal the most plausible origins of these subviral RNAs.

A New Class of SINEs with snRNA Gene-Derived Heads

Eukaryotic genomes are colonized by various transposons including short interspersed elements (SINEs). The 5' region (head) of the majority of SINEs is derived from one of the three types of RNA genes--7SL RNA, transfer RNA (tRNA), or 5S ribosomal RNA (rRNA)--and the internal promoter inside the head promotes the transcription of the entire SINEs. Here I report a new group of SINEs whose heads originate from either the U1 or U2 small nuclear RNA gene. These SINEs, named SINEU, are distributed among crocodilians and classified into three families. The structures of the SINEU-1 subfamilies indicate the recurrent addition of a U1- or U2-derived sequence onto the 5' end of SINEU-1 elements. SINEU-1 and SINEU-3 are ancient and shared among alligators, crocodiles, and gharials, while SINEU-2 is absent in the alligator genome. SINEU-2 is the only SINE family that was active after the split of crocodiles and gharials. All SINEU families, especially SINEU-3, are preferentially inserted into a family of Mariner DNA transposon, Mariner-N4_AMi. A group of Tx1 non-long terminal repeat retrotransposons designated Tx1-Mar also show target preference for Mariner-N4_AMi, indicating that SINEU was mobilized by Tx1-Mar.

Dead element replicating: degenerate R2 element replication and rDNA genomic turnover in the Bacillus rossius stick insect (Insecta: Phasmida)

R2 is an extensively investigated non-LTR retrotransposon that specifically inserts into the 28S rRNA gene sequences of a wide range of metazoans, disrupting its functionality. During R2 integration, first strand synthesis can be incomplete so that 5’ end deleted copies are occasionally inserted. While active R2 copies repopulate the locus by retrotransposing, the non-functional truncated elements should frequently be eliminated by molecular drive processes leading to the concerted evolution of the rDNA array(s). Although, multiple R2 lineages have been discovered in the genome of many animals, the rDNA of the stick insect Bacillus rossius exhibits a peculiar situation: it harbors both a canonical, functional R2 element (R2Br^fun) as well as a full-length but degenerate element (R2Br^deg). An intensive sequencing survey in the present study reveals that all truncated variants in stick insects are present in multiple copies suggesting they were duplicated by unequal recombination. Sequencing results also demonstrate that all R2Br^deg copies are full-length, i. e. they have no associated 5' end deletions, and functional assays indicate they have lost the active ribozyme necessary for R2 RNA maturation. Although it cannot be completely ruled out, it seems unlikely that the degenerate elements replicate via reverse transcription, exploiting the R2Br^fun element enzymatic machinery, but rather via genomic amplification of inserted 28S by unequal recombination. That inactive copies (both R2Br^deg or 5'-truncated elements) are not eliminated in a short term in stick insects contrasts with findings for the Drosophila R2, suggesting a widely different management of rDNA loci and a lower efficiency of the molecular drive while achieving the concerted evolution.

Deletion of Drosophila Nopp140 induces subcellular ribosomopathies

The nucleolar and Cajal body phosphoprotein of 140 kDa (Nopp140) is considered a ribosome assembly factor, but its precise functions remain unknown. To approach this problem, we deleted the Nopp140 gene in Drosophila using FLP-FRT recombination. Genomic PCR, reverse transcriptase-PCR (RT-PCR), and immunofluorescence microscopy confirmed the loss of Nopp140, its messenger RNA (mRNA), and protein products from all tissues examined. Nopp140-/- larvae arrested in the second instar stage and most died within 8 days. While nucleoli appeared intact in Nopp140-/- cells, the C/D small nucleolar ribonucleoprotein (snoRNP) methyltransferase, fibrillarin, redistributed to the nucleoplasm in variable amounts depending on the cell type; RT-PCRs showed that 2'-O-methylation of ribosomal RNA (rRNA) in Nopp140-/- cells was reduced at select sites within both the 18S and 28S rRNAs. Ultrastructural analysis showed that Nopp140-/- cells were deficient in cytoplasmic ribosomes, but instead contained abnormal electron-dense cytoplasmic granules. Immunoblot analysis showed a loss of RpL34, and metabolic labeling showed a significant drop in protein translation, supporting the loss of functional ribosomes. Northern blots showed that pre-RNA cleavage pathways were generally unaffected by the loss of Nopp140, but that R2 retrotransposons that naturally reside within the 28S region of normally silent heterochromatic Drosophila ribosomal DNA (rDNA) genes were selectively expressed in Nopp140-/- larvae. Unlike copia elements and the related R1 retrotransposon, R2 expression appeared to be preferentially dependent on the loss of Nopp140 and not on environmental stresses. We believe the phenotypes described here define novel intracellular ribosomopathies resulting from the loss of Nopp140.

Pr77 and L1TcRz: A dual system within the 5'-end of L1Tc retrotransposon, internal promoter and HDV-like ribozyme

The sequence corresponding to the first 77 nucleotides of the L1Tc and NARTc non-LTR retrotransposons from Trypanosoma cruzi is an internal promoter (Pr77) that generates abundant, although poorly translatable, un-spliced transcripts. It has been recently described that L1TcRz, an HDV-like ribozyme, resides within the 5'-end of the RNA from the L1Tc and NARTc retrotransposons. Remarkably, the same first 77 nucleotides of L1Tc/NARTc elements comprise both the Pr77 internal promoter and the HDV-like L1TcRz. The L1TcRz cleaves on the 5'-side of the +1 nucleotide of the L1Tc element insuring that the promoter and the ribozyme functions travel with the transposon during retrotransposition. The ribozyme activity would prevent the mobilization of upstream sequences and insure the individuality of the L1Tc/NARTc copies transcribed from associated tandems. The Pr77/L1TcRz sequence is also found in other trypanosomatid's non-LTR retrotransposons and degenerated retroposons. The possible conservation of the ribozyme activity in a widely degenerated retrotransposon, as the Leishmania SIDERs, could indicate that the presence of this element and the catalytic activity could play some favorable genetic regulation. The functional implications of the Pr77/L1TcRz dual system in the regulation of the L1Tc/NARTc retrotransposons and in the gene expression of trypanosomatids are also discussed in this paper.

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 wide expansion of hepatitis delta virus-like ribozymes throughout trypanosomatid genomes is linked to the spreading of L1Tc/ingi clade mobile elements

Background

Hepatitis Delta Virus (HDV)-like ribozymes have recently been found in many mobile elements in which they take part in a mechanism that releases intermediate RNAs from cellular co-transcripts. L1Tc in Trypanosoma cruzi is one of the elements in which such a ribozyme is located. It lies in the so-called Pr77-hallmark, a conserved region shared by retrotransposons belonging to the trypanosomatid L1Tc/ingi clade. The wide distribution of the Pr77-hallmark detected in trypanosomatid retrotransposons renders the potential catalytic activity of these elements worthy of study: their distribution might contribute to host genetic regulation at the mRNA level. Indeed, in Leishmania spp, the pervasive presence of these HDV-like ribozyme-containing mobile elements in certain 3′-untranslated regions of protein-coding genes has been linked to mRNA downregulation.

Results

Intensive screening of publicly available trypanosomatid genomes, combined with manual folding analyses, allowed the isolation of putatively Pr77-hallmarks with HDV-like ribozyme activity. This work describes the conservation of an HDV-like ribozyme structure in the Pr77 sequence of retrotransposons in a wide range of trypanosomatids, the catalytic function of which is maintained in the majority.

These results are consistent with the previously suggested common phylogenetic origin of the elements that belong to this clade, although in some cases loss of functionality appears to have occurred and/or perhaps molecular domestication by the host.

Conclusions

These HDV-like ribozymes are widely distributed within retrotransposons across trypanosomatid genomes. This type of ribozyme was once thought to be rare in nature, but in fact it would seem to be abundant in trypanosomatid transcripts. It can even form part of the pool of mRNA 3′-untranslated regions, particularly in Leishmania spp. Its putative regulatory role in host genetic expression is discussed.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-340) contains supplementary material, which is available to authorized users.