Formation of metastable RNA structures by sequential folding during transcription: time-resolved structural analysis of potato spindle tuber viroid (-)-stranded RNA by temperature-gradient gel electrophoresis

Viroids, Satellite RNAs and Prions: Folding of Nucleic Acids and Misfolding of Proteins

Theodor ("Ted") Otto Diener (* 28 February 1921 in Zürich, Switzerland; † 28 March 2023 in Beltsville, MD, USA) pioneered research on viroids while working at the Plant Virology Laboratory, Agricultural Research Service, USDA, in Beltsville. He coined the name viroid and defined viroids' important features like the infectivity of naked single-stranded RNA without protein-coding capacity. During scientific meetings in the 1970s and 1980s, viroids were often discussed at conferences together with other "subviral pathogens". This term includes what are now called satellite RNAs and prions. Satellite RNAs depend on a helper virus and have linear or, in the case of virusoids, circular RNA genomes. Prions, proteinaceous infectious particles, are the agents of scrapie, kuru and some other diseases. Many satellite RNAs, like viroids, are non-coding and exert their function by thermodynamically or kinetically controlled folding, while prions are solely host-encoded proteins that cause disease by misfolding, aggregation and transmission of their conformations into infectious prion isoforms. In this memorial, we will recall the work of Ted Diener on subviral pathogens.

Effect of pausing on the cotranscriptional folding kinetics of RNAs

The pausing event in RNA molecule folding that occurs during the cotranscription process plays a critical role in cellular RNA function. Based on Conformational Resampling through Kinetic Relaxation (CRKR), this paper investigates a method to reveal the specifics of pausing for RNA cotranscriptional folding of 117-nucleotide E. coli SRP RNA and 73-nucleotide HIV-1 TAR RNA. It can be inferred from the results that pausing events generate valid cotranscriptional conformational rearrangement to protect the function structures and influence the folding pathway, which is remarkably consistent with the experimental results. Additionally, different transcription speeds result in different levels of protection capability. The folding pathway or conformational rearrangement can also be affected by a change in transcription speed after pausing site. These findings illuminate that RNAs with functional structures undergo complex rearrangement folding in pausing events, and different RNAs prefer different speeds for specific folding pathways. Moreover, pausing in particular circumstances can be employed to regulate the population of final functional structures. In conclusion, this paper offers fresh perspectives on the pausing event in the cotranscriptional folding of RNAs.

Predicting the Structure of a Viroid : Structure, Structure Distribution, Consensus Structure, and Structure Drawing

Viroids are small non-coding RNAs that require a special sequence and structure to be replicated and transported by the host machinery. Many of these features can be predicted and later experimentally verified. Here, we will present workflows to predict viroid structures and draw the predicted structures in a pleasing and descriptive way using recently developed software.

Insights Into Potato Spindle Tuber Viroid Quasi-Species From Infection to Disease

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.

Current overview on viroid-host interactions

Viroids are one of the most enigmatic highly structured, circular, single-stranded RNA phytopathogens. Although they are not known to code for any peptide, viroids induce visible symptoms in susceptible host plants that resemble those associated with many plant viruses. It is known that viroids induce disease symptoms by direct interaction with host factors; however, the precise mechanism by which this occurs remains poorly understood. Studies on the host's responses to viroid infection, host susceptibility and nonhost resistance have been underway for several years, but much remains to be done in order to fully understand the complex nature of viroid-host interactions. Recent progress using molecular biology techniques combined with computational algorithms, in particular evidence of the role of viroid-derived small RNAs in the RNA silencing pathways of a disease network, has widened the knowledge of viroid pathogenicity. The complexity of viroid-host interactions has been revealed in the past decades to include, but not be limited to, the involvement of host factors, viroid structural complexity, and viroid-induced ribosomal stress, which is further boosted by the discovery of long noncoding RNAs (lncRNAs). In this review, the current understanding of the viroid-host interaction has been summarized with the goal of simplifying the complexity of viroid biology for future research. This article is categorized under: RNA in Disease and Development > RNA in Disease.

A genetic selection reveals functional metastable structures embedded in a toxin-encoding mRNA

Post-transcriptional regulation plays important roles to fine-tune gene expression in bacteria. In particular, regulation of type I toxin-antitoxin (TA) systems is achieved through sophisticated mechanisms involving toxin mRNA folding. Here, we set up a genetic approach to decipher the molecular underpinnings behind the regulation of a type I TA in Helicobacter pylori. We used the lethality induced by chromosomal inactivation of the antitoxin to select mutations that suppress toxicity. We found that single point mutations are sufficient to allow cell survival. Mutations located either in the 5’ untranslated region or within the open reading frame of the toxin hamper its translation by stabilizing stem-loop structures that sequester the Shine-Dalgarno sequence. We propose that these short hairpins correspond to metastable structures that are transiently formed during transcription to avoid premature toxin expression. This work uncovers the co-transcriptional inhibition of translation as an additional layer of TA regulation in bacteria.

Viroid research and its significance for RNA technology and basic biochemistry

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

Conformational Dynamics of the RNA G-Quadruplex and its Effect on Translation Efficiency

During translation, intracellular mRNA folds co-transcriptionally and must refold following the passage of ribosome. The mRNAs can be entrapped in metastable structures during these folding events. In the present study, we evaluated the conformational dynamics of the kinetically favored, metastable, and hairpin-like structure, which disturbs the thermodynamically favored G-quadruplex structure, and its effect on co-transcriptional translation in prokaryotic cells. We found that nascent mRNA forms a metastable hairpin-like structure during co-transcriptional folding instead of the G-quadruplex structure. When the translation progressed co-transcriptionally before the metastable hairpin-like structure transition to the G-quadruplex, function of the G-quadruplex as a roadblock of the ribosome was sequestered. This suggested that kinetically formed RNA structures had a dominant effect on gene expression in prokaryotes. The results of this study indicate that it is critical to consider the conformational dynamics of RNA-folding to understand the contributions of the mRNA structures in controlling gene expression.

