Integrative Analysis of Zika Virus Genome RNA Structure Reveals Critical Determinants of Viral Infectivity

Promotion of TLR7-MyD88-dependent inflammation and autoimmunity in mice through stem-loop changes in Lnc-Atg16l1

Uncontrolled TLR signaling can cause inflammatory immunopathology and trigger autoimmune diseases. For example, TLR7 promotes pathogenesis of systemic lupus erythematosus. However, whether RNA structural changes affect nucleic acids-sensing TLRs signaling and impact disease progression is unclear. Here by iCLIP-seq we identify a TLR7-binding long non-coding RNA, Lnc-Atg16l1, and find that it promotes TLR7 and other MyD88-dependent TLRs signaling in various types of immune cells. Depletion of Lnc-Atg16l1 attenuates development of TLR7-linked autoimmune phenotypes in the mouse SLE model. Mechanistically, we find that Lnc-Atg16l1 binds to TLR7 at bases near U84 and MyD88 at bases around A129. The analysis of Lnc-Atg16l1 in situ structures show that it strengthens the interaction between TIR domain of TLR7 and MyD88 through specific stem-loop structure changes as a molecular scaffold after TLR7 activation to promote TLR7 downstream signaling. Therefore, we discover a mechanism for host RNA regulation of innate signaling and autoimmune disease through its structural changes. These findings provide insights into the pro-inflammatory function of self RNA in a structure-dependent manner and suggest a potential target for TLR-related autoimmune disorders.

An accumulated mutation gained in mosquito cells enhances Zika virus virulence and fitness in mice

Zika virus (ZIKV) remains a significant public health threat worldwide. A number of adaptive mutations have accumulated within the genome of ZIKV during global transmission, some of which have been linked to specific phenotypes. ZIKV maintains an alternating cycle of replication between mosquitoes and vertebrate hosts, but the role of mosquito-specific adaptive mutations in ZIKV has not been well investigated. In this study, we demonstrated that serial passaging of ZIKV in mosquito Aag2 cells led to the emergence of critical amino acid substitutions, including A94V in the prM protein and V153D and H401Y in the E protein. Further characterization via reverse genetics revealed that the H401Y substitution in the E protein did not augment viral replication in mosquitoes but significantly enhanced neurovirulence and lethality compared with those of the wild-type (WT) virus in mice. More importantly, the H401Y mutant maintained its virulence phenotype in mice after propagation in mosquitoes in mosquito-mouse cycle model. In particular, recombinant ZIKV harboring the H401Y substitution showed enhanced competitive fitness over WT ZIKV in various mammalian cells and mouse brains, but not in mosquito cells. Notably, the H401Y substitution in the ZIKV E protein has been detected in recent isolates derived from both mosquitoes and humans in Asia and the Americas. In summary, our findings not only identify a novel virulence determinant of ZIKV but also highlight the complexity of the relationship between the evolution of vector-borne viruses and their clinical outcome in nature.

IMPORTANCE

Zika virus (ZIKV) is an important arbovirus with a global impact. Experimental evolution by serial passaging of ZIKV in susceptible cells has led to the identification of a panel of critical amino acid substitutions with specific functions. Herein, we identified a mosquito cell-derived substitution, H401Y, in the ZIKV E protein via experimental evolution. The H401Y substitution significantly enhanced viral virulence and fitness in mammal cells and mice. Notably, the H401Y substitution has been detected in recent mosquito and human isolates from regions spanning Asia to the Americas. Our work elucidates unrecognized virulence determinant in the ZIKV genome that warrants urgent attention. Moreover, the findings underscore the critical need for extensive molecular surveillance and rigorous clinical observation to establish the potential impact in natural circulation. These endeavors are crucial for unraveling the potential of mutation to act as a catalyst for future epidemics, thereby preempting the public health challenges it may pose.

Viral RNA Interactome: The Ultimate Researcher's Guide to RNA-Protein Interactions

RNA molecules in the cell are bound by a multitude of RNA-binding proteins (RBPs) with a variety of regulatory consequences. Often, interactions with these RNA-binding proteins are facilitated by the complex secondary and tertiary structures of RNA molecules. Viral RNAs especially are known to be heavily structured and interact with many RBPs, with roles including genome packaging, immune evasion, enhancing replication and transcription, and increasing translation efficiency. As such, the RNA-protein interactome represents a critical facet of the viral replication cycle. Characterization of these interactions is necessary for the development of novel therapeutics targeted at the disruption of essential replication cycle events. In this review, we aim to summarize the various roles of RNA structures in shaping the RNA-protein interactome, the regulatory roles of these interactions, as well as up-to-date methods developed for the characterization of the interactome and directions for novel, RNA-directed therapeutics.

RNA architecture of porcine deltacoronavirus genome inside virions detected by vRIC-seq

Porcine deltacoronavirus (PDCoV) is a newly emerging and special delta coronavirus, which infect mammals such as pigs, cattle and humans, as well as chickens and birds. Exploring RNA structures in the viral genome benefits the understanding of the role of RNA in the lifecycle of viruses. In this study, vRIC-seq is employed to analyze the RNA-RNA interaction in the whole genome structure of PDCoV in virions. About 12.87 and 13.52 million paired reads are obtained in two biological replicates, respectively, with 17.9% and 14.8% of them are identified as valid chimeric reads. These are employed to predict the RNA secondary structure, which is compact and highly structured. A twisted-cyclized conformation is observed in the RNA-RNA interaction map of PDCoV for the first time. 77 multi-way junctions are evenly distributed in the PDCoV genome. Our work provides fundamental structural insights that are essential for understanding the genomic structure and function, genetic evolution, and packaging characteristics of PDCoV.

Genome-Wide Identification of Stable RNA Secondary Structures Across Multiple Organisms Using Chemical Probing Data: Insights into Short Structural Motifs and RNA-Targeting Therapeutics

