Search

Viroids are circular, highly structured, single-stranded, non-coding RNA pathogens known to infect and cause disease in several plant species. They are known to trigger the host plant's RNA silencing machinery. The detection of viroid-derived small RNAs (vd-sRNA) in viroid-infected host plants opened a new avenue of study in host–viroid pathogenicity. Since then, several viroid research groups have studied the vd-sRNA retrieved from different host–viroid combinations. Such studies require the segregation of 21- to 24-nucleotide long small RNAs (sRNA) from a deep-sequencing databank, followed by separating the vd-sRNA from any sRNA within this group that showed sequence similarity with either the genomic or the antigenomic strands of the viroid. Such mapped vd-sRNAs are then profiled on both the viroid's genomic and antigenomic strands for visualization. Although several commercial interfaces are currently available for this purpose, they are all programmed for linear RNA molecules. Hence, viroid researchers must develop a computer program that accommodates the sRNAs derived from the circular viroid genome. This is a laborious process, and consequently, it often creates a bottleneck for biologists. In order to overcome this constraint, and to help the research community in general, in this study, a python-based pattern matching interface was developed so as to be able to both profile and map sRNAs on a circular genome. A "matching tolerance" feature has been included in the program, thus permitting the mapping of the sRNAs derived from the quasi-species. Additionally, the "topology" feature allows the researcher to profile sRNA derived from both linear and circular RNA molecules. The efficiency of the program was tested using previously reported deep-sequencing data obtained from two independent studies. Clearly, this novel software should be a key tool with which to both evaluate the production of sRNA and to profile them on their target RNA species, irrespective of the topology of the target RNA molecule. [ABSTRACT FROM AUTHOR]

Viroids are the smallest non-coding infectious RNAs (between 246 and 401 nucleotides) known to be highly structured and replicate autonomously in the host plants. Although they do not encode any peptides, viroids induce visible symptoms in susceptible host plants. This article provides an overview of their physical and biological properties, the diseases they cause and their significance for the plants. The mechanisms underlying the expression of symptoms in host plants, their detection and various strategies employed for diseases prevention are also developed.

Single nucleotide polymorphisms (SNPs) are common genetic variations that are present in over 1% of the population and can significantly modify the structures of both DNA and RNA. G-quadruplex structures (G4) are formed by the superposition of tetrads of guanines. To date, the impact of SNPs on both G4 ligands' binding efficacies and specificities has not been investigated. Here, using a bioinformatically predicted G4 and SNPs found in the α-synuclein gene as a proof-of-concept, it was demonstrated that SNPs can modulate both DNA and RNA G4s' responses to ligands. Specifically, six widely recognized ligands (Phen-DC3, PDS, 360A, RHPS4, BRACO19, and TMPyP4) were shown to differentially affect both the structure and the polymerase stalling of the different SNPs. This work highlights the importance of choosing the appropriate G4 ligand when dealing with an SNP identified in a G-rich gene.

Background: G-quadruplexes (G4s) are secondary structures in DNA and RNA that impact various cellular processes, such as transcription, splicing, and translation. Due to their numerous functions, G4s are involved in many diseases, making their study important. Yet, G4s evolution remains largely unknown, due to their low sequence similarity and the poor quality of their sequence alignments across several species. To address this, we designed a strategy that avoids direct G4s alignment to study G4s evolution in the 3 species kingdoms. We also explored the coevolution between RBPs and G4s., Methods: We retrieved one-to-one orthologous genes from the Ensembl Compara database and computed groups of one-to-one orthologous genes. For each group, we aligned gene sequences and identified G4 families as groups of overlapping G4s in the alignment. We analyzed these G4 families using Count, a tool to infer feature evolution into a gene or a species tree. Additionally, we utilized these G4 families to predict G4s by homology. To establish a control dataset, we performed mono-, di- and tri-nucleotide shuffling., Results: Only a few conserved G4s occur among all living kingdoms. In eukaryotes, G4s exhibit slight conservation among vertebrates, and few are conserved between plants. In archaea and bacteria, at most, only 2 G4s are common. The G4 homology-based prediction increases the number of conserved G4s in common ancestors. The coevolution between RNA-binding proteins and G4s was investigated and revealed a modest impact of RNA-binding proteins evolution on G4 evolution. However, the details of this relationship remain unclear., Conclusion: Even if G4 evolution still eludes us, the present study provides key information to compute groups of homologous G4 and to reveal the evolution history of G4 families., Competing Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2023.)

Hop latent viroid (HLVd) is the biggest concern for cannabis and hop growers worldwide. Although most HLVd-infected plants remain asymptomatic, research on hops has demonstrated a decrease in both the α-bitter acid and terpene content of hop cones, which affects their economic value. The HLVd-associated "dudding" or "duds" disease of cannabis was first reported in 2019 in California. Since then, the disease has become widespread in cannabis-growing facilities across North America. Although severe yield loss associated with duds disease has been recorded, little scientific information is available to growers in order to contain HLVd. Consequently, this review aims to summarise all of the scientific information available on HLVd so as to be able to understand the effect of HLVd on yield loss, cannabinoid content, terpene profile, disease management and inform crop protection strategies.

G-quadruplexes (G4) are 3D structures that are found in both DNA and RNA. Interest in this structure has grown over the past few years due to both its implication in diverse biological mechanisms and its potential use as a therapeutic target, to name two examples. G4s in humans have been widely studied; however, the level of their study in other species remains relatively minimal. That said, progress in this field has resulted in the prediction of G4s structures in various species, ranging from bacteria to eukaryotes. These predictions were analysed in a previous study which revealed that G4s are present in all living kingdoms. To date, eleven different databases have grouped the various G4s depending on either their structures, on the proteins that might bind them, or on their location in the various genomes. However, none of these databases contains information on their location in the transcriptome of many of the implicated species. The GAIA database was designed so as to make this data available online in a user-friendly manner. Through its web interface, users can query GAIA to filter G4s, which, we hope, will help the research in this field. GAIA is available at: https://gaia.cobius.usherbrooke.ca., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Mammalian orthoreovirus (MRV) is a double-stranded RNA virus from the Reoviridae family presenting a promising activity as an oncolytic virus. Recent studies have underlined MRV's ability to alter cellular alternative splicing (AS) during infection, with a limited understanding of the mechanisms at play. In this study, we investigated how MRV modulates AS. Using a combination of cell biology and reverse genetics experiments, we demonstrated that the M1 gene segment, encoding the μ2 protein, is the primary determinant of MRV's ability to alter AS, and that the amino acid at position 208 in μ2 is critical to induce these changes. Moreover, we showed that the expression of μ2 by itself is sufficient to trigger AS changes, and its ability to enter the nucleus is not required for all these changes. Moreover, we identified core components of the U5 snRNP (i.e. EFTUD2, PRPF8, and SNRNP200) as interactors of μ2 that are required for MRV modulation of AS. Finally, these U5 snRNP components are reduced at the protein level by both MRV infection and μ2 expression. Our findings identify the reduction of U5 snRNP components levels as a new mechanism by which viruses alter cellular AS., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