Predicting Cotranscriptional Folding Kinetics For Riboswitch

On the basis of a helix-based transition rate model, we developed a new method for sampling cotranscriptional RNA conformational ensemble and the prediction of cotranscriptional folding kinetics. Applications to E. coli. SRP RNA and pbuE riboswitch indicate that the model may provide reliable predictions for the cotranscriptional folding pathways and population kinetics. For E. coli. SRP RNA, the predicted population kinetics and the folding pathway are consistent with the SHAPE profiles in the recent cotranscriptional SHAPE-seq experiments. For the pbuE riboswitch, the model predicts the transcriptional termination efficiency as a function of the force. The theoretical results show (a) a force-induced transition from the aptamer (antiterminator) to the terminator structure and (b) the different folding pathways for the riboswitch with and without the ligand (adenine). More specifically, without adenine, the aptamer structure emerges as a short-lived kinetic transient state instead of a thermodynamically stable intermediate state. Furthermore, from the predicted extension-time curves, the model identifies a series of conformational switches in the pulling process, where the predicted relative residence times for the different structures are in accordance with the experimental data. The model may provide a new tool for quantitative predictions of cotranscriptional folding kinetics, and results can offer useful insights into cotranscriptional folding-related RNA functions such as regulation of gene expression with riboswitches.

Transcription elongation rate affects nascent histone pre-mRNA folding and 3' end processing

In this study, Saldi et al. investigated how transcription elongation rate influences cotranscriptional pre-mRNA maturation. Their findings show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

Transcription elongation rate influences cotranscriptional pre-mRNA maturation, but how such kinetic coupling works is poorly understood. The formation of nonadenylated histone mRNA 3′ ends requires recognition of an RNA structure by stem–loop-binding protein (SLBP). We report that slow transcription by mutant RNA polymerase II (Pol II) caused accumulation of polyadenylated histone mRNAs that extend past the stem–loop processing site. UV irradiation, which decelerates Pol II elongation, also induced long poly(A)⁺ histone transcripts. Inhibition of 3′ processing by slow Pol II correlates with failure to recruit SLBP to histone genes. Chemical probing of nascent RNA structure showed that the stem–loop fails to fold in transcripts made by slow Pol II, thereby explaining the absence of SLBP and failure to process 3′ ends. These results show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

RNA Regulations and Functions Decoded by Transcriptome-wide RNA Structure Probing

RNA folds into intricate structures that are crucial for its functions and regulations. To date, a multitude of approaches for probing structures of the whole transcriptome, i.e., RNA structuromes, have been developed. Applications of these approaches to different cell lines and tissues have generated a rich resource for the study of RNA structure–function relationships at a systems biology level. In this review, we first introduce the designs of these methods and their applications to study different RNA structuromes. We emphasize their technological differences especially their unique advantages and caveats. We then summarize the structural insights in RNA functions and regulations obtained from the studies of RNA structuromes. And finally, we propose potential directions for future improvements and studies.

Dissecting the secondary structure of the circular RNA of a nuclear viroid in vivo: A "naked" rod-like conformation similar but not identical to that observed in vitro

With a minimal (250-400 nt), non-protein-coding, circular RNA genome, viroids rely on sequence/structural motifs for replication and colonization of their host plants. These motifs are embedded in a compact secondary structure whose elucidation is crucial to understand how they function. Viroid RNA structure has been tackled in silico with algorithms searching for the conformation of minimal free energy, and in vitro by probing in solution with RNases, dimethyl sulphate and bisulphite, and with selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), which interrogates the RNA backbone at single-nucleotide resolution. However, in vivo approaches at that resolution have not been assayed. Here, after confirming by 3 termodynamics-based predictions and by in vitro SHAPE that the secondary structure adopted by the infectious monomeric circular (+) RNA of potato spindle tuber viroid (PSTVd) is a rod-like conformation with double-stranded segments flanked by loops, we have probed it in vivo with a SHAPE modification. We provide direct evidence that a similar, but not identical, rod-like conformation exists in PSTVd-infected leaves of Nicotiana benthamiana, verifying the long-standing view that this RNA accumulates in planta as a "naked" form rather than tightly associated with protecting host proteins. However, certain nucleotides of the central conserved region, including some of the loop E involved in key functions such as replication, are more SHAPE-reactive in vitro than in vivo. This difference is most likely due to interactions with proteins mediating some of these functions, or to structural changes promoted by other factors of the in vivo habitat.

Contact Statistics Highlight Distinct Organizing Principles of Proteins and RNA

Although both RNA and proteins have densely packed native structures, chain organizations of these two biopolymers are fundamentally different. Motivated by the recent discoveries in chromatin folding that interphase chromosomes have territorial organization with signatures pointing to metastability, we analyzed the biomolecular structures deposited in the Protein Data Bank and found that the intrachain contact probabilities, P(s) as a function of the arc length s, decay in power-law ∼s(-γ) over the intermediate range of s, 10 ≲ s ≲ 110. We found that the contact probability scaling exponent is γ ≈ 1.11 for large RNA (N > 110), γ ≈ 1.41 for small-sized RNA (N < 110), and γ ≈ 1.65 for proteins. Given that Gaussian statistics is expected for a fully equilibrated chain in polymer melts, the deviation of γ-value from γ = 1.5 for the subchains of large RNA in the native state suggests that the chain configuration of RNA is not fully equilibrated. It is visually clear that folded structures of large-sized RNA (N ≳ 110) adopt crumpled structures, partitioned into modular multidomains assembled by proximal sequences along the chain, whereas the polypeptide chain of folded proteins looks better mixed with the rest of the structure. Our finding of γ ≈ 1 for large RNA might be an ineluctable consequence of the hierarchical ordering of the secondary to tertiary elements in the folding process.

