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

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.

Subgenomic flaviviral RNAs (sfRNAs) are produced by pathogenic flaviviruses to counteract the host antiviral response. sfRNAs are products of incomplete degradation of viral RNA by the host exoribonuclease XRN1 which stalls on XRN1-resistant structural elements (xrRNAs) in the 3‟UTR. Here, Slonchak et al. identify xrRNAs for mosquito-specific flaviviruses, characterize their ability to stall XRN1 and produce sfRNAs, and apply structure-based phylogenetic analysis to provide evidence for evolutionary selection of these elements.

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.

A New Subclass of Exoribonuclease-Resistant RNA Found in Multiple Genera of Flaviviridae

Viruses have developed innovative strategies to exploit the cellular machinery and overcome the antiviral defenses of the host, often using specifically structured RNA elements. Examples are found in the Flavivirus genus (in the family Flaviviridae), where during flaviviral infection, pathogenic subgenomic flaviviral RNAs (sfRNAs) accumulate in the cell. These sfRNAs are formed when a host cell 5' to 3' exoribonuclease degrades the viral genomic RNA but is blocked by an exoribonuclease-resistant RNA structure (xrRNA) located in the viral genome's 3' untranslated region (UTR). Although known to exist in several Flaviviridae genera, the full distribution and diversity of xrRNAs in this family were unknown. Using the recently solved high-resolution structure of an xrRNA from the divergent flavivirus Tamana bat virus (TABV) as a reference, we used bioinformatic searches to identify xrRNAs in the remaining three genera of Flaviviridae: Pegivirus, Pestivirus, and Hepacivirus We biochemically and structurally characterized several examples, determining that they are genuine xrRNAs with a conserved fold. These new xrRNAs look superficially similar to the previously described xrRNAs but possess structural differences making them distinct from previous classes of xrRNAs. Overall, we have identified the presence of xrRNA in all four genera of Flaviviridae, but not in all species. Our findings thus require adjustments of previous xrRNA classification schemes and expand the previously known distribution of xrRNA in Flaviviridae.IMPORTANCE The members of the Flaviviridae comprise one of the largest families of positive-sense single-stranded RNA (+ssRNA) and are divided into the Flavivirus, Pestivirus, Pegivirus, and Hepacivirus genera. The genus Flavivirus contains many medically relevant viruses such as Zika virus, dengue virus, and Powassan virus. In these, a part of the RNA of the virus twists up into a distinct three-dimensional shape called an exoribonuclease-resistant RNA (xrRNA) that blocks the ability of the cell to "chew up" the viral RNA. Hence, part of the RNA of the virus remains intact, and this protected part is important for viral infection. These xrRNAs were known to occur in flaviviruses, but whether they existed in the other members of the family was not known. In this study, we identified a new subclass of xrRNA found not only in flaviviruses but also in the remaining three genera. The fact that these structured viral RNAs exist throughout the Flaviviridae family suggests they are important parts of the infection strategy of diverse pathogens, which could lead to new avenues of research.

A New Subclass of Exoribonuclease-Resistant RNA Found in Multiple Genera of Flaviviridae

The members of the Flaviviridae comprise one of the largest families of positive-sense single-stranded RNA (+ssRNA) and are divided into the Flavivirus, Pestivirus, Pegivirus, and Hepacivirus genera. The genus Flavivirus contains many medically relevant viruses such as Zika virus, dengue virus, and Powassan virus. In these, a part of the RNA of the virus twists up into a distinct three-dimensional shape called an exoribonuclease-resistant RNA (xrRNA) that blocks the ability of the cell to “chew up” the viral RNA. Hence, part of the RNA of the virus remains intact, and this protected part is important for viral infection. These xrRNAs were known to occur in flaviviruses, but whether they existed in the other members of the family was not known. In this study, we identified a new subclass of xrRNA found not only in flaviviruses but also in the remaining three genera. The fact that these structured viral RNAs exist throughout the Flaviviridae family suggests they are important parts of the infection strategy of diverse pathogens, which could lead to new avenues of research.

