RNAcanvas: interactive drawing and exploration of nucleic acid structures

R2DT: a comprehensive platform for visualizing RNA secondary structure

RNA secondary (2D) structure visualization is an essential tool for understanding RNA function. R2DT is a software package designed to visualize RNA 2D structures in consistent, recognizable, and reproducible layouts. The latest release, R2DT 2.0, introduces multiple significant features, including the ability to display position-specific information, such as single nucleotide polymorphisms or SHAPE reactivities. It also offers a new template-free mode allowing visualization of RNAs without pre-existing templates, alongside a constrained folding mode and support for animated visualizations. Users can interactively modify R2DT diagrams, either manually or using natural language prompts, to generate new templates or create publication-quality images. Additionally, R2DT features faster performance, an expanded template library, and a growing collection of compatible tools and utilities. Already integrated into multiple biological databases, R2DT has evolved into a comprehensive platform for RNA 2D visualization, accessible at https://r2dt.bio.

Selective deuteration of an RNA:RNA complex for structural analysis using small-angle scattering

The structures of RNA:RNA complexes regulate many biological processes. Despite their importance, protein-free RNA:RNA complexes represent a tiny fraction of experimentally determined structures. Here, we describe a joint small-angle X-ray and neutron scattering (SAXS/SANS) approach to structurally interrogate conformational changes in a model RNA:RNA complex. Using SAXS, we measured the solution structures of the individual RNAs and of the overall RNA:RNA complex. With SANS, we demonstrate, as a proof of principle, that isotope labeling and contrast matching (CM) can be combined to probe the bound state structure of an RNA within a selectively deuterated RNA:RNA complex. Furthermore, we show that experimental scattering data can validate and improve predicted AlphaFold 3 RNA:RNA complex structures to reflect its solution structure. Our work demonstrates that in silico modeling, SAXS, and CM-SANS can be used in concert to directly analyze conformational changes within RNAs when in complex, enhancing our understanding of RNA structure in functional assemblies.

A Cascade of Conformational Switches in SARS-CoV-2 Frameshifting: Coregulation by Upstream and Downstream Elements

Targeting ribosomal frameshifting has emerged as a potential therapeutic intervention strategy against COVID-19. In this process, a -1 shift in the ribosomal reading frame encodes alternative viral proteins. Any interference with this process profoundly affects viral replication and propagation. For SARS-CoV-2, two RNA sites associated with ribosomal frameshifting are positioned on the 5' and 3' of the frameshifting residues. Although much attention has been focused on the 3' frameshift element (FSE), the 5' stem-loop (attenuator hairpin, AH) can play a role. Yet the relationship between the two regions is unknown. In addition, multiple folds of the FSE and FSE-containing RNA regions have been discovered. To gain more insight into these RNA folds in the larger sequence context that includes AH, we apply our graph-theory-based modeling tools to represent RNA secondary structures, "RAG" (RNA-As-Graphs), to generate conformational landscapes that suggest length-dependent conformational distributions. We show that the AH region can coexist as a stem-loop with main and alternative 3-stem pseudoknots of the FSE (dual graphs 3_6 and 3_3 in our notation) but that an alternative stem 1 (AS1) can disrupt the FSE pseudoknots and trigger other folds. A critical length for AS1 of 10-bp regulates key folding transitions. Together with designed mutants and available experimental data, we present a sequential view of length-dependent folds during frameshifting and suggest their mechanistic roles. These structural and mutational insights into both ends of the FSE advance our understanding of the SARS-CoV-2 frameshifting mechanism by suggesting how alternative folds play a role in frameshifting and defining potential therapeutic intervention techniques that target specific folds.

Role of the Psi Packaging Signal and Dimerization Initiation Sequence in the Organization of Rous Sarcoma Virus Gag-gRNA Co-Condensates

Retroviral genome selection and virion assembly remain promising targets for novel therapeutic intervention. Recent studies have demonstrated that the Gag proteins of Rous sarcoma virus (RSV) and human immunodeficiency virus type-1 (HIV-1) undergo nuclear trafficking, colocalize with nascent genomic viral RNA (gRNA) at transcription sites, may interact with host transcription factors, and display biophysical properties characteristic of biomolecular condensates. In the present work, we utilized a controlled in vitro condensate assay and advanced imaging approaches to investigate the effects of interactions between RSV Gag condensates and viral and nonviral RNAs on condensate abundance and organization. We observed that the psi (Ψ) packaging signal and the dimerization initiation sequence (DIS) had stabilizing effects on RSV Gag condensates, while RNAs lacking these features promoted or antagonized condensation, depending on local protein concentration and condensate architecture. An RNA containing Ψ, DIS, and the dimerization linkage structure (DLS) that is capable of stable dimer formation was observed to act as a bridge between RSV Gag condensates. These observations suggest additional, condensate-related roles for Gag-Ψ binding, gRNA dimerization, and Gag dimerization/multimerization in gRNA selection and packaging, representing a significant step forward in our understanding of how these interactions collectively facilitate efficient genome packaging.