In silico methods for co-transcriptional RNA secondary structure prediction and for investigating alternative RNA structure expression

RNA transcripts are the primary products of active genes in any living organism, including many viruses. Their cellular destiny not only depends on primary sequence signals, but can also be determined by RNA structure. Recent experimental evidence shows that many transcripts can be assigned more than a single functional RNA structure throughout their cellular life and that structure formation happens co-transcriptionally, i.e. as the transcript is synthesised in the cell. Moreover, functional RNA structures are not limited to non-coding transcripts, but can also feature in coding transcripts. The picture that now emerges is that RNA structures constitute an additional layer of information that can be encoded in any RNA transcript (and on top of other layers of information such as protein-context) in order to exert a wide range of functional roles. Moreover, different encoded RNA structures can be expressed at different stages of a transcript's life in order to alter the transcript's behaviour depending on its actual cellular context. Similar to the concept of alternative splicing for protein-coding genes, where a single transcript can yield different proteins depending on cellular context, it is thus appropriate to propose the notion of alternative RNA structure expression for any given transcript. This review introduces several computational strategies that my group developed to detect different aspects of RNA structure expression in vivo. Two aspects are of particular interest to us: (1) RNA secondary structure features that emerge during co-transcriptional folding and (2) functional RNA structure features that are expressed at different times of a transcript's life and potentially mutually exclusive.

The transcription initiation sites of eggplant latent viroid strands map within distinct motifs in their in vivo RNA conformations

Eggplant latent viroid (ELVd), like other members of family Avsunviroidae, replicates in plastids through a symmetric rolling-circle mechanism in which elongation of RNA strands is most likely catalyzed by a nuclear-encoded polymerase (NEP) translocated to plastids. Here we have addressed where NEP initiates transcription of viroid strands. Because this step is presumably directed by sequence/structural motifs, we have previously determined the conformation of the monomeric linear (+) and (-) RNAs of ELVd resulting from hammerhead-mediated self-cleavage. In silico predictions with 3 softwares led to similar bifurcated conformations for both ELVd strands. In vitro examination by non-denaturing PAGE showed that they migrate as prominent single bands, with the ELVd (+) RNA displaying a more compact conformation as revealed by its faster electrophoretic mobility. In vitro SHAPE analysis corroborated the ELVd conformations derived from thermodynamics-based predictions in silico. Moreover, sequence analysis of 94 full-length natural ELVd variants disclosed co-variations, and mutations converting canonical into wobble pairs or vice versa, which confirmed in vivo most of the stems predicted in silico and in vitro, and additionally helped to introduce minor structural refinements. Therefore, results from the 3 experimental approaches were essentially consistent among themselves. Application to RNA preparations from ELVd-infected tissue of RNA ligase-mediated rapid amplification of cDNA ends, combined with pretreatments to modify the 5' ends of viroid strands, mapped the transcription initiation sites of ELVd (+) and (-) strands in vivo at different sequence/structural motifs, in contrast with the situation previously observed in 2 other members of the family Avsunviroidae.

Coupling of RNA Polymerase II Transcription Elongation with Pre-mRNA Splicing

Pre-mRNA maturation frequently occurs at the same time and place as transcription by RNA polymerase II. The co-transcriptionality of mRNA processing has permitted the evolution of mechanisms that functionally couple transcription elongation with diverse events that occur on the nascent RNA. This review summarizes the current understanding of the relationship between transcriptional elongation through a chromatin template and co-transcriptional splicing including alternative splicing decisions that affect the expression of most human genes.

Structure and Associated Biological Functions of Viroids

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.

Dual masking of specific negative splicing regulatory elements resulted in maximal exon 7 inclusion of SMN2 gene

Spinal muscular atrophy (SMA) is a fatal autosomal recessive disease caused by survival motor neuron (SMN) protein insufficiency due to SMN1 mutations. Boosting SMN2 expression is a potential therapy for SMA. SMN2 has identical coding sequence as SMN1 except for a silent C-to-T transition at the 6th nucleotide of exon 7, converting a splicing enhancer to a silencer motif. Consequently, most SMN2 transcripts lack exon 7. More than ten putative splicing regulatory elements (SREs) were reported to regulate exon 7 splicing. To investigate the relative strength of each negative SRE in inhibiting exon 7 inclusion, antisense oligonucleotides (AONs) were used to mask each element, and the fold increase of full-length SMN transcripts containing exon 7 were compared. The most potent negative SREs are at intron 7 (in descending order): ISS-N1, 3' splice site of exon 8 (ex8 3'ss) and ISS+100. Dual-targeting AONs were subsequently used to mask two nonadjacent SREs simultaneously. Notably, masking of both ISS-N1 and ex8 3'ss induced the highest fold increase of full-length SMN transcripts and proteins. Therefore, efforts should be directed towards the two elements simultaneously for the development of optimal AONs for SMA therapy.

On the importance of cotranscriptional RNA structure formation

The expression of genes, both coding and noncoding, can be significantly influenced by RNA structural features of their corresponding transcripts. There is by now mounting experimental and some theoretical evidence that structure formation in vivo starts during transcription and that this cotranscriptional folding determines the functional RNA structural features that are being formed. Several decades of research in bioinformatics have resulted in a wide range of computational methods for predicting RNA secondary structures. Almost all state-of-the-art methods in terms of prediction accuracy, however, completely ignore the process of structure formation and focus exclusively on the final RNA structure. This review hopes to bridge this gap. We summarize the existing evidence for cotranscriptional folding and then review the different, currently used strategies for RNA secondary-structure prediction. Finally, we propose a range of ideas on how state-of-the-art methods could be potentially improved by explicitly capturing the process of cotranscriptional structure formation.

