Functional long-range RNA-RNA interactions in positive-strand RNA viruses

Recruitment of the 40S ribosomal subunit by the West Nile virus 3' UTR promotes the cross-talk between the viral genomic ends for translation regulation

Flaviviral RNA genomes are composed of discrete RNA structural units arranged in an ordered fashion and grouped into complex folded domains that regulate essential viral functions, e.g. replication and translation. This is achieved by adjusting the overall structure of the RNA genome via the establishment of inter- and intramolecular interactions. Translation regulation is likely the main process controlling flaviviral gene expression. Although the genomic 3' UTR is a key player in this regulation, little is known about the molecular mechanisms underlying this role. The present work provides evidence for the specific recruitment of the 40S ribosomal subunit by the 3' UTR of the West Nile virus RNA genome, showing that the joint action of both genomic ends contributes the positioning of the 40S subunit at the 5' end. The combination of structural mapping techniques revealed specific conformational requirements at the 3' UTR for 40S binding, involving the highly conserved SL-III, 5'DB, 3'DB and 3'SL elements, all involved in the translation regulation. These results point to the 40S subunit as a bridge to ensure cross-talk between both genomic ends during viral translation and support a link between 40S recruitment by the 3' UTR and translation control.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Mapping of the influenza A virus genome RNA structure and interactions reveals essential elements of viral replication

Influenza A virus (IAV) represents a constant public health threat. The single-stranded, segmented RNA genome of IAV is replicated in host cell nuclei as a series of 8 ribonucleoprotein complexes (vRNPs) with RNA structures known to exert essential function to support viral replication. Here, we investigate RNA secondary structures and RNA interactions networks of the IAV genome and construct an in vivo structure model for each of the 8 IAV genome segments. Our analyses reveal an overall in vivo and in virio resemblance of the IAV genome conformation but also wide disparities among long-range and intersegment interactions. Moreover, we identify a long-range RNA interaction that exerts an essential role in genome packaging. Disrupting this structure displays reduced infectivity, attenuating virus pathogenicity in mice. Our findings characterize the in vivo RNA structural landscape of the IAV genome and reveal viral RNA structures that can be targeted to develop antiviral interventions.

Host-like RNA Elements Regulate Virus Translation

Viruses are obligate, intracellular parasites that co-opt host cell machineries for propagation. Critical among these machineries are those that translate RNA into protein and their mechanisms of control. Most regulatory mechanisms effectuate their activity by targeting sequence or structural features at the RNA termini, i.e., at the 5' or 3' ends, including the untranslated regions (UTRs). Translation of most eukaryotic mRNAs is initiated by 5' cap-dependent scanning. In contrast, many viruses initiate translation at internal RNA regions at internal ribosome entry sites (IRESs). Eukaryotic mRNAs often contain upstream open reading frames (uORFs) that permit condition-dependent control of downstream major ORFs. To offset genome compression and increase coding capacity, some viruses take advantage of out-of-frame overlapping uORFs (oORFs). Lacking the essential machinery of protein synthesis, for example, ribosomes and other translation factors, all viruses utilize the host apparatus to generate virus protein. In addition, some viruses exhibit RNA elements that bind host regulatory factors that are not essential components of the translation machinery. SARS-CoV-2 is a paradigm example of a virus taking advantage of multiple features of eukaryotic host translation control: the virus mimics the established human GAIT regulatory element and co-opts four host aminoacyl tRNA synthetases to form a stimulatory binding complex. Utilizing discontinuous transcription, the elements are present and identical in all SARS-CoV-2 subgenomic RNAs (and the genomic RNA). Thus, the virus exhibits a post-transcriptional regulon that improves upon analogous eukaryotic regulons, in which a family of functionally related mRNA targets contain elements that are structurally similar but lacking sequence identity. This "thrifty" virus strategy can be exploited against the virus since targeting the element can suppress the expression of all subgenomic RNAs as well as the genomic RNA. Other 3' end viral elements include 3'-cap-independent translation elements (3'-CITEs) and 3'-tRNA-like structures. Elucidation of virus translation control elements, their binding proteins, and their mechanisms can lead to novel therapeutic approaches to reduce virus replication and pathogenicity.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

Subversion of selective autophagy for the biogenesis of tombusvirus replication organelles inhibits autophagy

Elaborate viral replication organelles (VROs) are formed to support positive-strand RNA virus replication in infected cells. VRO formation requires subversion of intracellular membranes by viral replication proteins. Here, we showed that the key ATG8f autophagy protein and NBR1 selective autophagy receptor were co-opted by Tomato bushy stunt virus (TBSV) and the closely-related carnation Italian ringspot virus. Knockdown of ATG8f or NBR1 in plants led to reduced tombusvirus replication, suggesting pro-viral function for selective autophagy. BiFC and proximity-labeling experiments showed that the TBSV p33 replication protein interacted with ATG8f and NBR1 to recruit them to VROs. In addition, we observed that several core autophagy proteins, such as ATG1a, ATG4, ATG5, ATG101 and the plant-specific SH3P2 autophagy adaptor proteins were also re-localized to TBSV VROs, suggesting that TBSV hijacks the autophagy machinery in plant cells. We demonstrated that subversion of autophagy components facilitated the recruitment of VPS34 PI3 kinase and enrichment of phospholipids, such as phosphatidylethanolamine and PI3P phosphoinositide in the VRO membranes. Hijacking of autophagy components into TBSV VROs led to inhibition of autophagic flux. We also found that a fraction of the subverted ATG8f and NBR1 was sequestered in biomolecular condensates associated with VROs. We propose that the VRO-associated condensates trap those autophagy proteins, taking them away from the autophagy pathway. Overall, tombusviruses hijack selective autophagy to provide phospholipid-rich membranes for replication and to regulate the antiviral autophagic flux.

SARS-CoV-2 RNA stabilizes host mRNAs to elicit immunopathogenesis

SARS-CoV-2 RNA interacts with host factors to suppress interferon responses and simultaneously induces cytokine release to drive the development of severe coronavirus disease 2019 (COVID-19). However, how SARS-CoV-2 hijacks host RNAs to elicit such imbalanced immune responses remains elusive. Here, we analyzed SARS-CoV-2 RNA in situ structures and interactions in infected cells and patient lung samples using RIC-seq. We discovered that SARS-CoV-2 RNA forms 2,095 potential duplexes with the 3' UTRs of 205 host mRNAs to increase their stability by recruiting RNA-binding protein YBX3 in A549 cells. Disrupting the SARS-CoV-2-to-host RNA duplex or knocking down YBX3 decreased host mRNA stability and reduced viral replication. Among SARS-CoV-2-stabilized host targets, NFKBIZ was crucial for promoting cytokine production and reducing interferon responses, probably contributing to cytokine storm induction. Our study uncovers the crucial roles of RNA-RNA interactions in the immunopathogenesis of RNA viruses such as SARS-CoV-2 and provides valuable host targets for drug development.

