Mechanism and structural diversity of exoribonuclease-resistant RNA structures in flaviviral RNAs

Tick-borne flavivirus exoribonuclease-resistant RNAs contain a double loop structure

Viruses from the Flaviviridae family contain human relevant pathogens that generate subgenomic noncoding RNAs during infection using structured exoribonuclease resistant RNAs (xrRNAs). These xrRNAs block progression of host cell's 5' to 3' exoribonucleases. The structures of several xrRNAs from mosquito-borne and insect-specific flaviviruses reveal a conserved fold in which a ring-like motif encircles the 5' end of the xrRNA. However, the xrRNAs found in tick-borne and no known vector flaviviruses have distinct characteristics, and their 3-D fold was unsolved. Here, we verify the presence of xrRNAs in the encephalitis-causing tick-borne Powassan Virus. We characterize their secondary structure and obtain a mid-resolution map of one of these xrRNAs using cryo-EM, revealing a unique double-loop ring element. Integrating these results with covariation analysis, biochemical data, and existing high-resolution structural information yields a model in which the core of the fold matches the previously solved xrRNA fold, but the expanded double loop ring is remodeled upon encountering the exoribonuclease. These results are representative of a broad class of xrRNAs and reveal a conserved strategy of structure-based exoribonuclease resistance achieved through a unique topology across a viral family of importance to global health.

Identifying the interactions conferring functional mechanical rigidity on RNase-resistant RNA from Zika virus

Some viruses counter host-cell efforts to digest invading viral RNA by using special structures resistant to host RNases, known as exoribonuclease-resistant RNAs (xrRNAs). xrRNAs typically form an unusual fold with the 5'-end threaded through a ring consisting of a multihelix junction closed by a pseudoknot. By using single-molecule force spectroscopy (SMFS), we previously showed that a Zika virus xrRNA is extremely rigid mechanically, withstanding very high forces, and that this mechanical resistance-not simply the knot-like fold topology-is essential for RNase resistance. Here, we have determined which interactions are most important for generating mechanical rigidity in the Zika virus xrRNA, by systematically mutating tertiary contacts. We found that removing any of the tertiary contacts involving the threaded 5' end was sufficient to abrogate mechanical resistance. In contrast, breaking a single pseudoknot base pair was not sufficient to do so: Two broken pairs were needed. This hierarchy of interaction importance for mechanical rigidity was supported by simulations mapping how mechanical tension was distributed within the xrRNA. For all mutants, RNase resistance varied in lock-step with mechanical resistance, confirming the primary role of mechanical rigidity in xrRNA function. This work reveals which interactions are most important for Zika xrRNA function, with implications for targeting the xrRNA therapeutically.

The zinc finger protein ZFP36L2 inhibits flavivirus infection via the 5'-3' XRN1-mediated RNA decay pathway in the replication complexes

BACKGROUND

The zinc finger protein 36-like (ZFP36L) family is a CCCH-type group consisting of RNA-binding proteins, i.e., ZFP36L1 and ZFP36L2, which regulate cellular mRNA through the RNA decay pathway. ZFP36L1 combats flavivirus infections through the 5'-3' XRN1 and 3'-5' RNA exosome decay pathways. The present study clarified the role of human ZFP36L2 in the defense response of the host against flavivirus infection.

METHODS

Cell lines with overexpression or knockdown of ZFP36L2 were established using lentiviral vectors carrying genes for overexpression and short-hairpin RNA targeting specific genes, respectively. A plaque assay was employed to determine the viral titer. Immunofluorescence and real-time quantitative polymerase chain reaction were used to measure the viral RNA levels. The in vitro-transcribed RNA transcript derived from a replication-dead Japanese encephalitis virus (JEV) replicon containing the renilla luciferase reporter gene (J-R2A-NS5mt) was used to assess the stability of the flavivirus RNA. An RNA immunoprecipitation assay was used to detect the protein-RNA binding ability. Confocal microscopic images were captured to analyze protein colocalization.

RESULTS

ZFP36L2 served as an innate host defender against JEV and dengue virus. ZFP36L2 inhibited flavivirus infection solely through the 5'-3' XRN1 RNA decay pathway, whereas ZFP36L1 inhibited JEV infection via the 5'-3' XRN1 and 3'-5' RNA exosome RNA decay pathways. The direct binding between viral RNA and ZFP36L2 via its CCCH-type zinc finger motifs facilitated the degradation of flavivirus RNA mediated by 5'-3' XRN1. Furthermore, ZFP36L2 was localized in processing bodies (PBs), which participate in the 5'-3' XRN1-mediated RNA decay pathway. Nonetheless, the disruption of PBs did not affect the antiviral activity of ZFP36L2, suggesting that its localization is not essential for the function of the protein. Interestingly, the colocalization of ZFP36L2 and XRN1 with viral RNA and NS3 revealed that the antiviral activity of ZFP36L2 occurred within the replication complexes (RCs).

CONCLUSIONS

In summary, ZFP36L2 bound to and degraded viral RNA through the XRN1-mediated RNA decay pathway in the RCs, thereby inhibiting flavivirus replication. These findings provide valuable insights into the diverse antiviral mechanisms of the ZFP36-like family of proteins in the innate immune response against flavivirus infection.

Lifetime of ground conformational state determines the activity of structured RNA

Biomolecules continually sample alternative conformations. Consequently, even the most energetically favored ground conformational state has a finite lifetime. Here, we show that, in addition to the three-dimensional (3D) structure, the lifetime of a ground conformational state determines its biological activity. Using hydrogen-deuterium exchange nuclear magnetic resonance spectroscopy, we found that Zika virus exoribonuclease-resistant RNA (xrRNA) encodes a ground conformational state with a lifetime that is ~105-107 longer than that of canonical base pairs. Mutations that shorten the apparent lifetime of the ground state without affecting its 3D structure decreased exoribonuclease resistance in vitro and impaired virus replication in cells. Additionally, we observed this exceptionally long-lived ground state in xrRNAs from diverse infectious mosquito-borne flaviviruses. These results demonstrate the biological importance of the lifetime of a preorganized ground state and further suggest that elucidating the lifetimes of dominant 3D structures of biomolecules may be crucial for understanding their behaviors and functions.

Purification of Enzymatically Active Xrn1 for Removal of Non-capped mRNAs from In Vitro Transcription Reactions and Evaluation of mRNA Decapping Status In Vivo

The cap is a 7-methylguanosine attached to the first messenger RNA (mRNA) nucleotide with a 5'-5' triphosphate bridge. This conserved eukaryotic modification confers stability to the transcripts and is essential for translation initiation. The specific mechanisms that govern transcript cytoplasmic longevity and translatability were always of substantial interest. Multiple works aimed at modeling mRNA decay mechanisms, including the onset of decapping, which is the rate-limiting step of mRNA decay. Additionally, with the recent advances in RNA-based vaccines, the importance of efficient synthesis of fully functional mRNAs has increased. Non-capped mRNAs arising during in vitro transcription are highly immunogenic, and multiple approaches were developed to reduce their levels. Efficient and low-cost methods for elimination of non-capped mRNAs in vitro are therefore essential to basic sciences and to pharmaceutical applications. Here, we present a protocol for heterologous expression and purification of catalytically active recombinant Xrn1 from Thermothelomyces (Myceliophthora) thermophilus (Tt_Xrn1). We also describe protocols needed to verify the enzyme quality.

The pseudoknot structure of a viral RNA reveals a conserved mechanism for programmed exoribonuclease resistance

Exoribonuclease-resistant RNAs (xrRNAs) are viral RNA structures that block degradation by cellular 5'-3' exoribonucleases to produce subgenomic viral RNAs during infection. Initially discovered in flaviviruses, xrRNAs have since been identified in wide range of RNA viruses, including those that infect plants. High sequence variability among viral xrRNAs raises questions about the shared molecular features that characterize this functional RNA class. Here, we present the first structure of a plant-virus xrRNA in its active exoribonuclease-resistant conformation. The xrRNA forms a 9 base pair pseudoknot that creates a knot-like topology similar to that of flavivirus xrRNAs, despite lacking sequence similarity. Biophysical assays confirm a compact pseudoknot structure in solution, and functional studies validate its relevance both in vitro and during infection. Our study reveals how viral RNAs achieve a common functional outcome through highly divergent sequences and identifies the knot-like topology as a defining feature of xrRNAs.