G-quadruplexes are motifs found in DNA and RNA that can fold into tertiary structures. Until now, they have been studied experimentally mainly in humans and a few other species. Recently, predictions have been made with bacterial and archaeal genomes. Nevertheless, a global comparison of predicted G4s (pG4s) across and within the three living kingdoms has not been addressed. In this study, we aimed to predict G4s in genes and transcripts of all kingdoms of living organisms and investigated the differences in their distributions. The relation of the predictions with GC content was studied. It appears that GC content is not the only parameter impacting G4 predictions and abundance. The distribution of pG4 densities varies depending on the class of transcripts and the group of species. Indeed, we have observed that, in coding transcripts, there are more predicted G4s than expected for eukaryotes but not for archaea and bacteria, while in noncoding transcripts, there are as many or fewer predicted G4s in all species groups. We even noticed that some species with the same GC content presented different pG4 profiles. For instance, Leishmania major and Chlamydomonas reinhardtii both have 60% of GC content, but the former has a pG4 density of 0.07 and the latter 1.16., (© The Author(s) 2022. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

Quebec is the third largest wine grape producer in Canada in acreage, tonnage, and wine grape sales (Carisse et al. 2017; Ben Moussa et al. 2019). To evaluate the diversity of viruses infecting grapevine in Quebec, a total of 77 leaf tissue samples (cv. Vidal) were collected from July to October in 2020 in three different vineyards located in Frelighsburg, Hemmingford and Saint-Jacques-le-Mineur in Quebec, Canada. Double-stranded RNA was extracted from each sample and used for cDNA library preparation with the Nextera XT DNA Library Preparation Kit (Illumina) as described previously (Kesanakurti et al. 2016). High-throughput sequencing (HTS, 2x300 bp) was conducted on dual-indexed libraries in a v3 flow cell using the Illumina MiSeq platform (Adkar-Purushothama et al. 2020). The obtained raw FASTQ data was de-multiplexed into 154 separate sequence files, and the adapters and barcode sequences were trimmed. The quality of the sequences was verified using Trimmomatic V.0.32 and the "clean" sequences were analyzed using Virtool and VirFind virus detection pipelines described elsewhere (Ho and Tzanetakis 2014; Rott et al. 2017) to screen for all possible viruses in the databases. Over 100,000 reads per sample were obtained with a percentage of mapped viral reads ranging from 1.47 to 19.43% of total number of reads. Out of 77 samples, 16 revealed the sequence of grapevine yellow speckle viroid 1 (GYSVd-1), for which the length coverage ranged from 98.5 to 99.1%; the depth ranged from 2X to 856X. The GYSVd-1 positive sequence files were subjected to whole genome assembly on CLC genomics Workbench v20.0.4 with the isolate SY-BR from Brazil (KU880715) used as reference. Seven complete genomes of GYSVd-1 of 366-368 nucleotides (nt) in size were deposited (GenBank Acc. MW732682 to MW732688). BLASTN analysis of the sequences showed 98-100% nt identities with isolate SY-BR. Other viruses and viroids such as Grapevine fleck virus, Grapevine rupestris stem pitting-associated virus, Grapevine rupestris vein feathering virus and Hop stunt viroid were also detected. To confirm GYSVd-1 presence in Quebec vineyards, seven of the 16 HTS-positive grapevine leaf tissue samples were subjected to total RNA extraction, followed by RT-PCR assay as before (Adkar-Purushothama et al. 2015; Sahana et al. 2013); all were positive by RT-PCR. The PCR products were directly Sanger-sequenced, and they showed 100% nt identity to the HTS derived sequences. Three of the seven GYSVd-1 positive grapevines exhibited yellow leaf spots and flecks and tiny yellow leaves, but their mixed infection status makes definitive symptoms association difficult to determine. Previously, Hop stunt viroid was reported from grapevines in Canada (Xiao et al. 2019; Fall et al. 2020) but to the best of our knowledge, this is the first report of GYSVd-1 infecting grapevines in Canada, specifically in the province of Quebec. Further research is required to assess the GYSVd-1 related yield loss. Monitoring and testing for GYSVd-1 infection is necessary to prevent propagation of infected materials, spread, and potential negative impact for the Canadian grapevine industry.

The early 1970s marked two breakthroughs in the field of biology: (i) The development of nucleotide sequencing technology; and, (ii) the discovery of the viroids. The first DNA sequences were obtained by two-dimensional chromatography which was later replaced by sequencing using electrophoresis technique. The subsequent development of fluorescence-based sequencing method which made DNA sequencing not only easier, but many orders of magnitude faster. The knowledge of DNA sequences has become an indispensable tool for both basic and applied research. It has shed light biology of viroids, the highly structured, circular, single-stranded non-coding RNA molecules that infect numerous economically important plants. Our understanding of viroid molecular biology and biochemistry has been intimately associated with the evolution of nucleic acid sequencing technologies. With the development of the next-generation sequence method, viroid research exponentially progressed, notably in the areas of the molecular mechanisms of viroids and viroid diseases, viroid pathogenesis, viroid quasi-species, viroid adaptability, and viroid-host interactions, to name a few examples. In this review, the progress in the understanding of viroid biology in conjunction with the improvements in nucleotide sequencing technology is summarized. The future of viroid research with respect to the use of third-generation sequencing technology is also briefly envisaged.

Viroids are non-coding RNA plant pathogens that are characterized by their possession of a high mutation level. Although the sequence heterogeneity in viroid infected plants is well understood, shifts in viroid population dynamics due to mutations over the course of infection remain poorly understood. In this study, the ten most abundant sequence variants of potato spindle tuber viroid RG1 (PSTVd) expressed at different time intervals in PSTVd infected tomato plants were identified by high-throughput sequencing. The sequence variants, forming a quasi-species, were subjected to both the identification of the regions favoring mutations and the effect of the mutations on viroid secondary structure and viroid derived small RNAs (vd-sRNA). At week 1 of PSTVd infection, 25% of the sequence variants were similar to the "master" sequence (i.e., the sequence used for inoculation). The frequency of the master sequence within the population increased to 70% at week 2 after PSTVd infection, and then stabilized for the rest of the disease cycle (i.e., weeks 3 and 4). While some sequence variants were abundant at week 1 after PSTVd infection, they tended to decrease in frequency over time. For example, the variants with insertions at positions 253 or 254, positions that could affect the Loop E as well as the metastable hairpin I structure that has been shown important during replication and viroid infectivity, resulted in decreased frequency. Data obtained by in silico analysis of the viroid derived small RNAs (vd-sRNA) was also analyzed. A few mutants had the potential of positively affecting the viroid's accumulation by inducing the RNA silencing of the host's defense related genes. Variants with mutations that could negatively affect viroid abundance were also identified because their derived vd-sRNA were no longer capable of targeting any host mRNA or of changing its target sequence from a host defense gene to some other non-important host gene. Together, these findings open avenues into understanding the biological role of sequence variants, this viroid's interaction with host components, stable and metastable structures generated by mutants during the course of infection, and the influence of sequence variants on stabilizing viroid population dynamics., (Copyright © 2020 Adkar-Purushothama, Bolduc, Bru and Perreault.)

It has been demonstrated that RNA G-quadruplexes (G4) are structural motifs present in transcriptomes and play important regulatory roles in several post-transcriptional mechanisms. However, the full picture of RNA G4 locations and the extent of their implication remain elusive. Solely computational prediction analysis of the whole transcriptome may reveal all potential G4, since experimental identifications are always limited to specific conditions or specific cell lines. The present study reports the first in-depth computational prediction of potential G4 region across the complete human transcriptome. Although using a relatively stringent approach based on three prediction scores that accounts for the composition of G4 sequences, the composition of their neighboring sequences, and the various forms of G4, over 1.1 million of potential G4 (pG4) were predicted. The abundance of G4 was computationally confirmed in both 5' and 3'UTR as well as splicing junction of mRNA, appreciate for the first time in the long ncRNA, while almost absent of most of the small ncRNA families. The present results constitute an important step toward a full understanding of the roles of G4 in post-transcriptional mechanisms., (© The Author(s) 2019. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

G-Quadruplexes are non-canonical secondary structures that can be adopted under physiological conditions by guanine-rich DNA and RNA molecules. They have been reported to occur, and to perform multiple biological functions, in the genomes and transcriptomes of many species, including humans. This chapter focuses specifically on RNA G-quadruplexes and reviews the most recent discoveries in the field, as well as addresses the upcoming challenges researchers studying these structures face.

The anti-apoptotic BAG-1 protein isoforms are known to be overexpressed in colorectal tumors and are considered to be potential therapeutic targets. The isoforms are derived from alternative translation initiations occuring at four in-frame start codons of a single mRNA transcript. Its 5'UTR also contains an internal ribosome entry site (IRES) regulating the cap-independent translation of the transcript. An RNA G-quadruplex (rG4) is located at the 5'end of the BAG-1 5'UTR, upstream of the known cis-regulatory elements. Herein, we observed that the expression of BAG-1 isoforms is post-transcriptionally regulated in colorectal cancer cells and tumors, and that stabilisation of the rG4 by small molecules ligands reduces the expression of endogenous BAG-1 isoforms. We demonstrated a critical role for the rG4 in the control of both cap-dependent and independent translation of the BAG-1 mRNA in colorectal cancer cells. Additionally, we found an upstream ORF that also represses BAG-1 mRNA translation. The structural probing of the complete 5'UTR showed that the rG4 acts as a steric block which controls the initiation of translation at each start codon of the transcript and also maintains the global 5'UTR secondary structure required for IRES-dependent translation., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Viroids are naked RNAs that do not code for any known protein and yet are able to infect plants causing severe diseases. Because of their RNA nature, many studies have focused on the involvement of viroids in RNA-mediated gene silencing as being their pathogenesis mechanism. Here, the alterations caused by the Citrus exocortis viroid (CEVd) on the tomato translation machinery were studied as a new aspect of viroid pathogenesis. The presence of viroids in the ribosomal fractions of infected tomato plants was detected. More precisely, CEVd and its derived viroid small RNAs were found to co-sediment with tomato ribosomes in vivo, and to provoke changes in the global polysome profiles, particularly in the 40S ribosomal subunit accumulation. Additionally, the viroid caused alterations in ribosome biogenesis in the infected tomato plants, affecting the 18S rRNA maturation process. A higher expression level of the ribosomal stress mediator NAC082 was also detected in the CEVd-infected tomato leaves. Both the alterations in the rRNA processing and the induction of NAC082 correlate with the degree of viroid symptomatology. Taken together, these results suggest that CEVd is responsible for defective ribosome biogenesis in tomato, thereby interfering with the translation machinery and, therefore, causing ribosomal stress., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related deaths around the world. Recent advances in genomic technologies have allowed the identification of various molecular signatures in HCC tissues. For instance, differential gene expression levels of various cytochrome P450 genes (CYP450) have been reported in studies performed on limited numbers of HCC tissue samples, or focused on a small subset on CYP450s. In the present study, we monitored the expression landscape of all the members of the CYP450 family (57 genes) in more than 200 HCC tissues using RNA-Seq data from The Cancer Genome Atlas. Using stringent statistical filters and data from paired tissues, we identified significantly dysregulated CYP450 genes in HCC. Moreover, the expression level of selected CYP450s was validated by qPCR on cDNA samples from an independent cohort. Threshold values (sensitivity and specificity) based on dysregulated gene expression were also determined to allow for confident identification of HCC tissues. Finally, a global look at expression levels of the 57 members of the CYP450 family across ten different cancer types revealed specific expression signatures. Overall, this study provides useful information on the transcriptomic landscape of CYP450 genes in HCC and on new potential HCC biomarkers., Competing Interests: CONFLICTS OF INTEREST The authors declare that they have no competing interests.

RNA G-quadruplexes (rG4) are stable non-canonical secondary structures composed of G-rich sequences. Many rG4 structures located in the 5'UTRs of mRNAs act as translation repressors due to their high stability which is thought to impede ribosomal scanning. That said, it is not known if these are mRNA-specific examples, or if they are indicative of a global expression regulation mechanism of the mRNAs involved in a common pathway based on structure folding recognition. Gene-ontology analysis of mRNAs bearing a predicted rG4 motif in their 5'UTRs revealed an enrichment for mRNAs associated with the colorectal cancer pathway. Bioinformatic tools for rG4 prediction, and experimental in vitro validations were used to confirm and compare the folding of the predicted rG4s of the mRNAs associated with dysregulated pathways in colorectal cancer. The rG4 folding was confirmed for the first time for 9 mRNAs. A repressive effect of 3 rG4 candidates on the expression of a reporter gene was also measured in colorectal cancer cell lines. This work highlights the fact that rG4 prediction is not yet accurate, and that experimental characterization is still essential in order to identify the precise rG4 folding sequences and the possible common features shared between the rG4 overrepresented in important biological pathways., Competing Interests: The authors have declared that no competing interests exist.

Though RNA G-quadruplexes became a focus of study over a decade ago, the main challenge associated with the identification of new potential G-quadruplexes remains a bottleneck step. It slows the study of these non-canonical structures in nucleic acids, and thus the understanding of their significance. The G4RNA screener is an accurate tool for the prediction of RNA G-quadruplexes but its deployment has brought to light an issue with its accessibility to G-quadruplex experts and biologists. G4RNA screener web server is a platform that provides a much needed interface to manage the input, parameters and result display of the main command-line ready tool. It is accessible at http://scottgroup.med.usherbrooke.ca/G4RNA_screener/., (Copyright © 2018 Elsevier B.V. and Société Française de Biochimie et Biologie Moléculaire (SFBBM). All rights reserved.)

Understanding in intimate details how the viroid interaction with host's defense genes is a cornerstone for developing viroid resistant plants. In this present study, small RNAs (sRNA) derived from Potato spindle tuber viroid (PSTVd) were studied in silico in order to detect any interactions with the serine threonine kinase receptor, a transmembrane protein that plays a role in disease resistance in plants. Using molecular biology techniques, it was determined that PSTVd infection negatively affects at least three serine threonine kinase receptors as well as with three other genes that are known to be involved in the overall development of the tomato plants. The transient expression of these putative PSTVd-sRNAs, using the microRNA sequence as a backbone, in tomato plants induced phenotypes similar to viroid infection. Mutants created by altering the sequence of PSTVd in these regions failed to infect the tomato plant. The data presented here illustrates the importance of these regions in viroid survival, and suggests a possible avenue of exploration for the development of viroid resistant plants.

Motivation: G-quadruplex structures in RNA molecules are known to have regulatory impacts in cells but are difficult to locate in the genome. The minimal requirements for G-quadruplex folding in RNA (G≥3N1-7 G≥3N1-7 G≥3N1-7 G≥3) is being challenged by observations made on specific examples in recent years. The definition of potential G-quadruplex sequences has major repercussions on the observation of the structure since it introduces a bias. The canonical motif only describes a sub-population of the reported G-quadruplexes. To address these issues, we propose an RNA G-quadruplex prediction strategy that does not rely on a motif definition., Results: We trained an artificial neural network with sequences of experimentally validated G-quadruplexes from the G4RNA database encoded using an abstract definition of their sequence. This artificial neural network, G4NN, evaluates the similarity of a given sequence to known G-quadruplexes and reports it as a score. G4NN has a predictive power comparable to the reported G richness and G/C skewness evaluations that are the current state-of-the-art for the identification of potential RNA G-quadruplexes. We combined these approaches in the G4RNA screener, a program designed to manage and evaluate the sequences to identify potential G-quadruplexes., Availability and Implementation: G4RNA screener is available for download at http://gitlabscottgroup.med.usherbrooke.ca/J-Michel/g4rna_screener., Contact: jean-michel.garant@usherbrooke.ca or jean-pierre.perreault@usherbrooke.ca or michelle.scott@usherbrooke.ca., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2017. Published by Oxford University Press.)

Summary: Discovering function-related structural features, such as the cloverleaf shape of transfer RNA secondary structures, is essential to understand RNA function. With this aim, we have developed a platform, named Structurexplor, to facilitate the exploration of structural features in populations of RNA secondary structures. It has been designed and developed to help biologists interactively search for, evaluate and select interesting structural features that can potentially explain RNA functions., Availability and Implementation: Structurxplor is a web application available at http://structurexplor.dinf.usherbrooke.ca. The source code can be found at http://jpsglouzon.github.io/structurexplor/., Contact: shengrui.wang@usherbrooke.ca., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author (2017). Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com)

It is well established that viroid derived small RNA (vd-sRNA) induces RNA silencing of endogenous mRNA. However, it remains not clear how exactly viroid infections can lead to severe symptom induction given the fact that fewer vd-sRNAs binding the specific target mRNAs were recovered from the infected plants. To answer this question, the two least expressed (+) and (-) strand vd-sRNAs of potato spindle tuber viroid (PSTVd) binding to both the 3' UTR and the coding region of tomato mRNAs were analyzed by infecting tomato plants with two variants of PSTVd. As products of these putative target mRNAs are involved in plant phenotype, the effect of this viroid on these genes were analyzed by infecting tomato plants with two variants of PSTVd. The direct interaction between the vd-sRNAs and putative mRNAs was validated by artificial microRNA experiments in a transient expression system and by RNA ligase-mediated rapid amplification of cDNA ends. Parallel analysis of RNA ends of viroid infected plants revealed the widespread cleavage of the target mRNAs in locations other than the vd-sRNA binding site during the viroid infection implying the viroid-infection induced vd-sRNA independent degradation of endogenous mRNAs during viroid infection.

The simplest known plant pathogens are the viroids. Because of their non-coding single-stranded circular RNA genome, they depend on both their sequence and their structure for both a successful infection and their replication. In the recent years, important progress in the elucidation of their structures was achieved using an adaptation of the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) protocol in order to probe viroid structures in solution. Previously, SHAPE has been adapted to elucidate the structures of all of the members of the family Avsunviroidae, as well as those of a few members of the family Pospiviroidae. In this study, with the goal of providing an entire compendium of the secondary structures of the various viroid species, a total of thirteen new Pospiviroidae members were probed in solution using the SHAPE protocol. More specifically, the secondary structures of eleven species for which the genus was previously known were initially elucidated. At this point, considering all of the SHAPE elucidated secondary structures, a classification system for viroids in their respective genera was proposed. On the basis of the structural classification reported here, the probings of both the Grapevine latent viroid and the Dahlia latent viroid provide sound arguments for the determination of their respective genera, which appear to be Apscaviroid and Hostuviroid, respectively. More importantly, this study provides the complete repertoire of the secondary structures, mapped in solution, of all of the accepted viroid species reported thus far. In addition, a classification scheme based on structural hallmarks, an important tool for many biological studies, is proposed.

MicroRNAs (miRNAs) are small noncoding RNAs that repress the translation of their target genes. It has previously been shown that a target's availability to miRNA can be affected by its structure. G-quadruplexes (G4) are noncanonical structures adopted by G-rich nucleic acids that have been shown to have multiple biological functions. In this study, whether or not G4 structures' presence in the 3' UTRs of mRNAs can hinder miRNA binding was investigated. Putative G4 overlapping with predicted miRNAs' binding sites was searched for, and 44,294 hits were found in humans. The FADS2 mRNA/mir331-3p pair was selected as a model example. In-line probing and G4-specific fluorescent ligand experiments binding were performed and confirmed the presence of a G4 near the predicted miRNA binding site. Subsequent luciferase assays showed that the presence of the G4 prevents the binding of mir331-3p in cellulo. Together, these results served as proof of concept that a G4 structure present in a 3' UTR sequence should be taken into consideration when predicting miRNA binding sites., (© 2017 Rouleau et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Motivation: Comparing ribonucleic acid (RNA) secondary structures of arbitrary size uncovers structural patterns that can provide a better understanding of RNA functions. However, performing fast and accurate secondary structure comparisons is challenging when we take into account the RNA configuration (i.e. linear or circular), the presence of pseudoknot and G-quadruplex (G4) motifs and the increasing number of secondary structures generated by high-throughput probing techniques. To address this challenge, we propose the super-n-motifs model based on a latent analysis of enhanced motifs comprising not only basic motifs but also adjacency relations. The super-n-motifs model computes a vector representation of secondary structures as linear combinations of these motifs., Results: We demonstrate the accuracy of our model for comparison of secondary structures from linear and circular RNA while also considering pseudoknot and G4 motifs. We show that the super-n-motifs representation effectively captures the most important structural features of secondary structures, as compared to other representations such as ordered tree, arc-annotated and string representations. Finally, we demonstrate the time efficiency of our model, which is alignment free and capable of performing large-scale comparisons of 10 000 secondary structures with an efficiency up to 4 orders of magnitude faster than existing approaches., Availability and Implementation: The super-n-motifs model was implemented in C ++. Source code and Linux binary are freely available at http://jpsglouzon.github.io/supernmotifs/ ., Contact: Shengrui.Wang@Usherbrooke.ca., Supplementary Information: Supplementary data are available at Bioinformatics o nline., (© The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com)

DNA guanine (G)-rich 4-stranded helical nucleic acid structures called G-quadruplexes (G4), have been extensively studied during the last decades. However, emerging evidence reveals that 5'- and 3'-untranslated regions (5'- and 3'-UTRs) as well as open reading frames (ORFs) contain putative RNA G-quadruplexes. These stable secondary structures play key roles in telomere homeostasis and RNA metabolism including pre-mRNA splicing, polyadenylation, mRNA targeting and translation. Interestingly, multiple RNA binding proteins such as nucleolin, FMRP, DHX36, and Aven were identified to bind RNA G-quadruplexes. Moreover, accumulating reports suggest that RNA G-quadruplexes regulate translation in cap-dependent and -independent manner. Herein, we discuss potential roles of RNA G-quadruplexes and associated trans- acting factors in the regulation of mRNA translation.

Mature viroids consist of a noncoding, covalently closed circular RNA that is able to autonomously infect respective host plants. Thus, they must utilize proteins of the host for most biological functions such as replication, processing, transport, and pathogenesis. Therefore, viroids can be regarded as minimal parasites of the host machinery. They have to present to the host machinery the appropriate signals based on either their sequence or their structure. Here, we summarize such sequence and structural features critical for the biological functions of viroids., (© 2016 Elsevier Inc. All rights reserved.)

Antisense-based molecules targeting HIV-1 RNA have the potential to be used as part of gene or drug therapy to treat HIV-1 infection. In this study, HIV-1 RNA was screened to identify more conserved and accessible target sites for ribozymes based on the hepatitis delta virus motif. Using a quantitative screen for effects on HIV-1 production, we identified a ribozyme targeting a highly conserved site in the Gag coding sequence with improved inhibitory potential compared to our previously described candidates targeting the overlapping Tat/Rev coding sequence. We also demonstrate that this target site is highly accessible to short hairpin directed RNA interference, suggesting that it may be available for the binding of antisense RNAs with different modes of action. We provide evidence that this target site is structurally conserved in diverse viral strains and that it is sufficiently different from the human transcriptome to limit off-target effects from antisense therapies. We also show that the modified hepatitis delta virus ribozyme is more sensitive to a mismatch in its target site compared to the short hairpin RNA. Overall, our results validate the potential of a new target site in HIV-1 RNA to be used for the development of antisense therapies.

G-quadruplexes are widespread four-stranded structures that are adopted by G-rich regions of both DNA and RNA and are involved in essential biological processes such as mRNA translation. They are formed by the stacking of two or more G-quartets that are linked together by three loops. Although the maximal loop length is usually fixed to 7 nt in most G-quadruplex-predicting software, it has already been demonstrated that artificial DNA G-quadruplexes containing two distal loops that are limited to 1 nt each and a central loop up to 30 nt long are likely to form in vitro. This report demonstrates that such structures possessing a long central loop are actually found in the 5'-UTRs of human mRNAs. Firstly, 1453 potential G-quadruplex-forming sequences (PG4s) were identified through a bioinformatic survey that searched for sequences respecting the requirement for two 1-nt long distal loops and a long central loop of 2-90 nt in length. Secondly, in vitro in-line probing experiments confirmed and characterized the folding of eight candidates possessing central loops of 10-70 nt long. Finally, the biological effect of several G-quadruplexes with a long central loop on mRNA expression was studied in cellulo using a luciferase gene reporter assay. Clearly, the actual definition of G-quadruplex-forming sequences is too conservative and must be expanded to include the long central loop. This greatly expands the number of expected PG4s in the transcriptome. Consideration of these new candidates might aid in elucidating the potentially important biological implications of the G-quadruplex structure., (© 2014 Jodoin et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Viroids are small, circular, single stranded RNA molecules that infect plants. Since they are non-coding, their structures play a critical role in their life cycles. To date, little effort has been spend on elucidating viroid structures in solution due to both the experimental difficulties and the time-consuming nature of the methodologies implicated. Recently, the technique of high-throughput selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) was adapted for the probing of the members of family Avsunviroidae, all of whom replicate in the chloroplast and demonstrate ribozyme activity. In the present work, twelve viroid species belonging to four different genera of the family Pospiviroidae, whose members are characterized by the presence of a central conserved region (CCR) and who replicate in nucleus of the host, were probed. Given that the structures of five distinct viroid species from the family Pospiviroidae have been previously reported, an overview of the different structural characteristics for all genera and the beginning of a manual classification of the different viroids based on their structural features are presented here.

Nucleic acid therapies targeting HIV replication have the potential to be used in conjunction with or in place of the standard small-molecule therapies. Among the different classes of nucleic acid therapies, several ribozymes (Rzs, RNA enzymes) have been developed to target HIV RNA. The design of Rzs targeting HIV RNA is complicated by the sequence diversity of viral strains and the structural diversity of their target sites. Using the SOFA-HDV Rz as an example, this chapter describes methods that can be used to design Rzs for controlling HIV replication. We describe how to (1) identify highly conserved Rz target sites in HIV RNA; (2) generate a set of Rzs with the potential to be used as therapeutics; and (3) screen these Rzs for activity against HIV production.

G-quadruplexes (G4s) are non-canonical structures involved in many important cellular processes. To date, the prediction of potential G-quadruplex structures (PG4s) has been based almost exclusively on the sequence of interest agreeing with the algorithm Gx-N-1-7-Gx-N1-7-Gx-N1-7-Gx (where x ≥ 3 and N = A, U, G or C). However, many sequences agreeing with this algorithm do not form G4s and are considered false-positive predictions. Here we show the RNA PG4 candidate in the 3'-untranslated region (UTR) of the TTYH1 gene to be one such false positive. Specifically, G4 folding was observed to be inhibited by the presence of multiple-cytosine tracks, located in the candidate's genomic context, that adopted a Watson-Crick base-paired structure. Clearly, the neighbouring sequences of a PG4 may influence its folding. The secondary structure of 12 PG4 motifs along with either 15 or 50 nucleotides of their upstream and downstream genomic contexts were evaluated by in-line probing. Data permitted the development of a scoring system for the prediction of PG4s taking into account the effect of the neighbouring sequences. The accuracy of this scoring system was assessed by probing 14 other novel PG4 candidates retrieved in human 5'-UTRs. This new scoring system can be used, in combination with the standard algorithm, to better predict the folding of RNA G4s.

Guanine-rich RNA sequences can fold into non-canonical, four stranded helical structures called G-quadruplexes that have been shown to be widely distributed within the mammalian transcriptome, as well as being key regulatory elements in various biological mechanisms. That said, their role within the 3'-untranslated region (UTR) of mRNA remains to be elucidated and appreciated. A bioinformatic analysis of the 3'-UTRs of mRNAs revealed enrichment in G-quadruplexes. To shed light on the role(s) of these structures, those found in the LRP5 and FXR1 genes were characterized both in vitro and in cellulo. The 3'-UTR G-quadruplexes were found to increase the efficiencies of alternative polyadenylation sites, leading to the expression of shorter transcripts and to possess the ability to interfere with the miRNA regulatory network of a specific mRNA. Clearly, G-quadruplexes located in the 3'-UTRs of mRNAs are cis-regulatory elements that have a significant impact on gene expression.

Small cis-acting ribozymes have been converted into trans-acting ribozymes possessing the ability to cleave RNA substrates. The Hepatitis Delta Virus (HDV) ribozyme is one of the rare examples of these that is derived from an RNA species that is found in human cells. Consequently, it possesses the natural ability to function in the presence of human proteins in addition to an outstanding stability in human cells, two significant advantages in its use. The development of an additional specific on/off adaptor (SOFA) has led to the production of a new generation of HDV ribozymes with improved specificities that provide a tool with significant potential for future development in the fields of both functional genomics and gene -therapy. SOFA-HDV ribozyme-based gene inactivation systems have been reported in both prokaryotic and eukaryotic cells. Here, a step-by-step approach for the efficient design of highly specific SOFA-HDV ribozymes with a minimum investment of time and effort is described.

Given that greater than 90% of the human genome is expressed, it is logical to assume that post-transcriptional regulatory mechanisms must be the primary means of controlling the flow of information from mRNA to protein. This report describes a robust approach that includes in silico, in vitro and in cellulo experiments permitting an in-depth evaluation of the impact of G-quadruplexes as translational repressors. Sequences including potential G-quadruplexes were selected within nine distinct genes encoding proteins involved in various biological processes. Their abilities to fold into G-quadruplex structures in vitro were evaluated using circular dichroism, thermal denaturation and the novel use of in-line probing. Six sequences were observed to fold into G-quadruplex structures in vitro, all of which exhibited translational inhibition in cellulo when linked to a reporter gene. Sequence analysis, direct mutagenesis and subsequent experiments were performed in order to define the rules governing the folding of G-quadruplexes. In addition, the impact of single-nucleotide polymorphism was shown to be important in the formation of G-quadruplexes located within the 5'-untranslated region of an mRNA. In light of these results, clearly the 5'-UTR G-quadruplexes represent a class of translational repressors that is broadly distributed in the cell.

Decapping enzymes are required for the removal of the 5'-end (m7)GpppN cap of mRNAs to allow their decay in cells. While many cap-binding proteins recognize the cap structure via the stacking of the methylated guanosine ring between two aromatic residues, the precise mechanism of cap recognition by decapping enzymes has yet to be determined. In order to get insights into the interaction of decapping enzymes with the cap structure, we studied the vaccinia virus D10 decapping enzyme as a model to investigate the important features for substrate recognition by the enzyme. We demonstrate that a number of chemically modified purines can competitively inhibit the decapping reaction, highlighting the molecular features of the cap structure that are required for recognition by the enzyme, such as the nature of the moiety at positions 2 and 6 of the guanine base. A 3D structural model of the D10 protein was generated which suggests amino acids implicated in cap binding. Consequently, we expressed 17 mutant proteins with amino acid substitutions in the active site of D10 and found that eight are critical for the decapping activity. These data underscore the functional features involved in the non-canonical cap-recognition by the vaccinia virus D10 decapping enzyme.

The transcription factor Pax-5, is vital during B lymphocyte differentiation and is known to contribute to the oncogenesis of certain cancers. The Pax-5 locus generates multiple yet structurally related mRNA transcripts through the specific activation of alternative promoter regions and/or alternative splicing events which poses challenges in the study of specific isoform function. In this study, we investigated the function of a major Pax-5 transcript, Pax-5B using an enhanced version of the Hepatitis Delta Virus ribozyme (HDV Rz) suppression system that is specifically designed to recognize and cleave the human Pax-5B mRNA. The activity of these ribozymes resulted in the specific suppression of the Pax-5B transcripts without altering the transcript levels of other closely related Pax-5 isoforms mRNAs both in vitro and in an intracellular setting. Following stable transfection of the ribozymes into a model B cell line (REH), we showed that Pax-5B suppression led to an increase of CD19 mRNA and cell surface protein expression. In response to this Pax-5B specific deregulation, a marked increase in apoptotic activity compared to control cell lines was observed. These results suggest that Pax-5B has distinct roles in physiological processes in cell fate events during lymphocyte development.

An RNA guanylyltransferase activity is involved in the synthesis of the cap structure found at the 5' end of eukaryotic mRNAs. The RNA guanylyltransferase activity is a two-step ping-pong reaction in which the enzyme first reacts with GTP to produce the enzyme-GMP covalent intermediate with the concomitant release of pyrophosphate. In the second step of the reaction, the GMP moiety is then transferred to a diphosphorylated RNA. Both reactions were previously shown to be reversible. In this study, we report a biochemical and thermodynamic characterization of both steps of the reaction of the RNA guanylyltransferase from Paramecium bursaria Chlorella virus 1, the prototype of a family of viruses infecting green algae. Using a combination of real-time fluorescence spectroscopy, radioactive kinetic assays, and inhibition assays, the complete kinetic parameters of the RNA guanylyltransferase were determined. We produced a thermodynamic scheme for the progress of the reaction as a function of the energies involved in each step. We were able to demonstrate that the second step comprises the limiting steps for both the direct and reverse overall reactions. In both cases, the binding to the RNA substrates is the step requiring the highest energy and generating unstable intermediates that will promote the catalytic activites of the enzyme. This study reports the first thorough kinetic and thermodynamic characterization of the reaction catalyzed by an RNA capping enzyme.

The Chlorella virus RNA triphosphatase (cvRTPase) is involved in the formation of the RNA cap structure found at the 5'-end of the viral mRNAs and requires magnesium ions to mediate its catalytic activity. To extend our studies on the role of metal ions in phosphohydrolysis, we have used a combination of fluorescence spectroscopy, circular dichroism, denaturation studies and thermodynamic analyses to monitor the binding of magnesium ions to the cvRTPase. Using these techniques, the thermodynamic forces responsible for the interaction of metal ions with an RNA triphosphatase were also evaluated for the first time. Our thermodynamic analyses indicate that the initial association of magnesium with the cvRTPase is dominated by a favorable entropic effect and is accompanied by the release of eight water molecules from the enzyme. Moreover, both fluorescence spectroscopy and circular dichroism assays indicated that minor conformational changes were occurring upon magnesium binding. Mutational studies were also performed and confirmed the importance of three specific glutamate residues located in the active site of the enzyme for the binding of magnesium ions. Finally, in contrast to the yeast RNA triphosphatase, we demonstrate that the binding of magnesium ions to the cvRTPase does not lead to the stabilization of the ground state binding of the RNA substrate. Based on the results of the present study, we hypothesize that the binding of magnesium ions induces local conformational perturbations in the active site residues that ultimately positions the lateral chains of critical amino acids involved in catalysis. Our results highlight fundamental differences in the role of magnesium ions in the phosphohydrolase reactions catalyzed by the cvRTPase and the closely related yeast RNA triphosphatase.

The Y genes encode small noncoding RNAs whose functions remain elusive, whose numbers vary between species, and whose major property is to be bound by the Ro60 protein (or its ortholog in other species). To better understand the evolution of the Y gene family, we performed a homology search in 27 different genomes along with a structural search using Y RNA specific motifs. These searches confirmed that Y RNAs are well conserved in the animal kingdom and resulted in the detection of several new Y RNA genes, including the first Y RNAs in insects and a second Y RNA detected in Caenorhabditis elegans. Unexpectedly, Y5 genes were retrieved almost as frequently as Y1 and Y3 genes, and, consequently are not the result of a relatively recent apparition as is generally believed. Investigation of the organization of the Y genes demonstrated that the synteny was conserved among species. Interestingly, it revealed the presence of six putative "fossil" Y genes, all of which were Y4 and Y5 related. Sequence analysis led to inference of the ancestral sequences for all Y RNAs. In addition, the evolution of existing Y RNAs was deduced for many families, orders and classes. Moreover, a consensus sequence and secondary structure for each Y species was determined. Further evolutionary insight was obtained from the analysis of several thousand Y retropseudogenes among various species. Taken together, these results confirm the rich and diversified evolution history of Y RNAs.

A bioinformatic covariation analysis of a collection of 119 novel variants of the antigenomic, self-cleaving hepatitis delta virus (HDV) RNA motif supported the formation of all of the Watson-Crick base pairs (bp) of the catalytic centre except the C19-G81 pair located at the bottom of the P2 stem. In fact, a novel Watson-Crick bp between C19 and G80 is suggested by the data. Both chemical and enzymatic probing demonstrated that initially the C19-G81 pair is formed in the ribozyme (Rz), but upon substrate (S) binding and the formation of the P1.1 pseudoknot C19 switches its base-pairing partner from G81 to G80. As a result of this finding, the secondary structure of this ribozyme has been redrawn. The formation of the C19-G80 bp results in a J4/2 junction composed of four nucleotides, similar to that seen in the genomic counterpart, thereby increasing the similarities between these two catalytic RNAs. Additional mutagenesis, cleavage activity and probing experiments yield an original characterization of the structural features involving the residues of the J4/2 junction.

In order to revisit the architecture of the catalytic center of the antigenomic hepatitis delta virus (HDV) ribozyme we developed an unbiased in vitro selection procedure that efficiently selected novel variants from a relatively small set of sequences. Using this procedure we examined all possible variants from a pool of HDV ribozymes that had been randomized at 25 positions (4(25)). The isolated set of sequences shows more variability than do the natural variants. Nucleotide variations were found at all randomized positions, even at positions when the general belief was that the specific base was absolutely required for catalytic activity. Covariation analysis supports the presence of several base pairs, although it failed to propose any new tertiary contacts. HDV ribozyme appears to possess a greater number of constraints, in terms of sequences capable of supporting the catalysed cleavage, than do other catalytic RNAs. This supports the idea that the appearance of this catalytic RNA structure has a low probability (i.e. is a rare event), which may explain why to date it has been found in nature only in the HDV. These contrasts with the hammerhead self-cleaving motif that is proposed to have multiple origins, and that is widespread among different organisms. Thus, just because a self-cleaving RNA motif is small does not imply that it occurs easily.

The trans-acting antigenomic delta ribozyme, isolated from the human hepatitis delta virus, was shown to be highly stable and active in vitro, as well as in mammalian cell lines. However, the stability and gene-targeting competence of this small ribozyme have not been studied previously in bacterial cells. In this paper we describe the use of two variants of the trans-acting antigenomic delta ribozyme targeting the abundant EF-Tu mRNA in the industrially important gram-positive bacterium Lactococcus lactis. These two delta ribozyme variants were expressed at significant levels and were shown to be highly stable in vivo. The half-life of the EF-Tu mRNA was slightly but consistently reduced in the presence of the classical delta ribozymes (7 to 13%). In contrast, delta ribozymes harboring a specific on/off riboswitch (SOFA-delta ribozymes) targeting the same sites on the EF-Tu mRNA considerably reduced the half-life of this mRNA (22 to 47%). The rates of catalysis of the SOFA-delta ribozymes in L. lactis were similar to the rates determined in vitro, showing that this new generation of delta ribozymes was highly efficient in these bacterial cells. Clearly, SOFA-delta ribozymes appear to be an ideal means for development of gene inactivation systems in bacteria.

We report the characterization in the human genome of 966 pseudogenes derived from the four human Y (hY) RNAs, components of the Ro/SS-A autoantigen. About 95% of the Y RNA pseudogenes are found in corresponding locations on the chimpanzee and human chromosomes. On the contrary, Y pseudogenes in mice are both infrequent and found in different genomic regions. In addition to this rodent/primate discrepancy, the conservation of hY pseudogenes relative to hY genes suggests that they occurred after rodent/primate divergence. Flanking regions of hY pseudogenes contain convincing evidence for involvement of the L1 retrotransposition machinery. Although Alu elements are found in close proximity to most hY pseudogenes, these are not chimeric retrogenes. Point mutations in hY RNA transcripts specifically affecting binding of Ro60 protein likely contributed to their selection for direct trans retrotransposition. This represents a novel requirement for the selection of specific RNAs for their genomic integration by the L1 retrotransposition machinery. Over 40% of the hY pseudogenes are found in intronic regions of protein-coding genes. Considering the functions of proteins known to bind subsets of hY RNAs, hY pseudogenes constitute a new class of L1-dependent non-autonomous retroelements, potentially involved in post-transcriptional regulation of gene expression.

Ribozymes, RNA molecules that catalyze the cleavage of RNA substrates, provide an interesting alternative to the RNA interference (RNAi) approach to gene inactivation, especially given the fact that RNAi seems to trigger an immunological response. Unfortunately, the limited substrate specificity of ribozymes is considered to be a significant hurdle in their development as molecular tools. Here, we report the molecular engineering of a ribozyme possessing a new biosensor module that switches the cleavage activity from 'off' (a 'safety lock') to 'on' solely in the presence of the appropriate RNA target substrate. Both proof-of-concept and the mechanism of action of this man-made riboswitch are demonstrated using hepatitis delta virus ribozymes that cleave RNA transcripts derived from the hepatitis B and C viruses. To our knowledge, this is the first report of a ribozyme bearing a target-dependent module that is activated by its RNA substrate, an arrangement which greatly diminishes non-specific effects. This new approach provides a highly specific and improved tool with significant potential for application in the fields of both functional genomics and gene therapy.

The ability of ribozymes (i.e. RNA enzymes) to specifically recognize and subsequently catalyze the cleavage of an RNA substrate makes them attractive for the development of therapeutic tools for the inactivation of both viral RNAs and mRNAs associated with various diseases. Several applicable ribozyme models have been tested both in vitro and in a cellular environment, and have shown significant promise. However, several hurdles remain to be surpassed before we generate a useful gene-inactivation system based on a ribozyme. Among the most important requirements for further progress are a better understanding of the features that contribute to defining the substrate specificity for cleavage by a ribozyme, and the identification of the potential cleavage sites in a given target RNA. The goal of this review is to illustrate the importance of both of these factors at the RNA level in the development of any type of ribozyme based gene-therapy. This is achieved by reviewing the recent progress in both the structure-function relationships and the development of a gene-inactivation system of a model ribozyme, specifically delta ribozyme.

Hepatitis delta virus RNAs possess self-cleavage activities that produce 2',3'-cyclic phosphate and 5'-hydroxyl termini (i.e. cis-acting delta ribozyme). Trans-acting delta ribozymes have been engineered by removing a junction from the cis version, thereby producing one molecule possessing the substrate sequence and the other the catalytic domain. According to the pseudoknot model, the secondary structure of the delta ribozyme includes a pseudoknot (i.e. P1.1 stem) formed by two base pairs from residues of the L3 loop and J1/4 junction. A collection of 48 P1.1 stem mutants was synthesized in order to provide an original characterization of both the importance and the structure of this pseudoknot in a trans-acting version of the ribozyme. Several structural differences were noted compared to the results reported for cis-acting ribozymes. For example, a combination of two stable Watson-Crick base pairs composing the essential P1.1 stem was demonstrated to be crucial for a significant level of activity, while the cis version required only one base pair. In addition, we present the first physical evidences revealing that the composition of the P1.1 stem affects the substrate specificity for ribozyme cleavage. Depending on the residues forming the J1/4 junction, non-productive ribozyme-substrate complexes can be observed. This phenomenon is proposed to be important for further development of a gene-inactivation system based on delta ribozyme.

Delta ribozyme possesses several unique features related to the fact that it is the only catalytic RNA known to be naturally active in human cells. This makes it attractive as a therapeutic tool for the inactivation of clinically relevant RNAs. However, several hurdles must be overcome prior to the development of useful gene-inactivation systems based on delta ribozyme. We have developed three procedures for the selection of potential delta ribozyme target sites within the hepatitis B virus (HBV) pregenome: (i) the use of bioinformatic tools coupled to biochemical assays; (ii) RNase H hydrolysis with a pool of oligonucleotides; and (iii) cleavage assays with a pool of ribozymes. The results obtained with delta ribozyme show that these procedures are governed by several rules, some of which are different from those both for other catalytic RNAs and antisense oligonucleotides. Together, these procedures identified 12 sites in the HBV pregenome that can be cleaved by delta ribozymes, although with different efficiencies. Clearly, both target site accessibility and the ability to form an active ribozyme-substrate complex constitute interdependent factors that can best be addressed using a combinatorial library of either oligonucleotides or ribozymes.

We have identified ribose 2'-hydroxyl groups (2'-OHs) that are critical for the activity of a trans-cleaving delta ribozyme derived from the antigenomic strand of the hepatitis delta virus. Initially, an RNA-DNA mixed ribozyme composed of 26 deoxyribo- (specifically the nucleotides forming the P2 stem and the P4 stem-loop) and 31 ribonucleotides (those forming the catalytic center) was engineered. This mixed ribozyme catalyzed the cleavage of a small substrate with kinetic parameters virtually identical to those of the all-RNA ribozyme. The further substitution of deoxyribose for ribose residues permitted us to investigate the contribution of all 2'-OHs to catalysis. Determination of the kinetic parameters for the cleavage reaction of the resulting ribozymes revealed (i) 10 2'-OH groups appear to be important in supporting the formation of several hydrogen bonds within the catalytic core, (ii) none of the important 2'-OHs seem to coordinate a magnesium cation, and (iii) 1 of the tested RNA-DNA mixed polymers appeared to stabilize the ribozyme-substrate transition-state complex, resulting in an improvement over the all-RNA counterpart. The contribution of the 2'-OHs to the catalytic mechanism is discussed, and differences with the crystal structure of a genomic delta self-cleaved product are explained. Clearly, the 2'-OHs are essential components of the network of interactions involved in the formation of the catalytic center of the delta ribozyme.