KARR-seq reveals cellular higher-order RNA structures and RNA-RNA interactions

RNA fate and function are affected by their structures and interactomes. However, how RNA and RNA-binding proteins (RBPs) assemble into higher-order structures and how RNA molecules may interact with each other to facilitate functions remain largely unknown. Here we present KARR-seq, which uses N3-kethoxal labeling and multifunctional chemical crosslinkers to covalently trap and determine RNA-RNA interactions and higher-order RNA structures inside cells, independent of local protein binding to RNA. KARR-seq depicts higher-order RNA structure and detects widespread intermolecular RNA-RNA interactions with high sensitivity and accuracy. Using KARR-seq, we show that translation represses mRNA compaction under native and stress conditions. We determined the higher-order RNA structures of respiratory syncytial virus (RSV) and vesicular stomatitis virus (VSV) and identified RNA-RNA interactions between the viruses and the host RNAs that potentially regulate viral replication.

Cooperation and Competition of RNA Secondary Structure and RNA-Protein Interactions in the Regulation of Alternative Splicing

The regulation of alternative splicing in eukaryotic cells is carried out through the coordinated action of a large number of factors, including RNA-binding proteins and RNA structure. The RNA structure influences alternative splicing by blocking cis-regulatory elements, or bringing them closer or farther apart. In combination with RNA-binding proteins, it generates transcript conformations that help to achieve the necessary splicing outcome. However, the binding of regulatory proteins depends on RNA structure and, vice versa, the formation of RNA structure depends on the interaction with regulators. Therefore, RNA structure and RNA-binding proteins are inseparable components of common regulatory mechanisms. This review highlights examples of alternative splicing regulation by RNA-binding proteins, the regulation through local and long-range RNA structures, as well as how these elements work together, cooperate, and compete.

Dimerization of an umbravirus RNA genome activates subgenomic mRNA transcription

Many eukaryotic RNA viruses transcribe subgenomic (sg) mRNAs during infections to control expression of a subset of viral genes. Such transcriptional events are commonly regulated by local or long-range intragenomic interactions that form higher-order RNA structures within these viral genomes. In contrast, here we report that an umbravirus activates sg mRNA transcription via base pair-mediated dimerization of its plus-strand RNA genome. Compelling in vivo and in vitro evidence demonstrate that this viral genome dimerizes via a kissing-loop interaction involving an RNA stem-loop structure located just upstream from its transcriptional initiation site. Both specific and non-specific features of the palindromic kissing-loop complex were found to contribute to transcriptional activation. Structural and mechanistic aspects of the process in umbraviruses are discussed and compared with genome dimerization events in other RNA viruses. Notably, probable dimer-promoting RNA stem-loop structures were also identified in a diverse group of umbra-like viruses, suggesting broader utilization of this unconventional transcriptional strategy.

Exploring structural determinants and the role of nucleolin in formation of the long-range interactions between untranslated regions of p53 mRNA

p53 protein is a key regulator of cellular homeostasis by coordinating the framework of antiproliferative pathways as a response to various stress factors. Although the main mechanism of stress-dependent induction of p53 protein relies on post-translational modifications influencing its stability and activity, a growing amount of evidence suggests that complex regulation of p53 expression occurs also at the mRNA level. This study explores structural determinants of long-range RNA-RNA interactions in p53 mRNA, crucial for stress-dependent regulation of p53 protein translation. We demonstrate that the 8-nt bulge motif plays a key structural role in base-pairing of complementary sequences from the 5' and 3' untranslated regions of p53 mRNA. We also show that one of the p53 translation regulators, nucleolin, displays an RNA chaperone activity and facilitates the association of sequences involved in the formation of long-range interactions in p53 mRNA. Nucleolin promotes base-pairing of complementary sequences through the bulge motif, because mutations of this region reduce or inhibit pairing while compensatory mutations restore this interaction. Mutational analysis of nucleolin reveals that all four RNA recognition motifs are indispensable for optimal RNA chaperone activity of nucleolin. These observations help to decipher the unique mechanism of p53 protein translation regulation pointing to bulge motif and nucleolin as the critical factors during intramolecular RNA-RNA recognition in p53 mRNA.

RNA levers and switches controlling viral gene expression

RNA viruses are diverse and abundant pathogens that are responsible for numerous human diseases. RNA viruses possess relatively compact genomes and have therefore evolved multiple mechanisms to maximize their coding capacities, often by encoding overlapping reading frames. These reading frames are then decoded by mechanisms such as alternative splicing and ribosomal frameshifting to produce multiple distinct proteins. These solutions are enabled by the ability of the RNA genome to fold into 3D structures that can mimic cellular RNAs, hijack host proteins, and expose or occlude regulatory protein-binding motifs to ultimately control key process in the viral life cycle. We highlight recent findings focusing on less conventional mechanisms of gene expression and new discoveries on the role of RNA structures.

Measuring intramolecular connectivity in long RNA molecules using two-dimensional DNA patch-probe arrays

We describe a simple method to infer intramolecular connections in a population of long RNA molecules in vitro. First we add DNA oligonucleotide "patches" that perturb the RNA connections, then we use a microarray containing a complete set of DNA oligonucleotide "probes" to record where perturbations occur. The pattern of perturbations reveals couplings between different regions of the RNA sequence, from which we infer connections as well as their prevalences in the population. We validate this patch-probe method using the 1,058-nucleotide RNA genome of satellite tobacco mosaic virus (STMV), which has previously been shown to have multiple long-range connections. Our results not only indicate long duplexes that agree with previous structures but also reveal the prevalence of competing connections. Together, these results suggest that globally-folded and locally-folded structures coexist in solution. We show that the prevalence of connections changes when pseudouridine, an important component of natural and synthetic RNA molecules, is substituted for uridine in STMV RNA.

A potential long-range RNA-RNA interaction in the HIV-1 RNA

It is well-established that viral and cellular mRNAs alike harbour functional long-range intra-molecular RNA-RNA interactions. Despite the biological importance of such interactions, their identification and characterization remain challenging. Here we present a computational method for the identification of certain kinds of long-range intra-molecular RNA-RNA interactions involving the loop nucleotides of a hairpin loop. Using the computational method, we analysed 4272 HIV-1 genomic mRNAs. A potential long-range intra-molecular RNA-RNA interaction within the HIV-1 genomic RNA was identified. The long-range interaction is mediated by a kissing loop structure between two stem-loops of the previously reported SHAPE-based secondary structure of the entire HIV-1 genome. Structural modelling studies were carried out to show that the kissing loop structure not only is sterically feasible, but also contains a conserved RNA structural motif often found in compact RNA pseudoknots. The computational method should be generally applicable to the identification of potential long-range intra-molecular RNA-RNA interactions in any viral or cellular mRNA sequence.Communicated by Ramaswamy H. Sarma.