Subgenomic flaviviral RNAs and human proteins: in silico exploration of anti-host defense mechanisms

Flaviviruses pose significant global health threats, infecting over 300 million people annually. Among their evasion strategies, the production of subgenomic flaviviral RNAs (sfRNAs) from the 3' UTR of viral genomes is particularly notable. Utilizing a comprehensive in silico approach with the catRAPID algorithm, we analyzed over 300,000 interactions between sfRNAs and human proteins derived from more than 8000 flavivirus genomes, including Dengue, Zika, Yellow Fever, West Nile, and Japanese Encephalitis viruses. By providing the first extensive atlas of sfRNA interactions, we offer new insights into how flaviviruses can manipulate host cellular machinery to facilitate viral survival and persistence. Our study not only validated known interactions but also revealed novel human proteins that could be involved in sfRNA-mediated host defense evasion, including helicases, splicing factors, and chemokines. These findings significantly expand the known interactome of sfRNAs with human proteins, underscoring their role in modulating host cellular pathways. Intriguingly, we predict interaction with stress granules, a critical component of the cellular response to viral infection, suggesting a mechanism by which flaviviruses inhibit their formation to evade host defenses. Moreover, a set of highly-interacting proteins in common among the sfRNAs showed predictive power to identify sfRNA-forming regions, highlighting how protein signatures could be used to annotate viruses. This atlas not only serves as a resource for exploring therapeutic targets but also aids in the identification of sfRNA biomarkers for improved flavivirus diagnostics.

Pan-flavivirus analysis reveals sfRNA-independent, 3' UTR-biased siRNA production from an insect-specific flavivirus

RNA interference (RNAi) plays an essential role in mosquito antiviral immunity, but it is not known whether viral small interfering RNA (siRNA) profiles differ between mosquito-borne and mosquito-specific viruses. A pan-Orthoflavivirus analysis in Aedes albopictus cells revealed that viral siRNAs were evenly distributed across the viral genome of most representatives of the Flavivirus genus. In contrast, siRNA production was biased toward the 3' untranslated region (UTR) of the genomes of classical insect-specific flaviviruses (cISF), which was most pronounced for Kamiti River virus (KRV), a virus with a unique, 1.2 kb long 3' UTR. KRV-derived siRNAs were produced in high quantities and almost exclusively mapped to the 3' UTR. We mapped the 5' end of KRV subgenomic flavivirus RNAs (sfRNAs), products of the 5'-3' exoribonuclease XRN1/Pacman stalling on secondary RNA structures in the 3' UTR of the viral genome. We found that KRV produces high copy numbers of a long, 1,017 nt sfRNA1 and a short, 421 nt sfRNA2, corresponding to two predicted XRN1-resistant elements. Expression of both sfRNA1 and sfRNA2 was reduced in Pacman-deficient Aedes albopictus cells; however, this did not correlate with a shift in viral siRNA profiles. We suggest that cISFs, particularly KRV, developed a unique mechanism to produce high amounts of siRNAs as a decoy for the antiviral RNAi response in an sfRNA-independent manner.IMPORTANCEThe Flavivirus genus contains diverse mosquito viruses ranging from insect-specific viruses circulating exclusively in mosquito populations to mosquito-borne viruses that cause disease in humans and animals. Studying the mechanisms of virus replication and antiviral immunity in mosquitoes is important to understand arbovirus transmission and may inform the development of disease control strategies. In insects, RNA interference (RNAi) provides broad antiviral activity and constitutes a major immune response against viruses. Comparing diverse members of the Flavivirus genus, we found that all flaviviruses are targeted by RNAi. However, the insect-specific Kamiti River virus was unique in that small interfering RNAs are highly skewed toward its uniquely long 3' untranslated region. These results suggest that mosquito-specific viruses have evolved unique mechanisms for genome replication and immune evasion.

The role of structure in regulatory RNA elements

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e., RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Magnesium Regulates RNA Ring Dynamics and Folding in Subgenomic Flaviviral RNA

Mosquito-borne flaviviruses including dengue, Zika, yellow fever, and regional encephalitis produce a large amount of short subgenomic flaviviral RNAs during infection. A segment of these RNAs named as xrRNA1 features a multi-pseudoknot (PK)-associated structure, which resists the host cell enzyme (XRN1) from degrading the viral RNA. We investigate how this long-range RNA PK folds in the presence of counterions, specifically in a mix of monovalent (K+) and divalent (Mg2+) salts at physiological concentrations. In this study, we use extensive explicit solvent molecular dynamics (MD) simulations to characterize the RNA ion environment of the folded RNA conformation, as determined by the crystal structure. This allowed us to identify the precise locations of various coordinated RNA-Mg2+ interactions, including inner-sphere/chelated and outer-sphere coordinated Mg2+. Given that RNA folding involves large-scale conformational changes, making it challenging to explore through classical MD simulations, we investigate the folding mechanism of xrRNA1 using an all-atom structure-based RNA model with a hybrid implicit-explicit treatment of the ion environment via the dynamic counterion condensation model, both with and without physiological Mg2+ concentration. The study reveals potential folding pathways for this xrRNA1, which is consistent with the results obtained from optical tweezer experiments. The equilibrium and free energy simulations both capture a dynamic equilibrium between the ring-open and ring-close states of the RNA, driven by a long-range PK interaction. Free energy calculations reveal that with the addition of Mg2+ ions, the equilibrium shifts more toward the ring-close state. A detailed analysis of the free energy pathways and ion-mediated contact probability map highlights the critical role of Mg2+ in bridging G50 and A33. This Mg2+-mediated connection helps form the long-range PK which in turn controls the transition between the ring-open and ring-close states. The study underscores the critical role of Mg2+ in the RNA folding transition, highlighting specific locations of Mg2+ contributing to the stabilization of long-range PK connections likely to enhance the robustness of Xrn1 resistance of flaviviral xrRNAs.

A Comparative Overview of the Role of Human Ribonucleases in Nonsense-Mediated mRNA Decay

Eukaryotic cells possess surveillance mechanisms that detect and degrade defective transcripts. Aberrant transcripts include mRNAs with a premature termination codon (PTC), targeted by the nonsense-mediated decay (NMD) pathway, and mRNAs lacking a termination codon, targeted by the nonstop decay (NSD) pathway. The eukaryotic exosome, a ribonucleolytic complex, plays a crucial role in mRNA processing and turnover through its catalytic subunits PM/Scl100 (Rrp6 in yeast), DIS3 (Rrp44 in yeast), and DIS3L1. Additionally, eukaryotic cells have other ribonucleases, such as SMG6 and XRN1, that participate in RNA surveillance. However, the specific pathways through which ribonucleases recognize and degrade mRNAs remain elusive. In this study, we characterized the involvement of human ribonucleases, both nuclear and cytoplasmic, in the mRNA surveillance mechanisms of NMD and NSD. We performed knockdowns of SMG6, PM/Scl100, XRN1, DIS3, and DIS3L1, analyzing the resulting changes in mRNA levels of selected natural NMD targets by RT-qPCR. Additionally, we examined the levels of different human β-globin variants under the same conditions: wild-type, NMD-resistant, NMD-sensitive, and NSD-sensitive. Our results demonstrate that all the studied ribonucleases are involved in the decay of certain endogenous NMD targets. Furthermore, we observed that the ribonucleases SMG6 and DIS3 contribute to the degradation of all β-globin variants, with an exception for βNS in the former case. This is also the case for PM/Scl100, which affects all β-globin variants except the NMD-sensitive variants. In contrast, DIS3L1 and XRN1 show specificity for β-globin WT and NMD-resistant variants. These findings suggest that eukaryotic ribonucleases are target-specific rather than pathway-specific. In addition, our data suggest that ribonucleases play broader roles in mRNA surveillance and degradation mechanisms beyond just NMD and NSD.

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.

Sub-genomic flaviviral RNA elements increase the stability and abundance of recombinant AAV vector transcripts

