Identification of an hepatitis delta virus-like ribozyme at the mRNA 5'-end of the L1Tc retrotransposon from Trypanosoma cruzi

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.

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.

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.

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.

Variability of the Pr77 sequence of L1Tc retrotransposon among six T. cruzi strains belonging to different discrete typing units (DTUs)

All trypanosomatid genomes are colonized by non-LTR retrotransposons which exhibit a highly conserved 77-nt sequence at their 5' ends, known as the Pr77-hallmark (Pr77). The wide distribution of Pr77 is expected to be related to the gene regulation processes in these organisms as it has promoter and HDV-like ribozyme activities at the DNA and RNA levels, respectively. The identification of Pr77 hallmark-bearing retrotransposons and the study of the associations of mobile elements with relevant genes have been analyzed in the genomes of six strains of Trypanosoma cruzi belonging to different discrete typing units (DTUs) and with different geographical origins and host/vectors. The genomes have been sequenced, assembled and annotated. BUSCO analyses indicated a good quality for the assemblies that were used in comparative analyses. The results show differences among the six genomes in the copy number of genes related to virulence processes, the abundance of retrotransposons bearing the Pr77 sequence and the presence of the Pr77 hallmarks not associated with retroelements. The analyses also show frequent associations of Pr77-bearing retrotransposons and single Pr77 hallmarks with genes coding for trans-sialidases, RHS, MASP or hypothetical proteins, showing variable proportion depending on the type of retroelement, gene class and parasite strain. These differences in the genomic distribution of active retroelements and other Pr77-containing elements have shaped the genome architecture of these six strains and might be contributing to the phenotypic variability existing among them.

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.

TBP and SNAP50 transcription factors bind specifically to the Pr77 promoter sequence from trypanosomatid non-LTR retrotransposons

Background

Trypanosomatid genomes are colonized by active and inactive mobile DNA elements, such as LINE, SINE-like, SIDER and DIRE retrotransposons. These elements all share a 77-nucleotide-long sequence at their 5′ ends, known as Pr77, which activates transcription, thereby generating abundant unspliced and translatable transcripts. However, transcription factors that mediates this process have still not been reported.

Methods

TATA-binding protein (TBP) and small nuclear RNA-activating protein 50 kDa (SNAP50) recombinant proteins and specific antibodies raised against them were generated. Protein capture assay, electrophoretic mobility-shift assays (EMSA) and EMSA competition assays carried out using these proteins and nuclear proteins of the parasite together to specific DNA sequences used as probes allowed detecting direct interaction of these transcription factors to Pr77 sequence.

Results

This study identified TBP and SNAP50 as part of the DNA-protein complex formed by the Pr77 promoter sequence and nuclear proteins of Trypanosoma cruzi. TBP establishes direct and specific contact with the Pr77 sequence, where the DPE and DPE downstream regions are docking sites with preferential binding. TBP binds cooperatively (Hill coefficient = 1.67) to Pr77 and to both strands of the Pr77 sequence, while the conformation of this highly structured sequence is not involved in TBP binding. Direct binding of SNAP50 to the Pr77 sequence is weak and may be mediated by protein–protein interactions through other trypanosomatid nuclear proteins.

Conclusions

Identification of the transcription factors that mediate Pr77 transcription may help to elucidate how these retrotransposons are mobilized within the trypanosomatid genomes and their roles in gene regulation processes in this human parasite.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04803-5.

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.

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.

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.

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.

Biology of Trypanosoma cruzi Retrotransposons: From an Enzymatic to a Structural Point of View

INTRODUCTION

An important portion of the Trypanosoma cruzi genome is composed of mobile genetic elements, which are interspersed with genes on all chromosomes. The L1Tc non-LTR retrotransposon and its truncated version NARTc are the most highly represented and best studied of these elements. L1Tc is actively transcribed in all three forms of the Trypanosoma parasite and encodes the proteins that enable it to autonomously mobilize. This mini review discusses the enzymatic properties of L1Tc that enable its mobilization and possibly the mobilization of other non-autonomous retrotransposons in Trypanosoma. We also briefly review the Hepatitis Delta Virus-like autocatalytic and 2A self-cleaving viral-like sequences contained in L1Tc that regulate post-transcriptional properties such as relative protein abundance and mRNA stability. Special emphasis is placed on the Pr77 dual system, which is based on the RNA pol II-dependent internal promoter of L1Tc and NARTc and the HDV-like ribozyme activity encoded by the first 77 nucleotides of the element's DNA and RNA. The high degree of conservation of the Pr77 sequence, referred to as the "Pr77-hallmark", among different trypanosomatid retroelements suggests that these mobile elements are responsible for the distribution of regulatory sequences within the genome they inhabit.

CONCLUSION