Rotavirus RNA chaperone mediates global transcriptome-wide increase in RNA backbone flexibility

Due to genome segmentation, rotaviruses must co-package eleven distinct genomic RNAs. The packaging is mediated by virus-encoded RNA chaperones, such as the rotavirus NSP2 protein. While the activities of distinct RNA chaperones are well studied on smaller RNAs, little is known about their global effect on the entire viral transcriptome. Here, we used Selective 2'-hydroxyl Acylation Analyzed by Primer Extension and Mutational Profiling (SHAPE-MaP) to examine the secondary structure of the rotavirus transcriptome in the presence of increasing amounts of NSP2. SHAPE-MaP data reveals that despite the well-documented helix-unwinding activity of NSP2 in vitro, its incubation with cognate rotavirus transcripts does not induce a significant change in the SHAPE reactivities. However, a quantitative analysis of mutation rates measured by mutational profiling reveals a global 5-fold rate increase in the presence of NSP2. We demonstrate that the normalization procedure used in deriving SHAPE reactivities from mutation rates can mask an important global effect of an RNA chaperone. Analysis of the mutation rates reveals a larger effect on stems rather than loops. Together, these data provide the first experimentally derived secondary structure model of the rotavirus transcriptome and reveal that NSP2 acts by globally increasing RNA backbone flexibility in a concentration-dependent manner.

A Putative long-range RNA-RNA interaction between ORF8 and Spike of SARS-CoV-2

SARS-CoV-2 has affected people worldwide as the causative agent of COVID-19. The virus is related to the highly lethal SARS-CoV-1 responsible for the 2002-2003 SARS outbreak in Asia. Research is ongoing to understand why both viruses have different spreading capacities and mortality rates. Like other beta coronaviruses, RNA-RNA interactions occur between different parts of the viral genomic RNA, resulting in discontinuous transcription and production of various sub-genomic RNAs. These sub-genomic RNAs are then translated into other viral proteins. In this work, we performed a comparative analysis for novel long-range RNA-RNA interactions that may involve the Spike region. Comparing in-silico fragment-based predictions between reference sequences of SARS-CoV-1 and SARS-CoV-2 revealed several predictions amongst which a thermodynamically stable long-range RNA-RNA interaction between (23660-23703 Spike) and (28025-28060 ORF8) unique to SARS-CoV-2 was observed. The patterns of sequence variation using data gathered worldwide further supported the predicted stability of the sub-interacting region (23679-23690 Spike) and (28031-28042 ORF8). Such RNA-RNA interactions can potentially impact viral life cycle including sub-genomic RNA production rates.

Measurement of the specific and non-specific binding energies of Mg2+ to RNA

Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl2), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg2+ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl2) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals ΔG≃ 10 kcal/mol. Besides stabilizing the 3WJ, Mg2+ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

The helper oligonucleotides enable detection of folded single-stranded DNA by lateral flow immunoassay after HCR signal amplification

A combination of Hybridization Chain Reaction (HCR) and Lateral Flow Immunoassay (LFIA) is an attractive strategy for a simple signal amplification DNA/RNA detection. The present study aimed to report a strategy used to solve a problem encountered when the target DNA contained folded secondary structure during HCR, enabling HCR hairpin probes to easily access the target site. The 24-nt conserved sequence within 3'-UTR, present only in dengue virus genome but not in other species, is an ideal target to use as a probe binding site for pan-dengue virus detection. Thus, the 105-nt target containing the 24-nt target sequence was chosen as a target with secondary structures. The 24-nucleotide (nt) synthetic target DNA successfully induced HCR reaction within 5 min at room temperature. However, the HCR detection of the 105-nt synthetic target DNA with secondary structures was problematic. The probe hybridization was prevented by the secondary structures of the target, resulting in a failure to generate HCR product. To solve this problem, two helper oligonucleotides (helper1 and helper2) were designed to linearize the folded structure of the 105-nt target through strand-displacement mechanism, allowing the HCR hairpin probes to easily access the target site. The HCR product with the labeled helper oligonucleotides and the labeled probes were successfully detected by LFIA. With this strategy, the combination of the helper-enhanced HCR and LFIA exhibited a limit of detection (LOD) in a nanomolar range of the 105-nt DENV synthetic target DNA. Our study demonstrated that signal amplification by the combination of HCR and LFIA could successfully detect the target DNA with secondary structure, but not target RNA with secondary structure. In summary, this work provided a proof of concept of two main issues including probe hybridization enhancement by helper oligonucleotide for the target with complicated secondary structure and the advantage of a combination of labeled helper and HCR probes design for LFIA to overcome the false positive result from HCR probe leakage. Our findings on the use of helper oligonucleotides may be beneficial for the development of other isothermal amplification, since the secondary structure of the target is one of the major obstacles among hybridization-based methods.

Picornavirus translation strategies

The genome of viruses classified as picornaviruses consists of a single monocistronic positive strand RNA. The coding capacity of these RNA viruses is rather limited, and thus, they rely on the cellular machinery for their viral replication cycle. Upon the entry of the virus into susceptible cells, the viral RNA initially competes with cellular mRNAs for access to the protein synthesis machinery. Not surprisingly, picornaviruses have evolved specialized strategies that successfully allow the expression of viral gene products, which we outline in this review. The main feature of all picornavirus genomes is the presence of a heavily structured RNA element on the 5´UTR, referred to as an internal ribosome entry site (IRES) element, which directs viral protein synthesis as well and, consequently, triggers the subsequent steps required for viral replication. Here, we will summarize recent studies showing that picornavirus IRES elements consist of a modular structure, providing sites of interaction for ribosome subunits, eIFs, and a selective group of RNA‐binding proteins.

Identification of Novel 5' and 3' Translation Enhancers in Umbravirus-Like Coat Protein-Deficient RNA Replicons