Many viruses have evolved structured RNA elements that can influence transcript abundance and translational efficiency, and help evade host immune factors by hijacking cellular machinery during replication. Here, we evaluated the functional impact of sub-genomic flaviviral RNAs (sfRNAs) known to stall exoribonuclease activity, by incorporating these elements into recombinant adeno-associated viral (AAV) genome cassettes. Specifically, sfRNAs from Dengue, Zika, Japanese Encephalitis, Yellow Fever, Murray Valley Encephalitis, and West Nile viruses increased transcript stability and transgene expression compared to a conventional woodchuck hepatitis virus element (WPRE). Further dissection of engineered transcripts revealed that sfRNA elements (i) require incorporation in cis within the 3' untranslated region (UTR) of AAV genomes, (ii) require minimal dumbbell structures to exert the observed effects, and (iii) can stabilize AAV transcripts independent of 5'-3' exoribonuclease 1 (XRN1)-mediated decay. Additionally, preliminary in vivo assessment of AAV vectors bearing sfRNA elements in mice achieved increased transcript abundance and expression in cardiac tissue. Leveraging the functional versatility of engineered viral RNA elements may help improve the potency of AAV vector-based gene therapies.

IMPORTANCE

Viral RNA elements can hijack host cell machinery to control stability of transcripts and consequently, infection. Studies that help better understand such viral elements can provide insights into antiviral strategies and also potentially leverage these features for therapeutic applications. In this study, by incorporating structured flaviviral RNA elements into recombinant adeno-associated viral (AAV) vector genomes, we show improved AAV transcript stability and transgene expression can be achieved, with implications for gene transfer.

An enhanced Eco1 retron editor enables precision genome engineering in human cells from a single-copy integrated lentivirus

Retrons are a retroelement class found in diverse prokaryotes that can be adapted to augment CRISPR-Cas9 genome engineering technology to efficiently rewrite short stretches of genetic information in bacteria and yeast; however, efficiency in human cells has been limited by unknown factors. We identified non-coding RNA (ncRNA) instability and impaired Cas9 activity as major contributors to poor retron editor efficiency. We re-engineered the Eco1 ncRNA to incorporate an exoribonuclease-resistant RNA pseudoknot from the Zika virus 3' UTR and devised an RNA processing strategy using Csy4 ribonuclease to liberate the sgRNA and ncRNA. These modifications yielded a ncRNA with 5'- and 3'-end protection and an sgRNA with minimal 5' extension. This strategy increased steady-state ncRNA levels and rescued Cas9 activity leading to enhanced efficiency of the Eco1 retron editor in human cells. The enhanced Eco1 retron editor enabled the insertion of missense mutations in human cells from a single integrated lentivirus, thereby ensuring genotype-phenotype linkage over multiple cell divisions. This work reveals a previously unappreciated role for ncRNA stability in retron editor efficiency in human cells. Here we present an enhanced Eco1 retron editor that enables efficient introduction of missense mutations in human cells from a single heritable genome copy.

The West Nile virus genome harbors essential riboregulatory elements with conserved and host-specific functional roles

West Nile virus (WNV) is an arthropod-borne, positive-sense RNA virus that poses an increasing global threat due to warming climates and lack of effective therapeutics. Like other enzootic viruses, little is known about how host context affects the structure of the full-length RNA genome. Here, we report a complete secondary structure of the entire WNV genome within infected mammalian and arthropod cell lines. Our analysis affords structural insights into multiple, conserved aspects of flaviviral biology. We show that the WNV genome folds with minimal host dependence, and we prioritize well-folded regions for functional validation using structural homology between hosts as a guide. Using structure-disrupting, antisense locked nucleic acids, we then demonstrate that the WNV genome contains riboregulatory structures with conserved and host-specific functional roles. These results reveal promising RNA drug targets within flaviviral genomes, and they highlight the therapeutic potential of ASO-LNAs as both WNV-specific and pan-flaviviral therapeutic agents.

West Nile Virus Subgenomic RNAs Modulate Gene Expression in a Neuronal Cell Line

Subgenomic flaviviral RNAs (sfRNAs) are small non-coding products of the incomplete degradation of viral genomic RNA. They accumulate during flaviviral infection and have been associated with many functional roles inside the host cell. Studies so far have demonstrated that sfRNA plays a crucial role in determining West Nile virus (WNV) pathogenicity. However, its modulatory role on neuronal homeostasis has not been studied in depth. In this study, we investigated the mechanism of sfRNA biosynthesis and its importance for WNV replication in neuronal cells. We found that sfRNA1 is functionally redundant for both replication and translation of WNV. However, the concurrent absence of sfRNA1 and sfRNA2 species is detrimental for the survival of the virus. Differential expression analysis on RNA-seq data from WT and ΔsfRNA replicon cell lines revealed transcriptional changes induced by sfRNA and identified a number of putative targets. Overall, it was shown that sfRNA contributes to the viral evasion by suppressing the interferon-mediated antiviral response. An additional differential expression analysis among replicon and control Neuro2A cells also clarified the transcriptional changes that support WNV replication in neuronal cells. Increased levels of translation and oxidative phosphorylation, post-translational modification processes, and activated DNA repair pathways were observed in replicon cell lines, while developmental processes such as axonal growth were deficient.

Xinyang flavivirus, from Haemaphysalis flava ticks in Henan Province, China, defines a basal, likely tick-only Orthoflavivirus clade

Tick-borne orthoflaviviruses (TBFs) are classified into three conventional groups based on genetics and ecology: mammalian, seabird and probable-TBF group. Recently, a fourth basal group has been identified in Rhipicephalus ticks from Africa: Mpulungu flavivirus (MPFV) in Zambia and Ngoye virus (NGOV) in Senegal. Despite attempts, isolating these viruses in vertebrate and invertebrate cell lines or intracerebral injection of newborn mice with virus-containing homogenates has remained unsuccessful. In this study, we report the discovery of Xinyang flavivirus (XiFV) in Haemaphysalis flava ticks from Xìnyáng, Henan Province, China. Phylogenetic analysis shows that XiFV was most closely related to MPFV and NGOV, marking the first identification of this tick orthoflavivirus group in Asia. We developed a reverse transcriptase quantitative PCR assay to screen wild-collected ticks and egg clutches, with absolute infection rates of 20.75 % in adult females and 15.19 % in egg clutches, suggesting that XiFV could be potentially spread through transovarial transmission. To examine potential host range, dinucleotide composition analyses revealed that XiFV, MPFV and NGOV share a closer composition to classical insect-specific orthoflaviviruses than to vertebrate-infecting TBFs, suggesting that XiFV could be a tick-only orthoflavivirus. Additionally, both XiFV and MPFV lack a furin cleavage site in the prM protein, unlike other TBFs, suggesting these viruses might exist towards a biased immature particle state. To examine this, chimeric Binjari virus with XIFV-prME (bXiFV) was generated, purified and analysed by SDS-PAGE and negative-stain transmission electron microscopy, suggesting prototypical orthoflavivirus size (~50 nm) and bias towards uncleaved prM. In silico structural analyses of the 3'-untranslated regions show that XiFV forms up to five pseudo-knot-containing stem-loops and a prototypical orthoflavivirus dumbbell element, suggesting the potential for multiple exoribonuclease-resistant RNA structures.

Intermolecular base-pairing interactions, a unique topology and exoribonuclease-resistant noncoding RNAs drive formation of viral chimeric RNAs in plants

In plants, exoribonuclease-resistant RNAs (xrRNAs) are produced by many viruses. Whereas xrRNAs contribute to the pathogenicity of these viruses, the role of xrRNAs in the virus infectious cycle remains elusive. Here, we show that xrRNAs produced by a benyvirus (a multipartite RNA virus with four genomic segments) in plants are involved in the formation of monocistronic coat protein (CP)-encoding chimeric RNAs. Naturally occurring chimeric RNAs, we discovered, are composed of 5'-end of RNA 2 and 3'-end of either RNA 3 or RNA 4 bearing conservative exoribonuclease-resistant 'coremin' region. Using computational tools and site-directed mutagenesis, we show that de novo formation of chimeric RNAs requires intermolecular base-pairing interaction between 'coremin' and 3'-proximal part of the CP gene of RNA 2 as well as a stem-loop structure immediately adjacent to the CP gene. Moreover, knockdown of the expression of the XRN4 gene, encoding 5'→3' exoribonuclease, inhibits biogenesis of both xrRNAs and chimeric RNAs. Our findings suggest a novel mechanism involving a unique tropology of the intermolecular base-pairing complex between xrRNAs and RNA2 to promote formation of chimeric RNAs in plants. XrRNAs, essential for chimeric RNA biogenesis, are generated through the action of cytoplasmic Xrn 4 5'→3' exoribonuclease conserved in all plant species.