Small molecules targeting specific RNA binding sites, including stable and transient RNA structures, are emerging as effective pharmacological approaches for modulating gene expression. However, little is understood about how stable RNA secondary structures are shared across organisms, an important factor in controlling drug selectivity. In this study, I provide an analytical pipeline named RNA Secondary Structure Finder (R2S-Finder) to discover short, stable RNA structural motifs for humans, Escherichia coli ( E. coli ), SARS-CoV-2, and Zika virus by leveraging existing in vivo and in vitro genome-wide chemical RNA-probing datasets. I found several common features across organisms. For example, apart from the well-documented tetraloops, AU-rich tetraloops are widely present in different organisms. I also found that the 5' untranslated region (UTR) contains a higher proportion of stable structures than the coding sequences in humans, SARS-CoV-2, and Zika virus. In general, stable structures predicted from in vitro (protein-free) and in vivo datasets are consistent in humans, E. coli , and SARS-CoV-2, indicating that most stable structure formation were driven by RNA folding alone, while a larger variation was found between in vitro and in vivo data with certain RNA types, such as human long intergenic non-coding RNAs (lincRNAs). Finally, I predicted stable three- and four-way RNA junctions that exist both in vivo and in vitro conditions, which can potentially serve as drug targets. All results of stable sequences, stem-loops, internal loops, bulges, and three- and four-way junctions have been collated in the R2S-Finder database ( https://github.com/JingxinWangLab/R2S-Finder ), which is coded in hyperlinked HTML pages and tabulated in CSV files.

Mapping the structural landscape of the yeast Ty3 retrotransposon RNA genome

Long terminal repeat (LTR)-retrotransposons are significant contributors to the evolution and diversity of eukaryotic genomes. Their RNA genomes (gRNA) serve as a template for protein synthesis and reverse transcription to a DNA copy, which can integrate into the host genome. Here, we used the SHAPE-MaP strategy to explore Ty3 retrotransposon gRNA structure in yeast and under cell-free conditions. Our study reveals the structural dynamics of Ty3 gRNA and the well-folded core, formed independently of the cellular environment. Based on the detailed map of Ty3 gRNA structure, we characterized the structural context of cis-acting sequences involved in reverse transcription and frameshifting. We also identified a novel functional sequence as a potential initiator for Ty3 gRNA dimerization. Our data indicate that the dimer is maintained by direct interaction between short palindromic sequences at the 5' ends of the two Ty3 gRNAs, resembling the model characteristic for other retroelements like HIV-1 and Ty1. This work points out a range of cell-dependent and -independent Ty3 gRNA structural changes that provide a solid background for studies on RNA structure-function relationships important for retroelement biology.

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.

SHAPEwarp-web: sequence-agnostic search for structurally homologous RNA regions across databases of chemical probing data

RNA molecules perform a variety of functions in cells, many of which rely on their secondary and tertiary structures. Chemical probing methods coupled with high-throughput sequencing have significantly accelerated the mapping of RNA structures, and increasingly large datasets of transcriptome-wide RNA chemical probing data are becoming available. Analogously to what has been done for decades in the protein world, this RNA structural information can be leveraged to aid the discovery of structural similarity to a known RNA (or RNA family), which, in turn, can inform about the function of transcripts. We have previously developed SHAPEwarp, a sequence-agnostic method for the search of structurally homologous RNA segments in a database of reactivity profiles derived from chemical probing experiments. In its original implementation, however, SHAPEwarp required substantial computational resources, even for moderately sized databases, as well as significant Linux command line know-how. To address these limitations, we introduce here SHAPEwarp-web, a user-friendly web interface to rapidly query large databases of RNA chemical probing data for structurally similar RNAs. Aside from featuring a completely rewritten core, which speeds up by orders of magnitude the search inside large databases, the web server hosts several high-quality chemical probing databases across multiple species. SHAPEwarp-web is available from https://shapewarp.incarnatolab.com.

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

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

The 3' terminal region of Zika virus RNA contains a conserved G-quadruplex and is unfolded by human DDX17

Zika virus (ZIKV) infection remains a worldwide concern, and currently no effective treatments or vaccines are available. Novel therapeutics are an avenue of interest that could probe viral RNA-human protein communication to stop viral replication. One specific RNA structure, G-quadruplexes (G4s), possess various roles in viruses and all domains of life, including transcription and translation regulation and genome stability, and serves as nucleation points for RNA liquid-liquid phase separation. Previous G4 studies on ZIKV using a quadruplex forming G-rich sequences Mapper located a potential G-quadruplex sequence in the 3' terminal region (TR) and was validated structurally using a 25-mer oligo. It is currently unknown if this structure is conserved and maintained in a large ZIKV RNA transcript and its specific roles in viral replication. Using bioinformatic analysis and biochemical assays, we demonstrate that the ZIKV 3' TR G4 is conserved across all ZIKV isolates and maintains its structure in a 3' TR full-length transcript. We further established the G4 formation using pyridostatin and the BG4 G4-recognizing antibody binding assays. Our study also demonstrates that the human DEAD-box helicases, DDX3X132-607 and DDX17135-555, bind to the 3' TR and that DDX17135-555 unfolds the G4 present in the 3' TR. These findings provide a path forward in potential therapeutic targeting of DDX3X or DDX17's binding to the 3' TR G4 region for novel treatments against ZIKV.

Experimental and computational methods for studying the dynamics of RNA-RNA interactions in SARS-COV2 genomes

Long-range ribonucleic acid (RNA)-RNA interactions (RRI) are prevalent in positive-strand RNA viruses, including Beta-coronaviruses, and these take part in regulatory roles, including the regulation of sub-genomic RNA production rates. Crosslinking of interacting RNAs and short read-based deep sequencing of resulting RNA-RNA hybrids have shown that these long-range structures exist in severe acute respiratory syndrome coronavirus (SARS-CoV)-2 on both genomic and sub-genomic levels and in dynamic topologies. Furthermore, co-evolution of coronaviruses with their hosts is navigated by genetic variations made possible by its large genome, high recombination frequency and a high mutation rate. SARS-CoV-2's mutations are known to occur spontaneously during replication, and thousands of aggregate mutations have been reported since the emergence of the virus. Although many long-range RRIs have been experimentally identified using high-throughput methods for the wild-type SARS-CoV-2 strain, evolutionary trajectory of these RRIs across variants, impact of mutations on RRIs and interaction of SARS-CoV-2 RNAs with the host have been largely open questions in the field. In this review, we summarize recent computational tools and experimental methods that have been enabling the mapping of RRIs in viral genomes, with a specific focus on SARS-CoV-2. We also present available informatics resources to navigate the RRI maps and shed light on the impact of mutations on the RRI space in viral genomes. Investigating the evolution of long-range RNA interactions and that of virus-host interactions can contribute to the understanding of new and emerging variants as well as aid in developing improved RNA therapeutics critical for combating future outbreaks.