Translation of plant plus-strand RNA viral genomes that lack a 5' cap frequently requires the use of cap-independent translation enhancers (CITEs) located in or near the 3' untranslated region (UTR). 3'CITEs are grouped based on secondary structure and ability to interact with different translation initiation factors or ribosomal subunits, which assemble a complex at the 3' end that is nearly always transferred to the 5' end via a long-distance kissing-loop interaction between sequences in the 3'CITE and 5' hairpins. We report here the identification of a novel 3'CITE in coat protein-deficient RNA replicons that are related to umbraviruses. Umbra-like associated RNAs (ulaRNAs), such as citrus yellow vein-associated virus (CYVaV), are a new type of subviral RNA that do not encode movement proteins, coat proteins, or silencing suppressors but can independently replicate using their encoded RNA-dependent RNA polymerase. An extended hairpin structure containing multiple internal loops in the 3' UTR of CYVaV is strongly conserved in the most closely related ulaRNAs and structurally resembles an I-shaped structure (ISS) 3'CITE. However, unlike ISS, the CYVaV structure binds to eIF4G and no long-distance interaction is discernible between the CYVaV ISS-like structure and sequences at or near the 5' end. We also report that the ∼30-nucleotide (nt) 5' terminal hairpin of CYVaV and related ulaRNAs can enhance translation of reporter constructs when associated with either the CYVaV 3'CITE or the 3'CITEs of umbravirus pea enation mosaic virus (PEMV2) and even independent of a 3'CITE. These findings introduce a new type of 3'CITE and provide the first information on translation of ulaRNAs. IMPORTANCE Umbra-like associated RNAs (ulaRNAs) are a recently discovered type of subviral RNA that use their encoded RNA-dependent RNA polymerase for replication but do not encode any coat proteins, movement proteins, or silencing suppressors yet can be found in plants in the absence of any discernible helper virus. We report the first analysis of their translation using class 2 ulaRNA citrus yellow vein-associated virus (CYVaV). CYVaV uses a novel eIF4G-binding I-shaped structure as its 3' cap-independent translation enhancer (3'CITE), which does not connect with the 5' end by a long-distance RNA:RNA interaction that is typical of 3'CITEs. ulaRNA 5' terminal hairpins can also enhance translation in association with cognate 3'CITEs or those of related ulaRNAs and, to a lesser extent, with 3'CITEs of umbraviruses, or even independent of a 3'CITE. These findings introduce a new type of 3'CITE and provide the first information on translation of ulaRNAs.

OUP accepted manuscript

Subdomain 5BSL3.2 of hepatitis C virus RNA lies at the core of a network of distal RNA–RNA contacts that connect the 5′ and 3′ regions of the viral genome and regulate the translation and replication stages of the viral cycle. Using small-angle X-ray scattering and NMR spectroscopy experiments, we have determined at low resolution the structural models of this subdomain and its distal complex with domain 3′X, located at the 3′-terminus of the viral RNA chain. 5BSL3.2 adopts a characteristic ‘L’ shape in solution, whereas the 5BSL3.2–3′X distal complex forms a highly unusual ‘Y’-shaped kissing junction that blocks the dimer linkage sequence of domain 3′X and promotes translation. The structure of this complex may impede an effective association of the viral polymerase with 5BSL3.2 and 3′X to start negative-strand RNA synthesis, contributing to explain the likely mechanism used by these sequences to regulate viral replication and translation. In addition, sequence and shape features of 5BSL3.2 are present in functional RNA motifs of flaviviruses, suggesting conserved regulatory processes within the Flaviviridae family.

Structures of artificially designed discrete RNA nanoarchitectures at near-atomic resolution

Although advances in nanotechnology have enabled the construction of complex and functional synthetic nucleic acid–based nanoarchitectures, high-resolution discrete structures are lacking because of the difficulty in obtaining good diffracting crystals. Here, we report the design and construction of RNA nanostructures based on homooligomerizable one-stranded tiles for x-ray crystallographic determination. We solved three structures to near-atomic resolution: a 2D parallelogram, a 3D nanobracelet unexpectedly formed from an RNA designed for a nanocage, and, eventually, a bona fide 3D nanocage designed with the guidance of the two previous structures. Structural details of their constituent motifs, such as kissing loops, branched kissing loops, and T-junctions, that resemble natural RNA motifs and resisted x-ray determination are revealed, providing insights into those natural motifs. This work unveils the largely unexplored potential of crystallography in gaining high-resolution feedback for nanoarchitectural design and suggests a route to investigate RNA motif structures by configuring them into nanoarchitectures.

Cotranslational prolyl hydroxylation is essential for flavivirus biogenesis

Viral pathogens are an ongoing threat to public health worldwide. Analysing their dependence on host biosynthetic pathways could lead to effective antiviral therapies1. Here we integrate proteomic analyses of polysomes with functional genomics and pharmacological interventions to define how enteroviruses and flaviviruses remodel host polysomes to synthesize viral proteins and disable host protein production. We find that infection with polio, dengue or Zika virus markedly modifies polysome composition, without major changes to core ribosome stoichiometry. These viruses use different strategies to evict a common set of translation initiation and RNA surveillance factors from polysomes while recruiting host machineries that are specifically required for viral biogenesis. Targeting these specialized viral polysomes could provide a new approach for antiviral interventions. For example, we find that both Zika and dengue use the collagen proline hydroxylation machinery to mediate cotranslational modification of conserved proline residues in the viral polyprotein. Genetic or pharmacological inhibition of proline hydroxylation impairs nascent viral polyprotein folding and induces its aggregation and degradation. Notably, such interventions prevent viral polysome remodelling and lower virus production. Our findings delineate the modular nature of polysome specialization at the virus-host interface and establish a powerful strategy to identify targets for selective antiviral interventions.

RNA structure probing uncovers RNA structure-dependent biological functions

RNA molecules fold into complex structures that enable their diverse functions in cells. Recent revolutionary innovations in transcriptome-wide RNA structural probing of living cells have ushered in a new era in understanding RNA functions. Here, we summarize the latest technological advances for probing RNA secondary structures and discuss striking discoveries that have linked RNA regulation and biological processes through interrogation of RNA structures. In particular, we highlight how different long noncoding RNAs form into distinct secondary structures that determine their modes of interactions with protein partners to realize their unique functions. These dynamic structures mediate RNA regulatory functions through altering interactions with proteins and other RNAs. We also outline current methodological hurdles and speculate about future directions for development of the next generation of RNA structure-probing technologies of higher sensitivity and resolution, which could then be applied in increasingly physiologically relevant studies.

The architecture of the SARS-CoV-2 RNA genome inside virion

SARS-CoV-2 carries the largest single-stranded RNA genome and is the causal pathogen of the ongoing COVID-19 pandemic. How the SARS-CoV-2 RNA genome is folded in the virion remains unknown. To fill the knowledge gap and facilitate structure-based drug development, we develop a virion RNA in situ conformation sequencing technology, named vRIC-seq, for probing viral RNA genome structure unbiasedly. Using vRIC-seq data, we reconstruct the tertiary structure of the SARS-CoV-2 genome and reveal a surprisingly "unentangled globule" conformation. We uncover many long-range duplexes and higher-order junctions, both of which are under purifying selections and contribute to the sequential package of the SARS-CoV-2 genome. Unexpectedly, the D614G and the other two accompanying mutations may remodel duplexes into more stable forms. Lastly, the structure-guided design of potent small interfering RNAs can obliterate the SARS-CoV-2 in Vero cells. Overall, our work provides a framework for studying the genome structure, function, and dynamics of emerging deadly RNA viruses.