We also discuss how the involvement of L1Tc and NARTc in the gene regulatory processes of these parasites could justify their domestication and long-term coexistence in these ancient organisms.

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.

RNA secondary structure and nucleotide composition of the conserved hallmark sequence of Leishmania SIDER2 retroposons are essential for endonucleolytic cleavage and mRNA degradation

We have reported previously that Short Interspersed Degenerate Retroposons of the SIDER2 subfamily, largely located within 3'UTRs of Leishmania transcripts, promote rapid turnover of mRNAs through endonucleolytic cleavage within the highly conserved second tandem 79-nt hallmark sequence (79-nt SII). Here, we used site-directed mutagenesis and in silico RNA structural studies to delineate the cis-acting requirements within 79-nt SII for cleavage and mRNA degradation. The putative cleavage site(s) and other nucleotides predicted to alter the RNA secondary structure of 79-nt SII were either deleted or mutated and their effect on mRNA turnover was monitored using a gene reporter system. We found that short deletions of 8-nt spanning the two predicted cleavage sites block degradation of SIDER2-containing transcripts, leading to mRNA accumulation. Furthermore, single or double substitutions of the dinucleotides targeted for cleavage as well as mutations altering the predicted RNA secondary structure encompassing both cleavage sites also prevent mRNA degradation, confirming that these dinucleotides are the bona fide cleavage sites. In line with these results, we show that stage-regulated SIDER2 inactivation correlates with the absence of endonucleolytic cleavage. Overall, these data demonstrate that both cleavage sites within the conserved 79-nt SII as well as RNA folding in this region are essential for SIDER2-mediated mRNA decay, and further support that SIDER2-harboring transcripts are targeted for degradation by endonucleolytic cleavage.

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.

Genomic cartography and proposal of nomenclature for the repeated, interspersed elements of the Leishmania major SIDER2 family and identification of SIDER2-containing transcripts

The genomes of most eukaryotic organisms contain a large number of transposable elements that are able to move from one genomic site to another either by transferring of DNA mobile elements (transposons) or transpose via reverse transcription of an RNA intermediate (retroposons). An exception to this rule is found in protists of the subgenus Leishmania, in which active retroposons degenerated after a flourishing era, leaving only retroposon remains; these have been classified into two families: SIDER1 and SIDER2. In this work, we have re-examined the elements belonging to the family SIDER2 present in the genome of Leishmania major with the aim of providing a nomenclature that will facilitate a future reference to particular elements. According to sequence conservation, the 1100 SIDER2 elements have been grouped into subfamilies, and the inferred taxonomic relationships have also been incorporated into the nomenclature. Additionally, we are providing detailed data regarding the genomic distribution of these elements and their association with specific transcripts, based on the recently established transcriptome for L. major. Thus, the presented data can help to study and better understand the roles played by these degenerated retroposons in both regulation of gene expression and genome plasticity.

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.

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.

The Trypanosomatid Pr77-hallmark contains a downstream core promoter element essential for transcription activity of the Trypanosoma cruzi L1Tc retrotransposon

Background

Trypanosomatid genomes are highly colonized by non-LTR retroelements that make up to 5 % of the nuclear genome. These elements are mainly accumulated in the strand switch regions (SSRs) where polycistronic transcription is initiated and have a 77 nt-long sequence - Pr77 - at their 5′ ends. L1Tc is the best represented retrotransposon in the Trypanosoma cruzi genome and is a potentially functional autonomous element that encodes its own retrotransposition machinery. The Pr77 of the T. cruzi L1Tc element activates gene transcription via RNA polymerase II, generating abundant, unspliced transcripts which are translated.

Results

The present manuscript describes the identification of a downstream core promoter element (DPE) in the L1Tc Pr77 sequence. Just four nucleotides long (CGTG), it covers in Pr77 positions +25 to +28 of the described L1Tc transcription start site. The Pr77-DPE motif is conserved in terms of sequence composition and position in the Pr77 of most trypanosomatid non-LTR retrotransposons, independent of the coding or non-coding capacity of these retroelements. Transcription assays in T. cruzi stable transfectants with vector containing point mutations at 17 locations of the Pr77 nucleotide sequence evidence that the DPE motif is essential for the promoter function of Pr77. Furthermore, the obtained data show that other nucleotides also contributed to the promoter function of Pr77. In addition, the presented results indicate that parasite nuclear proteins specifically bind to different regions of the Pr77 sequence although the strongest binding is to the DPE motif. Moreover, it is shown that the DPE sense single-stranded sequence is being required in DNA-protein recognition of nuclear factors.

Conclusions