De novo and inherited dominant variants in U4 and U6 snRNAs cause retinitis pigmentosa

The U4 small nuclear RNA (snRNA) forms a duplex with the U6 snRNA and, together with U5 and ~30 proteins, is part of the U4/U6.U5 tri-snRNP complex, located at the core of the major spliceosome. Recently, recurrent de novo variants in the U4 RNA, transcribed from the RNU4-2 gene, and in at least two other RNU genes were discovered to cause neurodevelopmental disorder. We detected inherited and de novo heterozygous variants in RNU4-2 (n.18_19insA and n.56T>C) and in four out of the five RNU6 paralogues (n.55_56insG and n.56_57insG) in 135 individuals from 62 families with non-syndromic retinitis pigmentosa (RP), a rare form of hereditary blindness. We show that these variants are recurrent among RP families and invariably cluster in close proximity within the three-way junction (between stem-I, the 5' stem-loop and stem-II) of the U4/U6 duplex, affecting its natural conformation. Interestingly, this region binds to numerous splicing factors of the tri-snRNP complex including PRPF3, PRPF8 and PRPF31, previously associated with RP as well. The U4 and U6 variants identified seem to affect snRNP biogenesis, namely the U4/U6 di-snRNP, which is an assembly intermediate of the tri-snRNP. Based on the number of positive cases observed, deleterious variants in RNU4-2 and in RNU6 paralogues could be a significant cause of isolated or dominant RP, accounting for up to 1.2% of all undiagnosed RP cases. This study highlights the role of non-coding genes in rare Mendelian disorders and uncovers pleiotropy in RNU4-2, where different variants underlie neurodevelopmental disorder and RP.

Structural impact of 3-methylcytosine modification on the anticodon stem of a neuronally-enriched arginine tRNA

All tRNAs undergo a series of chemical modifications to fold and function correctly. In mammals, the C32 nucleotide in the anticodon loop of tRNA-Arg-CCU and UCU is methylated to form 3-methylcytosine (m3C). Deficiency of m3C in arginine tRNAs has been linked to human neurodevelopmental disorders, indicating a critical biological role for m3C modification. However, the structural repercussions of m3C modification are not well understood. Here, we examine the structural effects of m3C32 modification on the anticodon stem loop (ASL) of human tRNA-Arg-UCU-4-1, a unique tRNA with enriched expression in the central nervous system. Optical melting experiments demonstrate that m3C modification can locally disrupt nearby base pairing within the ASL while simultaneously stabilizing the ASL electrostatically, resulting in little net change thermodynamically. The isoenergetic nature of the C32 - A38 pair vs the m3C32 - A38 pair may help discriminate against structures not adopting canonical C32 - A38 pairings, as most other m3C pairings are unfavorable. Furthermore, multidimensional NMR reveals that after m3C modification there are changes in hairpin loop structure and dynamics, the structure of A37, and the neighboring A31 - U39 base pair. However, these structural changes after modification are made while maintaining the shape of the C32 - A38 pairing, which is essential for efficient tRNA function in translation. These findings suggests that m3C32 modification could alter interactions of tRNA-Arg isodecoders with one or more binding partners while simultaneously maintaining the tRNA's ability to function in translation.

Electrochemical Aptamer-Based Biosensor for Detecting Pap31, a Biomarker for Carrion's Disease

Carrion's disease, caused by infection with the bacterium Bartonella bacilliformis (B. bacilliformis), is effectively treated with antibiotics, but reaches fatality rates of ~90% if untreated. Current diagnostic methods are limited, insufficiently sensitive, or require laboratory technology unavailable in endemic areas. Electrochemical aptamer-based (E-AB) biosensors provide a potential solution for this unmet need, as these biosensors are portable, sensitive, and can rapidly report the detection of small molecule targets. Here, we developed an E-AB biosensor to detect Pap31, a biomarker of Carrion's disease and an outer membrane protein in B. bacilliformis. We identified an aptamer with Pap31-specific binding affinity using a magnetic pull-down assay with magnetic bead-bound Pap31 and an aptamer library followed by electrophoretic mobility shift assays. We incorporated the Pap31-binding aptamer into a DNA oligonucleotide that changes conformation upon binding Pap31. The resultant Pap31 E-AB biosensor produced a rapid, significant, and target-specific electrical current readout in the buffer, demonstrating an apparent KD of 0.95 nM with a limit of detection of 0.1 nM, and no significant signal change when challenged with off-target proteins. This proof-of-concept Pap31 biosensor design is a first step toward the development of more rapid, sensitive, and portable diagnostic tools for detecting Carrion's disease.

Functional Validation of SAM Riboswitch Element A from Listeria monocytogenes

SreA is one of seven candidate S-adenosyl methionine (SAM) class I riboswitches identified in Listeria monocytogenes, a saprophyte and opportunistic foodborne pathogen. SreA precedes genes encoding a methionine ATP-binding cassette (ABC) transporter, which imports methionine and is presumed to regulate transcription of its downstream genes in a SAM-dependent manner. The proposed role of SreA in controlling the transcription of genes encoding an ABC transporter complex may have important implications for how the bacteria senses and responds to the availability of the metabolite SAM in the diverse environments in which L. monocytogenes persists. Here we validate SreA as a functional SAM-I riboswitch through ligand binding studies, structure characterization, and transcription termination assays. We determined that SreA has both a structure and SAM binding properties similar to those of other well-characterized SAM-I riboswitches. Despite the apparent structural similarities to previously described SAM-I riboswitches, SreA induces transcription termination in response to comparatively lower (nanomolar) ligand concentrations. Furthermore, SreA is a leaky riboswitch that permits some transcription of the downstream gene even in the presence of millimolar SAM, suggesting that L. monocytogenes may "dampen" the expression of genes for methionine import but likely does not turn them "OFF".

Stick-slip unfolding favors self-association of expanded HTT mRNA

In Huntington's Disease (HD) and related disorders, expansion of CAG trinucleotide repeats produces a toxic gain of function in affected neurons. Expanded huntingtin (expHTT) mRNA forms aggregates that sequester essential RNA binding proteins, dysregulating mRNA processing and translation. The physical basis of RNA aggregation has been difficult to disentangle owing to the heterogeneous structure of the CAG repeats. Here, we probe the folding and unfolding pathways of expHTT mRNA using single-molecule force spectroscopy. Whereas normal HTT mRNAs unfold reversibly and cooperatively, expHTT mRNAs with 20 or 40 CAG repeats slip and unravel non-cooperatively at low tension. Slippage of CAG base pairs is punctuated by concerted rearrangement of adjacent CCG trinucleotides, trapping partially folded structures that readily base pair with another RNA strand. We suggest that the conformational entropy of the CAG repeats, combined with stable CCG base pairs, creates a stick-slip behavior that explains the aggregation propensity of expHTT mRNA.

De novo gene synthesis by an antiviral reverse transcriptase

Defense-associated reverse transcriptase (DRT) systems perform DNA synthesis to protect bacteria against viral infection, but the identities and functions of their DNA products remain largely unknown. We show that DRT2 systems encode an unprecedented immune pathway that involves de novo gene synthesis through rolling circle reverse transcription of a noncoding RNA (ncRNA). Programmed template jumping on the ncRNA generates a concatemeric cDNA, which becomes double-stranded upon viral infection. This DNA product constitutes a protein-coding, nearly endless open reading frame (neo) gene whose expression leads to potent cell growth arrest, restricting the viral infection. Our work highlights an elegant expansion of genome coding potential through RNA-templated gene creation and challenges conventional paradigms of genetic information encoded along the one-dimensional axis of genomic DNA.

A specific domain within the 3' untranslated region of Usutu virus confers resistance to the exonuclease ISG20

Usutu virus (USUV) and West Nile virus (WNV) are two closely related emerging mosquito-borne flaviviruses. Their natural hosts are wild birds, but they can also cause severe neurological disorders in humans. Both viruses are efficiently suppressed by type I interferon (IFN), which interferes with viral replication, dissemination, pathogenesis and transmission. Here, we show that the replication of USUV and WNV are inhibited through a common set of IFN-induced genes (ISGs), with the notable exception of ISG20, which USUV is resistant to. Strikingly, USUV was the only virus among all the other tested mosquito-borne flaviviruses that demonstrated resistance to the 3'-5' exonuclease activity of ISG20. Our findings highlight that the intrinsic resistance of the USUV genome, irrespective of the presence of cellular or viral proteins or protective post-transcriptional modifications, relies on a unique sequence present in its 3' untranslated region. Importantly, this genomic region alone can confer ISG20 resistance to a susceptible flavivirus, without compromising its infectivity, suggesting that it could be acquired by other flaviviruses. This study provides new insights into the strategy employed by emerging flaviviruses to overcome host defense mechanisms.

R2DT: A COMPREHENSIVE PLATFORM FOR VISUALISING RNA SECONDARY STRUCTURE

RNA secondary (2D) structure visualisation is an essential tool for understanding RNA function. R2DT is a software package designed to visualise RNA 2D structures in consistent, recognisable, and reproducible layouts. The latest release, R2DT 2.0, introduces multiple significant features, including the ability to display position-specific information, such as single nucleotide polymorphisms (SNPs) or SHAPE reactivities. It also offers a new template-free mode allowing visualisation of RNAs without pre-existing templates, alongside a constrained folding mode and support for animated visualisations. Users can interactively modify R2DT diagrams, either manually or using natural language prompts, to generate new templates or create publication-quality images. Additionally, R2DT features faster performance, an expanded template library, and a growing collection of compatible tools and utilities. Already integrated into multiple biological databases, R2DT has evolved into a comprehensive platform for RNA 2D visualisation, accessible at https://r2dt.bio.

Requirement of the N-terminal region of nonstructural protein 1 in cis for SARS-CoV-2 defective RNA replication

SARS-CoV-2 belongs to the family Coronaviridae and carries a single-stranded positive-sense RNA genome. During coronavirus (CoV) replication, defective or defective interfering RNAs that lack a large portion of the genome often emerge. These defective RNAs typically carry the necessary RNA elements that are required for replication and packaging. We identified the minimum requirement of the 5' proximal region necessary for viral RNA replication by using artificially generated SARS-CoV-2 minigenomes. The minigenomes consist of the 5'-proximal region, an open reading frame (ORF) that encodes a fusion protein consisting of the N-terminal of viral NSP1 and a reporter gene, and the 3' untranslated region of the SARS-CoV-2 genome. We used a modified SARS-CoV-2 variant to support replication of the minigenomes. A minigenome carrying the 5' proximal 634 nucleotides replicated, whereas those carrying shorter than 634 nucleotides did not, demonstrating that the entire 265 nt-long 5' untranslated region and N-terminal portion of the NSP1 coding region are required for the minigenome replication. Minigenome RNAs carrying a specific amino acid substitution or frame shift insertions in the partial NSP1 coding sequence abrogated minigenome replication. Introduction of synonymous mutations in the minigenome RNAs also affected the replication efficiency of the minigenomes. These data suggest that the expression of the N-terminal portion of NSP1 and the primary sequence of the 5' proximal 634 nucleotides are important for minigenome replication.IMPORTANCESARS-CoV-2, the causative agent of COVID-19, is highly transmissible and continues to have a significant impact on public health and the global economy. While several vaccines mitigate the severe consequences of SARS-CoV-2 infection, mutant viruses with reduced reactivity to current vaccines continue to emerge and circulate. This study aimed to identify the minimal 5' proximal region of SARS-CoV-2 genomic RNA required for SARS-CoV-2 defective RNA replication and investigate the importance of an ORF encoded in these defective RNAs. Identifying cis-acting replication signals of SARS-CoV-2 genomic RNA is critical for the development of antivirals that target these signals. Additionally, replication-competent defective RNAs can serve as therapeutic reagents to interfere with SARS-CoV-2 replication. Our findings provide valuable insights into the mechanisms of SARS-CoV-2 RNA replication and the development of reagents that suppress SARS-CoV-2 replication.

Pseudomonas aeruginosa modulates both Caenorhabditis elegans attraction and pathogenesis by regulating nitrogen assimilation

Detecting chemical signals is important for identifying food sources and avoiding harmful agents. Like many animals, C. elegans use olfaction to chemotax towards their main food source, bacteria. However, little is known about the bacterial compounds governing C. elegans attraction to bacteria and the physiological importance of these compounds to bacteria. Here, we address these questions by investigating the function of a small RNA, P11, in the pathogen, Pseudomonas aeruginosa, that was previously shown to mediate learned pathogen avoidance. We discovered that this RNA also affects the attraction of untrained C. elegans to P. aeruginosa and does so by controlling production of ammonia, a volatile odorant produced during nitrogen assimilation. We describe the complex regulation of P. aeruginosa nitrogen assimilation, which is mediated by a partner-switching mechanism involving environmental nitrates, sensor proteins, and P11. In addition to mediating C. elegans attraction, we demonstrate that nitrogen assimilation mutants perturb bacterial fitness and pathogenesis during C. elegans infection by P. aeruginosa. These studies define ammonia as a major mediator of trans-kingdom signaling, implicate nitrogen assimilation as important for both bacteria and host organisms, and highlight how a bacterial metabolic pathway can either benefit or harm a host in different contexts.

Selective deuteration of an RNA:RNA complex for structural analysis using small-angle scattering

The structures of RNA:RNA complexes regulate many biological processes. Despite their importance, protein-free RNA:RNA complexes represent a tiny fraction of experimentally-determined structures. Here, we describe a joint small-angle X-ray and neutron scattering (SAXS/SANS) approach to structurally interrogate conformational changes in a model RNA:RNA complex. Using SAXS, we measured the solution structures of the individual RNAs in their free state and of the overall RNA:RNA complex. With SANS, we demonstrate, as a proof-of-principle, that isotope labeling and contrast matching (CM) can be combined to probe the bound state structure of an RNA within a selectively deuterated RNA:RNA complex. Furthermore, we show that experimental scattering data can validate and improve predicted AlphaFold 3 RNA:RNA complex structures to reflect its solution structure. Our work demonstrates that in silico modeling, SAXS, and CM-SANS can be used in concert to directly analyze conformational changes within RNAs when in complex, enhancing our understanding of RNA structure in functional assemblies.

RiboVision2: A Web Server for Advanced Visualization of Ribosomal RNAs

RiboVision2 is a web server designed to visualize phylogenetic, structural, and evolutionary properties of ribosomal RNAs simultaneously at the levels of primary, secondary, and three-dimensional structure and in the context of full ribosomal complexes. RiboVision2 instantly computes and displays a broad variety of data; it has no login requirements, is open-source, free for all users, and available at https://ribovision2.chemistry.gatech.edu.

Nucleocapsids of the Rift Valley fever virus ambisense S segment contain an exposed RNA element in the center that overlaps with the intergenic region

Rift Valley fever virus (RVFV) is a mosquito-borne zoonotic pathogen. Its RNA genome consists of two negative-sense segments (L and M) with one gene each, and one ambisense segment (S) with two opposing genes separated by the noncoding "intergenic region" (IGR). These vRNAs and the complementary cRNAs are encapsidated by nucleoprotein (N). Using iCLIP2 (individual-nucleotide resolution UV crosslinking and immunoprecipitation) to map all N-vRNA and N-cRNA interactions, we detect N coverage along the L and M segments. However, the S segment vRNA and cRNA each contain approximately 100 non-encapsidated nucleotides stretching from the IGR into the 5'-adjacent reading frame. These exposed regions are RNase-sensitive and predicted to form stem-loop structures with the mRNA transcription termination motif positioned near the top. Moreover, optimal S segment transcription and replication requires the entire exposed region rather than only the IGR. Thus, the RVFV S segment contains a central, non-encapsidated RNA region with a functional role.

Direct nanopore RNA sequencing of umbra-like virus-infected plants reveals long non-coding RNAs, specific cleavage sites, D-RNAs, foldback RNAs, and temporal- and tissue-specific profiles

The traditional view of plus (+)-strand RNA virus transcriptomes is that infected cells contain a limited variety of viral RNAs, such as full-length (+)-strand genomic RNA(s), (-)-strand replication intermediate(s), 3' co-terminal subgenomic RNA(s), and viral recombinant defective (D)-RNAs. To ascertain the full complement of viral RNAs associated with the simplest plant viruses, long-read direct RNA nanopore sequencing was used to perform transcriptomic analyses of two related umbra-like viruses: citrus yellow vein-associated virus (CY1) from citrus and CY2 from hemp. Analysis of different timepoints/tissues in CY1- and CY2-infected Nicotiana benthamiana plants and CY2-infected hemp revealed: (i) three 5' co-terminal RNAs of 281 nt, 442 nt and 671 nt, each generated by a different mechanism; (ii) D-RNA populations containing the 671 fragment at their 5'ends; (iii) many full-length genomic RNAs and D-RNAs with identical 3'end 61 nt truncations; (iv) virtually all (-)-strand reads missing 3 nt at their 3' termini; (v) (±) foldback RNAs comprising about one-third of all (-)-strand reads and (vi) a higher proportion of full-length gRNAs in roots than in leaves, suggesting that roots may be functioning as a gRNA reservoir. These findings suggest that viral transcriptomes are much more complex than previously thought.

Determinants of pegivirus persistence, cross-species infection, and adaptation in the laboratory mouse

Viruses capable of causing persistent infection have developed sophisticated mechanisms for evading host immunity, and understanding these processes can reveal novel features of the host immune system. One such virus, human pegivirus (HPgV), infects ~15% of the global human population, but little is known about its biology beyond the fact that it does not cause overt disease. We passaged a pegivirus isolate of feral brown rats (RPgV) in immunodeficient laboratory mice to develop a mouse-adapted virus (maPgV) that established persistent high-titer infection in a majority of wild-type laboratory mice. maRPgV viremia was detected in the blood of mice for >300 days without apparent disease, closely recapitulating the hallmarks of HPgV infection in humans. We found a pro-viral role for type-I interferon in chronic infection; a lack of PD-1-mediated tolerance to PgV infection; and multiple mechanisms by which PgV immunity can be achieved by an immunocompetent host. These data indicate that the PgV immune evasion strategy has aspects that are both common and unique among persistent viral infections. The creation of maPgV represents the first PgV infection model in wild-type mice, thus opening the entire toolkit of the mouse host to enable further investigation of this persistent RNA virus infections.

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

Antagonistic conflict between transposon-encoded introns and guide RNAs

TnpB nucleases represent the evolutionary precursors to CRISPR-Cas12 and are widespread in all domains of life. IS605-family TnpB homologs function as programmable RNA-guided homing endonucleases in bacteria, driving transposon maintenance through DNA double-strand break-stimulated homologous recombination. In this work, we uncovered molecular mechanisms of the transposition life cycle of IS607-family elements that, notably, also encode group I introns. We identified specific features for a candidate "IStron" from Clostridium botulinum that allow the element to carefully control the relative levels of spliced products versus functional guide RNAs. Our results suggest that IStron transcripts evolved an ability to balance competing and mutually exclusive activities that promote selfish transposon spread while limiting adverse fitness costs on the host. Collectively, this work highlights molecular innovation in the multifunctional utility of transposon-encoded noncoding RNAs.

RNAscape: geometric mapping and customizable visualization of RNA structure

Analyzing and visualizing the tertiary structure and complex interactions of RNA is essential for being able to mechanistically decipher their molecular functions in vivo. Secondary structure visualization software can portray many aspects of RNA; however, these layouts are often unable to preserve topological correspondence since they do not consider tertiary interactions between different regions of an RNA molecule. Likewise, quaternary interactions between two or more interacting RNA molecules are not considered in secondary structure visualization tools. The RNAscape webserver produces visualizations that can preserve topological correspondence while remaining both visually intuitive and structurally insightful. RNAscape achieves this by designing a mathematical structural mapping algorithm which prioritizes the helical segments, reflecting their tertiary organization. Non-helical segments are mapped in a way that minimizes structural clutter. RNAscape runs a plotting script that is designed to generate publication-quality images. RNAscape natively supports non-standard nucleotides, multiple base-pairing annotation styles and requires no programming experience. RNAscape can also be used to analyze RNA/DNA hybrid structures and DNA topologies, including G-quadruplexes. Users can upload their own three-dimensional structures or enter a Protein Data Bank (PDB) ID of an existing structure. The RNAscape webserver allows users to customize visualizations through various settings as desired. URL: https://rnascape.usc.edu/.

An intricate balancing act: Upstream and downstream frameshift co-regulatory elements

Targeting ribosomal frameshifting has emerged as a potential therapeutic intervention strategy against Covid-19. During ribosomal translation, a fraction of elongating ribosomes slips by one base in the 5' direction and enters a new reading frame for viral protein synthesis. Any interference with this process profoundly affects viral replication and propagation. For Covid-19, two RNA sites associated with ribosomal frameshifting for SARS-CoV-2 are positioned on the 5' and 3' of the frameshifting residues. Although much attention has been on the 3' frameshift element (FSE), the 5' stem-loop (attenuator hairpin, AH) can play a role. The formation of AH has been suggested to occur as refolding of the 3' RNA structure is triggered by ribosomal unwinding. However, the attenuation activity and the relationship between the two regions are unknown. To gain more insight into these two related viral RNAs and to further enrich our understanding of ribosomal frameshifting for SARS-CoV-2, we explore the RNA folding of both 5' and 3' regions associated with frameshifting. Using our graph-theory-based modeling tools to represent RNA secondary structures, "RAG" (RNA- As-Graphs), and conformational landscapes to analyze length-dependent conformational distributions, we show that AH coexists with the 3-stem pseudoknot of the 3' FSE (graph 3_6 in our dual graph notation) and alternative pseudoknot (graph 3_3) but less likely with other 3' FSE alternative folds (such as 3-way junction 3_5). This is because an alternative length-dependent Stem 1 (AS1) can disrupt the FSE pseudoknots and trigger other folds. In addition, we design four mutants for long lengths that stabilize or disrupt AH, AS1 or FSE pseudoknot to illustrate the deduced AH/AS1 roles and favor the 3_5, 3_6 or stem-loop. These mutants further show how a strengthened pseudoknot can result from a weakened AS1, while a dominant stem-loop occurs with a strengthened AS1. These structural and mutational insights into both ends of the FSE in SARS-CoV-2 advance our understanding of the SARS-CoV-2 frameshifting mechanism by suggesting a sequence of length-dependent folds, which in turn define potential therapeutic intervention techniques involving both elements. Our work also highlights the complexity of viral landscapes with length-dependent folds, and challenges in analyzing these multiple conformations.

De novo gene synthesis by an antiviral reverse transcriptase

Bacteria defend themselves from viral infection using diverse immune systems, many of which sense and target foreign nucleic acids. Defense-associated reverse transcriptase (DRT) systems provide an intriguing counterpoint to this immune strategy by instead leveraging DNA synthesis, but the identities and functions of their DNA products remain largely unknown. Here we show that DRT2 systems execute an unprecedented immunity mechanism that involves de novo gene synthesis via rolling-circle reverse transcription of a non-coding RNA (ncRNA). Unbiased profiling of RT-associated RNA and DNA ligands in DRT2-expressing cells revealed that reverse transcription generates concatenated cDNA repeats through programmed template jumping on the ncRNA. The presence of phage then triggers second-strand cDNA synthesis, leading to the production of long double-stranded DNA. Remarkably, this DNA product is efficiently transcribed, generating messenger RNAs that encode a stop codon-less, never-ending ORF (neo) whose translation causes potent growth arrest. Phylogenetic analyses and screening of diverse DRT2 homologs further revealed broad conservation of rolling-circle reverse transcription and Neo protein function. Our work highlights an elegant expansion of genome coding potential through RNA-templated gene creation, and challenges conventional paradigms of genetic information encoded along the one-dimensional axis of genomic DNA.

Analysis of Telomerase RNA Structure in Physcomitrium patens Indicates Functionally Relevant Transitions Between OPEN and CLOSED Conformations

Telomerase RNA (TR) conformation determines its function as a template for telomere synthesis and as a scaffold for the assembly of the telomerase nucleoprotein complex. Experimental analyses of TR secondary structure using DMS-Map Seq and SHAPE-Map Seq techniques show its CLOSED conformation as the consensus structure where the template region cannot perform its function. Our data show that the apparent discrepancy between experimental results and predicted TR functional conformation, mostly ignored in published studies, can be explained using data analysis based on single-molecule structure prediction from individual sequencing reads by the recently established DaVinci method. This method results in several clusters of secondary structures reflecting the structural dynamics of TR, possibly related to its multiple functional states. Interestingly, the presumed active (OPEN) conformation of TR corresponds to a minor fraction of TR under in vivo conditions. Therefore, structural polymorphism and dynamic TR transitions between CLOSED and OPEN conformations may be involved in telomerase activity regulation as a switch that functions independently of total TR transcript levels.

Advances in Understanding the Mechanism of Cap-Independent Cucurbit Aphid-Borne Yellows Virus Protein Synthesis

Non-canonical translation mechanisms have been described for many viral RNAs. In the case of several plant viruses, their protein synthesis is controlled by RNA elements in their genomic 3'-ends that are able to enhance cap-independent translation (3'-CITE). The proposed general mechanism of 3'-CITEs includes their binding to eukaryotic translation initiation factors (eIFs) that reach the 5'-end and AUG start codon through 5'-3'-UTR-interactions. It was previously shown that cucurbit aphid-borne yellows virus (CABYV) has a 3'-CITE, which varies in sequence and structure depending on the phylogenetic group to which the isolate belongs, possibly as a result of adaptation to the different geographical regions. In this work, the cap-independent translation mechanisms of two CABYV 3'-CITEs belonging to the Mediterranean (CMTE) and Asian (CXTE) groups, respectively, were studied. In vivo cap-independent translation assays show that these 3'-CITEs require the presence of the CABYV short genomic 5'-UTR with at least 40% adenines in cis and an accessible 5'-end for its activity. Additionally, they suggest that the eIF4E-independent CABYV 3'-CITE activities may not require either eIF4A or the eIF4F complex, but may depend on eIF4G and PABP. By pulling down host proteins using RNA baits containing both 5'- and 3'-CABYV-UTRs, 80 RNA binding proteins were identified. These interacted preferentially with either CMTE, CXTE, or both. One of these proteins, specifically interacting with the RNA containing CMTE, was HSP70.2. Preliminary results suggested that HSP70.2 may be involved in CMTE- but not CXTE-mediated cap-independent translation activity.

