A genomewide search for ribozymes reveals an HDV-like sequence in the human CPEB3 gene

The distinct spatiotemporal evolutionary landscape of HBV and HDV largely determines the unique epidemic features of HDV globally

Chronic infection of hepatitis B virus (HBV) and hepatitis D virus (HDV) causes the most severe form of viral hepatitis. Due to the dependence on HBV, HDV was deemed to co-evolve and co-migrate with HBV. However, we previously found that the naturally occurred HDV/HBV combinations do not always reflect the most efficient virological adaptation (Wang et al., 2021). Moreover, regions with heavy HBV burden do not always correlate with high HDV prevalence (e.g., East Asia), and vice versa (e.g., Central Asia). Herein, we systematically elucidated the spatiotemporal evolutionary landscape of HDV to understand the unique epidemic features of HDV. We found that the MRCA of HDV was from South America around the late 13th century, was globally dispersed mainly via Central Asia, and evolved into eight genotypes from the 19th to 20th century. In contrast, the MRCA of HBV was from Europe ∼23.7 thousand years ago (Kya), globally dispersed mainly via Africa and East Asia, and evolved into eight genotypes ∼1100 years ago. When HDV stepped in, all present-day HBV genotypes had already formed and its global genotypic distribution had stayed stable geographically. Nevertheless, regionalized HDV adapted to local HBV genotypes and human lineages, contributing to the global geographical separation of HDV genotypes. Additionally, a sharp increase in HDV infections was observed after the 20th century. In conclusion, HDV exhibited a distinct spatiotemporal distribution path compared with HBV. This unique evolutionary relationship largely fostered the unique epidemic features we observe nowadays. Moreover, HDV infections may continue to ramp up globally, thus more efforts are urgently needed to combat this disease.

Direct testing of natural twister ribozymes from over a thousand organisms reveals a broad tolerance for structural imperfections

Twister ribozymes are an extensively studied class of nucleolytic RNAs. Thousands of natural twisters have been proposed using sequence homology and structural descriptors. Yet, most of these candidates have not been validated experimentally. To address this gap, we developed CHiTA (Cleavage High-Throughput Assay), a high-throughput pipeline utilizing massively parallel oligonucleotide synthesis and next-generation sequencing to test putative ribozymes en masse in a scarless fashion. As proof of principle, we applied CHiTA to a small set of known active and mutant ribozymes. We then used CHiTA to test two large sets of naturally occurring twister ribozymes: over 1, 600 previously reported putative twisters and ∼1, 000 new candidate twisters. The new candidates were identified computationally in ∼1, 000 organisms, representing a massive increase in the number of ribozyme-harboring organisms. Approximately 94% of the twisters we tested were active and cleaved site-specifically. Analysis of their structural features revealed that many substitutions and helical imperfections can be tolerated. We repeated our computational search with structural descriptors updated from this analysis, whereupon we identified and confirmed the first intrinsically active twister ribozyme in mammals. CHiTA broadly expands the number of active twister ribozymes found in nature and provides a powerful method for functional analyses of other RNAs.

GRAPHICAL ABSTRACT

New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn2+ and Mg2

CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.

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.

Flanking Sequence Cotranscriptionally Regulates Twister Ribozyme Activity

Small nucleolytic ribozymes are RNAs that cleave their own phosphodiester backbone. While proteinaceous enzymes are regulated by a variety of known mechanisms, methods of regulation for ribozymes remain unclear. Twister is one ribozyme class for which many structural and catalytic properties have been elucidated. However, few studies have analyzed the activity of twister ribozymes in the context of a native flanking sequence, even though ribozymes as transcribed in nature do not exist in isolation. Interactions between the ribozyme and its neighboring sequences can induce conformational changes that inhibit self-cleavage, providing a regulatory mechanism that could naturally determine ribozyme activity in vivo and in synthetic applications. To date, eight twister ribozymes have been identified within the staple crop rice (Oryza sativa). Herein, we select several twister ribozymes from rice and show that they are differentially regulated by their flanking sequence using published RNA-seq data sets, structure probing, and cotranscriptional cleavage assays. We found that the Osa 1-2 ribozyme does not interact with its flanking sequences. However, sequences flanking the Osa 1-3 and Osa 1-8 ribozymes form inactive conformations, referred to here as "ribozymogens", that attenuate ribozyme self-cleavage activity. For the Osa 1-3 ribozyme, we show that activity can be rescued upon addition of a complementary antisense oligonucleotide, suggesting ribozymogens can be controlled via external signals. In all, our data provide a plausible mechanism wherein flanking sequence differentially regulates ribozyme activity in vivo. More broadly, the ability to regulate ribozyme behavior locally has potential applications in control of gene expression and synthetic biology.

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

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

Assessment of three-dimensional RNA structure prediction in CASP15

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

RNA target highlights in CASP15: Evaluation of predicted models by structure providers