The Pr77 sequence present in most of non-LTR retrotransposons of trypanosomatids contains a downstream core promoter element (DPE) which is conserved in terms of nucleotide composition and location. The Pr77-DPE motif is essential for the transcriptional activity of Pr77 although other nucleotides are also involved. DPE has a high affinity binding for nuclear proteins in T. cruzi. The wide retroelement-mediated distribution of Pr77 suggests that it may represent an important tool for regulating gene expression in trypanosomatids.

Electronic supplementary material

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

Mobile Element Evolution Playing Jigsaw - SINEs in Gastropod and Bivalve Mollusks

SINEs (Short INterspersed Elements) are widely distributed among eukaryotes. Some SINE families are organized in superfamilies characterized by a shared central domain. These central domains are conserved across species, classes, and even phyla. Here we report the identification of two novel such superfamilies in the genomes of gastropod and bivalve mollusks. The central conserved domain of the first superfamily is present in SINEs in Caenogastropoda and Vetigastropoda as well as in all four subclasses of Bivalvia. We designated the domain MESC (Romanian for MElc-snail and SCoica-mussel) because it appears to be restricted to snails and mussels. The second superfamily is restricted to Caenogastropoda. Its central conserved domain-Snail-is related to the Nin-DC domain. Furthermore, we provide evidence that a 40-bp subdomain of the SINE V-domain is conserved in SINEs in mollusks and arthropods. It is predicted to form a stable stem-loop structure that is preserved in the context of the overall SINE RNA secondary structure in invertebrates. Our analysis also recovered short retrotransposons with a Long INterspersed Element (LINE)-derived 5' end. These share the body and/or the tail with transfer RNA (tRNA)-derived SINEs within and across species. Finally, we identified CORE SINEs in gastropods and bivalves-extending the distribution range of this superfamily.

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.

New classes of self-cleaving ribozymes revealed by comparative genomics analysis

Enzymes made of RNA catalyze reactions that are essential for protein synthesis and RNA processing. However, such natural ribozymes are exceedingly rare, as evident by the fact that the discovery rate for new classes has dropped to one per decade from about one per year during the 1980s. Indeed, only 11 distinct ribozyme classes have been experimentally validated to date. Recently, we recognized that self-cleaving ribozymes frequently associate with certain types of genes from bacteria. Herein this synteny was exploited to identify divergent architectures for two previously known ribozyme classes and to discover additional noncoding RNA motifs that are self-cleaving RNA candidates. Three new self-cleaving classes, named twister sister, pistol and hatchet, have been identified from this collection, suggesting that even more ribozymes remain hidden in modern cells.

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.

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.

Identification of minimal HDV-like ribozymes with unique divalent metal ion dependence in the human microbiome

HDV-like self-cleaving ribozymes have been found in a wide variety of organisms, implicated in diverse biological processes, and their activity typically shows a strong divalent metal dependence, but little metal specificity. Recent studies suggested that very short variants of these ribozymes exist in nature, but their distribution and biochemical properties have not been established. To map out the distribution of small HDV-like ribozymes, the drz-Spur-3 sequence was minimized to yield a core construct for structure-based bioinformatic searches. These searches revealed several microbial ribozymes, particularly in the human microbiome. Kinetic profile of the smallest ribozyme revealed two distinct metal binding sites, only one of which promotes fast catalysis. Furthermore, this ribozyme showed markedly reduced activity in Ca(2+), even in the presence of physiological Mg(2+) concentrations. Our study substantially expands the number of microbial HDV-like ribozymes and provides an example of cleavage regulation by divalent metals.

A widespread self-cleaving ribozyme class is revealed by bioinformatics

Ribozymes are noncoding RNAs that promote chemical transformations with rate enhancements approaching those of protein enzymes. Although ribozymes are likely to have been abundant during the RNA world era, only ten classes are known to exist among contemporary organisms. We report the discovery and analysis of an additional self-cleaving ribozyme class, called twister, which is present in many species of bacteria and eukarya. Nearly 2,700 twister ribozymes were identified that conform to a secondary structure consensus that is small yet complex, with three stems conjoined by internal and terminal loops. Two pseudoknots provide tertiary structure contacts that are critical for catalytic activity. The twister ribozyme motif provides another example of a natural RNA catalyst and calls attention to the potentially varied biological roles of this and other classes of widely distributed self-cleaving RNAs.

Evolution of the R2 retrotransposon ribozyme and its self-cleavage site