The SARS-CoV-2 subgenome landscape and its novel regulatory features

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

Structural Analysis and Whole Genome Mapping of a New Type of Plant Virus Subviral RNA: Umbravirus-Like Associated RNAs

We report the biological and structural characterization of umbravirus-like associated RNAs (ulaRNAs), a new category of coat-protein dependent subviral RNA replicons that infect plants. These RNAs encode an RNA-dependent RNA polymerase (RdRp) following a -1 ribosomal frameshift event, are 2.7-4.6 kb in length, and are related to umbraviruses, unlike similar RNA replicons that are related to tombusviruses. Three classes of ulaRNAs are proposed, with citrus yellow vein associated virus (CYVaV) placed in Class 2. With the exception of CYVaV, Class 2 and Class 3 ulaRNAs encode an additional open reading frame (ORF) with movement protein-like motifs made possible by additional sequences just past the RdRp termination codon. The full-length secondary structure of CYVaV was determined using Selective 2' Hydroxyl Acylation analyzed by Primer Extension (SHAPE) structure probing and phylogenic comparisons, which was used as a template for determining the putative structures of the other Class 2 ulaRNAs, revealing a number of distinctive structural features. The ribosome recoding sites of nearly all ulaRNAs, which differ significantly from those of umbraviruses, may exist in two conformations and are highly efficient. The 3' regions of Class 2 and Class 3 ulaRNAs have structural elements similar to those of nearly all umbraviruses, and all Class 2 ulaRNAs have a unique, conserved 3' cap-independent translation enhancer. CYVaV replicates independently in protoplasts, demonstrating that the reported sequence is full-length. Additionally, CYVaV contains a sequence in its 3' UTR that confers protection to nonsense mediated decay (NMD), thus likely obviating the need for umbravirus ORF3, a known suppressor of NMD. This initial characterization lays down a road map for future investigations into these novel virus-like RNAs.

In vivo structural characterization of the SARS-CoV-2 RNA genome identifies host proteins vulnerable to repurposed drugs

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic. Understanding of the RNA virus and its interactions with host proteins could improve therapeutic interventions for COVID-19. By using icSHAPE, we determined the structural landscape of SARS-CoV-2 RNA in infected human cells and from refolded RNAs, as well as the regulatory untranslated regions of SARS-CoV-2 and six other coronaviruses. We validated several structural elements predicted in silico and discovered structural features that affect the translation and abundance of subgenomic viral RNAs in cells. The structural data informed a deep-learning tool to predict 42 host proteins that bind to SARS-CoV-2 RNA. Strikingly, antisense oligonucleotides targeting the structural elements and FDA-approved drugs inhibiting the SARS-CoV-2 RNA binding proteins dramatically reduced SARS-CoV-2 infection in cells derived from human liver and lung tumors. Our findings thus shed light on coronavirus and reveal multiple candidate therapeutics for COVID-19 treatment.

RanDeL-Seq: a High-Throughput Method to Map Viral cis- and trans-Acting Elements

Recent studies have renewed interest in developing novel antiviral therapeutics and vaccines based on defective interfering particles (DIPs)—a subset of viral deletion mutants that conditionally replicate. Identifying and engineering DIPs require that viral cis- and trans-acting elements be accurately mapped.

It has long been known that noncoding genomic regions can be obligate cis elements acted upon in trans by gene products. In viruses, cis elements regulate gene expression, encapsidation, and other maturation processes, but mapping these elements relies on targeted iterative deletion or laborious prospecting for rare spontaneously occurring mutants. Here, we introduce a method to comprehensively map viral cis and trans elements at single-nucleotide resolution by high-throughput random deletion. Variable-size deletions are randomly generated by transposon integration, excision, and exonuclease chewback and then barcoded for tracking via sequencing (i.e., random deletion library sequencing [RanDeL-seq]). Using RanDeL-seq, we generated and screened >23,000 HIV-1 variants to generate a single-base resolution map of HIV-1’s cis and trans elements. The resulting landscape recapitulated HIV-1’s known cis-acting elements (i.e., long terminal repeat [LTR], Ψ, and Rev response element [RRE]) and, surprisingly, indicated that HIV-1’s central DNA flap (i.e., central polypurine tract [cPPT] to central termination sequence [CTS]) is as critical as the LTR, Ψ, and RRE for long-term passage. Strikingly, RanDeL-seq identified a previously unreported ∼300-bp region downstream of RRE extending to splice acceptor 7 that is equally critical for sustained viral passage. RanDeL-seq was also used to construct and screen a library of >90,000 variants of Zika virus (ZIKV). Unexpectedly, RanDeL-seq indicated that ZIKV’s cis-acting regions are larger than the untranscribed (UTR) termini, encompassing a large fraction of the nonstructural genes. Collectively, RanDeL-seq provides a versatile framework for generating viral deletion mutants, enabling discovery of replication mechanisms and development of novel antiviral therapeutics, particularly for emerging viral infections.

Development of Tomato bushy stunt virus-based vectors for fusion and non-fusion expression of heterologous proteins in an alternative host Nicotiana excelsiana

Plant virus-based expression systems are an alternative expression platform for the production of clinically and industrially useful recombinant proteins. Nonetheless, due to a lack of viral vector with the commercial potentials, it is urgent to design and develop new, versatile, and efficient plant virus vectors. The genome of Tomato bushy stunt virus (TBSV) offers an attractive alternative to being modified as a vector for producing heterologous proteins in plants. Here, we developed a set of novel fusion and non-fusion TBSV-CP replacement vectors, which provide more flexible and efficient tools for expressing proteins of interest in plants. An alternative tobacco plant, Nicotiana excelsiana, was used in this study as a host for newly constructed TBSV vectors because the unwanted necrotic effects were reported on the commonly used Nicotiana benthamiana host associated with expression of TBSV-encoded P19 protein. The data showed that TBSV vectors caused a symptomless infection and overexpressed reporter gene in N. excelsiana leaves, demonstrating that N. excelsiana is an ideal host plant for TBSV-mediated heterologous gene expression. Moreover, a TBSV non-fusion vector, dAUG, shows the similar accumulation level of reporter proteins to that of TMV- and PVX-based vectors in side-by-side comparison and provides more flexible aspects than the previously developed TBSV vectors. Collectively, our newly developed TBSV expression system adds a new member to the family of plant viral expression vectors and meanwhile offers a flexible and highly effective approach for producing proteins of interest in plants. KEY POINTS: • The TBSV-based transient expression system has been significantly improved. • The necrotic effects caused by viral P19 protein were avoided by the usage of N. excelsiana as a host plant. • The expression level of the non-fusion vector was similar to the most effective virus vectors reported so far.