The first RNA category of the Critical Assessment of Techniques for Structure Prediction competition was only made possible because of the scientists who provided experimental structures to challenge the predictors. In this article, these scientists offer a unique and valuable analysis of both the successes and areas for improvement in the predicted models. All 10 RNA-only targets yielded predictions topologically similar to experimentally determined structures. For one target, experimentalists were able to phase their x-ray diffraction data by molecular replacement, showing a potential application of structure predictions for RNA structural biologists. Recommended areas for improvement include: enhancing the accuracy in local interaction predictions and increased consideration of the experimental conditions such as multimerization, structure determination method, and time along folding pathways. The prediction of RNA-protein complexes remains the most significant challenge. Finally, given the intrinsic flexibility of many RNAs, we propose the consideration of ensemble models.

Assessment of three-dimensional RNA structure prediction in CASP15

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and X-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as non-canonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Discovery of the major 15-30 nt mammalian small RNAs, their biogenesis and function

Small RNAs (sRNAs) within 15-30 nt such as miRNA, tsRNA, srRNA with 3'-OH have been identified. However, whether these sRNAs are the major 15-30 nt sRNAs is still unknown. Here we show about 90% mammalian sRNAs within 15-30 nt end with 2',3'-cyclic phosphate (3'-cP). TANT-seq was developed to simultaneously profile sRNAs with 3'-cP (sRNA-cPs) and sRNA-OHs, and huge amount of sRNA-cPs were detected. Surprisingly, sRNA-cPs and sRNA-OHs usually have distinct sequences. The data from TANT-seq were validated by a novel method termed TE-qPCR, and Northern blot. Furthermore, we found that Angiogenin and RNase 4 contribute to the biogenesis of sRNA-cPs. Moreover, much more sRNA-cPs than sRNA-OHs bind to Ago2, and can regulate gene expression. Particularly, snR-2-cP regulates Bcl2 by targeting to its 3'UTR dependent on Ago2, and subsequently regulates apoptosis. In addition, sRNA-cPs can guide the cleavage of target RNAs in Ago2 complex as miRNAs without the requirement of 3'-cP. Our discovery greatly expands the repertoire of mammalian sRNAs, and provides strategies and powerful tools towards further investigation of sRNA-cPs.

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.

Mining metatranscriptomes reveals a vast world of viroid-like circular RNAs

Viroids and viroid-like covalently closed circular (ccc) RNAs are minimal replicators that typically encode no proteins and hijack cellular enzymes for replication. The extent and diversity of viroid-like agents are poorly understood. We developed a computational pipeline to identify viroid-like cccRNAs and applied it to 5,131 metatranscriptomes and 1,344 plant transcriptomes. The search yielded 11,378 viroid-like cccRNAs spanning 4,409 species-level clusters, a 5-fold increase compared to the previously identified viroid-like elements. Within this diverse collection, we discovered numerous putative viroids, satellite RNAs, retrozymes, and ribozy-like viruses. Diverse ribozyme combinations and unusual ribozymes within the cccRNAs were identified. Self-cleaving ribozymes were identified in ambiviruses, some mito-like viruses and capsid-encoding satellite virus-like cccRNAs. The broad presence of viroid-like cccRNAs in diverse transcriptomes and ecosystems implies that their host range is far broader than currently known, and matches to CRISPR spacers suggest that some cccRNAs replicate in prokaryotes.

RNA sequence to structure analysis from comprehensive pairwise mutagenesis of multiple self-cleaving ribozymes

Self-cleaving ribozymes are RNA molecules that catalyze the cleavage of their own phosphodiester backbones. These ribozymes are found in all domains of life and are also a tool for biotechnical and synthetic biology applications. Self-cleaving ribozymes are also an important model of sequence-to-function relationships for RNA because their small size simplifies synthesis of genetic variants and self-cleaving activity is an accessible readout of the functional consequence of the mutation. Here, we used a high-throughput experimental approach to determine the relative activity for every possible single and double mutant of five self-cleaving ribozymes. From this data, we comprehensively identified non-additive effects between pairs of mutations (epistasis) for all five ribozymes. We analyzed how changes in activity and trends in epistasis map to the ribozyme structures. The variety of structures studied provided opportunities to observe several examples of common structural elements, and the data was collected under identical experimental conditions to enable direct comparison. Heatmap-based visualization of the data revealed patterns indicating structural features of the ribozymes including paired regions, unpaired loops, non-canonical structures, and tertiary structural contacts. The data also revealed signatures of functionally critical nucleotides involved in catalysis. The results demonstrate that the data sets provide structural information similar to chemical or enzymatic probing experiments, but with additional quantitative functional information. The large-scale data sets can be used for models predicting structure and function and for efforts to engineer self-cleaving ribozymes.

Ribocentre: a database of ribozymes

Ribozymes are excellent systems in which to study 'sequence - structure - function' relationships in RNA molecules. Understanding these relationships may greatly help structural modeling and design of functional RNA structures and some functional structural modules could be repurposed in molecular design. At present, there is no comprehensive database summarising all the natural ribozyme families. We have therefore created Ribocentre, a database that collects together sequence, structure and mechanistic data on 21 ribozyme families. This includes available information on timelines, sequence families, secondary and tertiary structures, catalytic mechanisms, applications of the ribozymes together with key publications. The database is publicly available at https://www.ribocentre.org.