Transient RNA structure features are evolutionarily conserved and can be computationally predicted

Functional RNA structures tend to be conserved during evolution. This finding is, for example, exploited by comparative methods for RNA secondary structure prediction that currently provide the state-of-art in terms of prediction accuracy. We here provide strong evidence that homologous RNA genes not only fold into similar final RNA structures, but that their folding pathways also share common transient structural features that have been evolutionarily conserved. For this, we compile and investigate a non-redundant data set of 32 sequences with known transient and final RNA secondary structures and devise a dedicated computational analysis pipeline.

Maximum expected accuracy structural neighbors of an RNA secondary structure

Background

Since RNA molecules regulate genes and control alternative splicing by allostery, it is important to develop algorithms to predict RNA conformational switches. Some tools, such as paRNAss, RNAshapes and RNAbor, can be used to predict potential conformational switches; nevertheless, no existent tool can detect general (i.e., not family specific) entire riboswitches (both aptamer and expression platform) with accuracy. Thus, the development of additional algorithms to detect conformational switches seems important, especially since the difference in free energy between the two metastable secondary structures may be as large as 15-20 kcal/mol. It has recently emerged that RNA secondary structure can be more accurately predicted by computing the maximum expected accuracy (MEA) structure, rather than the minimum free energy (MFE) structure.

Results

Given an arbitrary RNA secondary structure S₀ for an RNA nucleotide sequence a = a₁,..., a_n, we say that another secondary structure S of a is a k-neighbor of S₀, if the base pair distance between S₀ and S is k. In this paper, we prove that the Boltzmann probability of all k-neighbors of the minimum free energy structure S₀ can be approximated with accuracy ε and confidence 1 - p, simultaneously for all 0 ≤ k < K, by a relative frequency count over N sampled structures, provided that N>N(ε,p,K)=Φ-1p2K24ε2, where Φ(z) is the cumulative distribution function (CDF) for the standard normal distribution. We go on to describe the algorithm RNAborMEA, which for an arbitrary initial structure S₀ and for all values 0 ≤ k < K, computes the secondary structure MEA(k), having maximum expected accuracy over all k-neighbors of S₀. Computation time is O(n³ · K²), and memory requirements are O(n² · K). We analyze a sample TPP riboswitch, and apply our algorithm to the class of purine riboswitches.

Conclusions

The approximation of RNAbor by sampling, with rigorous bound on accuracy, together with the computation of maximum expected accuracy k-neighbors by RNAborMEA, provide additional tools toward conformational switch detection. Results from RNAborMEA are quite distinct from other tools, such as RNAbor, RNAshapes and paRNAss, hence may provide orthogonal information when looking for suboptimal structures or conformational switches. Source code for RNAborMEA can be downloaded from http://sourceforge.net/projects/rnabormea/ or http://bioinformatics.bc.edu/clotelab/RNAborMEA/.

Characterization of plant miRNAs and small RNAs derived from potato spindle tuber viroid (PSTVd) in infected tomato

To defend against invading pathogens, plants possess RNA silencing mechanisms involving small RNAs (miRNAs, siRNAs). Also viroids - plant infectious, non-coding, unencapsidated RNA - cause the production of viroid-specific small RNAs (vsRNA), but viroids do escape the cytoplasmic silencing mechanism. Viroids with minor sequence variations can produce different symptoms in infected plants, suggesting an involvement of vsRNAs in symptom production. We analyzed by deep sequencing the spectrum of vsRNAs induced by the PSTVd strain AS1, which causes strong symptoms such as dwarfing and necrosis upon infection of tomato plants cv Rutgers. Indeed, vsRNAs found with highest frequency mapped to the pathogenicity-modulating domain of PSTVd, supporting an involvement of vsRNAs in symptom production. Furthermore, in PSTVd AS1-infected plants the accumulation of some endogenous miRNAs, which are involved in leaf development via regulation of transcription factors, is suppressed. The latter finding supports the hypothesis that a miRNA-dependent (mis)regulation of transcription factors causes the viroid symptoms.

TRANSAT-- method for detecting the conserved helices of functional RNA structures, including transient, pseudo-knotted and alternative structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment.

Many non-coding genes exert their function via an RNA structure which starts emerging while the RNA sequence is being transcribed from the genome. The resulting folding pathway is known to depend on a variety of features such as the transcription speed, the concentration of various ions and the binding of proteins and other molecules. Not all of these influences can be adequately captured by the existing computational methods which try to replicate what happens in vivo. So far, it has been challenging to experimentally investigate co-transcriptional folding pathways in vivo and only little data from in vitro experiments exists. In order to investigate if functionally similar RNA sequences from different organisms fold in a similar way, we have developed a new computational method, called Transat, which does not require the detailed computational modeling of the cellular environment. We show in a comprehensive analysis that our method is capable of detecting known structural features and provide evidence that structural features of the in vivo folding pathways have been conserved for several biologically interesting classes of RNA sequences.

The RNA strands of the plus and minus polarities of peach latent mosaic viroid fold into different structures

It is believed that peach latent mosaic viroid (PLMVd) strands of both the plus and minus polarities fold into similar secondary and tertiary structures. In order to verify this hypothesis, the behavior of both strands in three biophysical assays was examined. PLMVd transcripts of plus and minus polarity were found to exhibit distinct electrophoretic mobility properties under native conditions, to precipitate differently in the presence of lithium chloride, and to possess variable thermal denaturation profiles. Subsequently, the structure of PLMVd transcripts of minus polarity was elucidated by biochemical methods, thereby permitting comparison to the known structure of the plus polarity. Specifically, enzymatic probing, electrophoretic mobility shift assay, and ribonuclease H hydrolysis were performed in order to resolve the secondary structure of the minus polarity. The left domains of the strands of both polarities appear to be similar, while the right domain exhibited several differences even though they both adopted a branched structure. The pseudoknot P8 formed in the plus strand seemed not formed in the minus strands. The structural differences between the two polarities might have important implications in various steps of the PLMVd life cycle.