A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

Subgenomic Flaviviral RNAs of Dengue Viruses

Subgenomic flaviviral RNAs (sfRNAs) are produced during flavivirus infections in both arthropod and vertebrate cells. They are undegraded products originating from the viral 3' untranslated region (3' UTR), a result of the action of the host 5'-3' exoribonuclease, Xrn1, when it encounters specific RNA structures known as Xrn1-resistant RNAs (xrRNAs) within the viral 3' UTR. Dengue viruses generate three to four distinct species of sfRNAs through the presence of two xrRNAs and two dumbbell structures (DBs). The tertiary structures of xrRNAs have been characterized to form a ringlike structure around the 5' end of the viral RNA, effectively inhibiting the activity of Xrn1. The most important role of DENV sfRNAs is to inhibit host antiviral responses by interacting with viral and host proteins, thereby influencing viral pathogenicity, replicative fitness, epidemiological fitness, and transmission. In this review, we aimed to summarize the biogenesis, structures, and functions of DENV sfRNAs, exploring their implications for viral interference.

SKI complex: A multifaceted cytoplasmic RNA exosome cofactor in mRNA metabolism with links to disease, developmental processes, and antiviral responses

RNA stability and quality control are integral parts of gene expression regulation. A key factor shaping eukaryotic transcriptomes, mainly via 3'-5' exoribonucleolytic trimming or degradation of diverse transcripts in nuclear and cytoplasmic compartments, is the RNA exosome. Precise exosome targeting to various RNA molecules requires strict collaboration with specialized auxiliary factors, which facilitate interactions with its substrates. The predominant class of cytoplasmic RNA targeted by the exosome are protein-coding transcripts, which are carefully scrutinized for errors during translation. Normal, functional mRNAs are turned over following protein synthesis by the exosome or by Xrn1 5'-3'-exonuclease, acting in concert with Dcp1/2 decapping complex. In turn, aberrant transcripts are eliminated by dedicated surveillance pathways, triggered whenever ribosome translocation is impaired. Cytoplasmic 3'-5' mRNA decay and surveillance are dependent on the tight cooperation between the exosome and its evolutionary conserved co-factor-the SKI (superkiller) complex (SKIc). Here, we summarize recent findings from structural, biochemical, and functional studies of SKIc roles in controlling cytoplasmic RNA metabolism, including links to various cellular processes. Mechanism of SKIc action is illuminated by presentation of its spatial structure and details of its interactions with exosome and ribosome. Furthermore, contribution of SKIc and exosome to various mRNA decay pathways, usually converging on recycling of ribosomal subunits, is delineated. A crucial physiological role of SKIc is emphasized by describing association between its dysfunction and devastating human disease-a trichohepatoenteric syndrome (THES). Eventually, we discuss SKIc functions in the regulation of antiviral defense systems, cell signaling and developmental transitions, emerging from interdisciplinary investigations. This article is categorized under: RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Turnover and Surveillance > Regulation of RNA Stability RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Xrn1-resistant RNA motifs are disseminated throughout the RNA virome and are able to block scanning ribosomes

RNAs that are able to prevent degradation by the 5'-3' exoribonuclease Xrn1 have emerged as crucial structures during infection by an increasing number of RNA viruses. Several plant viruses employ the so-called coremin motif, an Xrn1-resistant RNA that is usually located in 3' untranslated regions. Investigation of its structural and sequence requirements has led to its identification in plant virus families beyond those in which the coremin motif was initially discovered. In this study, we identified coremin-like motifs that deviate from the original in the number of nucleotides present in the loop region of the 5' proximal hairpin. They are present in a number of viral families that previously did not have an Xrn1-resistant RNA identified yet, including the double-stranded RNA virus families Hypoviridae and Chrysoviridae. Through systematic mutational analysis, we demonstrated that a coremin motif carrying a 6-nucleotide loop in the 5' proximal hairpin generally requires a YGNNAD consensus for stalling Xrn1, similar to the previously determined YGAD consensus required for Xrn1 resistance of the original coremin motif. Furthermore, we determined the minimal requirements for the 3' proximal hairpin. Since some putative coremin motifs were found in intergenic regions or coding sequences, we demonstrated their capacity for inhibiting translation through an in vitro ribosomal scanning inhibition assay. Consequently, this study provides a further expansion on the number of viral families with known Xrn1-resistant elements, while adding a novel, potentially regulatory function for this structure.

RNA editing: Expanding the potential of RNA therapeutics

RNA therapeutics have had a tremendous impact on medicine, recently exemplified by the rapid development and deployment of mRNA vaccines to combat the COVID-19 pandemic. In addition, RNA-targeting drugs have been developed for diseases with significant unmet medical needs through selective mRNA knockdown or modulation of pre-mRNA splicing. Recently, RNA editing, particularly antisense RNA-guided adenosine deaminase acting on RNA (ADAR)-based programmable A-to-I editing, has emerged as a powerful tool to manipulate RNA to enable correction of disease-causing mutations and modulate gene expression and protein function. Beyond correcting pathogenic mutations, the technology is particularly well suited for therapeutic applications that require a transient pharmacodynamic effect, such as the treatment of acute pain, obesity, viral infection, and inflammation, where it would be undesirable to introduce permanent alterations to the genome. Furthermore, transient modulation of protein function, such as altering the active sites of enzymes or the interface of protein-protein interactions, opens the door to therapeutic avenues ranging from regenerative medicine to oncology. These emerging RNA-editing-based toolsets are poised to broadly impact biotechnology and therapeutic applications. Here, we review the emerging field of therapeutic RNA editing, highlight recent laboratory advancements, and discuss the key challenges on the path to clinical development.

Nanopore sequencing of internal 2'-PO4 modifications installed by RNA repair

Ligation by plant and fungal RNA ligases yields an internal 2'-phosphate group on each RNA ligation product. In budding yeast, this covalent mark occurs at the splice junction of two targets of ligation: intron-containing tRNAs and the messenger RNA HAC1 The repertoire of RNA molecules repaired by RNA ligation has not been explored due to a lack of unbiased approaches for identifying RNA ligation products. Here, we define several unique signals produced by 2'-phosphorylated RNAs during nanopore sequencing. A 2'-phosphate at the splice junction of HAC1 mRNA inhibits 5' → 3' degradation, enabling detection of decay intermediates in yeast RNA repair mutants by nanopore sequencing. During direct RNA sequencing, intact 2'-phosphorylated RNAs on HAC1 and tRNAs produce diagnostic changes in nanopore current properties and base calling features, including stalls produced as the modified RNA translocates through the nanopore motor protein. These approaches enable directed studies to identify novel RNA repair events in other contexts.

Isolation and Genetic Characterization of Japanese Encephalitis Virus Two Decades after Its Elimination in Singapore

Japanese encephalitis virus (JEV) is an important arbovirus in Asia that can cause serious neurological disease. JEV is transmitted by mosquitoes in an enzootic cycle involving porcine and avian reservoirs, in which humans are accidental, dead-end hosts. JEV is currently not endemic in Singapore, after pig farming was abolished in 1992; the last known human case was reported in 2005. However, due to its location along the East-Asian Australasian Flyway (EAAF), Singapore is vulnerable to JEV re-introduction from the endemic regions. Serological and genetic evidence in the last decade suggests JEV's presence in the local fauna. In the present study, we report the genetic characterization and the first isolation of JEV from 3214 mosquito pools consisting of 41,843 Culex mosquitoes, which were trapped from April 2014 to May 2021. The findings demonstrated the presence of genotype I of JEV (n = 10), in contrast to the previous reports of the presence of genotype II of JEV in Singapore. The genetic analyses also suggested that JEV has entered Singapore on several occasions and has potentially established an enzootic cycle in the local fauna. These observations have important implications in the risk assessment and the control of Japanese encephalitis in non-endemic countries, such as Singapore, that are at risk for JEV transmission.