The myriad roles of RNA structure in the flavivirus life cycle

As positive-sense RNA viruses, the genomes of flaviviruses serve as the template for all stages of the viral life cycle, including translation, replication, and infectious particle production. Yet, they encode just 10 proteins, suggesting that the structure and dynamics of the viral RNA itself helps shepherd the viral genome through these stages. Herein, we highlight advances in our understanding of flavivirus RNA structural elements through the lens of their impact on the viral life cycle. We highlight how RNA structures impact translation, the switch from translation to replication, negative- and positive-strand RNA synthesis, and virion assembly. Consequently, we describe three major themes regarding the roles of RNA structure in flavivirus infections: 1) providing a layer of specificity; 2) increasing the functional capacity; and 3) providing a mechanism to support genome compaction. While the interactions described herein are specific to flaviviruses, these themes appear to extend more broadly across RNA viruses.

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

Small Molecules Targeting Viral RNA

The majority of antivirals available target viral proteins; however, RNA is emerging as a new and promising antiviral target due to the presence of highly structured RNA in viral genomes fundamental for their replication cycle. Here, we discuss methods for the identification of RNA-targeting compounds, starting from the determination of RNA structures either from purified RNA or in living cells, followed by in silico screening on RNA and phenotypic assays to evaluate viral inhibition. Moreover, we review the small molecules known to target the programmed ribosomal frameshifting element of SARS-CoV-2, the internal ribosomal entry site of different viruses, and RNA elements of HIV.

Zika virus dumbbell-1 structure is critical for sfRNA presence and cytopathic effect during infection

All flaviviruses contain conserved RNA structures in the 3' untranslated region (3' UTR) that are important for flavivirus RNA replication, translation, and pathogenesis. Flaviviruses like Zika virus (ZIKV) contain multiple conserved RNA structures in the viral 3' UTR, including the structure known as dumbbell-1 (DB-1). Previous research has shown that the DB-1 structure is important for flavivirus positive-strand genome replication, but the functional role of the flavivirus DB-1 structure and the mechanism by which it contributes to viral pathogenesis are not known. Using the recently solved flavivirus DB RNA structural data, we designed two DB-1 mutant ZIKV infectious clones, termed ZIKV-TL.PK and ZIKV-p.2.5', which disrupt DB-1 tertiary folding. We found that viral positive-strand genome replication of both ZIKV DB-1 mutant clones is similar to wild-type (WT) ZIKV, but ZIKV DB-1 mutants exhibit significantly decreased cytopathic effect due to reduced caspase-3 activation. We next show that ZIKV DB-1 mutants exhibit decreased levels of sfRNA species compared to ZIKV-WT during infection. However, ZIKV DB-1 mutant 3' UTRs exhibit unchanged sfRNA biogenesis following XRN1 degradation in vitro. We also found that ZIKV DB-1 mutant virus (ZIKV-p.2.5') exhibited enhanced sensitivity to type I interferon treatment, and both ZIKV-DB-1 mutants exhibit reduced morbidity and mortality due to tissue-specific attenuated viral replication in brain tissue of interferon type I/II receptor knockout mice. We propose that the flavivirus DB-1 RNA structure maintains sfRNA levels during infection despite maintained sfRNA biogenesis, and these results indicate that ZIKV DB-dependent maintenance of sfRNA levels support caspase-3-dependent, cytopathic effect, type I interferon resistance, and viral pathogenesis in mammalian cells and in a ZIKV murine model of disease. IMPORTANCE The group of viruses termed flaviviruses cause important disease throughout the world and include dengue virus, Zika virus, Japanese encephalitis virus, and many more. All of these flaviviruses have highly conserved RNA structures in the untranslated regions of the virus genome. One of the shared RNA structures, termed the dumbbell region, is not well studied, but mutations in this region are important for vaccine development. In this study, we made structure-informed targeted mutations in the Zika virus dumbbell region and studied the effect on the virus. We found that Zika virus dumbbell mutants are significantly weakened or attenuated due to a decreased ability to produce non-coding RNA that is needed to support infection, support virus-induced cell death, and support escape from the host immune system. These data show that targeted mutations in the flavivirus dumbbell RNA structure may be an important approach to develop future vaccine candidates.

5'/3' RACE method for sequencing the 5' and 3' untranslated regions of Zika virus

The spread of Zika virus (ZIKV) from the African continent to the Americas promoted its molecular evolution, as reflected by mutations in its RNA genome. Most of the ZIKV genome sequences in the GenBank database have incomplete 5' and 3' UTR sequences, reflecting the deficiency of whole-genome sequencing technologies to resolve the sequences of the genome ends. We modified a protocol for rapid amplification of cDNA ends (RACE) to determine the complete sequences of the 5' and 3' UTRs of a previously reported ZIKV isolate (GenBank no. MH544701.1). This strategy is useful for determining 5' and 3' UTR sequences of ZIKV isolates and will be useful for comparative genomics applications.

ASFV transcription reporter screening system identifies ailanthone as a broad antiviral compound

African swine fever (ASF) is an acute, highly contagious and deadly viral disease in swine that jeopardizes the worldwide pig industry. Unfortunately, there are no authoritative vaccine and antiviral drug available for ASF control. African swine fever virus (ASFV) is the etiological agent of ASF. Among the ASFV proteins, p72 is the most abundant component in the virions and thus a potential target for anti-ASFV drug design. Here, we constructed a luciferase reporter system driven by the promoter of p72, which is transcribed by the co-transfected ASFV RNA polymerase complex. Using this system, we screened over 3200 natural product compounds and obtained three potent candidates against ASFV. We further evaluated the anti-ASFV effects and proved that among the three candidates, ailanthone (AIL) inhibits the replication of ASFV at the nanomolar concentration (IC50 ​= ​15 ​nmol/L). Our in vitro experiments indicated that the antiviral effect of AIL is associated with its inhibition of the HSP90-p23 cochaperone. Finally, we showed the antiviral activity of AIL on Zika virus and hepatitis B virus (HBV), which supports that AIL is a potential broad-spectrum antiviral agent.