Dynamic extended folding: modeling the RNA secondary structures during co-transcriptional folding

For RNA secondary structure prediction, it is an important issue that how to deal with co-transcriptional folding during the RNA synthesis in the cell. On one hand, co-transcriptional folding, leads to the correct final structure of the whole RNA molecule. On the other hand, it may form the recognition sites for the progress of the transcription. Considering the hurdles in the experimental determination of RNA folding structures, we proposed a so-called "dynamic extended folding simulation" approach. We used two human pre-mRNA samples, the first functional alpha-gene HBZ and the fifth beta-gene HBB, to "display" the co-transcriptional folding images in detail. The modeling process starts from the prediction of a 30-nucleotide (nt) sequence, then in each update 30 nts was extended, say, 1-30, 1-60, 1-90, 1-120,..., 1-1651 nts (for HBB, 1-1606 nts). We selected the RNAstructure program to predict the folding secondary structures of all the segments. We defined "hairpin" as the unit of the secondary structure and analyzed the states of such unit during the sequential dynamic extended folding processes. We found that some hairpins are "conserved", i.e., after its appearance, it always is there in the followed foldings. Some hairpins present partially in the folding segments, and some hairpins appear for only once or twice. This phenomenon vividly depicts the generation and adjusting of the temporal structural units during the co-transcriptional folding process. It is these "hairpins" that support the thermodynamically stable structure at the end of the RNA synthesis. They may also play a role in RNA splicing process and even in the folding structure of the synthesized protein.

The jerky and knotty dynamics of RNA

RNA is known to exhibit a jerky dynamics, as intramolecular thermal motion, on <0.1 micros time scales, is punctuated by infrequent structural rearrangements on much longer time scales, i.e. from >10 micros up to a few minutes or even hours. These rare stochastic events correspond to the formation or dissociation of entire stems through cooperative base pairing/unpairing transitions. Such a clear separation of time scales in RNA dynamics has made it possible to implement coarse grained RNA simulations, which predict RNA folding and unfolding pathways including kinetically trapped structures on biologically relevant time scales of seconds to minutes. RNA folding simulations also enable to predict the formation of pseudoknots, that is, helices interior to loops, which mechanically restrain the relative orientations of other non-nested helices. But beyond static structural constraints, pseudoknots can also strongly affect the folding and unfolding dynamics of RNA, as the order by which successive helices are formed and dissociated can lead to topologically blocked transition intermediates. The resulting knotty dynamics can enhance the stability of RNA switches, improve the efficacy of co-transcriptional folding pathways and lead to unusual self-assembly properties of RNA.

The fraction of RNA that folds into the correct branched secondary structure determines hepatitis delta virus type 3 RNA editing levels

RNA editing by the host RNA adenosine deaminase ADAR1 at the amber/W site of hepatitis delta virus RNA plays a central role in the viral replication cycle by affecting the balance between viral RNA synthesis and packaging. Previously, we found that HDV genotype III (HDV-3) RNA can form two secondary structures following transcription: an unbranched rod structure, which is characteristic of HDV, and a metastable branched structure that serves as the substrate for editing. The unstable nature of the branched editing substrate structure raised the possibility that structural dynamics of the RNA following transcription could determine the rate at which editing occurs. Here, editing and its control are examined in two HDV-3 isolates, from Peru and Ecuador. Analysis of editing in vitro by ADAR1 indicated that the branched structure formed by RNA derived from the Peruvian isolate is edited more efficiently than that from the Ecuadorian isolate. In contrast, in the context of replication, Peruvian RNA is edited less efficiently than RNA containing Ecuadorian sequences. Computational analyses of RNA folding using the massively parallel genetic algorithm (MPGAfold) indicated that the Peruvian RNA is less likely to form the branched structure required for editing than the Ecuadorian isolate. This difference was confirmed by in vitro transcription of these RNAs. Overall, our data indicate that HDV-3 controls RNA editing levels via (1) the fraction of the RNA that folds, during transcription, into the metastable branched structure required for editing and (2) the efficiency with which ADAR1 edits this branched substrate RNA.

RNA folding during transcription: protocols and studies

RNA folds during transcription in the cell. Compared to most in vitro studies where the focus is generally on Mg(2+)-initiated refolding of fully synthesized transcripts, cotranscriptional RNA folding studies better replicate how RNA folds in a cellular environment. Unique aspects of cotranscriptional folding include the 5'- to 3'-polarity of RNA, the transcriptional speed, pausing properties of the RNA polymerase, the effect of the transcriptional complex and associated factors, and the effect of RNA-binding proteins. Identifying strategic pause sites can reveal insights on the folding pathway of the nascent transcript. Structural mapping of the paused transcription complexes identifies important folding intermediates along these pathways. Oligohybridization assays and the appearance of the catalytic activity of a ribozyme either in trans or in cis can be used to monitor cotranscriptional folding under a wide range of conditions. In our laboratory, these methodologies have been applied to study the folding of three highly conserved RNAs (RNase P, SRP, and tmRNA), several circularly permuted forms of a bacterial RNase P RNA, a riboswitch (thiM), and an aptamer-activated ribozyme (glmS).

Nuclear targeting by fragmentation of the potato spindle tuber viroid genome