Characterization of dengue virus 3'UTR RNA binding proteins in mosquitoes reveals that AeStaufen reduces subgenomic flaviviral RNA in saliva

Dengue viruses (DENV) are expanding global pathogens that are transmitted through the bite of mosquitoes, mostly Aedes aegypti. As RNA viruses, DENV rely on RNA-binding proteins (RBPs) to complete their life cycle. Alternatively, RBPs can act as restriction factors that prevent DENV multiplication. While the importance of RBPs is well-supported in humans, there is a dearth of information about their influence on DENV transmission by mosquitoes. Such knowledge could be harnessed to design novel, effective interventions against DENV. Here, we successfully adapted RNA-affinity chromatography coupled with mass spectrometry-a technique initially developed in mammalian cells-to identify RBPs in Ae. aegypti cells. We identified fourteen RBPs interacting with DENV serotype 2 3'UTR, which is involved in the viral multiplication and produces subgenomic flaviviral RNA (sfRNA). We validated the RNA affinity results for two RBPs by confirming that AePur binds the 3'UTR, whereas AeStaufen interacts with both 3'UTR and sfRNA. Using in vivo functional evaluation, we determined that RBPs like AeRan, AeExoRNase, and AeRNase have pro-viral functions, whereas AeGTPase, AeAtu, and AePur have anti-viral functions in mosquitoes. Furthermore, we showed that human and mosquito Pur homologs have a shared affinity to DENV2 RNA, although the anti-viral effect is specific to the mosquito protein. Importantly, we revealed that AeStaufen mediates a reduction of gRNA and sfRNA copies in several mosquito tissues, including the salivary glands and that AeStaufen-mediated sfRNA reduction diminishes the concentration of transmission-enhancing sfRNA in saliva, thereby revealing AeStaufen's role in DENV transmission. By characterizing the first RBPs that associate with DENV2 3'UTR in mosquitoes, our study unravels new pro- and anti-viral targets for the design of novel therapeutic interventions as well as provides foundation for studying the role of RBPs in virus-vector interactions.

Evolutionary traits of Tick-borne encephalitis virus: Pervasive non-coding RNA structure conservation and molecular epidemiology

Tick-borne encephalitis virus (TBEV) is the aetiological agent of tick-borne encephalitis, an infectious disease of the central nervous system that is often associated with severe sequelae in humans. While TBEV is typically classified into three subtypes, recent evidence suggests a more varied range of TBEV subtypes and lineages that differ substantially in the architecture of their 3ʹ untranslated region (3ʹUTR). Building on comparative genomic approaches and thermodynamic modelling, we characterize the TBEV UTR structureome diversity and propose a unified picture of pervasive non-coding RNA structure conservation. Moreover, we provide an updated phylogeny of TBEV, building on more than 220 publicly available complete genomes, and investigate the molecular epidemiology and phylodynamics with Nextstrain, a web-based visualization framework for real-time pathogen evolution.

The In Vivo and In Vitro Architecture of the Hepatitis C Virus RNA Genome Uncovers Functional RNA Secondary and Tertiary Structures

Hepatitis C virus (HCV) is a positive-strand RNA virus that remains one of the main contributors to chronic liver disease worldwide. Studies over the last 30 years have demonstrated that HCV contains a highly structured RNA genome and many of these structures play essential roles in the HCV life cycle. Despite the importance of riboregulation in this virus, most of the HCV RNA genome remains functionally unstudied. Here, we report a complete secondary structure map of the HCV RNA genome in vivo, which was studied in parallel with the secondary structure of the same RNA obtained in vitro. Our results show that HCV is folded extensively in the cellular context. By performing comprehensive structural analyses on both in vivo data and in vitro data, we identify compact and conserved secondary and tertiary structures throughout the genome. Genetic and evolutionary functional analyses demonstrate that many of these elements play important roles in the virus life cycle. In addition to providing a comprehensive map of RNA structures and riboregulatory elements in HCV, this work provides a resource for future studies aimed at identifying therapeutic targets and conducting further mechanistic studies on this important human pathogen. IMPORTANCE HCV has one of the most highly structured RNA genomes studied to date, and it is a valuable model system for studying the role of RNA structure in protein-coding genes. While previous studies have identified individual cases of regulatory RNA structures within the HCV genome, the full-length structure of the HCV genome has not been determined in vivo. Here, we present the complete secondary structure map of HCV determined both in cells and from corresponding transcripts generated in vitro. In addition to providing a comprehensive atlas of functional secondary structural elements throughout the genomic RNA, we identified a novel set of tertiary interactions and demonstrated their functional importance. In terms of broader implications, the pipeline developed in this study can be applied to other long RNAs, such as long noncoding RNAs. In addition, the RNA structural motifs characterized in this study broaden the repertoire of known riboregulatory elements.

Enhancement of prime editing via xrRNA motif-joined pegRNA

The prime editors (PEs) have shown great promise for precise genome modification. However, their suboptimal efficiencies present a significant technical challenge. Here, by appending a viral exoribonuclease-resistant RNA motif (xrRNA) to the 3'-extended portion of pegRNAs for their increased resistance against degradation, we develop an upgraded PE platform (xrPE) with substantially enhanced editing efficiencies in multiple cell lines. A pan-target average enhancement of up to 3.1-, 4.5- and 2.5-fold in given cell types is observed for base conversions, small deletions, and small insertions, respectively. Additionally, xrPE exhibits comparable edit:indel ratios and similarly minimal off-target editing as the canonical PE3. Of note, parallel comparison of xrPE to the most recently developed epegRNA-based PE system shows their largely equivalent editing performances. Our study establishes a highly adaptable platform of improved PE that shall have broad implications.

Structural analysis of 3'UTRs in insect flaviviruses reveals novel determinants of sfRNA biogenesis and provides new insights into flavivirus evolution

Subgenomic flaviviral RNAs (sfRNAs) are virus-derived noncoding RNAs produced by pathogenic mosquito-borne flaviviruses (MBF) to counteract the host antiviral response. To date, the ability of non-pathogenic flaviviruses to produce and utilise sfRNAs remains largely unexplored, and it is unclear what role XRN1 resistance plays in flavivirus evolution and host adaptation. Herein the production of sfRNAs by several insect-specific flaviviruses (ISFs) that replicate exclusively in mosquitoes is shown, and the secondary structures of their complete 3'UTRs are determined. The xrRNAs responsible for the biogenesis of ISF sfRNAs are also identified, and the role of these sfRNAs in virus replication is demonstrated. We demonstrate that 3'UTRs of all classical ISFs, except Anopheles spp-asscoaited viruses, and of the dual-host associated ISF Binjari virus contain duplicated xrRNAs. We also reveal novel structural elements in the 3'UTRs of dual host-associated and Anopheles-associated classical ISFs. Structure-based phylogenetic analysis demonstrates that xrRNAs identified in Anopheles spp-associated ISF are likely ancestral to xrRNAs of ISFs and MBFs. In addition, our data provide evidence that duplicated xrRNAs are selected in the evolution of flaviviruses to provide functional redundancy, which preserves the production of sfRNAs if one of the structures is disabled by mutations or misfolding.

Intron-encoded cistronic transcripts for minimally invasive monitoring of coding and non-coding RNAs

Despite their fundamental role in assessing (patho)physiological cell states, conventional gene reporters can follow gene expression but leave scars on the proteins or substantially alter the mature messenger RNA. Multi-time-point measurements of non-coding RNAs are currently impossible without modifying their nucleotide sequence, which can alter their native function, half-life and localization. Thus, we developed the intron-encoded scarless programmable extranuclear cistronic transcript (INSPECT) as a minimally invasive transcriptional reporter embedded within an intron of a gene of interest. Post-transcriptional excision of INSPECT results in the mature endogenous RNA without sequence alterations and an additional engineered transcript that leaves the nucleus by hijacking the nuclear export machinery for subsequent translation into a reporter or effector protein. We showcase its use in monitoring interleukin-2 (IL2) after T cell activation and tracking the transcriptional dynamics of the long non-coding RNA (lncRNA) NEAT1 during CRISPR interference-mediated perturbation. INSPECT is a method for monitoring gene transcription without altering the mature lncRNA or messenger RNA of the target of interest.