PrismNet: predicting protein-RNA interaction using in vivo RNA structural information

Fundamental to post-transcriptional regulation, the in vivo binding of RNA binding proteins (RBPs) on their RNA targets heavily depends on RNA structures. To date, most methods for RBP-RNA interaction prediction are based on RNA structures predicted from sequences, which do not consider the various intracellular environments and thus cannot predict cell type-specific RBP-RNA interactions. Here, we present a web server PrismNet that uses a deep learning tool to integrate in vivo RNA secondary structures measured by icSHAPE experiments with RBP binding site information from UV cross-linking and immunoprecipitation in the same cell lines to predict cell type-specific RBP-RNA interactions. Taking an RBP and an RNA region with sequential and structural information as input ('Sequence & Structure' mode), PrismNet outputs the binding probability of the RBP and this RNA region, together with a saliency map and a sequence-structure integrative motif. The web server is freely available at http://prismnetweb.zhanglab.net.

Dengue Virus Capsid Protein Facilitates Genome Compaction and Packaging

Dengue virus (DENV) is a single-stranded (+)-sense RNA virus that infects humans and mosquitoes, posing a significant health risk in tropical and subtropical regions. Mature virions are composed of an icosahedral shell of envelope (E) and membrane (M) proteins circumscribing a lipid bilayer, which in turn contains a complex of the approximately 11 kb genomic RNA with capsid (C) proteins. Whereas the structure of the envelope is clearly defined, the structure of the packaged genome in complex with C proteins remains elusive. Here, we investigated the interactions of C proteins with viral RNA, in solution and inside mature virions, via footprinting and cross-linking experiments. We demonstrated that C protein interaction with DENV genomes saturates at an RNA:C protein ratio below 1:250. Moreover, we also showed that the length of the RNA genome interaction sites varies, in a multimodal distribution, consistent with the C protein binding to each RNA site mostly in singlets or pairs (and, in some instances, higher numbers). We showed that interaction sites are preferentially sites with low base pairing, as previously measured by 2'-acetylation analyzed by primer extension (SHAPE) reactivity indicating structuredness. We found a clear association pattern emerged: RNA-C protein binding sites are strongly associated with long-range RNA-RNA interaction sites, particularly inside virions. This, in turn, explains the need for C protein in viral genome packaging: the protein has a chief role in coordinating these key interactions, promoting proper packaging of viral RNA. Such sites are, thus, highly consequential for viral assembly, and, as such, may be targeted in future drug development strategies against these and related viruses.

Specialized cis-Acting RNA Elements Balance Genome Cyclization to Ensure Efficient Replication of Yellow Fever Virus

Genome cyclization is essential for viral RNA (vRNA) replication of the vertebrate-infecting flaviviruses, and yet its regulatory mechanisms are not fully understood. Yellow fever virus (YFV) is a notorious pathogenic flavivirus. Here, we demonstrated that a group of cis-acting RNA elements in YFV balance genome cyclization to govern efficient vRNA replication. It was shown that the downstream of the 5'-cyclization sequence hairpin (DCS-HP) is conserved in the YFV clade and is important for efficient YFV propagation. By using two different replicon systems, we found that the function of the DCS-HP is determined primarily by its secondary structure and, to a lesser extent, by its base-pair composition. By combining in vitro RNA binding and chemical probing assays, we found that the DCS-HP orchestrates the balance of genome cyclization through two different mechanisms, as follows: the DCS-HP assists the correct folding of the 5' end in a linear vRNA to promote genome cyclization, and it also limits the overstabilization of the circular form through a potential crowding effect, which is influenced by the size and shape of the DCS-HP structure. We also provided evidence that an A-rich sequence downstream of the DCS-HP enhances vRNA replication and contributes to the regulation of genome cyclization. Interestingly, diversified regulatory mechanisms of genome cyclization, involving both the downstream of the 5'-cyclization sequence (CS) and the upstream of the 3'-CS elements, were identified among different subgroups of the mosquito-borne flaviviruses. In summary, our work highlighted how YFV precisely controls the balance of genome cyclization to ensure viral replication. IMPORTANCE Yellow fever virus (YFV), the prototype of the Flavivirus genus, can cause devastating yellow fever disease. Although it is preventable by vaccination, there are still tens of thousands of yellow fever cases per year, and no approved antiviral medicine is available. However, the understandings about the regulatory mechanisms of YFV replication are obscure. In this study, by a combination of bioinformatics, reverse genetics, and biochemical approaches, it was shown that the downstream of the 5'-cyclization sequence hairpin (DCS-HP) promotes efficient YFV replication by modulating the conformational balance of viral RNA. Interestingly, we found specialized combinations for the downstream of the 5'-cyclization sequence (CS) and upstream of the 3'-CS elements in different groups of the mosquito-borne flaviviruses. Moreover, possible evolutionary relationships among the various downstream of the 5'-CS elements were implied. This work highlighted the complexity of RNA-based regulatory mechanisms in the flaviviruses and will facilitate the design of RNA structure-targeted antiviral therapies.

RNA levers and switches controlling viral gene expression

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

Probing the dynamic RNA structurome and its functions

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

Chemical RNA Cross-Linking: Mechanisms, Computational Analysis, and Biological Applications

In recent years, RNA has emerged as a multifaceted biomolecule that is involved in virtually every function of the cell and is critical for human health. This has led to a substantial increase in research efforts to uncover the many chemical and biological aspects of RNA and target RNA for therapeutic purposes. In particular, analysis of RNA structures and interactions in cells has been critical for understanding their diverse functions and druggability. In the last 5 years, several chemical methods have been developed to achieve this goal, using chemical cross-linking combined with high-throughput sequencing and computational analysis. Applications of these methods resulted in important new insights into RNA functions in a variety of biological contexts. Given the rapid development of new chemical technologies, a thorough perspective on the past and future of this field is provided. In particular, the various RNA cross-linkers and their mechanisms, the computational analysis and challenges, and illustrative examples from recent literature are discussed.