Transient expression of engineered reporter RNAs encoding an intron-containing green fluorescent protein (GFP) from a Potato virus X-based expression vector previously demonstrated the nuclear targeting capability of the 359 nucleotide Potato spindle tuber viroid (PSTVd) RNA genome. To further delimit the putative nuclear-targeting signal, PSTVd subgenomic fragments were embedded within the intron, and recombinant reporter RNAs were inoculated onto Nicotiana benthamiana plants. Appearance of green fluorescence in leaf tissue inoculated with PSTVd-fragment-containing constructs indicated shuttling of the RNA into the nucleus by fragments as short as 80 nucleotides in length. Plant-to-plant variation in the timing of intron removal and subsequent GFP fluorescence was observed; however, earliest and most abundant GFP expression was obtained with constructs containing the conserved hairpin I palindrome structure and embedded upper central conserved region. Our results suggest that this conserved sequence and/or the stem-loop structure it forms is sufficient for import of PSTVd into the nucleus.

Dynamics of co-transcriptional pre-mRNA folding influences the induction of dystrophin exon skipping by antisense oligonucleotides

Antisense oligonucleotides (AONs) mediated exon skipping offers potential therapy for Duchenne muscular dystrophy. However, the identification of effective AON target sites remains unsatisfactory for lack of a precise method to predict their binding accessibility. This study demonstrates the importance of co-transcriptional pre-mRNA folding in determining the accessibility of AON target sites for AON induction of selective exon skipping in DMD. Because transcription and splicing occur in tandem, AONs must bind to their target sites before splicing factors. Furthermore, co-transcriptional pre-mRNA folding forms transient secondary structures, which redistributes accessible binding sites. In our analysis, to approximate transcription elongation, a “window of analysis” that included the entire targeted exon was shifted one nucleotide at a time along the pre-mRNA. Possible co-transcriptional secondary structures were predicted using the sequence in each step of transcriptional analysis. A nucleotide was considered “engaged” if it formed a complementary base pairing in all predicted secondary structures of a particular step. Correlation of frequency and localisation of engaged nucleotides in AON target sites accounted for the performance (efficacy and efficiency) of 94% of 176 previously reported AONs. Four novel insights are inferred: (1) the lowest frequencies of engaged nucleotides are associated with the most efficient AONs; (2) engaged nucleotides at 3′ or 5′ ends of the target site attenuate AON performance more than at other sites; (3) the performance of longer AONs is less attenuated by engaged nucleotides at 3′ or 5′ ends of the target site compared to shorter AONs; (4) engaged nucleotides at 3′ end of a short target site attenuates AON efficiency more than at 5′ end.

A genomic map of viroid RNA motifs critical for replication and systemic trafficking

RNA replication and systemic trafficking play significant roles in developmental regulation and host-pathogen interactions. Viroids are the simplest noncoding eukaryotic RNA pathogens and genetic units that are capable of autonomous replication and systemic trafficking and offer excellent models to investigate the role of RNA structures in these processes. Like other RNAs, the predicted secondary structure of a viroid RNA contains many loops and bulges flanked by double-stranded helices, the biological functions of which are mostly unknown. Using Potato spindle tuber viroid infection of Nicotiana benthamiana as the experimental system, we tested the hypothesis that these loops/bulges are functional motifs that regulate replication in single cells or trafficking in a plant. Through a genome-wide mutational analysis, we identified multiple loops/bulges essential or important for each of these biological processes. Our results led to a genomic map of viroid RNA motifs that mediate single-cell replication and systemic trafficking, respectively. This map provides a framework to enable high-throughput studies on the tertiary structures and functional mechanisms of RNA motifs that regulate viroid replication and trafficking. Our model and approach should also be valuable for comprehensive investigations of the replication and trafficking motifs in other RNAs.

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

Encoding folding paths of RNA switches

RNA co-transcriptional folding has long been suspected to play an active role in helping proper native folding of ribozymes and structured regulatory motifs in mRNA untranslated regions (UTRs). Yet, the underlying mechanisms and coding requirements for efficient co-transcriptional folding remain unclear. Traditional approaches have intrinsic limitations to dissect RNA folding paths, as they rely on sequence mutations or circular permutations that typically perturb both RNA folding paths and equilibrium structures. Here, we show that exploiting sequence symmetries instead of mutations can circumvent this problem by essentially decoupling folding paths from equilibrium structures of designed RNA sequences. Using bistable RNA switches with symmetrical helices conserved under sequence reversal, we demonstrate experimentally that native and transiently formed helices can guide efficient co-transcriptional folding into either long-lived structure of these RNA switches. Their folding path is controlled by the order of helix nucleations and subsequent exchanges during transcription, and may also be redirected by transient antisense interactions. Hence, transient intra- and inter-molecular base pair interactions can effectively regulate the folding of nascent RNA molecules into different native structures, provided limited coding requirements, as discussed from an information theory perspective. This constitutive coupling between RNA synthesis and RNA folding regulation may have enabled the early emergence of autonomous RNA-based regulation networks.

Transcription of potato spindle tuber viroid by RNA polymerase II starts in the left terminal loop

Viroids are single-stranded, circular RNAs of 250 to 400 bases, that replicate autonomously in their host plants but do not code for a protein. Viroids of the family Pospiviroidae, of which potato spindle tuber viroid (PSTVd) is the type strain, are replicated by the host's DNA-dependent RNA polymerase II in the nucleus. To analyze the initiation site of transcription from the (+)-stranded circles into (-)-stranded replication intermediates, we used a nuclear extract from a non-infected cell culture of the host plant S. tuberosum. The (-)-strands, which were de novo-synthesized in the extract upon addition of circular (+)-PSTVd, were purified by affinity chromatography. This purification avoided contamination by host nucleic acids that had resulted in a misassignment of the start site in an earlier study. Primer-extension analysis of the de novo-synthesized (-)-strands revealed a single start site located in the hairpin loop of the left terminal region in circular PSTVd's secondary structure. This start site is supported further by analysis of the infectivity and replication behavior of site-directed mutants in planta.