Viruses have developed innovative strategies to exploit the cellular machinery and overcome the antiviral defenses of the host, often using specifically structured RNA elements. Examples are found in the Flavivirus genus (in the family Flaviviridae), where during flaviviral infection, pathogenic subgenomic flaviviral RNAs (sfRNAs) accumulate in the cell. These sfRNAs are formed when a host cell 5′ to 3′ exoribonuclease degrades the viral genomic RNA but is blocked by an exoribonuclease-resistant RNA structure (xrRNA) located in the viral genome’s 3′ untranslated region (UTR). Although known to exist in several Flaviviridae genera, the full distribution and diversity of xrRNAs in this family were unknown. Using the recently solved high-resolution structure of an xrRNA from the divergent flavivirus Tamana bat virus (TABV) as a reference, we used bioinformatic searches to identify xrRNAs in the remaining three genera of Flaviviridae: Pegivirus, Pestivirus, and Hepacivirus. We biochemically and structurally characterized several examples, determining that they are genuine xrRNAs with a conserved fold. These new xrRNAs look superficially similar to the previously described xrRNAs but possess structural differences making them distinct from previous classes of xrRNAs. Overall, we have identified the presence of xrRNA in all four genera of Flaviviridae, but not in all species. Our findings thus require adjustments of previous xrRNA classification schemes and expand the previously known distribution of xrRNA in Flaviviridae.

Discoveries of Exoribonuclease-Resistant Structures of Insect-Specific Flaviviruses Isolated in Zambia