Antagonistic conflict between transposon-encoded introns and guide RNAs

TnpB nucleases represent the evolutionary precursors to CRISPR-Cas12 and are widespread in all domains of life, presumably due to the critical roles they play in transposon proliferation. IS605family TnpB homologs function in bacteria as programmable homing endonucleases by exploiting transposon-encoded guide RNAs to cleave vacant genomic sites, thereby driving transposon maintenance through DSB-stimulated homologous recombination. Whether this pathway is conserved in other genetic contexts, and in association with other transposases, is unknown. Here we uncover molecular mechanisms of transposition and RNA-guided DNA cleavage by IS607-family elements that, remarkably, also encode catalytic, self-splicing group I introns. After reconstituting and systematically investigating each of these biochemical activities for a candidate 'IStron' derived from Clostridium botulinum, we discovered sequence and structural features of the transposon-encoded RNA that satisfy molecular requirements of a group I intron and TnpB guide RNA, while still retaining the ability to be faithfully mobilized at the DNA level by the TnpA transposase. Strikingly, intron splicing was strongly repressed not only by TnpB, but also by the secondary structure of ωRNA alone, allowing the element to carefully control the relative levels of spliced products versus functional guide RNAs. Our results suggest that IStron transcripts have evolved a sensitive equilibrium to balance competing and mutually exclusive activities that promote transposon maintenance while limiting adverse fitness costs on the host. Collectively, this work explains how diverse enzymatic activities emerged during the selfish spread of IS607-family elements and highlights molecular innovation in the multi-functional utility of transposon-encoded noncoding RNAs.