Loop-Mediated Isothermal Amplification on Paper Microfluidic Chips for Highly Sensitive and Specific Zika Virus Detection Using Smartphone

Zika virus (ZIKV) infection may cause serious birth defects and is a critical concern for women of child-bearing age in affected regions. A simple, portable, and easy-to-use ZIKV detection method would enable point-of-care testing, which may aid in prevention of the spread of the virus. Herein, we describe a reverse transcription isothermal loop-mediated amplification (RT-LAMP) method that detects the presence of ZIKV RNA in complex samples (e.g., blood, urine, and tap water). Phenol red is the colorimetric indicator of successful amplification. Color changes based on the amplified RT-LAMP product from the presence of viral target are monitored using a smartphone camera under ambient light conditions. A single viral RNA molecule per μL can be detected in as quickly as 15 min using this method with 100% sensitivity and 100% specificity in blood and tap water, while 100% sensitivity and 67% specificity in urine. This platform can also be used to identify other viruses including SARS-CoV-2 and improve the current state of field-based diagnostics.

ScanFold 2.0: a rapid approach for identifying potential structured RNA targets in genomes and transcriptomes

A major limiting factor in target discovery for both basic research and therapeutic intervention is the identification of structural and/or functional RNA elements in genomes and transcriptomes. This was the impetus for the original ScanFold algorithm, which provides maps of local RNA structural stability, evidence of sequence-ordered (potentially evolved) structure, and unique model structures comprised of recurring base pairs with the greatest structural bias. A key step in quantifying this propensity for ordered structure is the prediction of secondary structural stability for randomized sequences which, in the original implementation of ScanFold, is explicitly evaluated. This slow process has limited the rapid identification of ordered structures in large genomes/transcriptomes, which we seek to overcome in this current work introducing ScanFold 2.0. In this revised version of ScanFold, we no longer explicitly evaluate randomized sequence folding energy, but rather estimate it using a machine learning approach. For high randomization numbers, this can increase prediction speeds over 100-fold compared to ScanFold 1.0, allowing for the analysis of large sequences, as well as the use of additional folding algorithms that may be computationally expensive. In the testing of ScanFold 2.0, we re-evaluate the Zika, HIV, and SARS-CoV-2 genomes and compare both the consistency of results and the time of each run to ScanFold 1.0. We also re-evaluate the SARS-CoV-2 genome to assess the quality of ScanFold 2.0 predictions vs several biochemical structure probing datasets and compare the results to those of the original ScanFold program.

Advances and opportunities in RNA structure experimental determination and computational modeling

Beyond transferring genetic information, RNAs are molecules with diverse functions that include catalyzing biochemical reactions and regulating gene expression. Most of these activities depend on RNAs' specific structures. Therefore, accurately determining RNA structure is integral to advancing our understanding of RNA functions. Here, we summarize the state-of-the-art experimental and computational technologies developed to evaluate RNA secondary and tertiary structures. We also highlight how the rapid increase of experimental data facilitates the integrative modeling approaches for better resolving RNA structures. Finally, we provide our thoughts on the latest advances and challenges in RNA structure determination methods, as well as on future directions for both experimental approaches and artificial intelligence-based computational tools to model RNA structure. Ultimately, we hope the technological advances will deepen our understanding of RNA biology and facilitate RNA structure-based biomedical research such as designing specific RNA structures for therapeutics and deploying RNA-targeting small-molecule drugs.

The in vivo RNA structurome of the malaria parasite Plasmodium falciparum, a protozoan with an A/U-rich transcriptome

Plasmodium falciparum, a protozoan parasite and causative agent of human malaria, has one of the most A/T-biased genomes sequenced to date. This may give the genome and the transcriptome unusual structural features. Recent progress in sequencing techniques has made it possible to study the secondary structures of RNA molecules at the transcriptomic level. Thus, in this study we produced the in vivo RNA structurome of a protozoan parasite with a highly A/U-biased transcriptome. We showed that it is possible to probe the secondary structures of P. falciparum RNA molecules in vivo using two different chemical probes, and obtained structures for more than half of all transcripts in the transcriptome. These showed greater stability (lower free energy) than the same structures modelled in silico, and structural features appeared to influence translation efficiency and RNA decay. Finally, we compared the P. falciparum RNA structurome with the predicted RNA structurome of an A/U-balanced species, P. knowlesi, finding a bias towards lower overall transcript stability and more hairpins and multi-stem loops in P. falciparum. This unusual protozoan RNA structurome will provide a basis for similar studies in other protozoans and also in other unusual genomes.

Computational and Experimental Approaches to Study the RNA Secondary Structures of RNA Viruses

Most pandemics of recent decades can be traced to RNA viruses, including HIV, SARS, influenza, dengue, Zika, and SARS-CoV-2. These RNA viruses impose considerable social and economic burdens on our society, resulting in a high number of deaths and high treatment costs. As these RNA viruses utilize an RNA genome, which is important for different stages of the viral life cycle, including replication, translation, and packaging, studying how the genome folds is important to understand virus function. In this review, we summarize recent advances in computational and high-throughput RNA structure-mapping approaches and their use in understanding structures within RNA virus genomes. In particular, we focus on the genome structures of the dengue, Zika, and SARS-CoV-2 viruses due to recent significant outbreaks of these viruses around the world.

Recent advances in RNA structurome

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

RNAvigator: A Pipeline to Identify Candidates for Functional RNA Structure Elements

Identifying structural elements in long and complex RNAs, such as long non-coding and RNA viruses, can shed light on the functionality and mechanisms of such RNAs. Here we present RNAvigator, a tool able to identify elements of structural importance by using experimental SHAPE data or SHAPE-like predictions in conjunction with stability and entropy assessments. RNAvigator recognizes regions that are the most stable, unambiguous, and structured on RNA molecules, and thus potentially functional. When relying on predictions, RNAvigator uses the CROSS algorithm, a neural network trained on experimental data that achieved an AUC of 0.74 on hepatitis C virus SHAPE-MaP data and which was able to improve the predictive power of Superfold. By using RNAvigator, we can identify known functional regions on the complete hepatitis C virus genome, including the regulatory regions CRE and IRES, and the 3’ UTR of dengue virus, a region known for the presence of structural elements essential for its replication, and functional regions of long non-coding RNAs such as XIST and HOTAIR. We envision that RNAvigator will be a useful tool for studying long and complex RNA molecules using known chemical probing data or, if they are not available, by employing predicted profiles.