Effects of the noncoding subgenomic RNA of red clover necrotic mosaic virus in virus infection

In recent years, a new class of viral noncoding subgenomic RNA (ncsgRNA) has been identified. This RNA is generated as a stable degradation product via an exoribonuclease-resistant RNA (xrRNA) structure, which blocks the progression of 5′→3′ exoribonuclease on viral RNAs in infected cells. Here, we assess the effects of the ncsgRNA of red clover necrotic mosaic virus (RCNMV), called SR1f, in infected plants. We demonstrate the following: (i) the absence of SR1f reduces symptoms and decreases viral RNA accumulation in Nicotiana benthamiana and Arabidopsis thaliana plants; (ii) SR1f has an essential function other than suppression of RNA silencing; and (iii) the cytoplasmic exoribonuclease involved in mRNA turnover, XRN4, is not required for SR1f production or virus infection. A comparative transcriptomic analysis in N. benthamiana infected with wild-type RCNMV or an SR1f-deficient mutant RCNMV revealed that wild-type RCNMV infection, which produces SR1f and much higher levels of virus, has a greater and more significant impact on cellular gene expression than the SR1f-deficient mutant. Upregulated pathways include plant hormone signaling, plant-pathogen interaction, MAPK signaling, and several metabolic pathways, while photosynthesis-related genes were downregulated. We compare this to host genes known to participate in infection by other tombusvirids. Viral reads revealed a 10- to 100-fold ratio of positive to negative strand, and the abundance of reads of both strands mapping to the 3′ region of RCNMV RNA1 support the premature transcription termination mechanism of synthesis for the coding sgRNA. These results provide a framework for future studies of the interactions and functions of noncoding RNAs of plant viruses.

IMPORTANCE Knowledge of how RNA viruses manipulate host and viral gene expression is crucial to our understanding of infection and disease. Unlike viral protein-host interactions, little is known about the control of gene expression by viral RNA. Here, we begin to address this question by investigating the noncoding subgenomic RNA (ncsgRNA) of red clover necrotic mosaic virus (RCNMV), called SR1f. Similar exoribonuclease-resistant RNAs of flaviviruses are well studied, but the roles of plant viral ncsgRNAs, and how they arise, are poorly understood. Surprisingly, we find the likely exonuclease candidate, XRN4, is not required to generate SR1f, and we assess the effects of SR1f on virus accumulation and symptom development. Finally, we compare the effects of infection by wild-type RCNMV versus an SR1f-deficient mutant on host gene expression in Nicotiana benthamiana, which reveals that ncsgRNAs such as SR1f are key players in virus-host interactions to facilitate productive infection.

Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs

To understand how RNA dynamics is regulated and connected to its function, we investigate the folding, conformational dynamics and robustness of Xrn1 resistance of a set of flaviviral xrRNAs using SAXS, smFRET and in vitro enzymatic assays. Flaviviral xrRNAs form discrete ring-like 3D structures, in which the length of a conserved long-range pseudoknot (PK2) ranges from 2 bp to 7 bp. We find that xrRNAs' folding, conformational dynamics and Xrn1 resistance are strongly correlated and highly Mg²⁺-dependent, furthermore, the Mg²⁺-dependence is modulated by PK2 length variations. xrRNAs with long PK2 require less Mg²⁺ to stabilize their folding, exhibit reduced conformational dynamics and strong Xrn1 resistance even at low Mg²⁺, and tolerate mutations at key tertiary motifs at high Mg²⁺, which generally are destructive to xrRNAs with short PK2. These results demonstrate an unusual regulatory mechanism of RNA dynamics providing insights into the functions and future biomedical applications of xrRNAs.

Chimeric Flaviviral RNA-siRNA Molecules Resist Degradation by The Exoribonuclease Xrn1 and Trigger Gene Silencing in Mammalian Cells

RNA is an emerging platform for drug delivery, but the susceptibility of RNA to nuclease degradation remains a major barrier to its implementation in vivo. Here, we engineered flaviviral Xrn1-resistant RNA (xrRNA) motifs to host small interfering RNA (siRNA) duplexes. The xrRNA-siRNA molecules self-assemble in vitro, resist degradation by the conserved eukaryotic 5' to 3' exoribonuclease Xrn1, and trigger gene silencing in 293T cells. The resistance of the molecules to Xrn1 does not translate to stability in blood serum. Nevertheless, our results demonstrate that flavivirus-derived xrRNA motifs can confer Xrn1 resistance on a model therapeutic payload and set the stage for further investigations into using the motifs as building blocks in RNA nanotechnology.

Replication-Dependent Biogenesis of Turnip Crinkle Virus Long Noncoding RNAs

Long noncoding RNAs (lncRNAs) of virus origin accumulate in cells infected by many positive-strand (+) RNA viruses to bolster viral infectivity. Their biogenesis mostly utilizes exoribonucleases of host cells that degrade viral genomic or subgenomic RNAs in the 5'-to-3' direction until being stalled by well-defined RNA structures. Here, we report a viral lncRNA that is produced by a novel replication-dependent mechanism. This lncRNA corresponds to the last 283 nucleotides of the turnip crinkle virus (TCV) genome and hence is designated tiny TCV subgenomic RNA (ttsgR). ttsgR accumulated to high levels in TCV-infected Nicotiana benthamiana cells when the TCV-encoded RNA-dependent RNA polymerase (RdRp), also known as p88, was overexpressed. Both (+) and (-) strand forms of ttsgR were produced in a manner dependent on the RdRp functionality. Strikingly, templates as short as ttsgR itself were sufficient to program ttsgR amplification, as long as the TCV-encoded replication proteins p28 and p88 were provided in trans. Consistent with its replicational origin, ttsgR accumulation required a 5' terminal carmovirus consensus sequence (CCS), a sequence motif shared by genomic and subgenomic RNAs of many viruses phylogenetically related to TCV. More importantly, introducing a new CCS motif elsewhere in the TCV genome was alone sufficient to cause the emergence of another lncRNA. Finally, abolishing ttsgR by mutating its 5' CCS gave rise to a TCV mutant that failed to compete with wild-type TCV in Arabidopsis. Collectively, our results unveil a replication-dependent mechanism for the biogenesis of viral lncRNAs, thus suggesting that multiple mechanisms, individually or in combination, may be responsible for viral lncRNA production. IMPORTANCE Many positive-strand (+) RNA viruses produce long noncoding RNAs (lncRNAs) during the process of cellular infections and mobilize these lncRNAs to counteract antiviral defenses, as well as coordinate the translation of viral proteins. Most viral lncRNAs arise from 5'-to-3' degradation of longer viral RNAs being stalled at stable secondary structures. Here, we report a viral lncRNA that is produced by the replication machinery of turnip crinkle virus (TCV). This lncRNA, designated ttsgR, shares the terminal characteristics with TCV genomic and subgenomic RNAs and overaccumulates in the presence of moderately overexpressed TCV RNA-dependent RNA polymerase (RdRp). Furthermore, templates that are of similar sizes as ttsgR are readily replicated by TCV replication proteins (p28 and RdRp) provided from nonviral sources. In summary, this study establishes an approach for uncovering low abundance viral lncRNAs, and characterizes a replicating TCV lncRNA. Similar investigations on human-pathogenic (+) RNA viruses could yield novel therapeutic targets.

Shared properties and singularities of exoribonuclease-resistant RNAs in viruses

What viral RNA genomes lack in size, they make up for in intricacy. Elaborate RNA structures embedded in viral genomes can hijack essential cellular mechanisms aiding virus propagation. Exoribonuclease-resistant RNAs (xrRNAs) are an emerging class of viral elements, which resist degradation by host cellular exoribonucleases to produce viral RNAs with diverse roles during infection. Detailed three-dimensional structural studies of xrRNAs from flaviviruses and a subset of plant viruses led to a mechanistic model in which xrRNAs block enzymatic digestion using a ring-like structure that encircles the 5' end of the resistant structure. In this mini-review, we describe the state of our understanding of the phylogenetic distribution of xrRNAs, their structures, and their conformational dynamics. Because xrRNAs have now been found in several major superfamilies of RNA viruses, they may represent a more widely used strategy than currently appreciated. Could xrRNAs represent a 'molecular clock' that would help us understand virus evolution and pathogenicity? The more we study xrRNAs in viruses, the closer we get to finding xrRNAs within cellular RNAs.