Memory efficient folding algorithms for circular RNA secondary structures

BACKGROUND

A small class of RNA molecules, in particular the tiny genomes of viroids, are circular. Yet most structure prediction algorithms handle only linear RNAs. The most straightforward approach is to compute circular structures from 'internal' and 'external' substructures separated by a base pair. This is incompatible, however, with the memory-saving approach of the Vienna RNA Package which builds a linear RNA structure from shorter (internal) structures only.

RESULT

Here we describe how circular secondary structures can be obtained without additional memory requirements as a kind of 'post-processing' of the linear structures.

AVAILABILITY

The circular folding algorithm is implemented in the current version of the of RNAfold program of the Vienna RNA Package, which can be downloaded from http://www.tbi.univie.ac.at/RNA/

Kinetics and thermodynamics make different contributions to RNA folding in vitro and in yeast

RNAs somehow adopt specific functional structures despite the capacity to form alternative nonfunctional structures with similar stabilities. We analyzed RNA assembly during transcription in vitro and in yeast using hairpin ribozyme self-cleavage to assess partitioning between functional ribozyme structures and nonfunctional stem loops. Complementary insertions located upstream of the ribozyme inhibited ribozyme assembly more than downstream insertions during transcription in vitro, consistent with a sequential folding model in which the outcome is determined by the structure that forms first. In contrast, both upstream and downstream insertions strongly inhibited assembly of the same ribozyme variants when expressed as chimeric mRNAs in yeast, indicating that inhibitory stem loops can form even after the entire ribozyme sequence has been transcribed. Evidently, some feature unique to the intracellular environment modulates the influence of transcription polarity and enhances the contribution of thermodynamic stability to RNA folding in vivo.

The robustness of naturally and artificially selected nucleic acid secondary structures

Thermodynamic stability and mutational robustness of secondary structure are critical to the function and evolutionary longevity of RNA molecules. We hypothesize that natural and artificial selection for functional molecules favors the formation of structures that are stable to both thermal and mutational perturbation. There is little direct evidence, however, that functional RNA molecules have been selected for their stability. Here we use thermodynamic secondary structure prediction algorithms to compare the thermal and mutational robustness of over 1000 naturally and artificially evolved molecules. Although we find evidence for the evolution of both types of stability in both sets of molecules, the naturally evolved functional RNA molecules were significantly more stable than those selected in vitro, and artificially evolved catalysts (ribozymes) were more stable than artificially evolved binding species (aptamers). The thermostability of RNA molecules bred in the laboratory is probably not constrained by a lack of suitable variation in the sequence pool but, rather, by intrinsic biases in the selection process.

Co-transcriptional folding is encoded within RNA genes

Background

Most of the existing RNA structure prediction programs fold a completely synthesized RNA molecule. However, within the cell, RNA molecules emerge sequentially during the directed process of transcription. Dedicated experiments with individual RNA molecules have shown that RNA folds while it is being transcribed and that its correct folding can also depend on the proper speed of transcription.

Methods

The main aim of this work is to study if and how co-transcriptional folding is encoded within the primary and secondary structure of RNA genes. In order to achieve this, we study the known primary and secondary structures of a comprehensive data set of 361 RNA genes as well as a set of 48 RNA sequences that are known to differ from the originally transcribed sequence units. We detect co-transcriptional folding by defining two measures of directedness which quantify the extend of asymmetry between alternative helices that lie 5' and those that lie 3' of the known helices with which they compete.

Results

We show with statistical significance that co-transcriptional folding strongly influences RNA sequences in two ways: (1) alternative helices that would compete with the formation of the functional structure during co-transcriptional folding are suppressed and (2) the formation of transient structures which may serve as guidelines for the co-transcriptional folding pathway is encouraged.

Conclusions

These findings have a number of implications for RNA secondary structure prediction methods and the detection of RNA genes.

Cis and trans requirements for rolling circle replication of a satellite RNA

Satellite RNAs usurp the replication machinery of their helper viruses, even though they bear little or no sequence similarity to the helper virus RNA. In Cereal yellow dwarf polerovirus serotype RPV (CYDV-RPV), the 322-nucleotide satellite RNA (satRPV RNA) accumulates to high levels in the presence of the CYDV-RPV helper virus. Rolling circle replication generates multimeric satRPV RNAs that self-cleave via a double-hammerhead ribozyme structure. Alternative folding inhibits formation of a hammerhead in monomeric satRPV RNA. Here we determine helper virus requirements and the effects of mutations and deletions in satRPV RNA on its replication in oat cells. Using in vivo selection of a satRPV RNA pool randomized at specific bases, we found that disruption of the base pairing necessary to form the non-self-cleaving conformation reduced satRPV RNA accumulation. Unlike other satellite RNAs, both the plus and minus strands proved to be equally infectious. Accordingly, very similar essential replication structures were identified in each strand. A different region is required only for encapsidation. The CYDV-RPV RNA-dependent RNA polymerase (open reading frames 1 and 2), when expressed from the nonhelper Barley yellow dwarf luteovirus, was capable of replicating satRPV RNA. Thus, the helper virus's polymerase is the sole determinant of the ability of a virus to replicate a rolling circle satellite RNA. We present a framework for functional domains in satRPV RNA with three types of function: (i) conformational control elements comprising an RNA switch, (ii) self-functional elements (hammerhead ribozymes), and (iii) cis-acting elements that interact with viral proteins.

Structural characterization of the 69 nucleotide potato spindle tuber viroid left-terminal domain by NMR and thermodynamic analysis