Predicting higher-order mutational effects in an RNA enzyme by machine learning of high-throughput experimental data

Ribozymes are RNA molecules that catalyze biochemical reactions. Self-cleaving ribozymes are a common naturally occurring class of ribozymes that catalyze site-specific cleavage of their own phosphodiester backbone. In addition to their natural functions, self-cleaving ribozymes have been used to engineer control of gene expression because they can be designed to alter RNA processing and stability. However, the rational design of ribozyme activity remains challenging, and many ribozyme-based systems are engineered or improved by random mutagenesis and selection (in vitro evolution). Improving a ribozyme-based system often requires several mutations to achieve the desired function, but extensive pairwise and higher-order epistasis prevent a simple prediction of the effect of multiple mutations that is needed for rational design. Recently, high-throughput sequencing-based approaches have produced data sets on the effects of numerous mutations in different ribozymes (RNA fitness landscapes). Here we used such high-throughput experimental data from variants of the CPEB3 self-cleaving ribozyme to train a predictive model through machine learning approaches. We trained models using either a random forest or long short-term memory (LSTM) recurrent neural network approach. We found that models trained on a comprehensive set of pairwise mutant data could predict active sequences at higher mutational distances, but the correlation between predicted and experimentally observed self-cleavage activity decreased with increasing mutational distance. Adding sequences with increasingly higher numbers of mutations to the training data improved the correlation at increasing mutational distances. Systematically reducing the size of the training data set suggests that a wide distribution of ribozyme activity may be the key to accurate predictions. Because the model predictions are based only on sequence and activity data, the results demonstrate that this machine learning approach allows readily obtainable experimental data to be used for RNA design efforts even for RNA molecules with unknown structures. The accurate prediction of RNA functions will enable a more comprehensive understanding of RNA fitness landscapes for studying evolution and for guiding RNA-based engineering efforts.

Conformational dynamics in the disordered region of human CPEB3 linked to memory consolidation

Background

Current understanding of the molecular basis of memory consolidation points to an important function of amyloid formation by neuronal-specific isoforms of the cytoplasmic polyadenylation element binding (CPEB) protein family. In particular, CPEB is thought to promote memory persistence through formation of self-sustaining prion-like amyloid assemblies at synapses, mediated by its intrinsically disordered region (IDR) and leading to permanent physical alterations at the basis of memory persistence. Although the molecular mechanisms by which amyloid formation takes place in CPEB have been described in invertebrates, the way amyloid formation occurs in the human homolog CPEB3 (hCPEB3) remains unclear. Here, we characterize by NMR spectroscopy the atomic level conformation and ps-ms dynamics of the 426-residue IDR of hCPEB3, which has been associated with episodic memory in humans.

Results

We show that the 426-residue N-terminal region of hCPEB3 is a dynamic, intrinsically disordered region (IDR) which lacks stable folded structures. The first 29 residues, M₁QDDLLMDKSKTQPQPQQQQRQQQQPQP₂₉, adopt a helical + disordered motif, and residues 86–93: P₈₃QQPPPP₉₃, and 166–175: P₁₆₆PPPAPAPQP₁₇₅ form polyproline II (PPII) helices. The (VG)₅ repeat motif is completely disordered, and residues 200–250 adopt three partially populated α-helices. Residues 345–355, which comprise the nuclear localization signal (NLS), form a modestly populated α-helix which may mediate STAT5B binding. These findings allow us to suggest a model for nascent hCPEB3 structural transitions at single residue resolution, advancing that amyloid breaker residues, like proline, are a key difference between functional versus pathological amyloids.

Conclusion

Our NMR spectroscopic analysis of hCPEB3 provides insights into the first structural transitions involved in protein–protein and protein-mRNA interactions. The atomic level understanding of these structural transitions involved in hCPEB3 aggregation is a key first step toward understanding memory persistence in humans, as well as sequence features that differentiate beneficial amyloids from pathological ones.

Areas

Biophysics, Structural Biology, Biochemistry & Neurosciences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01310-6.

Application of RtcB ligase to monitor self-cleaving ribozyme activity by RNA-seq

Self-cleaving ribozymes are catalytic RNAs and can be found in all domains of life. They catalyze a site-specific cleavage that results in a 5' fragment with a 2',3' cyclic phosphate (2',3' cP) and a 3' fragment with a 5' hydroxyl (5' OH) end. Recently, several strategies to enrich self-cleaving ribozymes by targeted biochemical methods have been introduced by us and others. Here, we develop an alternative strategy in which 5' OH RNAs are specifically ligated by RtcB ligase, which first guanylates the 3' phosphate of the adapter and then ligates it directly to RNAs with 5' OH ends. Our results demonstrate that adapter ligation to highly structured ribozyme fragments is much more efficient using the thermostable RtcB ligase from Pyrococcus horikoshii than the broadly applied Escherichia coli enzyme. Moreover, we investigated DNA, RNA and modified RNA adapters for their suitability in RtcB ligation reactions. We used the optimized RtcB-mediated ligation to produce RNA-seq libraries and captured a spiked 3' twister ribozyme fragment from E. coli total RNA. This RNA-seq-based method is applicable to detect ribozyme fragments as well as other cellular RNAs with 5' OH termini from total RNA.