Cellular roadmaps of viroid infection

Viroids are single-stranded circular noncoding RNAs that infect plants. According to the International Committee on Taxonomy of Viruses, there are 44 viroids known to date. Notably, more than 20 000 distinct viroid-like RNA sequences have recently been identified in existing sequencing datasets, suggesting an unprecedented complexity in biological roles of viroids and viroid-like RNAs. Interestingly, a human pathogen, hepatitis delta virus (HDV), also replicates via a rolling circle mechanism like viroids. Therefore, knowledge of viroid infection is informative for research on HDV and other viroid-like RNAs reported from various organisms. Here, we summarize recent advancements in understanding viroid shuttling among subcellular compartments for completing replication cycles, emphasizing regulatory roles of RNA motifs and structural dynamics in diverse biological processes. We also compare the knowledge of viroid intracellular trafficking with known pathways governing cellular RNA movement in cells. Future investigations on regulatory RNA structures and cognate factors in regulating viroid subcellular trafficking and replication will likely provide new insights into RNA structure-function relationships and facilitate the development of strategies controlling RNA localization and function in cells.

-1 Programmed ribosomal frameshifting in Class 2 umbravirus-like RNAs uses multiple long-distance interactions to shift between active and inactive structures and destabilize the frameshift stimulating element

Plus-strand RNA viruses frequently employ -1 programmed ribosomal frameshifting (-1 PRF) to maximize their coding capacity. Ribosomes can frameshift at a slippery sequence if progression is impeded by a frameshift stimulating element (FSE), which is generally a stable, complex, dynamic structure with multiple conformations that contribute to the efficiency of -1 PRF. As FSE are usually analyzed separate from the viral genome, little is known about cis-acting long-distance interactions. Using full-length genomic RNA of umbravirus-like (ula)RNA citrus yellow vein associated virus (CY1) and translation in wheat germ extracts, six tertiary interactions were found associated with the CY1 FSE that span nearly three-quarters of the 2.7 kb genomic RNA. All six tertiary interactions are conserved in other Class 2 ulaRNAs and two are conserved in all ulaRNAs. Two sets of interactions comprise local and distal pseudoknots that involve overlapping FSE nucleotides and thus are structurally incompatible, suggesting that Class 2 FSEs assume multiple conformations. Importantly, two long-distance interactions connect with sequences on opposite sides of the critical FSE central stem, which would unzip the stem and destabilize the FSE. These latter interactions could allow a frameshifting ribosome to translate through a structurally disrupted upstream FSE that no longer blocks ribosome progression.