Activation of viral transcription by stepwise largescale folding of an RNA virus genome

The genomes of RNA viruses contain regulatory elements of varying complexity. Many plus-strand RNA viruses employ largescale intra-genomic RNA-RNA interactions as a means to control viral processes. Here, we describe an elaborate RNA structure formed by multiple distant regions in a tombusvirus genome that activates transcription of a viral subgenomic mRNA. The initial step in assembly of this intramolecular RNA complex involves the folding of a large viral RNA domain, which generates a discontinuous binding pocket. Next, a distally-located protracted stem-loop RNA structure docks, via base-pairing, into the binding site and acts as a linchpin that stabilizes the RNA complex and activates transcription. A multi-step RNA folding pathway is proposed in which rate-limiting steps contribute to a delay in transcription of the capsid protein-encoding viral subgenomic mRNA. This study provides an exceptional example of the complexity of genome-scale viral regulation and offers new insights into the assembly schemes utilized by large intra-genomic RNA structures.

Site-specific covalent labeling of large RNAs with nanoparticles empowered by expanded genetic alphabet transcription

Conjugation of RNAs with nanoparticles (NPs) is of significant importance because of numerous applications in biology and medicine, which, however, remains challenging especially for large ones. So far, the majority of RNA labeling relies on solid-phase chemical synthesis, which is generally limited to RNAs smaller than 100 nucleotides (nts). We, here, present an efficient and generally applicable labeling strategy for site-specific covalent conjugation of large RNAs with a gold nanoparticle (Nanogold) empowered by transcription of an expanded genetic alphabet containing the A-T/U and G-C natural base pairs (bps) and the TPT3-NaM unnatural base pair (UBP). We synthesize an amine-derivatized TPT3 (TPT3^A), which is site specifically incorporated into a 97-nt 3'SL RNA and a 719-nt minigenomic RNA (DENV-mini) from Dengue virus serotype 2 (DENV2) by in vitro T7 transcription. The TPT3^A-modified RNAs are covalently conjugated with mono-Sulfo-N-hydroxysuccinimidyl (NHS)-Nanogold NPs via an amine and NHS ester reaction and further purified under nondenaturing conditions. TPT3 modification and Nanogold labeling cause minimal structural perturbations to the RNAs by circular dichroism, small angle X-ray scattering (SAXS), and binding activity assay. We demonstrate the application of the Nanogold-RNA conjugates in large RNA structural biology by an emerging molecular ruler, X-ray scattering interferometry (XSI). The internanoparticle distance distributions in the 3'SL and DENV-mini RNAs derived from XSI measurements support the hypothetical model of flavivirus genome circularization, thus, validate the applicability of this labeling strategy. The presented strategy overcomes the size constraints in conventional RNA labeling strategies and is expected to have wide applications in large RNA structural biology and RNA nanotechnology.

Structures and Functions of the 3' Untranslated Regions of Positive-Sense Single-Stranded RNA Viruses Infecting Humans and Animals

The 3′ untranslated region (3′ UTR) of positive-sense single-stranded RNA [ssRNA(+)] viruses is highly structured. Multiple elements in the region interact with other nucleotides and proteins of viral and cellular origin to regulate various aspects of the virus life cycle such as replication, translation, and the host-cell response. This review attempts to summarize the primary and higher order structures identified in the 3′UTR of ssRNA(+) viruses and their functional roles.

Prolonged (6-Month) Shedding of Middle East Respiratory Syndrome Coronavirus RNA in the Sputum of a Lymphoma Patient

During the 2015 Korea Middle East Respiratory Syndrome coronavirus (MERS-CoV) outbreak, a lymphoma patient developed MERS pneumonia. His pneumonia improved by 45 days after illness onset, but the polymerase chain reaction tests remained (+) for 6 months. However, replication-competent virus was detected by 60 days after illness onset.

G2-quadruplex in the 3'UTR of IE180 regulates Pseudorabies virus replication by enhancing gene expression

RNA secondary structure elements in the mRNA 3'-untranslated regions (3'UTR) play important roles in post-transcriptional regulation. RNA structure elements in the viral RNA provide valuable model for studying diverse regulation mechanisms. Herpesvirus genomes are double-stranded DNA with GC-rich sequences, which can be transcribed into abundant GC-rich RNAs. It is valuable to explore the structures and function of those GC-rich RNAs. We identified a G₂-quadruplex-forming sequence named PQS18-1 in the 3'UTR of the unique immediate early gene of Pseudorabies virus (PRV), an important member of Alphaherpesvirinae subfamily. The RNA PQS18-1 was folded into parallel G-quadruplex structure, enhancing gene expression. Both non-G-quadruplex mutant and G₃-quadruplex mutant in the 3'UTR showed lower gene expression level than the wildtype G₂-quadruplex. TMPyP4 destroyed PQS18-1 G₂-quadruplex and suppressed gene expression, accordingly reducing PRV replication by one titre in the PK15 cells at 24 h post infection. Our findings indicated that the RNA G₂-quadruplex in 3'UTR was essential for high expression of IE180 gene, and it could be a specific post-transcription regulation element in response to small molecules or other macromolecules. This study discovers a novel RNA G₂-quadruplex in the 3'UTR of an immediate early gene of alphaherpesvirus and provides a new nucleic acid target for anti-virus drug design.

Different Degrees of 5'-to-3' DAR Interactions Modulate Zika Virus Genome Cyclization and Host-Specific Replication

Mosquito-borne flaviviruses, which include many important human pathogens, such as West Nile virus (WNV), dengue virus (DENV), and Zika virus (ZIKV), have caused numerous emerging epidemics in recent years. Details of the viral genome functions necessary for effective viral replication in mosquito and vertebrate hosts remain obscure. Here, using ZIKV as a model, we found that the conserved "downstream of AUG region" (DAR), which is known to be an essential element for genome cyclization, is involved in viral replication in a host-specific manner. Mutational analysis of the DAR element showed that a single-nucleotide mismatch between the 5' DAR and the 3' DAR had little effect on ZIKV replication in mammalian cells but dramatically impaired viral propagation in mosquito cells. The revertant viruses passaged in mosquito cells generated compensatory mutations restoring the base pairing of the DAR, further confirming the importance of the complementarity of the DAR in mosquito cells. We demonstrate that a single-nucleotide mutation in the DAR is sufficient to destroy long-range RNA interaction of the ZIKV genome and affects de novo RNA synthesis at 28°C instead of 37°C, resulting in the different replication efficiencies of the mutant viruses in mosquito and mammalian cells. Our results reveal a novel function of the circular form of the flavivirus genome in host-specific viral replication, providing new ideas to further explore the functions of the viral genome during host adaptation.IMPORTANCE Flaviviruses naturally cycle between the mosquito vector and vertebrate hosts. The disparate hosts provide selective pressures that drive virus genome evolution to maintain efficient replication during host alteration. Host adaptation may occur at different stages of the viral life cycle, since host-specific viral protein processing and virion conformations have been reported in the individual hosts. However, the viral determinants and the underlying mechanisms associated with host-specific functions remain obscure. In this study, using Zika virus, we found that the DAR-mediated genome cyclization regulates viral replication differently and is under different selection pressures in mammalian and mosquito cells. A more constrained complementarity of the DAR is required in mosquito cells than in mammalian cells. Since the DAR element is stably maintained among mosquito-borne flaviviruses, our findings could provide new information for understanding the role of flavivirus genome cyclization in viral adaptation and RNA evolution in the two hosts.