Diversity and evolution of the animal virome

The COVID-19 pandemic has given the study of virus evolution and ecology new relevance. Although viruses were first identified more than a century ago, we likely know less about their diversity than that of any other biological entity. Most documented animal viruses have been sampled from just two phyla - the Chordata and the Arthropoda - with a strong bias towards viruses that infect humans or animals of economic and social importance, often in association with strong disease phenotypes. Fortunately, the recent development of unbiased metagenomic next-generation sequencing is providing a richer view of the animal virome and shedding new light on virus evolution. In this Review, we explore our changing understanding of the diversity, composition and evolution of the animal virome. We outline the factors that determine the phylogenetic diversity and genomic structure of animal viruses on evolutionary timescales and show how this impacts assessment of the risk of disease emergence in the short term. We also describe the ongoing challenges in metagenomic analysis and outline key themes for future research. A central question is how major events in the evolutionary history of animals, such as the origin of the vertebrates and periodic mass extinction events, have shaped the diversity and evolution of the viruses they carry.

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.

[Animal delta-like viruses (Kolmioviridae: Deltavirus) and the origin of the human hepatitis D virus (HDV)]

Hepatitis D (delta, δ) virus (HDV) was discovered more than 40 years ago, but the understanding of its origin and evolution is poor. This is mainly due to the lack, until recently, of data on the existence of any viruses similar to HDV. The discovery in recent years of sequences of new delta-like agents in a wide range of vertebrate (Vertebrata) and invertebrate (Invertebrata) species has facilitated a revision of views on the origin of HDV and contributed to understanding the place of this unique virus among other animals' viral agents. The purpose of this review is to analyze the latest published data on new delta-like agents and their biological characteristics.

Self-cleaving ribozymes: substrate specificity and synthetic biology applications

Various self-cleaving ribozymes appearing in nature catalyze the sequence-specific intramolecular cleavage of RNA and can be engineered to catalyze cleavage of appropriate substrates in an intermolecular fashion, thus acting as true catalysts. The mechanisms of the small, self-cleaving ribozymes have been extensively studied and reviewed previously. Self-cleaving ribozymes can possess high catalytic activity and high substrate specificity; however, substrate specificity is also engineerable within the constraints of the ribozyme structure. While these ribozymes share a common fundamental catalytic mechanism, each ribozyme family has a unique overall architecture and active site organization, indicating that several distinct structures yield this chemical activity. The multitude of catalytic structures, combined with some flexibility in substrate specificity within each family, suggests that such catalytic RNAs, taken together, could access a wide variety of substrates. Here, we give an overview of 10 classes of self-cleaving ribozymes and capture what is understood about their substrate specificity and synthetic applications. Evolution of these ribozymes in an RNA world might be characterized by the emergence of a new ribozyme family followed by rapid adaptation or diversification for specific substrates.

Self-cleaving ribozymes have become important tools of synthetic biology. Here we summarize the substrate specificity and applications of the main classes of these ribozymes.

The potential versatility of RNA catalysis

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

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.

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

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

Inverse RNA Folding Workflow to Design and Test Ribozymes that Include Pseudoknots

Ribozymes are RNAs that catalyze reactions. They occur in nature, and can also be evolved in vitro to catalyze novel reactions. This chapter provides detailed protocols for using inverse folding software to design a ribozyme sequence that will fold to a known ribozyme secondary structure and for testing the catalytic activity of the sequence experimentally. This protocol is able to design sequences that include pseudoknots, which is important as all naturally occurring full-length ribozymes have pseudoknots. The starting point is the known pseudoknot-containing secondary structure of the ribozyme and knowledge of any nucleotides whose identity is required for function. The output of the protocol is a set of sequences that have been tested for function. Using this protocol, we were previously successful at designing highly active double-pseudoknotted HDV ribozymes.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Experimental Resurrection of Ancestral Mammalian CPEB3 Ribozymes Reveals Deep Functional Conservation

Self-cleaving ribozymes are genetic elements found in all domains of life, but their evolution remains poorly understood. A ribozyme located in the second intron of the cytoplasmic polyadenylation binding protein 3 gene (CPEB3) shows high sequence conservation in mammals, but little is known about the functional conservation of self-cleaving ribozyme activity across the mammalian tree of life or during the course of mammalian evolution. Here, we use a phylogenetic approach to design a mutational library and a deep sequencing assay to evaluate the in vitro self-cleavage activity of numerous extant and resurrected CPEB3 ribozymes that span over 100 My of mammalian evolution. We found that the predicted sequence at the divergence of placentals and marsupials is highly active, and this activity has been conserved in most lineages. A reduction in ribozyme activity appears to have occurred multiple different times throughout the mammalian tree of life. The in vitro activity data allow an evaluation of the predicted mutational pathways leading to extant ribozyme as well as the mutational landscape surrounding these ribozymes. The results demonstrate that in addition to sequence conservation, the self-cleavage activity of the CPEB3 ribozyme has persisted over millions of years of mammalian evolution.