A highly sensitive strand-specific multiplex RT-qPCR assay for quantitation of Zika virus replication

Reverse-transcription quantitative polymerase chain reaction (RT-qPCR) is widely used to quantify viral RNA genomes for diagnostics and research, yet conventional RT-qPCR protocols are unable to accurately distinguish between the different viral RNA species that exist during infection. Here we show that false-priming and self-priming occur during reverse transcription with several published Zika virus (ZIKV) primer sets. We developed a RT-qPCR assay using tagged primers and thermostable reverse transcriptase, which greatly reduced the occurrence of nonspecific cDNA products. Furthermore, we optimized the assay for use in multiplex qPCR which allows for simultaneous quantitative detection of positive-strand and negative-strand ZIKV RNA along with an internal control from both human and mosquito cells. Importantly, this assay is sensitive enough to study early stages of virus infection in vitro. Strikingly, using this assay, we detected ZIKV negative-strand RNA as early as 3 h post-infection in mammalian cell culture, at a time point prior to the onset of positive-strand RNA synthesis. Overall, the strand-specific RT-qPCR assay developed herein is a valuable tool to quantify ZIKV RNA and to study viral replication dynamics during infection. The application of these findings has the potential to increase accuracy of RNA detection methods for a variety of viral pathogens.

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.

SHAPE-guided RNA structure homology search and motif discovery

The rapidly growing popularity of RNA structure probing methods is leading to increasingly large amounts of available RNA structure information. This demands the development of efficient tools for the identification of RNAs sharing regions of structural similarity by direct comparison of their reactivity profiles, hence enabling the discovery of conserved structural features. We here introduce SHAPEwarp, a largely sequence-agnostic SHAPE-guided algorithm for the identification of structurally-similar regions in RNA molecules. Analysis of Dengue, Zika and coronavirus genomes recapitulates known regulatory RNA structures and identifies novel highly-conserved structural elements. This work represents a preliminary step towards the model-free search and identification of shared and conserved RNA structural features within transcriptomes.

Emerging roles of RNA-RNA interactions in transcriptional regulation

Pervasive transcription of the human genome generates a massive amount of noncoding RNAs (ncRNAs) that lack protein-coding potential but play crucial roles in development, differentiation, and tumorigenesis. To achieve these biological functions, ncRNAs must first fold into intricate structures via intramolecular RNA-RNA interactions (RRIs) and then interact with different RNA substrates via intermolecular RRIs. RRIs are usually facilitated, stabilized, or mediated by RNA-binding proteins. With this guiding principle, several protein-based high-throughput methods have been developed for unbiased mapping of defined or all RNA-binding protein-mediated RRIs in various species and cell lines. In addition, some chemical-based approaches are also powerful to detect RRIs globally based on the fact that RNA duplex can be cross-linked by psoralen or its derivative 4'-aminomethyltrioxsalen. These efforts have significantly expanded our understanding of RRIs in determining the specificity and variability of gene regulation. Here, we review the current knowledge of the regulatory roles of RRI, focusing on their emerging roles in transcriptional regulation and nuclear body formation. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Architecture of RNA-RNA interactions

RNA molecules tend to form intricate tertiary structures via intramolecular RNA-RNA interactions (RRIs) to regulate transcription, RNA processing, and translation processes. In these biological processes, RNAs, especially noncoding RNAs, usually achieve their regulatory specificity through intermolecular RNA-RNA base pairing and execute their regulatory outcomes via associated RNA-binding proteins. Decoding intramolecular and intermolecular RRIs is a prerequisite for understanding the architecture of various RNA molecules and their regulatory roles in development, differentiation, and disease. Many sequencing-based methods have recently been invented and have revealed extraordinarily complicated RRIs in mammalian cells. Here, we discuss the technical advances and limitations of various methodologies developed for studying cellular RRIs, with a focus on the emerging architectural roles of RRIs in gene regulation.

Evaluating RNA Structural Flexibility: Viruses Lead the Way

Our understanding of RNA structure has lagged behind that of proteins and most other biological polymers, largely because of its ability to adopt multiple, and often very different, functional conformations within a single molecule. Flexibility and multifunctionality appear to be its hallmarks. Conventional biochemical and biophysical techniques all have limitations in solving RNA structure and to address this in recent years we have seen the emergence of a wide diversity of techniques applied to RNA structural analysis and an accompanying appreciation of its ubiquity and versatility. Viral RNA is a particularly productive area to study in that this economy of function within a single molecule admirably suits the minimalist lifestyle of viruses. Here, we review the major techniques that are being used to elucidate RNA conformational flexibility and exemplify how the structure and function are, as in all biology, tightly linked.

In vivo structure and dynamics of the SARS-CoV-2 RNA genome

The dynamics of SARS-CoV-2 RNA structure and their functional relevance are largely unknown. Here we develop a simplified SPLASH assay and comprehensively map the in vivo RNA-RNA interactome of SARS-CoV-2 genome across viral life cycle. We report canonical and alternative structures including 5'-UTR and 3'-UTR, frameshifting element (FSE) pseudoknot and genome cyclization in both cells and virions. We provide direct evidence of interactions between Transcription Regulating Sequences, which facilitate discontinuous transcription. In addition, we reveal alternative short and long distance arches around FSE. More importantly, we find that within virions, while SARS-CoV-2 genome RNA undergoes intensive compaction, genome domains remain stable but with strengthened demarcation of local domains and weakened global cyclization. Taken together, our analysis reveals the structural basis for the regulation of replication, discontinuous transcription and translational frameshifting, the alternative conformations and the maintenance of global genome organization during the whole life cycle of SARS-CoV-2, which we anticipate will help develop better antiviral strategies.

Rhinovirus Inhibitors: Including a New Target, the Viral RNA