The evolution of Zika virus from Asia to the Americas

Zika virus (ZIKV) was once considered an obscure member of the large and diverse family of mosquito-borne flaviviruses, and human infections with ZIKV were thought to be sporadic, with mild and self-limiting symptoms. The large-scale ZIKV epidemics in the Americas and the unexpected uncovering of a link to congenital birth defects escalated ZIKV infections to the status of a global public health emergency. Recent studies that combined reverse genetics with modelling in multiple systems have provided evidence that ZIKV has acquired additional amino acid substitutions at the same time as congenital Zika syndrome and other birth defects were detected. In this Progress article, we summarize the evolution of ZIKV during its spread from Asia to the Americas and discuss potential links to pathogenesis.

Communication Is Key: 5'-3' Interactions that Regulate mRNA Translation and Turnover

Most eukaryotic mRNAs maintain a 5' cap structure and 3' poly(A) tail, cis-acting elements that are often separated by thousands of nucleotides. Nevertheless, multiple paradigms exist where mRNA 5' and 3' termini interact with each other in order to regulate mRNA translation and turnover. mRNAs recruit translation initiation factors to their termini, which in turn physically interact with each other. This physical bridging of the mRNA termini is known as the "closed loop" model, with years of genetic and biochemical evidence supporting the functional synergy between the 5' cap and 3' poly(A) tail to enhance mRNA translation initiation. However, a number of examples exist of "non-canonical" 5'-3' communication for cellular and viral RNAs that lack 5' cap structures and/or poly(A) tails. Moreover, in several contexts, mRNA 5'-3' communication can function to repress translation. Overall, we detail how various mRNA 5'-3' interactions play important roles in posttranscriptional regulation, wherein depending on the protein factors involved can result in translational stimulation or repression.

Highly efficacious antiviral protection of plants by small interfering RNAs identified in vitro

In response to a viral infection, the plant’s RNA silencing machinery processes viral RNAs into a huge number of small interfering RNAs (siRNAs). However, a very few of these siRNAs actually interfere with viral replication. A reliable approach to identify these immunologically effective siRNAs (esiRNAs) and to define the characteristics underlying their activity has not been available so far. Here, we develop a novel screening approach that enables a rapid functional identification of antiviral esiRNAs. Tests on the efficacy of such identified esiRNAs of a model virus achieved a virtual full protection of plants against a massive subsequent infection in transient applications. We find that the functionality of esiRNAs depends crucially on two properties: the binding affinity to Argonaute proteins and the ability to access the target RNA. The ability to rapidly identify functional esiRNAs could be of great benefit for all RNA silencing-based plant protection measures against viruses and other pathogens.

RNA2Drawer: geometrically strict drawing of nucleic acid structures with graphical structure editing and highlighting of complementary subsequences

RNA structure prediction programs remain imperfect and many substructures are still identified by manual exploration, which is most efficiently conducted within an RNA structure drawing program. However, most nucleic acid structure drawing programs have limited capability for structure modification (i.e., breaking and forming new bonds between bases), often requiring that the structure notation be textually edited. RNA2Drawer was developed to allow for graphical structure editing while maintaining the geometry of a drawing (e.g., ellipsoid loops, stems with evenly stacked base pairs) throughout structural changes and manual adjustments to the layout by the user. In addition, the program allows for annotations such as colouring and circling of bases and drawing of tertiary interactions (e.g., pseudoknots). RNA2Drawer can also draw commonly desired elements such as an optionally flattened outermost loop and assists structure editing by automatically highlighting complementary subsequences, which facilitates the discovery of potentially new and alternative pairings, particularly tertiary pairings over long-distances, which are biologically critical in the genomes of many RNA viruses and cannot be accurately predicted by current structure prediction programs. Additionally, RNA2Drawer outputs drawings either as PNG files, or as PPTX and SVG files, such that every object of a drawing (e.g., bases, bonds) is an individual PPTX or SVG object, allowing for further manipulation in Microsoft PowerPoint or a vector graphics editor such as Adobe Illustrator. PowerPoint is the standard for presentations and is often used to create figures for publications, and RNA2Drawer is the first program to export drawings as PPTX files.

An unexpected RNA distal interaction mode found in an essential region of the hepatitis C virus genome

The 3’X tail is a functionally essential 98-nt sequence located at the 3′-end of the hepatitis C virus (HCV) RNA genome. The domain contains two absolutely conserved dimer linkage sequence (DLS) and k nucleotide segments involved in viral RNA dimerization and in a distal base-pairing interaction with stem-loop 5BSL3.2, respectively. We have previously shown that domain 3’X forms an elongated structure comprising two coaxially stacked SL1’ and SL2’ stem-loops. This conformation favors RNA dimerization by exposing a palindromic DLS segment in an apical loop, but buries in the upper stem of hairpin SL2’ the k nucleotides involved in the distal contact with 5BSL3.2. Using nuclear magnetic resonance spectroscopy and gel electrophoresis experiments, here we show that the establishment of the complex between domain 3’X and stem-loop 5BSL3.2 only requires a rearrangement of the nucleotides forming the upper region of subdomain SL2’. The results indicate that the interaction does not occur through a canonical kissing loop mechanism involving the unpaired nucleotides of two terminal loops, but rather involves a base-paired stem and an apical loop and may result in a kissing three-way junction. On the basis of this information we suggest how the 3’X tail switches between monomer, homodimer and heterodimer states to regulate the HCV viral cycle.

Physical and Functional Analysis of Viral RNA Genomes by SHAPE

RNA viruses encode the information required to usurp cellular metabolism and gene regulation and to enable their own replication in two ways: in the linear sequence of their RNA genomes and in higher-order structures that form when the genomic RNA strand folds back on itself. Application of high-resolution SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) structure probing to viral RNA genomes has identified numerous new regulatory elements, defined new principles by which viral RNAs interact with the cellular host and evade host immune responses, and revealed relationships between virus evolution and RNA structure. This review summarizes our current understanding of genome structure-function interrelationships for RNA viruses, as informed by SHAPE structure probing, and outlines opportunities for future studies.