Improved RNA Secondary Structure and Tertiary Base-pairing Prediction Using Evolutionary Profile, Mutational Coupling and Two-dimensional Transfer Learning

MOTIVATION

The recent discovery of numerous non-coding RNAs (long non-coding RNAs, in particular) has transformed our perception about the roles of RNAs in living organisms. Our ability to understand them, however, is hampered by our inability to solve their secondary and tertiary structures in high resolution efficiently by existing experimental techniques. Computational prediction of RNA secondary structure, on the other hand, has received much-needed improvement, recently, through deep learning of a large approximate data, followed by transfer learning with gold-standard base-pairing structures from high-resolution 3-D structures. Here, we expand this single-sequence-based learning to the use of evolutionary profiles and mutational coupling.

RESULTS

The new method allows large improvement not only in canonical base-pairs (RNA secondary structures) but more so in base-pairing associated with tertiary interactions such as pseudoknots, noncanonical and lone base-pairs. In particular, it is highly accurate for those RNAs of more than 1000 homologous sequences by achieving >0.8 F1-score (harmonic mean of sensitivity and precision) for 14/16 RNAs tested. The method can also significantly improve base-pairing prediction by incorporating artificial but functional homologous sequences generated from deep mutational scanning without any modification. The fully automatic method (publicly available as server and standalone software) should provide the scientific community a new powerful tool to capture not only the secondary structure but also tertiary base-pairing information for building three-dimensional models. It also highlights the future of accurately solving the base-pairing structure by using a large number of natural and/or artificial homologous sequences.

AVAILABILITY

Standalone-version of SPOT-RNA2 is available at https://github.com/jaswindersingh2/SPOT-RNA2. Direct prediction can also be made at https://sparks-lab.org/server/spot-rna2/. The datasets used in this research can also be downloaded from the GITHUB and the webserver mentioned above.

Diversification of mammalian deltaviruses by host shifting

Hepatitis delta virus (HDV) is an unusual RNA agent that replicates using host machinery but exploits hepatitis B virus (HBV) to mobilize its spread within and between hosts. In doing so, HDV enhances the virulence of HBV. How this seemingly improbable hyperparasitic lifestyle emerged is unknown, but it underpins the likelihood that HDV and related deltaviruses may alter other host-virus interactions. Here, we show that deltaviruses diversify by transmitting between mammalian species. Among 96,695 RNA sequence datasets, deltaviruses infected bats, rodents, and an artiodactyl from the Americas but were absent from geographically overrepresented Old World representatives of each mammalian order, suggesting a relatively recent diversification within the Americas. Consistent with diversification by host shifting, both bat and rodent-infecting deltaviruses were paraphyletic, and coevolutionary modeling rejected cospeciation with mammalian hosts. In addition, a 2-y field study showed common vampire bats in Peru were infected by two divergent deltaviruses, indicating multiple introductions to a single host species. One vampire bat-associated deltavirus was detected in the saliva of up to 35% of individuals, formed phylogeographically compartmentalized clades, and infected a sympatric bat, illustrating horizontal transmission within and between species on ecological timescales. Consistent absence of HBV-like viruses in two deltavirus-infected bat species indicated acquisitions of novel viral associations during the divergence of bat and human-infecting deltaviruses. Our analyses support an American zoonotic origin of HDV and reveal prospects for future cross-species emergence of deltaviruses. Given their peculiar life history, deltavirus host shifts will have different constraints and disease outcomes compared to ordinary animal pathogens.

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.

The evolution and clinical impact of hepatitis B virus genome diversity

The global burden of hepatitis B virus (HBV) is enormous, with 257 million persons chronically infected, resulting in more than 880,000 deaths per year worldwide. HBV exists as nine different genotypes, which differ in disease progression, natural history and response to therapy. HBV is an ancient virus, with the latest reports greatly expanding the host range of the Hepadnaviridae (to include fish and reptiles) and casting new light on the origins and evolution of this viral family. Although there is an effective preventive vaccine, there is no cure for chronic hepatitis B, largely owing to the persistence of a viral minichromosome that is not targeted by current therapies. HBV persistence is also facilitated through aberrant host immune responses, possibly due to the diverse intra-host viral populations that can respond to host-mounted and therapeutic selection pressures. This Review summarizes current knowledge on the influence of HBV diversity on disease progression and treatment response and the potential effect on new HBV therapies in the pipeline. The mechanisms by which HBV diversity can occur both within the individual host and at a population level are also discussed.

RNA Drugs and RNA Targets for Small Molecules: Principles, Progress, and Challenges