R2 is a non-long terminal repeat retrotransposon that inserts site-specifically in the tandem 28S rRNA genes of many animals. Previously, R2 RNA from various species of Drosophila was shown to self-cleave from the 28S rRNA/R2 co-transcript by a hepatitis D virus (HDV)-like ribozyme encoded at its 5' end. RNA cleavage was at the precise 5' junction of the element with the 28S gene. Here we report that RNAs encompassing the 5' ends of R2 elements from throughout its species range fold into HDV-like ribozymes. In vitro assays of RNA self-cleavage conducted in many R2 lineages confirmed activity. For many R2s, RNA self-cleavage was not at the 5' end of the element but at 28S rRNA sequences up to 36 nucleotides upstream of the junction. The location of cleavage correlated well with the types of endogenous R2 5' junctions from different species. R2 5' junctions were uniform for most R2s in which RNA cleavage was upstream in the rRNA sequences. The 28S sequences remaining on the first DNA strand synthesized during retrotransposition are postulated to anneal to the target site and uniformly prime second strand DNA synthesis. In species where RNA cleavage occurred at the R2 5' end, the 5' junctions were variable. This junction variation is postulated to result from the priming of second strand DNA synthesis by chance microhomologies between the target site and the first DNA strand. Finally, features of R2 ribozyme evolution, especially changes in cleavage site and convergence on the same active site sequences, are discussed.

Characterization of the trans Watson-Crick GU base pair located in the catalytic core of the antigenomic HDV ribozyme

The HDV ribozyme’s folding pathway is, by far, the most complex folding pathway elucidated to date for a small ribozyme. It includes 6 different steps that have been shown to occur before the chemical cleavage. It is likely that other steps remain to be discovered. One of the most critical of these unknown steps is the formation of the trans Watson-Crick GU base pair within loop III. The U₂₃ and G₂₈ nucleotides that form this base pair are perfectly conserved in all natural variants of the HDV ribozyme, and therefore are considered as being part of the signature of HDV-like ribozymes. Both the formation and the transformation of this base pair have been studied mainly by crystal structure and by molecular dynamic simulations. In order to obtain physical support for the formation of this base pair in solution, a set of experiments, including direct mutagenesis, the site-specific substitution of chemical groups, kinetic studies, chemical probing and magnesium-induced cleavage, were performed with the specific goal of characterizing this trans Watson-Crick GU base pair in an antigenomic HDV ribozyme. Both U₂₃ and G₂₈ can be substituted for nucleotides that likely preserve some of the H-bond interactions present before and after the cleavage step. The formation of the more stable trans Watson-Crick base pair is shown to be a post-cleavage event, while a possibly weaker trans Watson-Crick/Hoogsteen interaction seems to form before the cleavage step. The formation of this unusually stable post-cleavage base pair may act as a driving force on the chemical cleavage by favouring the formation of a more stable ground state of the product-ribozyme complex. To our knowledge, this represents the first demonstration of a potential stabilising role of a post-cleavage conformational switch event in a ribozyme-catalyzed reaction.

R2 and R2/R1 hybrid non-autonomous retrotransposons derived by internal deletions of full-length elements

BACKGROUND

R2 is a non-long terminal repeat (non-LTR) retrotransposable element that inserts site specifically into the 28S genes of the ribosomal (r)RNA gene loci. Encoded at the 5' end is a ribozyme that generates the precise 5' end by self-cleavage of a 28S gene cotranscript. Sequences at the 3' end are necessary for the R2 protein to bind RNA and initiate the target primed reverse transcription (TPRT) reaction. These minimal RNA requirements suggested that if recombination/DNA repair conjoined the 5' and 3' ends of R2, the result would be a non-autonomous element that could survive as long as autonomous R2 elements supplied the TPRT activity.

RESULTS

A PCR-based survey of 39 Drosophila species aided by genomic sequences from 12 of these species revealed two types of non-autonomous elements. We call these elements SIDEs (for 'Short Internally Deleted Elements'). The first consisted of a 5' ribozyme and a 3' end of an R2 element as predicted. Variation at the 5' junctions of the R2 SIDE copies was typical for R2 insertions suggesting their propagation by TPRT. The second class of SIDE contained sequences from R1 elements, another non-LTR retrotransposon that inserts into rRNA gene loci. These insertions had an R2 ribozyme immediately upstream of R1 3' end sequences. These hybrid SIDEs were inserted at the R1 site with 14 bp target site duplications typical of R1 insertions suggesting they used the R1 machinery for retrotransposition. Finally, the survey revealed examples of U12 small nuclear (sn)RNA and tRNA sequences at the 5' end of R2 elements suggesting the R2 reverse transcriptase can template jump from the R2 transcript to a second RNA during TPRT.

CONCLUSIONS

The R2 SIDE and R2/R1 hybrid SIDEs are rare examples of non-autonomous retrotransposons in the Drosophila genome. Associated non-autonomous elements and in vivo template jumps are two additional characteristics R2 shares with other non-LTR retrotransposons such as mammalian L1s. Analysis of the hybrid SIDEs provides supporting evidence that R1 elements, like R2 elements, recognize their 3' untranslated region (UTR) sequences and, thus, belong to the stringent class of non-LTR elements.

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

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