The Third Annual Meeting of the European Virus Bioinformatics Center

The Third Annual Meeting of the European Virus Bioinformatics Center (EVBC) took place in Glasgow, United Kingdom, 28-29 March 2019. Virus bioinformatics has become central to virology research, and advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks, being successfully used to detect, control, and treat infections of humans and animals. This active field of research has attracted approximately 110 experts in virology and bioinformatics/computational biology from Europe and other parts of the world to attend the two-day meeting in Glasgow to increase scientific exchange between laboratory- and computer-based researchers. The meeting was held at the McIntyre Building of the University of Glasgow; a perfect location, as it was originally built to be a place for "rubbing your brains with those of other people", as Rector Stanley Baldwin described it. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The meeting featured eight invited and twelve contributed talks, on the four main topics: (1) systems virology, (2) virus-host interactions and the virome, (3) virus classification and evolution and (4) epidemiology, surveillance and evolution. Further, the meeting featured 34 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

Structure mapping of dengue and Zika viruses reveals functional long-range interactions

Dengue (DENV) and Zika (ZIKV) viruses are clinically important members of the Flaviviridae family with an 11 kb positive strand RNA genome that folds to enable virus function. Here, we perform structure and interaction mapping on four DENV and ZIKV strains inside virions and in infected cells. Comparative analysis of SHAPE reactivities across serotypes nominates potentially functional regions that are highly structured, conserved, and contain low synonymous mutation rates. Interaction mapping by SPLASH identifies many pair-wise interactions, 40% of which form alternative structures, suggesting extensive structural heterogeneity. Analysis of shared interactions between serotypes reveals a conserved macro-organization whereby interactions can be preserved at physical locations beyond sequence identities. We further observe that longer-range interactions are preferentially disrupted inside cells, and show the importance of new interactions in virus fitness. These findings deepen our understanding of Flavivirus genome organization and serve as a resource for designing therapeutics in targeting RNA viruses.

Translational Control in Virus-Infected Cells

As obligate intracellular parasites, virus reproduction requires host cell functions. Despite variations in genome size and configuration, nucleic acid composition, and their repertoire of encoded functions, all viruses remain unconditionally dependent on the protein synthesis machinery resident within their cellular hosts to translate viral messenger RNAs (mRNAs). A complex signaling network responsive to physiological stress, including infection, regulates host translation factors and ribosome availability. Furthermore, access to the translation apparatus is patrolled by powerful host immune defenses programmed to restrict viral invaders. Here, we review the tactics and mechanisms used by viruses to appropriate control over host ribosomes, subvert host defenses, and dominate the infected cell translational landscape. These not only define aspects of infection biology paramount for virus reproduction, but continue to drive fundamental discoveries into how cellular protein synthesis is controlled in health and disease.

SilentMutations (SIM): A tool for analyzing long-range RNA-RNA interactions in viral genomes and structured RNAs

•

We developed a tool to analyze the effect of multiple point mutations on the secondary structures of two interacting viral RNAs.

•

Our tool simulates destructive and compensatory mutants of two key regions from a single-stranded RNA.

•

The simulated mutants can be utilized for the combinatorial in vitro analysis of RNA–RNA interactions.

•

We predicted potential mutants for in vitro validation experiments of influenza A virus and hepatitis C virus interactions.

A single nucleotide change in the coding region can alter the amino acid sequence of a protein. In consequence, natural or artificial sequence changes in viral RNAs may have various effects not only on protein stability, function and structure but also on viral replication.

In recent decades, several tools have been developed to predict the effect of mutations in structured RNAs such as viral genomes or non-coding RNAs. Some tools use multiple point mutations and also take coding regions into account. However, none of these tools was designed to specifically simulate the effect of mutations on viral long-range interactions.

Here, we developed SilentMutations (SIM), an easy-to-use tool to analyze the effect of multiple point mutations on the secondary structures of two interacting viral RNAs. The tool can simulate disruptive and compensatory mutants of two interacting single-stranded RNAs. This allows a fast and accurate assessment of key regions potentially involved in functional long-range RNA–RNA interactions and will eventually help virologists and RNA-experts to design appropriate experiments.

SIM only requires two interacting single-stranded RNA regions as input. The output is a plain text file containing the most promising mutants and a graphical representation of all interactions.

We applied our tool on two experimentally validated influenza A virus and hepatitis C virus interactions and we were able to predict potential double mutants for in vitro validation experiments.

The source code and documentation of SIM are freely available at github.com/desiro/silentMutations.

ScanFold: an approach for genome-wide discovery of local RNA structural elements-applications to Zika virus and HIV

In addition to encoding RNA primary structures, genomes also encode RNA secondary and tertiary structures that play roles in gene regulation and, in the case of RNA viruses, genome replication. Methods for the identification of functional RNA structures in genomes typically rely on scanning analysis windows, where multiple partially-overlapping windows are used to predict RNA structures and folding metrics to deduce regions likely to form functional structure. Separate structural models are produced for each window, where the step size can greatly affect the returned model. This makes deducing unique local structures challenging, as the same nucleotides in each window can be alternatively base paired. We are presenting here a new approach where all base pairs from analysis windows are considered and weighted by favorable folding. This results in unique base pairing throughout the genome and the generation of local regions/structures that can be ranked by their propensity to form unusually thermodynamically stable folds. We applied this approach to the Zika virus (ZIKV) and HIV-1 genomes. ZIKV is linked to a variety of neurological ailments including microcephaly and Guillain-Barré syndrome and its (+)-sense RNA genome encodes two, previously described, functionally essential structured RNA regions. HIV, the cause of AIDS, contains multiple functional RNA motifs in its genome, which have been extensively studied. Our approach is able to successfully identify and model the structures of known functional motifs in both viruses, while also finding additional regions likely to form functional structures. All data have been archived at the RNAStructuromeDB (www.structurome.bb.iastate.edu), a repository of RNA folding data for humans and their pathogens.

Revisiting the Closed-Loop Model and the Nature of mRNA 5'-3' Communication

Communication between the 5' and 3' ends of mature eukaryotic mRNAs lies at the heart of gene regulation, likely arising at the same time as the eukaryotic lineage itself. Our view of how and why it occurs has been shaped by elegant experiments that led to nearly universal acceptance of the "closed-loop model." However, new observations suggest that this classic model needs to be reexamined, revised, and expanded. Here, we address fundamental questions about the closed-loop model and discuss how a growing understanding of mRNA structure, dynamics, and intermolecular interactions presents new experimental opportunities. We anticipate that the application of emerging methods will lead to expanded models that include the role of intrinsic mRNA structure and quantitative dynamic descriptions of 5'-3' proximity linked to the functional status of an mRNA and will better reflect the messy realities of the crowded and rapidly changing cellular environment.