RNA-based therapies, including RNA molecules as drugs and RNA-targeted small molecules, offer unique opportunities to expand the range of therapeutic targets. Various forms of RNAs may be used to selectively act on proteins, transcripts, and genes that cannot be targeted by conventional small molecules or proteins. Although development of RNA drugs faces unparalleled challenges, many strategies have been developed to improve RNA metabolic stability and intracellular delivery. A number of RNA drugs have been approved for medical use, including aptamers (e.g., pegaptanib) that mechanistically act on protein target and small interfering RNAs (e.g., patisiran and givosiran) and antisense oligonucleotides (e.g., inotersen and golodirsen) that directly interfere with RNA targets. Furthermore, guide RNAs are essential components of novel gene editing modalities, and mRNA therapeutics are under development for protein replacement therapy or vaccination, including those against unprecedented severe acute respiratory syndrome coronavirus pandemic. Moreover, functional RNAs or RNA motifs are highly structured to form binding pockets or clefts that are accessible by small molecules. Many natural, semisynthetic, or synthetic antibiotics (e.g., aminoglycosides, tetracyclines, macrolides, oxazolidinones, and phenicols) can directly bind to ribosomal RNAs to achieve the inhibition of bacterial infections. Therefore, there is growing interest in developing RNA-targeted small-molecule drugs amenable to oral administration, and some (e.g., risdiplam and branaplam) have entered clinical trials. Here, we review the pharmacology of novel RNA drugs and RNA-targeted small-molecule medications, with a focus on recent progresses and strategies. Challenges in the development of novel druggable RNA entities and identification of viable RNA targets and selective small-molecule binders are discussed. SIGNIFICANCE STATEMENT: With the understanding of RNA functions and critical roles in diseases, as well as the development of RNA-related technologies, there is growing interest in developing novel RNA-based therapeutics. This comprehensive review presents pharmacology of both RNA drugs and RNA-targeted small-molecule medications, focusing on novel mechanisms of action, the most recent progress, and existing challenges.

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

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

Deciphering the Self-Cleavage Reaction Mechanism of Hairpin Ribozyme

Hairpin ribozyme catalyzes the reversible self-cleavage of phosphodiester bonds which plays prominent roles in key biological processes involving RNAs. Despite impressive advances on ribozymatic self-cleavage, critical aspects of its molecular reaction mechanism remain controversially debated. Here, we generate and analyze the multidimensional free energy landscape that underlies the reaction using extensive QM/MM metadynamics simulations to investigate in detail the full self-cleavage mechanism. This allows us to answer several pertinent yet controversial questions concerning activation of the 2'-OH group, the mechanistic role of water molecules present in the active site, and the full reaction pathway including the structures of transition states and intermediates. Importantly, we find that a sufficiently unrestricted reaction subspace must be mapped using accelerated sampling methods in order to compute the underlying free energy landscape. It is shown that lower-dimensional sampling where the bond formation and cleavage steps are coupled does not allow the system to sufficiently explore the landscape. On the basis of a three-dimensional free energy surface spanned by flexible generalized coordinates, we find that 2'-OH is indirectly activated by adjacent G8 nucleobase in conjunction with stabilizing H-bonding involving water. This allows the proton of the 2'-OH group to directly migrate toward the 5'-leaving group via a nonbridging oxygen of the phosphodiester link. At variance with similar enzymatic processes where water wires connected to protonable side chains of the protein matrix act as transient proton shuttles, no such de/reprotonation events of water molecules are found to be involved in this ribozymatic transesterification. Overall, our results support an acid-catalyzed reaction mechanism where A38 nucleobase directly acts as an acid whereas G8, in stark contrast, participates only indirectly via stabilizing the nascent nucleophile for subsequent attack. Moreover, we conclude that self-cleavage of hairpin ribozyme follows an A_N + D_N two-step associative pathway where the rate-determining step is the cleavage of the phosphodiester bond. These results provide a major advancement in our understanding of the unique catalytic mechanism of hairpin ribozyme which will fruitfully impact on the design of synthetic ribozymes.

Current knowledge on Hepatitis Delta Virus replication

Hepatitis B Virus (HBV) that infects liver parenchymal cells is responsible for severe liver diseases and co-infection with Hepatitis Delta Virus (HDV) leads to the most aggressive form of viral hepatitis. Even tough being different for their viral genome (relaxed circular partially double stranded DNA for HBV and circular RNA for HDV), HBV and HDV are both maintained as episomes in the nucleus of infected cells and use the cellular machinery for the transcription of their viral RNAs. We propose here an update on the current knowledge on HDV replication cycle that may eventually help to identify new antiviral targets.

Transcription polymerase-catalyzed emergence of novel RNA replicons