Rhinoviruses (RVs) are the main cause of recurrent infections with rather mild symptoms characteristic of the common cold. Nevertheless, RVs give rise to enormous numbers of absences from work and school and may become life-threatening in particular settings. Vaccination is jeopardised by the large number of serotypes eliciting only poorly cross-neutralising antibodies. Conversely, antivirals developed over the years failed FDA approval because of a low efficacy and/or side effects. RV species A, B, and C are now included in the fifteen species of the genus Enteroviruses based upon the high similarity of their genome sequences. As a result of their comparably low pathogenicity, RVs have become a handy model for other, more dangerous members of this genus, e.g., poliovirus and enterovirus 71. We provide a short overview of viral proteins that are considered potential drug targets and their corresponding drug candidates. We briefly mention more recently identified cellular enzymes whose inhibition impacts on RVs and comment novel approaches to interfere with infection via aggregation, virus trapping, or preventing viral access to the cell receptor. Finally, we devote a large part of this article to adding the viral RNA genome to the list of potential drug targets by dwelling on its structure, folding, and the still debated way of its exit from the capsid. Finally, we discuss the recent finding that G-quadruplex stabilising compounds impact on RNA egress possibly via obfuscating the unravelling of stable secondary structural elements.

Flavivirus enzymes and their inhibitors

Flaviviruses such as dengue, Japanese encephalitis, West Nile, Yellow Fever and Zika virus, cause viral hemorrhagic fever and encephalitis in humans. However, antiviral therapeutics to treat or prevent flavivirus infections are not yet available. Thus, there is pressing need to develop therapeutics and vaccines that target flavivirus infections. All flaviviruses carry a positive-sense single-stranded RNA genome, which encodes ten proteins; three structural proteins form the virus shell, and seven nonstructural (NS) proteins are involved in replication of the viral genome. While all NS proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) are part of a functional membrane-bound replication complex, enzymatic activities required for flaviviral replication reside in only two NS proteins, NS3 and NS5. NS3 functions as a protease, helicase, and triphosphatase, and NS5 as a capping enzyme, methyltransferase, and RNA-dependent RNA polymerase. In this chapter, we provide an overview of viral replication focusing on the structure and function of NS3 and NS5 replicases. We further describe strategies and examples of current efforts to identify potential flavivirus inhibitors against NS3 and NS5 enzymatic activities that can be developed as therapeutic agents to combat flavivirus infections.

Comprehensive mapping of SARS-CoV-2 interactions in vivo reveals functional virus-host interactions

SARS-CoV-2 is a major threat to global health. Here, we investigate the RNA structure and RNA-RNA interactions of wildtype (WT) and a mutant (Δ382) SARS-CoV-2 in cells using Illumina and Nanopore platforms. We identify twelve potentially functional structural elements within the SARS-CoV-2 genome, observe that subgenomic RNAs can form different structures, and that WT and Δ382 virus genomes fold differently. Proximity ligation sequencing identify hundreds of RNA-RNA interactions within the virus genome and between the virus and host RNAs. SARS-CoV-2 genome binds strongly to mitochondrial and small nucleolar RNAs and is extensively 2'-O-methylated. 2'-O-methylation sites are enriched in viral untranslated regions, associated with increased virus pair-wise interactions, and are decreased in host mRNAs upon virus infection, suggesting that the virus sequesters methylation machinery from host RNAs towards its genome. These studies deepen our understanding of the molecular and cellular basis of SARS-CoV-2 pathogenicity and provide a platform for targeted therapy.

RNA structure probing uncovers RNA structure-dependent biological functions

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

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

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

Global in situ profiling of RNA-RNA spatial interactions with RIC-seq

Emerging evidence has demonstrated that RNA-RNA interactions are vital in controlling diverse biological processes, including transcription, RNA splicing and protein translation. RNA in situ conformation sequencing (RIC-seq) is a technique for capturing protein-mediated RNA-RNA proximal interactions globally in living cells at single-base resolution. Cells are first treated with formaldehyde to fix all the protein-mediated RNA-RNA interactions in situ. After cell permeabilization and micrococcal nuclease digestion, the proximally interacting RNAs are 3' end-labeled with pCp-biotin and subsequently ligated using T4 RNA ligase. The chimeric RNAs are then enriched and converted into libraries for paired-end sequencing. After deep sequencing, computational analysis yields interaction strength scores for every base on proximally interacting RNAs in the starting populations. The whole experimental procedure is designed to be completed within 6 d, followed by an additional 8 d for computational analysis. RIC-seq technology can unbiasedly detect intra- and intermolecular RNA-RNA interactions, thereby rendering it useful for reconstructing RNA higher-order structures and identifying direct noncoding RNA targets.

Molecular Insights into the Flavivirus Replication Complex

Flaviviruses are vector-borne RNA viruses, many of which are clinically relevant human viral pathogens, such as dengue, Zika, Japanese encephalitis, West Nile and yellow fever viruses. Millions of people are infected with these viruses around the world each year. Vaccines are only available for some members of this large virus family, and there are no effective antiviral drugs to treat flavivirus infections. The unmet need for vaccines and therapies against these flaviviral infections drives research towards a better understanding of the epidemiology, biology and immunology of flaviviruses. In this review, we discuss the basic biology of the flavivirus replication process and focus on the molecular aspects of viral genome replication. Within the virus-induced intracellular membranous compartments, flaviviral RNA genome replication takes place, starting from viral poly protein expression and processing to the assembly of the virus RNA replication complex, followed by the delivery of the progeny viral RNA to the viral particle assembly sites. We attempt to update the latest understanding of the key molecular events during this process and highlight knowledge gaps for future studies.

Evaluation of Conserved RNA Secondary Structures within and between Geographic Lineages of Zika Virus

Zika virus (ZIKV), without a vaccine or an effective treatment approved to date, has globally spread in the last century. The infection caused by ZIKV in humans has changed progressively from mild to subclinical in recent years, causing epidemics with greater infectivity, tropism towards new tissues and other related symptoms as a product of various emergent ZIKV–host cell interactions. However, it is still unknown why or how the RNA genome structure impacts those interactions in differential evolutionary origin strains. Moreover, the genomic comparison of ZIKV strains from the sequence-based phylogenetic analysis is well known, but differences from RNA structure comparisons have barely been studied. Thus, in order to understand the RNA genome variability of lineages of various geographic distributions better, 410 complete genomes in a phylogenomic scanning were used to study the conservation of structured RNAs. Our results show the contemporary landscape of conserved structured regions with unique conserved structured regions in clades or in lineages within circulating ZIKV strains. We propose these structures as candidates for further experimental validation to establish their potential role in vital functions of the viral cycle of ZIKV and their possible associations with the singularities of different outbreaks that lead to ZIKV populations to acquire nucleotide substitutions, which is evidence of the local structure genome differentiation.