To monitor the arthropod-borne virus transmission in mosquitoes, we have attempted both to detect and isolate viruses from 3304 wild-caught female mosquitoes in the Livingstone (Southern Province) and Mongu (Western Province) regions in Zambia in 2017. A pan-flavivirus RT-PCR assay was performed to identify flavivirus genomes in total RNA extracted from mosquito lysates, followed by virus isolation and full genome sequence analysis using next-generation sequencing and rapid amplification of cDNA ends. We isolated a newly identified Barkedji virus (BJV Zambia) (10,899 nt) and a novel flavivirus, tentatively termed Barkedji-like virus (BJLV) (10,885 nt) from Culex spp. mosquitoes which shared 96% and 75% nucleotide identity with BJV which has been isolated in Israel, respectively. These viruses could replicate in C6/36 cells but not in mammalian and avian cell lines. In parallel, a comparative genomics screening was conducted to study evolutionary traits of the 5'- and 3'-untranslated regions (UTRs) of isolated viruses. Bioinformatic analyses of the secondary structures in the UTRs of both viruses revealed that the 5'-UTRs exhibit canonical stem-loop structures, while the 3'-UTRs contain structural homologs to exoribonuclease-resistant RNAs (xrRNAs), SL-III, dumbbell, and terminal stem-loop (3'SL) structures. The function of predicted xrRNA structures to stop RNA degradation by Xrn1 exoribonuclease was further proved by the in vitro Xrn1 resistance assay.

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.

Spatio-temporal characterization of the antiviral activity of the XRN1-DCP1/2 aggregation against cytoplasmic RNA viruses to prevent cell death

Host nucleases are implicated in antiviral response through the processing of pathogen-derived nucleic acids. Among many host RNases, decapping enzymes DCP1 and 2, and 5'→3' exonuclease XRN1, which are components of the RNA decay machinery, have been extensively studied in prokaryotes, plants, and invertebrates but less so in mammalian systems. As a result, the implication of XRN1 and DCPs in viral replication, in particular, the spatio-temporal dynamics during RNA viral infections remains elusive. Here, we highlight that XRN1 and DCPs play a critical role in limiting several groups of RNA viral infections. This antiviral activity was not obvious in wild-type cells but clearly observed in type I interferon (IFN-I)-deficient cells. Mechanistically, infection with RNA viruses induced the enrichment of XRN1 and DCPs in viral replication complexes (vRCs), hence forming distinct cytoplasmic aggregates. These aggregates served as sites for direct interaction between XRN1, DCP1/2, and viral ribonucleoprotein that contains viral RNA (vRNA). Although these XRN1-DCP1/2-vRC-containing foci resemble antiviral stress granules (SGs) or P-body (PB), they did not colocalize with known SG markers and did not correlate with critical PB functions. Furthermore, the presence of 5' mono- and 5' triphosphate structures on vRNA was not required for the formation of XRN1-DCP1/2-vRC-containing foci. On the other hand, single-, double-stranded, and higher-ordered vRNA species play a role but are not deterministic for efficient formation of XRN1-DCP1/2 foci and consequent antiviral activity in a manner proportional to RNA length. These results highlight the mechanism behind the antiviral function of XRN1-DCP1/2 in RNA viral infections independent of IFN-I response, protein kinase R and PB function.

Directional translocation resistance of Zika xrRNA

xrRNAs from flaviviruses survive in host cells because of their exceptional dichotomic response to the unfolding action of different enzymes. They can be unwound, and hence copied, by replicases, and yet can resist degradation by exonucleases. How the same stretch of xrRNA can encode such diverse responses is an open question. Here, by using atomistic models and translocation simulations, we uncover an elaborate and directional mechanism for how stress propagates when the two xrRNA ends, [Formula: see text] and [Formula: see text], are driven through a pore. Pulling the [Formula: see text] end, as done by replicases, elicits a progressive unfolding; pulling the [Formula: see text] end, as done by exonucleases, triggers a counterintuitive molecular tightening. Thus, in what appears to be a remarkable instance of intra-molecular tensegrity, the very pulling of the [Formula: see text] end is what boosts resistance to translocation and consequently to degradation. The uncovered mechanistic principle might be co-opted to design molecular meta-materials.

Impact of alphavirus 3'UTR plasticity on mosquito transmission

Alphaviruses such as chikungunya and western equine encephalitis viruses are important human pathogens transmitted by mosquitoes that have recently caused large epidemic and epizootic outbreaks. The epidemic potential of alphaviruses is often related to enhanced mosquito transmission. Tissue barriers and antiviral responses impose bottlenecks to viral populations in mosquitoes. Substitutions in the envelope proteins and the presence of repeated sequence elements (RSEs) in the 3'UTR of epidemic viruses were proposed to be specifically associated to efficient replication in mosquito vectors. Here, we discuss the molecular mechanisms that originated RSEs, the evolutionary forces that shape the 3'UTR of alphaviruses, and the significance of RSEs for mosquito transmission. Finally, the presence of RSEs in the 3'UTR of viral genomes appears as evolutionary trait associated to mosquito adaptation and emerges as a common feature among viruses from the alphavirus and flavivirus genera.

Comprehensive in-vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms

SARS-CoV-2 is the positive-sense RNA virus that causes COVID-19, a disease that has triggered a major human health and economic crisis. The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form stable RNA structures and yet, as much as 97% of its 30 kilobases have not been structurally explored in the context of a viral infection. Our limited knowledge of SARS-CoV-2 genomic architecture is a fundamental limitation to both our mechanistic understanding of coronavirus life cycle and the development of COVID-19 RNA-based therapeutics. Here, we apply a novel long amplicon strategy to determine for the first time the secondary structure of the SARS-CoV-2 RNA genome probed in infected cells. In addition to the conserved structural motifs at the viral termini, we report new structural features like a conformationally flexible programmed ribosomal frameshifting pseudoknot, and a host of novel RNA structures, each of which highlights the importance of studying viral structures in their native genomic context. Our in-depth structural analysis reveals extensive networks of well-folded RNA structures throughout Orf1ab and reveals new aspects of SARS-CoV-2 genome architecture that distinguish it from other single-stranded, positive-sense RNA viruses. Evolutionary analysis of RNA structures in SARS-CoV-2 shows that several features of its genomic structure are conserved across beta coronaviruses and we pinpoint individual regions of well-folded RNA structure that merit downstream functional analysis. The native, complete secondary structure of SAR-CoV-2 presented here is a roadmap that will facilitate focused studies on mechanisms of replication, translation and packaging, and guide the identification of new RNA drug targets against COVID-19.

VADR: validation and annotation of virus sequence submissions to GenBank

Background

GenBank contains over 3 million viral sequences. The National Center for Biotechnology Information (NCBI) previously made available a tool for validating and annotating influenza virus sequences that is used to check submissions to GenBank. Before this project, there was no analogous tool in use for non-influenza viral sequence submissions.

Results

We developed a system called VADR (Viral Annotation DefineR) that validates and annotates viral sequences in GenBank submissions. The annotation system is based on the analysis of the input nucleotide sequence using models built from curated RefSeqs. Hidden Markov models are used to classify sequences by determining the RefSeq they are most similar to, and feature annotation from the RefSeq is mapped based on a nucleotide alignment of the full sequence to a covariance model. Predicted proteins encoded by the sequence are validated with nucleotide-to-protein alignments using BLAST. The system identifies 43 types of “alerts” that (unlike the previous BLAST-based system) provide deterministic and rigorous feedback to researchers who submit sequences with unexpected characteristics. VADR has been integrated into GenBank’s submission processing pipeline allowing for viral submissions passing all tests to be accepted and annotated automatically, without the need for any human (GenBank indexer) intervention. Unlike the previous submission-checking system, VADR is freely available (https://github.com/nawrockie/vadr) for local installation and use. VADR has been used for Norovirus submissions since May 2018 and for Dengue virus submissions since January 2019. Since March 2020, VADR has also been used to check SARS-CoV-2 sequence submissions. Other viruses with high numbers of submissions will be added incrementally.

Conclusion

VADR improves the speed with which non-flu virus submissions to GenBank can be checked and improves the content and quality of the GenBank annotations. The availability and portability of the software allow researchers to run the GenBank checks prior to submitting their viral sequences, and thereby gain confidence that their submissions will be accepted immediately without the need to correspond with GenBank staff. Reciprocally, the adoption of VADR frees GenBank staff to spend more time on services other than checking routine viral sequence submissions.

Positive epistasis between viral polymerase and the 3' untranslated region of its genome reveals the epidemiologic fitness of dengue virus

Dengue virus (DENV) is a global health threat, causing repeated epidemics throughout the tropical world. While low herd immunity levels to any one of the four antigenic types of DENV predispose populations to outbreaks, viral genetic determinants that confer greater fitness for epidemic spread is an important but poorly understood contributor of dengue outbreaks. Here we report that positive epistasis between the coding and noncoding regions of the viral genome combined to elicit an epidemiologic fitness phenotype associated with the 1994 DENV2 outbreak in Puerto Rico. We found that five amino acid substitutions in the NS5 protein reduced viral genomic RNA (gRNA) replication rate to achieve a more favorable and relatively more abundant subgenomic flavivirus RNA (sfRNA), a byproduct of host 5'-3' exoribonuclease activity. The resulting increase in sfRNA relative to gRNA levels not only inhibited type I interferon (IFN) expression in infected cells through a previously described mechanism, but also enabled sfRNA to compete with gRNA for packaging into infectious particles. We suggest that delivery of sfRNA to new susceptible cells to inhibit type I IFN induction before gRNA replication and without the need for further de novo sfRNA synthesis could form a "preemptive strike" strategy against DENV.

Zika virus noncoding RNA suppresses apoptosis and is required for virus transmission by mosquitoes

Flaviviruses, including Zika virus (ZIKV), utilise host mRNA degradation machinery to produce subgenomic flaviviral RNA (sfRNA). In mammalian hosts, this noncoding RNA facilitates replication and pathogenesis of flaviviruses by inhibiting IFN-signalling, whereas the function of sfRNA in mosquitoes remains largely elusive. Herein, we conduct a series of in vitro and in vivo experiments to define the role of ZIKV sfRNA in infected Aedes aegypti employing viruses deficient in production of sfRNA. We show that sfRNA-deficient viruses have reduced ability to disseminate and reach saliva, thus implicating the role for sfRNA in productive infection and transmission. We also demonstrate that production of sfRNA alters the expression of mosquito genes related to cell death pathways, and prevents apoptosis in mosquito tissues. Inhibition of apoptosis restored replication and transmission of sfRNA-deficient mutants. Hence, we propose anti-apoptotic activity of sfRNA as the mechanism defining its role in ZIKV transmission.

The function on subgenomic flaviviral RNA (sfRNA) in the mosquito vector is not well understood. Here, Slonchak et al. show that sfRNA affects virus-induced apoptosis and dissemination of ZIKV in Aedes aegypti mosquitoes, suggesting a role of sfRNA in Zika virus replication and transmission.

Role of RNA-binding proteins during the late stages of Flavivirus replication cycle

The genus Flavivirus encompasses several worldwide-distributed arthropod-borne viruses including, dengue virus, Japanese encephalitis virus, West Nile virus, yellow fever virus, Zika virus, and tick-borne encephalitis virus. Infection with these viruses manifest with symptoms ranging from febrile illness to life- threatening hypotensive shock and encephalitis. Therefore, flaviviruses pose a great risk to public health. Currently, preventive measures are falling short to control epidemics and there are no antivirals against any Flavivirus.Flaviviruses carry a single stranded positive-sense RNA genome that plays multiple roles in infected cells: it is translated into viral proteins, used as template for genome replication, it is the precursor of the subgenomic flaviviral RNA and it is assembled into new virions. Furthermore, viral RNA genomes are also packaged into extracellular vesicles, e.g. exosomes, which represent an alternate mode of virus dissemination.Because RNA molecules are at the center of Flavivirus replication cycle, viral and host RNA-binding proteins (RBPs) are critical determinants of infection. Numerous studies have revealed the function of RBPs during Flavivirus infection, particularly at the level of RNA translation and replication. These proteins, however, are also critical participants at the late stages of the replication cycle. Here we revise the function of host RBPs and the viral proteins capsid, NS2A and NS3, during the packaging of viral RNA and the assembly of new virus particles. Furthermore, we go through the evidence pointing towards the importance of host RBPs in mediating cellular RNA export with the idea that the biogenesis of exosomes harboring Flavivirus RNA would follow an analogous pathway.

The Development of RNA-KISS, a Mammalian Three-Hybrid Method to Detect RNA-Protein Interactions in Living Mammalian Cells

RNA-protein interactions are essential for the regulation of mRNA and noncoding RNA functions and are implicated in many diseases, such as cancer and neurodegenerative disorders. A method that can detect RNA-protein interactions in living mammalian cells on a proteome-wide scale will be an important asset to identify and study these interactions. Here we show that a combination of the mammalian two-hybrid protein-protein detection method KISS (kinase substrate sensor) and the yeast RNA three-hybrid method, utilizing the specific interaction between the MS2 RNA and MS2 coat protein, is capable of detecting RNA-protein interactions in living mammalian cells. For conceptional proof we used the subgenomic flavivirus RNA (sfRNA) of the dengue virus (DENV), a highly structured noncoding RNA derived from the DENV genome known to target host cell proteins involved in innate immunity and antiviral defense, as bait. Using RNA-KISS, we could confirm the previously established interaction between the RNA-binding domain of DDX6 and the DENV sfRNA. Finally, we performed a human proteome-wide screen for DENV sfRNA-binding host factors, identifying several known flavivirus host factors such as DDX6 and PACT, further validating the RNA-KISS method as a robust and high-throughput cell-based RNA-protein interaction screening tool.

Universal RNA Secondary Structure Insight Into Mosquito-Borne Flavivirus (MBFV) cis-Acting RNA Biology

Mosquito-borne flaviviruses (MBFVs) spread between vertebrate (mammals and birds) and invertebrate (mosquitoes) hosts. The cis-acting RNAs of MBFV share common evolutionary origins and contain frequent alterations, which control the balance of linear and circular genome conformations and allow effective replication. Importantly, multiple cis-acting RNAs interact with trans-acting regulatory RNA-binding proteins (RBPs) and affect the MBFV lifecycle process, including viral replicase binding, viral RNA translation-cyclisation-synthesis and nucleocapsid assembly. Considering that extensive structural probing analyses have been performed on MBFV cis-acting RNAs, herein the homologous RNA structures are online folded and consensus structures are constructed by sort. The specific traits and underlying biology of MBFV cis-acting RNA are illuminated accordingly in a review of RNA structure. These findings deepen our understanding of MBFV cis-acting RNA biology and serve as a resource for designing therapeutics in targeting protein-viral RNA interaction or viral RNA secondary structures.