Transcription polymerases can exhibit an unusual mode of regenerating certain RNA templates from RNA, yielding systems that can replicate and evolve with RNA as the information carrier. Two classes of pathogenic RNAs (hepatitis delta virus in animals and viroids in plants) are copied by host transcription polymerases. Using in vitro RNA replication by the transcription polymerase of T7 bacteriophage as an experimental model, we identify hundreds of new replicating RNAs, define three mechanistic hallmarks of replication (subterminal de novo initiation, RNA shape-shifting, and interrupted rolling-circle synthesis), and describe emergence from DNA seeds as a mechanism for the origin of novel RNA replicons. These results inform models for the origins and replication of naturally occurring RNA genetic elements and suggest a means by which diverse RNA populations could be propagated as hereditary material in cellular contexts.

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

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

rPredictorDB: a predictive database of individual secondary structures of RNAs and their formatted plots

Secondary data structure of RNA molecules provides insights into the identity and function of RNAs. With RNAs readily sequenced, the question of their structural characterization is increasingly important. However, RNA structure is difficult to acquire. Its experimental identification is extremely technically demanding, while computational prediction is not accurate enough, especially for large structures of long sequences. We address this difficult situation with rPredictorDB, a predictive database of RNA secondary structures that aims to form a middle ground between experimentally identified structures in PDB and predicted consensus secondary structures in Rfam. The database contains individual secondary structures predicted using a tool for template-based prediction of RNA secondary structure for the homologs of the RNA families with at least one homolog with experimentally solved structure. Experimentally identified structures are used as the structural templates and thus the prediction has higher reliability than de novo predictions in Rfam. The sequences are downloaded from public resources. So far rPredictorDB covers 7365 RNAs with their secondary structures. Plots of the secondary structures use the Traveler package for readable display of RNAs with long sequences and complex structures, such as ribosomal RNAs. The RNAs in the output of rPredictorDB are extensively annotated and can be viewed, browsed, searched and downloaded according to taxonomic, sequence and structure data. Additionally, structure of user-provided sequences can be predicted using the templates stored in rPredictorDB.

B2 and ALU retrotransposons are self-cleaving ribozymes whose activity is enhanced by EZH2

Transposable elements make up half of the mammalian genome. One of the most abundant is the short interspersed nuclear element (SINE). Among their million copies, B2 accounts for ∼350,000 in the mouse genome and has garnered special interest because of emerging roles in epigenetic regulation. Our recent work demonstrated that B2 RNA binds stress genes to retard transcription elongation. Although epigenetically silenced, B2s become massively up-regulated during thermal and other types of stress. Specifically, an interaction between B2 RNA and the Polycomb protein, EZH2, results in cleavage of B2 RNA, release of B2 RNA from chromatin, and activation of thermal stress genes. Although an established RNA-binding protein and histone methyltransferase, EZH2 is not known to be a nuclease. Here, we provide evidence for the surprising conclusion that B2 is a self-cleaving ribozyme. Ribozyme activity depends on Mg⁺² and monovalent cations but is resistant to protease treatment. However, contact with EZH2 accelerates cleavage rate by >100-fold, suggesting that EZH2 promotes a cleavage-competent RNA conformation. B2 modification-interference analysis demonstrates that phosphorothioate changes at A and C nucleotides can substitute for EZH2. B2 nucleotides 45 to 55 and 100 to 101 are essential for activity. Finally, another family of SINEs, the human ALU element, also produces a self-cleaving RNA and is cleaved during T-cell activation as well as thermal and endoplasmic reticulum (ER) stress. Thus, B2/ALU SINEs may be classified as "epigenetic ribozymes" that function as transcriptional switches during stress. Given their high copy numbers, B2 and ALU may represent the predominant ribozyme activity in mammalian cells.

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.

Cellular Concentrations of Nucleotide Diphosphate-Chelated Magnesium Ions Accelerate Catalysis by RNA and DNA Enzymes

RNAs are involved in myriad cellular events. In general, RNA function is affected by cellular conditions. For instance, molecular crowding promotes RNA folding through compaction of the RNA. Metabolites generally destabilize RNA secondary structure, which improves RNA folding cooperativity and increases the fraction of functional RNA. Our recent studies demonstrate that cellular concentrations of amino acid-chelated magnesium (aaCM) stimulate RNA folding and catalysis while protecting RNAs from magnesium ion-induced degradation. However, effects of other cellular magnesium ion chelators on RNA function have not been tested. Herein, we report that nucleotide diphosphate-chelated magnesium, which is of intermediate strength, promotes RNA catalysis much like aaCM. Nucleotides are some of the major metabolites in cells and have one to three phosphates, which have increasingly tight binding of magnesium. On the basis of binding calculations, ∼85% ATP, ∼40% ADP, and only 5% AMP are estimated to possess a magnesium ion under cellular conditions of 0.50 mM Mg²⁺_free. We tested the self-cleaving activity of the hammerhead ribozyme in the presence of these chelated magnesium species. Our results indicate that NTP-chelated magnesium and NMP-chelated magnesium do not appreciably stimulate RNA catalysis, whereas NDP-chelated magnesium promotes RNA catalysis up to 6.5-fold. Inspired by NDP, we observed similar stimulatory effects for several other Mg²⁺ diphosphate-containing metabolites, including UDP-GlcNAC and UDP-Glc; in addition, we found similar effects for a DNAzyme. Thus, rate stimulatory effects are general with respect to the diphosphate and nucleic acid enzyme. These results implicate magnesium-chelated diphosphate metabolites as general facilitators of RNA function inside cells.