In vivo structure of the Ty1 retrotransposon RNA genome

Long terminal repeat (LTR)-retrotransposons constitute a significant part of eukaryotic genomes and influence their function and evolution. Like other RNA viruses, LTR-retrotransposons efficiently utilize their RNA genome to interact with host cell machinery during replication. Here, we provide the first genome-wide RNA secondary structure model for a LTR-retrotransposon in living cells. Using SHAPE probing, we explore the secondary structure of the yeast Ty1 retrotransposon RNA genome in its native in vivo state and under defined in vitro conditions. Comparative analyses reveal the strong impact of the cellular environment on folding of Ty1 RNA. In vivo, Ty1 genome RNA is significantly less structured and more dynamic but retains specific well-structured regions harboring functional cis-acting sequences. Ribosomes participate in the unfolding and remodeling of Ty1 RNA, and inhibition of translation initiation stabilizes Ty1 RNA structure. Together, our findings support the dual role of Ty1 genomic RNA as a template for protein synthesis and reverse transcription. This study also contributes to understanding how a complex multifunctional RNA genome folds in vivo, and strengthens the need for studying RNA structure in its natural cellular context.

The Pseudo-Circular Genomes of Flaviviruses: Structures, Mechanisms, and Functions of Circularization

The circularization of viral genomes fulfills various functions, from evading host defense mechanisms to promoting specific replication and translation patterns supporting viral proliferation. Here, we describe the genomic structures and associated host factors important for flaviviruses genome circularization and summarize their functional roles. Flaviviruses are relatively small, single-stranded, positive-sense RNA viruses with genomes of approximately 11 kb in length. These genomes contain motifs at their 5′ and 3′ ends, as well as in other regions, that are involved in circularization. These motifs are highly conserved throughout the Flavivirus genus and occur both in mature virions and within infected cells. We provide an overview of these sequence motifs and RNA structures involved in circularization, describe their linear and circularized structures, and discuss the proteins that interact with these circular structures and that promote and regulate their formation, aiming to clarify the key features of genome circularization and understand how these affect the flaviviruses life cycle.

The global and local distribution of RNA structure throughout the SARS-CoV-2 genome

SARS-CoV-2 is the causative viral agent of COVID-19, the disease at the center of the current global pandemic. While knowledge of highly structured regions is integral for mechanistic insights into the viral infection cycle, very little is known about the location and folding stability of functional elements within the massive, ∼30kb SARS-CoV-2 RNA genome. In this study, we analyze the folding stability of this RNA genome relative to the structural landscape of other well-known viral RNAs. We present an in-silico pipeline to predict regions of high base pair content across long genomes and to pinpoint hotspots of well-defined RNA structures, a method that allows for direct comparisons of RNA structural complexity within the several domains in SARS-CoV-2 genome. We report that the SARS-CoV-2 genomic propensity for stable RNA folding is exceptional among RNA viruses, superseding even that of HCV, one of the most structured viral RNAs in nature. Furthermore, our analysis suggests varying levels of RNA structure across genomic functional regions, with accessory and structural ORFs containing the highest structural density in the viral genome. Finally, we take a step further to examine how individual RNA structures formed by these ORFs are affected by the differences in genomic and subgenomic contexts, which given the technical difficulty of experimentally separating cellular mixtures of sgRNA from gRNA, is a unique advantage of our in-silico pipeline. The resulting findings provide a useful roadmap for planning focused empirical studies of SARS-CoV-2 RNA biology, and a preliminary guide for exploring potential SARS-CoV-2 RNA drug targets.Importance The RNA genome of SARS-CoV-2 is among the largest and most complex viral genomes, and yet its RNA structural features remain relatively unexplored. Since RNA elements guide function in most RNA viruses, and they represent potential drug targets, it is essential to chart the architectural features of SARS-CoV-2 and pinpoint regions that merit focused study. Here we show that RNA folding stability of SARS-CoV-2 genome is exceptional among viral genomes and we develop a method to directly compare levels of predicted secondary structure across SARS-CoV-2 domains. Remarkably, we find that coding regions display the highest structural propensity in the genome, forming motifs that differ between the genomic and subgenomic contexts. Our approach provides an attractive strategy to rapidly screen for candidate structured regions based on base pairing potential and provides a readily interpretable roadmap to guide functional studies of RNA viruses and other pharmacologically relevant RNA transcripts.

Accurate prediction of genome-wide RNA secondary structure profile based on extreme gradient boosting

MOTIVATION

RNA secondary structure plays a vital role in fundamental cellular processes, and identification of RNA secondary structure is a key step to understand RNA functions. Recently, a few experimental methods were developed to profile genome-wide RNA secondary structure, i.e. the pairing probability of each nucleotide, through high-throughput sequencing techniques. However, these high-throughput methods have low precision and cannot cover all nucleotides due to limited sequencing coverage.

RESULTS

Here, we have developed a new method for the prediction of genome-wide RNA secondary structure profile from RNA sequence based on the extreme gradient boosting technique. The method achieves predictions with areas under the receiver operating characteristic curve (AUC) >0.9 on three different datasets, and AUC of 0.888 by another independent test on the recently released Zika virus data. These AUCs are consistently >5% greater than those by the CROSS method recently developed based on a shallow neural network. Further analysis on the 1000 Genome Project data showed that our predicted unpaired probabilities are highly correlated (>0.8) with the minor allele frequencies at synonymous, non-synonymous mutations, and mutations in untranslated regions, which were higher than those generated by RNAplfold. Moreover, the prediction over all human mRNA indicated a consistent result with previous observation that there is a periodic distribution of unpaired probability on codons. The accurate predictions by our method indicate that such model trained on genome-wide experimental data might be an alternative for analytical methods.

AVAILABILITY AND IMPLEMENTATION

The GRASP is available for academic use at https://github.com/sysu-yanglab/GRASP.

SUPPLEMENTARY INFORMATION

Supplementary data are available online.