Mechanical strength of RNA knot in Zika virus protects against cellular defenses

Unusual knot-like structures recently discovered in viral exoribonuclease-resistant RNAs (xrRNAs) prevent digestion by host RNases to create subgenomic RNAs enhancing infection and pathogenicity. xrRNAs are proposed to prevent digestion through mechanical resistance to unfolding. However, their unfolding force has not been measured, and the factors determining RNase resistance are unclear. Furthermore, how these knots fold remains unknown. Unfolding a Zika virus xrRNA with optical tweezers revealed that it was the most mechanically stable RNA yet observed. The knot formed by threading the 5' end into a three-helix junction before pseudoknot interactions closed a ring around it. The pseudoknot and tertiary contacts stabilizing the threaded 5' end were both required to generate extreme force resistance, whereas removing a 5'-end contact produced a low-force knot lacking RNase resistance. These results indicate mechanical resistance plays a central functional role, with the fraction of molecules forming extremely high-force knots determining the RNase resistance level.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions

Motivation

The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3′ untranslated regions (3′UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies.

Results

We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3′UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3′UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses.

Availability and implementation

The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS.

Supplementary information

Supplementary data are available at Bioinformatics online.

High-Performance Allosteric Conditional Guide RNAs for Mammalian Cell-Selective Regulation of CRISPR/Cas

The activity of a conditional guide RNA (cgRNA) is dependent on the presence or absence of an RNA trigger, enabling cell-selective regulation of CRISPR/Cas function. cgRNAs are programmable at two levels, with the target-binding sequence controlling the target of Cas activity (edit, silence, or induce a gene of choice) and the trigger-binding sequence controlling the scope of Cas activity (subset of cells expressing the trigger RNA). Allosteric cgRNA mechanisms enable independent design of the target and trigger sequences, providing the flexibility to select the regulatory target and scope independently. Building on prior advances in dynamic RNA nanotechnology that demonstrated the cgRNA concept, here we set the goal of engineering high-performance allosteric cgRNA mechanisms for the mammalian setting, pursuing both ON → OFF logic (conditional inactivation by an RNA trigger) and OFF → ON logic (conditional activation by an RNA trigger). For each mechanism, libraries of orthogonal cgRNA/trigger pairs were designed using NUPACK. In HEK 293T cells expressing cgRNAs, triggers, and inducing dCas9: (1) a library of four ON → OFF "terminator switch" cgRNAs exhibit a median fold-change of ≈50×, a median fractional dynamic range of ≈20%, and a median crosstalk modulus of ≈9%; (2) a library of three OFF → ON "split-terminator switch" cgRNAs exhibit a median fold-change of ≈150×, a median fractional dynamic range of ≈50%, and a median crosstalk modulus of ≈4%. Further, we demonstrate that xrRNA elements that protect viral RNAs from degradation by exoribonucleases can dramatically enhance the performance of RNA synthetic biology. The high-performance allosteric cgRNAs demonstrated here for ON → OFF and OFF → ON logic in mammalian cells provide a foundation for pursuing applications of programmable cell-selective regulation.

Identification of host factors binding to dengue and Zika virus subgenomic RNA by efficient yeast three-hybrid screens of the human ORFeome

Flaviviruses such as the dengue (DENV) and the Zika virus (ZIKV) are important human pathogens causing around 100 million symptomatic infections each year. During infection, small subgenomic flavivirus RNAs (sfRNAs) are formed inside the infected host cell as a result of incomplete degradation of the viral RNA genome by cellular exoribonuclease XRN1. Although the full extent of sfRNA functions is to be revealed, these non-coding RNAs are key virulence factors and their detrimental effects on multiple cellular processes seem to consistently involve molecular interactions with RNA-binding proteins (RBPs). Discovery of such sfRNA-binding host-factors has followed established biochemical pull-down approaches skewed towards highly abundant proteins hampering proteome-wide coverage. Yeast three-hybrid (Y3H) systems represent an attractive alternative approach. To facilitate proteome-wide screens for RBP, we revisited and improved existing RNA-Y3H methodology by (1) implementing full-length ORF libraries in combination with (2) efficient yeast mating to increase screening depth and sensitivity, and (3) stringent negative controls to eliminate over-representation of non-specific RNA-binders. These improvements were validated employing the well-characterized interaction between DDX6 (DEAD-box helicase 6) and sfRNA of DENV as paradigm. Our advanced Y3H system was used to screen for human proteins binding to DENV and ZIKV sfRNA, resulting in a list of 69 putative sfRNA-binders, including several previously reported as well as numerous novel RBP host factors. Our methodology requiring no sophisticated infrastructure or analytic pipeline may be employed for the discovery of meaningful RNA-protein interactions at large scale in other fields.

Xrn1-resistant RNA structures are well-conserved within the genus flavivirus

Subgenomic RNAs are produced by several RNA viruses through incomplete degradation of their genomic RNA by the exoribonuclease Xrn1, and have been shown to be essential for viral growth and pathogenicity. Within the flavivirus genus of the Flaviviridae family, two distinct classes of Xrn1-resistant RNA motifs have been proposed; one for mosquito-borne and insect-specific flaviviruses, and one for tick-borne flaviviruses and no-known-vector flaviviruses. We investigated tick-borne and no-known-vector flavivirus Xrn1-resistant RNA motifs through systematic in vitro mutational analysis and showed that both classes actually possess very similar structural configurations, including a double pseudoknot and a base-triple at identical, conserved locations. For the no-known-vector flavivirus Modoc virus, we show that in vivo generation of subgenomic flaviviral RNA was affected by mutations targeted at nucleotides involved in the structural features of flaviviral Xrn1-resistant RNA motifs that were defined in this work. Our results suggest that throughout the genus flavivirus Xrn1-resistant RNA motifs adopt the same topologically conserved structure.

All genera of Flaviviridae host a conserved Xrn1-resistant RNA motif

After infection by flaviviruses like Zika and West Nile virus, eukaryotic hosts employ the well-conserved endoribonuclease Xrn1 to degrade the viral genomic RNA. Within the 3' untranslated regions, this enzyme encounters intricate Xrn1-resistant structures. This results in the accumulation of subgenomic flaviviral RNAs, an event that improves viral growth and aggravates viral pathogenicity. Xrn1-resistant RNAs have been established throughout the flaviviral genus, but not yet throughout the entire Flaviviridae family. In this work, we use previously determined characteristics of these structures to identify homologous sequences in many members of the genera pegivirus, hepacivirus and pestivirus. We used structural alignment and mutational analyses to establish that these sequences indeed represent Xrn1-resistant RNA and that they employ the general features of the flaviviral xrRNAs, consisting of a double pseudoknot formed by five base-paired regions stitched together by a crucial triple base interaction. Furthermore, we demonstrate that the pestivirus Bungowannah virus produces subgenomic RNA in vivo. Altogether, these results indicate that viruses make use of a universal Xrn1-resistant RNA throughout the Flaviviridae family.

Information Encoded by the Flavivirus Genomes beyond the Nucleotide Sequence

The genus Flavivirus comprises numerous, small, single positive-stranded RNA viruses, many of which are important human pathogens. To store all the information required for their successful propagation, flaviviruses use discrete structural genomic RNA elements to code for functional information by the establishment of dynamic networks of long-range RNA-RNA interactions that promote specific folding. These structural elements behave as true cis-acting, non-coding RNAs (ncRNAs) and have essential regulatory roles in the viral cycle. These include the control of the formation of subgenomic RNAs, known as sfRNAs, via the prevention of the complete degradation of the RNA genome. These sfRNAs are important in ensuring viral fitness. This work summarizes our current knowledge of the functions performed by the genome conformations and the role of RNA-RNA interactions in these functions. It also reviews the role of RNA structure in the production of sfRNAs across the genus Flavivirus, and their existence in related viruses.

An African tick flavivirus forming an independent clade exhibits unique exoribonuclease-resistant RNA structures in the genomic 3'-untranslated region

Tick-borne flaviviruses (TBFVs) infect mammalian hosts through tick bites and can cause various serious illnesses, such as encephalitis and hemorrhagic fevers, both in humans and animals. Despite their importance to public health, there is limited epidemiological information on TBFV infection in Africa. Herein, we report that a novel flavivirus, Mpulungu flavivirus (MPFV), was discovered in a Rhipicephalus muhsamae tick in Zambia. MPFV was found to be genetically related to Ngoye virus detected in ticks in Senegal, and these viruses formed a unique lineage in the genus Flavivirus. Analyses of dinucleotide contents of flaviviruses indicated that MPFV was similar to those of other TBFVs with a typical vertebrate genome signature, suggesting that MPFV may infect vertebrate hosts. Bioinformatic analyses of the secondary structures in the 3′-untranslated regions (UTRs) revealed that MPFV exhibited unique exoribonuclease-resistant RNA (xrRNA) structures. Utilizing biochemical approaches, we clarified that two xrRNA structures of MPFV in the 3′-UTR could prevent exoribonuclease activity. In summary, our findings provide new information regarding the geographical distribution of TBFV and xrRNA structures in the 3′-UTR of flaviviruses.

RNA Pore Translocation with Static and Periodic Forces: Effect of Secondary and Tertiary Elements on Process Activation and Duration

We use MD simulations to study the pore translocation properties of a pseudoknotted viral RNA. We consider the 71-nucleotide-long xrRNA from the Zika virus and establish how it responds when driven through a narrow pore by static or periodic forces applied to either of the two termini. Unlike the case of fluctuating homopolymers, the onset of translocation is significantly delayed with respect to the application of static driving forces. Because of the peculiar xrRNA architecture, activation times can differ by orders of magnitude at the two ends. Instead, translocation duration is much smaller than activation times and occurs on time scales comparable at the two ends. Periodic forces amplify significantly the differences at the two ends, for both activation times and translocation duration. Finally, we use a waiting-times analysis to examine the systematic slowing downs in xrRNA translocations and associate them to the hindrance of specific secondary and tertiary elements of xrRNA. The findings provide a useful reference to interpret and design future theoretical and experimental studies of RNA translocation.

Different tertiary interactions create the same important 3D features in a distinct flavivirus xrRNA

During infection by a flavivirus (FV), cells accumulate noncoding subgenomic flavivirus RNAs (sfRNAs) that interfere with several antiviral pathways. These sfRNAs are formed by structured RNA elements in the 3' untranslated region (UTR) of the viral genomic RNA, which block the progression of host cell exoribonucleases that have targeted the viral RNA. Previous work on these exoribonuclease-resistant RNAs (xrRNAs) from mosquito-borne FVs revealed a specific three-dimensional fold with a unique topology in which a ring-like structure protectively encircles the 5' end of the xrRNA. Conserved nucleotides make specific tertiary interactions that support this fold. Examination of more divergent FVs reveals differences in their 3' UTR sequences, raising the question of whether they contain xrRNAs and if so, how they fold. To answer this, we demonstrated the presence of an authentic xrRNA in the 3' UTR of the Tamana bat virus (TABV) and solved its structure by X-ray crystallography. The structure reveals conserved features from previously characterized xrRNAs, but in the TABV version these features are created through a novel set of tertiary interactions not previously seen in xrRNAs. This includes two important A-C interactions, four distinct backbone kinks, several ordered Mg²⁺ ions, and a C⁺-G-C base triple. The discovery that the same overall architecture can be achieved by very different sequences and interactions in distantly related flaviviruses provides insight into the diversity of this type of RNA and will inform searches for undiscovered xrRNAs in viruses and beyond.

The crystal structure of a Polerovirus exoribonuclease-resistant RNA shows how diverse sequences are integrated into a conserved fold

Exoribonuclease-resistant RNAs (xrRNAs) are discrete elements that block the progression of 5' to 3' exoribonucleases using specifically folded RNA structures. A recently discovered class of xrRNA is widespread in several genera of plant-infecting viruses, within both noncoding and protein-coding subgenomic RNAs. The structure of one such xrRNA from a dianthovirus revealed three-dimensional details of the resistant fold but did not answer all questions regarding the conservation and diversity of this xrRNA class. Here, we present the crystal structure of a representative polerovirus xrRNA that contains sequence elements that diverge from the previously solved structure. This new structure rationalizes previously unexplained sequence conservation patterns and shows interactions not present in the first structure. Together, the structures of these xrRNAs from dianthovirus and polerovirus genera support the idea that these plant virus xrRNAs fold through a defined pathway that includes a programmed intermediate conformation. This work deepens our knowledge of the structure-function relationship of xrRNAs and shows how evolution can craft similar RNA folds from divergent sequences.

The crystal structure of a Polerovirus exoribonuclease-resistant RNA shows how diverse sequences are integrated into a conserved fold

Exoribonuclease-resistant RNAs (xrRNAs) are discrete elements that block the progression of 5' to 3' exoribonucleases using specifically folded RNA structures. A recently discovered class of xrRNA is widespread in several genera of plant-infecting viruses, within both noncoding and protein-coding subgenomic RNAs. The structure of one such xrRNA from a dianthovirus revealed three-dimensional details of the resistant fold but did not answer all questions regarding the conservation and diversity of this xrRNA class. Here, we present the crystal structure of a representative polerovirus xrRNA that contains sequence elements that diverge from the previously solved structure. This new structure rationalizes previously unexplained sequence conservation patterns and shows interactions not present in the first structure. Together, the structures of these xrRNAs from dianthovirus and polerovirus genera support the idea that these plant virus xrRNAs fold through a defined pathway that includes a programmed intermediate conformation. This work deepens our knowledge of the structure-function relationship of xrRNAs and shows how evolution can craft similar RNA folds from divergent sequences.

Molecular mechanisms underlying the extreme mechanical anisotropy of the flaviviral exoribonuclease-resistant RNAs (xrRNAs)

Mechanical anisotropy is an essential property for many biomolecules to assume their structures, functions and applications, however, the mechanisms for their direction-dependent mechanical responses remain elusive. Herein, by using a single-molecule nanopore sensing technique, we explore the mechanisms of directional mechanical stability of the xrRNA1 RNA from ZIKA virus (ZIKV), which forms a complex ring-like architecture. We reveal extreme mechanical anisotropy in ZIKV xrRNA1 which highly depends on Mg²⁺ and the key tertiary interactions. The absence of Mg²⁺ and disruption of the key tertiary interactions strongly affect the structural integrity and attenuate mechanical anisotropy. The significance of ring structures in RNA mechanical anisotropy is further supported by steered molecular dynamics simulations in combination with force distribution analysis. We anticipate the ring structures can be used as key elements to build RNA-based nanostructures with controllable mechanical anisotropy for biomaterial and biomedical applications.

Positive selection within the genomes of SARS-CoV-2 and other Coronaviruses independent of impact on protein function

BACKGROUND

The emergence of a novel coronavirus (SARS-CoV-2) associated with severe acute respiratory disease (COVID-19) has prompted efforts to understand the genetic basis for its unique characteristics and its jump from non-primate hosts to humans. Tests for positive selection can identify apparently nonrandom patterns of mutation accumulation within genomes, highlighting regions where molecular function may have changed during the origin of a species. Several recent studies of the SARS-CoV-2 genome have identified signals of conservation and positive selection within the gene encoding Spike protein based on the ratio of synonymous to nonsynonymous substitution. Such tests cannot, however, detect changes in the function of RNA molecules.

METHODS

Here we apply a test for branch-specific oversubstitution of mutations within narrow windows of the genome without reference to the genetic code.

RESULTS

We recapitulate the finding that the gene encoding Spike protein has been a target of both purifying and positive selection. In addition, we find other likely targets of positive selection within the genome of SARS-CoV-2, specifically within the genes encoding Nsp4 and Nsp16. Homology-directed modeling indicates no change in either Nsp4 or Nsp16 protein structure relative to the most recent common ancestor. These SARS-CoV-2-specific mutations may affect molecular processes mediated by the positive or negative RNA molecules, including transcription, translation, RNA stability, and evasion of the host innate immune system. Our results highlight the importance of considering mutations in viral genomes not only from the perspective of their impact on protein structure, but also how they may impact other molecular processes critical to the viral life cycle.