Design of highly active double-pseudoknotted ribozymes: a combined computational and experimental study

Design of RNA sequences that adopt functional folds establishes principles of RNA folding and applications in biotechnology. Inverse folding for RNAs, which allows computational design of sequences that adopt specific structures, can be utilized for unveiling RNA functions and developing genetic tools in synthetic biology. Although many algorithms for inverse RNA folding have been developed, the pseudoknot, which plays a key role in folding of ribozymes and riboswitches, is not addressed in most algorithms. For the few algorithms that attempt to predict pseudoknot-containing ribozymes, self-cleavage activity has not been tested. Herein, we design double-pseudoknot HDV ribozymes using an inverse RNA folding algorithm and test their kinetic mechanisms experimentally. More than 90% of the positively designed ribozymes possess self-cleaving activity, whereas more than 70% of negative control ribozymes, which are predicted to fold to the necessary structure but with low fidelity, do not possess it. Kinetic and mutation analyses reveal that these RNAs cleave site-specifically and with the same mechanism as the WT ribozyme. Most ribozymes react just 50- to 80-fold slower than the WT ribozyme, and this rate can be improved to near WT by modification of a junction. Thus, fast-cleaving functional ribozymes with multiple pseudoknots can be designed computationally.

Identification and validation of a prognostic four-genes signature for hepatocellular carcinoma: integrated ceRNA network analysis

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most aggressive malignant tumors, with a poor long-term prognosis worldwide. The functional deregulations of global transcriptome were associated with the genesis and development of HCC, but lacks systematic research and validation.

METHODS

A total of 519 postoperative HCC patients were included. We built an interactive and visual competing endogenous RNA network. The prognostic signature was established with the least absolute shrinkage and selection operator algorithm. Multivariate Cox regression analysis was used to screen for independent prognostic factors for HCC overall survival.

RESULTS

In the training set, we identified a four-gene signature (PBK, CBX2, CLSPN, and CPEB3) and effectively predicted the overall survival. The survival times of patients in the high-score group were worse than those in the low-score group (p = 0.0004), and death was also more likely in the high-score group (HR 2.444, p < 0.001). The results were validated in internal validation set (p = 0.0057) and two external validation cohorts (HR 2.467 and 2.6). The signature (AUCs of 1, 2, 3 years were 0.716, 0.726, 0.714, respectively) showed high prognostic accuracy in the complete TCGA cohort.

CONCLUSIONS

In conclusion, we successfully built a more extensive ceRNA network for HCC and then identified a four-gene-based signature, enabling prediction of the overall survival of patients with HCC.

Novel hepatitis D-like agents in vertebrates and invertebrates

Hepatitis delta virus (HDV) is the smallest known RNA virus, encoding a single protein. Until recently, HDV had only been identified in humans, where it is strongly associated with co-infection with hepatitis B virus (HBV). However, the recent discovery of HDV-like viruses in metagenomic samples from birds and snakes suggests that this virus has a far longer evolutionary history. Herein, using additional meta-transcriptomic data, we show that highly divergent HDV-like viruses are also present in fish, amphibians, and invertebrates, with PCR and Sanger sequencing confirming the presence of the invertebrate HDV-like viruses. Notably, the novel viruses identified here share genomic features characteristic of HDV, such as a circular genome of only approximately 1.7 kb in length, and self-complementary, unbranched rod-like structures. Coiled-coil domains, leucine zippers, conserved residues with essential biological functions, and isoelectronic points similar to those in the human hepatitis delta virus antigens (HDAgs) were also identified in the putative non-human viruses. Importantly, none of these novel HDV-like viruses were associated with hepadnavirus infection, supporting the idea that the HDV–HBV association may be specific to humans. Collectively, these data not only broaden our understanding of the diversity and host range of HDV, but also shed light on its origin and evolutionary history.

Enveloped viruses distinct from HBV induce dissemination of hepatitis D virus in vivo

Hepatitis D virus (HDV) doesn't encode envelope proteins for packaging of its ribonucleoprotein (RNP) and typically relies on the surface glycoproteins (GPs) from hepatitis B virus (HBV) for virion assembly, envelopment and cellular transmission. HDV RNA genome can efficiently replicate in different tissues and species, raising the possibility that it evolved, and/or is still able to transmit, independently of HBV. Here we show that alternative, HBV-unrelated viruses can act as helper viruses for HDV. In vitro, envelope GPs from several virus genera, including vesiculovirus, flavivirus and hepacivirus, can package HDV RNPs, allowing efficient egress of HDV particles in the extracellular milieu of co-infected cells and subsequent entry into cells expressing the relevant receptors. Furthermore, HCV can propagate HDV infection in the liver of co-infected humanized mice for several months. Further work is necessary to evaluate whether HDV is currently transmitted by HBV-unrelated viruses in humans.