The 69 nucleotide left-terminal domain (T(L)) of the potato spindle tuber RNA viroid (PSTVd) constitutes one of its five structural elements. Due to a twofold complementary sequence repeat, two possible conformations are proposed for the T(L) secondary structure; an elongated-rod and a bifurcated form. In the present study, two T(L) mutants were designed that remove the symmetry of the sequence repeats and ensure that either the bifurcated or the elongated-rod conformation is thermodynamically favored. Imino 1H and 15N resonances were assigned for both mutants and the native T(L) domain based on 1H-1H NOESY and heteronuclear 1H-15N HSQC high-resolution NMR spectra. The NMR secondary structure analysis of all constructs establishes unambiguously the elongated-rod form as the secondary structure of the native T(L) domain. Temperature-gradient gel electrophoresis and UV melting experiments corroborate these results. A combined secondary structure and sequence analysis of T(L) domains of other Pospiviroidae family members indicates that the elongated-rod form is thermodynamically favored for the vast majority of these viroids.

Effect of transcription on folding of the Tetrahymena ribozyme

Sequential formation of RNA interactions during transcription can bias the folding pathway and ultimately determine the functional state of a transcript. The kinetics of cotranscriptional folding of the Tetrahymena L-21 ribozyme was compared with refolding of full-length transcripts under the same conditions. Sequential folding after transcription by phage T7 or Escherichia coli polymerase is only twice as fast as refolding, and the yield of native RNA is the same. By contrast, a greater fraction of circularly permuted variants folded correctly at early times during transcription than during refolding. Hybridization of complementary oligonucleotides suggests that cotranscriptional folding enables a permuted RNA beginning at G303 to escape non-native interactions in P3 and P9. We propose that base pairing of upstream sequences during transcription elongation favors branched secondary structures that increase the probability of forming the native ribozyme structure.

A mini-RNA containing the tetraloop, wobble-pair and loop E motifs of the central conserved region of potato spindle tuber viroid is processed into a minicircle

A Mini-RNA from potato spindle tuber viroid (PSTVd) was constructed specifically for cleavage and ligation to circles in vitro. It contains the C-domain with the so-called central conserved region (CCR) of PSTVd with a 17 nt duplication in the upper strand and hairpin structures at the left and rights ends of the secondary structure. The CCR was previously shown to be essential for processing of in vitro transcripts. When folded under conditions which favor formation of a kinetically controlled conformation and incubated in a potato nuclear extract, the Mini-RNA is cleaved correctly at the 5'- and the 3'-end and ligated to a circle. Thus, the CCR obviously contains all structural and functional requirements for correct processing and therefore may be regarded as 'processing domain' of PSTVd. Using the Mini-RNA as a model substrate, the structural and functional relevance of its conserved non-canonical motifs GAAA tetraloop, loop E and G:U wobble base pair were studied by mutational analysis. It was found that (i) the conserved GAAA tetraloop is essential for processing by favoring the kinetically controlled conformation, (ii) a G:U wobble base pair at the 5'-cleavage site contributes to its correct recognition and (iii) an unpaired nucleotide in loop E, which is different from the corresponding nucleotide in the conserved loop E motif, is essential for ligation of the 5'- with the 3'-end. Hence all three structural motifs are functional elements for processing in a potato nuclear extract.

Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy

We investigate 25-34 nucleotide RNA sequences, that have been rationally designed to adopt two different secondary structures that are in thermodynamic equilibrium. Experimental evidence for the co-existence of the two conformers results from the NH...N 1H NMR spectra. When compared to the NH...N 1H NMR spectra of appropriate reference sequences the equilibrium position is easily quantifiable even without the assignment of the individual NH resonances. The reference sequences represent several Watson-Crick base-paired double helical segments, each encountered in either of the two conformers of the bistable target sequence. In addition, we rationalize the influence of nucleotide mutations on the equilibrium position of one of the bistable RNA sequences. The approach further allows a detailed thermodynamic analysis and the evaluation of secondary structure predictions for multistable RNAs obtained by computational methods.

Replication of Potato spindle tuber viroid in cultured cells of tobacco and Nicotiana benthamiana: the role of specific nucleotides in determining replication levels for host adaptation

We have developed an electroporation protocol to inoculate cultured cells of tobacco and Nicotiana benthamiana with in vitro transcripts of Potato spindle tuber viroid (PSTVd) to characterize viroid structural features that determine replication efficiency at the cellular level. Both (+)- and (-)-strands of PSTVd were detected by Northern blots as early as 6 h postinoculation (h.p.i.). Accumulation of the (+)-circular PSTVd increased very rapidly starting at 24 h.p.i. and continued beyond 6 days postinoculation. Viroid accumulation in individual cells was visualized by in situ hybridization, which showed that 60-70% of the cells were infected. Previous work showed that C259 --> U substitution converted tomato isolate PSTVd(KF440-2) into a strain that is infectious on tobacco (M. Wassenegger, R. L. Spieker, S. Thalmeir, F.-U. Gast, L. Riedel, and H. L. Sänger, 1996. Virology 226, 191-197). Similarly, C259 --> U or U257 --> A substitution in the Intermediate strain (PSTVd(Int)) conferred infectivity in tobacco (Y. Zhu, Y. Qi, Y. Xun, R. Owens, and B. Ding, 2002. Plant Physiol. 130, 138-146). Our replication assays in tobacco-cultured cells demonstrated that U257 --> A and C259 --> U substitutions each enhanced PSTVd replication by 5- to 10-fold. Replacement of U257 with C, but not with G, also led to enhanced replication in tobacco cells. Replacement of C259 with nucleotide A or G did not enhance replication. Elevated accumulation of the (-)- and (+)-strands of these mutants was in part due to enhanced transcription. Interestingly, all of the nucleotide changes did not alter PSTVd replication levels in N. benthamiana cells. These results provide insights about PSTVd structures that modulate replication efficiency in adapting to a specific host.