COMRADES determines in vivo RNA structures and interactions

bpRNA-CosMoS: a robust and efficient RNA structural comparison method using k-mer based cosine similarity

MOTIVATION

RNA secondary structure is often essential to function. Recent work has led to the development of high-throughput experimental probing methods for structure determination. Although structure is more conserved than primary sequence, much of the bioinformatics pipelines to connect RNA structure to function rely on nucleotide sequence alignments rather than structural similarity. There is a need to develop methods for secondary structure comparisons that are also fast and efficient to navigate the vast amounts of structural data. K-mer based similarity approaches are valued for their computational efficiency and have been applied for protein, DNA, and RNA primary sequences. However, these approaches have yet to be implemented for RNA secondary structure.

RESULTS

Our method, bpRNA-CosMoS, fills this gap by using k-mers and length-weighted cosine similarity to compute similarity scores between RNA structures. bpRNA-CosMoS is built upon the bpRNA structure array, which represents the structural category of each nucleotide as a single-character structural code (e.g. hairpin=H, etc.). A structural comparison score is calculated through cosine similarity of the k-mer count vectors, generated from structure arrays. A major challenge with k-mer based methods is that they often ignore the length of the sequences being compared. We have overcome this with a length-weighted penalty that addresses cases of two RNAs of vastly different lengths. In addition, the use of "fuzzy counting" has added some optional flexibility to decrease the negative impact that small structural variations have on the similarity score. This results in a robust and efficient way to identify structural comparisons across large datasets.

AVAILABILITY AND IMPLEMENTATION

The code and application guidelines of bpRNA-CosMoS are made available at github (https://github.com/BLasher113/bpRNA-CosMoS) and Zenodo (10.5281/zenodo.14715285).

IPANEMAP Suite: a pipeline for probing-informed RNA structure modeling

In addition to their sequence, multiple functions of RNAs are encoded within their structure, which is often difficult to solve using physico-chemical methods. Incorporating low-resolution experimental data such as chemical probing into computational prediction significantly enhances RNA structure modeling accuracy. While medium- and high-throughput RNA structure probing techniques are widely accessible, the subsequent analysis process can be cumbersome, involving multiple software and manual data manipulation. In addition, the relevant interpretation of the data requires proper parameterization of the software and a strict consistency in the analysis pipeline. To streamline such workflows, we introduce IPANEMAP Suite, a comprehensive platform that guides users from chemically probing raw data to visually informative secondary structure models. IPANEMAP Suite seamlessly integrates various experimental datasets and facilitates comparative analysis of RNA structures under different conditions (footprinting), aiding in the study of protein or small molecule interactions with RNA. Here, we show that the unique ability of IPANEMAP Suite to perform integrative modeling using several chemical probing datasets with phylogenetic data can be instrumental in obtaining accurate secondary structure models. The platform's project-based approach ensures full traceability and generates publication-quality outputs, simplifying the entire RNA structure analysis process. IPANEMAP Suite is freely available at https://github.com/Sargueil-CiTCoM/ipasuite under a GPL-3.0 license.

Nuclear RNA-binding proteins meet cytoplasmic viruses

Cytoplasmic viruses interact intricately with the nuclear pore complex and nuclear import/export machineries, affecting nuclear-cytoplasmic trafficking. This can lead to the selective accumulation of nuclear RNA-binding proteins (RBPs) in the cytoplasm. Pioneering research has shown that relocated RBPs serve as an intrinsic defense mechanism against viruses, which involves RNA export, splicing, and nucleolar factors. For instance, the U2 small nuclear ribonucleoprotein (snRNP) relocates to the cytoplasm in infected cells and uses U2 snRNA to interact with viral genomes, repressing viral replication and gene expression. Here, we describe these emerging host-virus interactions and discuss the remaining questions to elucidate their antiviral mechanisms.

Structural features within the NORAD long noncoding RNA underlie efficient repression of Pumilio activity

Long noncoding RNAs (lncRNAs) are increasingly appreciated for their important functions in mammalian cells. However, how their functional capacities are encoded in their sequences and manifested in their structures remains largely unknown. Some lncRNAs bind to and modulate the availability of RNA-binding proteins, but the structural principles that underlie this mode of regulation are unknown. The NORAD lncRNA is a known decoy for Pumilio proteins, which modulate the translation and stability of hundreds of messenger RNAs and, consequently, a regulator of genomic stability and aging. Here we probed the RNA structure and long-range RNA-RNA interactions formed by human NORAD inside cells under different stressful conditions. We discovered a highly modular structure consisting of well-defined domains that contribute independently to NORAD function. Following arsenite stress, most structural domains undergo relaxation and form interactions with other RNAs that are targeted to stress granules. We further revealed a unique structural organization that spatially clusters the multiple Pumilio binding sites along NORAD and consequently contributes to the derepression of Pumilio targets. We then applied these structural principles to design an effective artificial decoy for the let-7 microRNA. Our work demonstrates how the sequence of a lncRNA spatially clusters its function into separated domains and how structural principles can be employed for the rational design of lncRNAs with desired activities.

snoRNA-facilitated protein secretion revealed by transcriptome-wide snoRNA target identification

Small nucleolar RNAs (snoRNAs) are non-coding RNAs known for guiding RNA modifications, including 2'-O-methylation (Nm) and pseudouridine (Ψ). While snoRNAs may also interact with other RNAs, such as mRNA, the full repertoire of RNAs targeted by snoRNA remains elusive due to the lack of effective technologies that identify snoRNA targets transcriptome wide. Here, we develop a chemical crosslinking-based approach that comprehensively detects cellular RNA targets of snoRNAs, yielding thousands of previously unrecognized snoRNA-mRNA interactions in human cells and mouse brain tissues. Many interactions occur outside of snoRNA-guided RNA modification sites, hinting at non-canonical functions beyond RNA modification. We find that one of these snoRNAs, SNORA73, targets mRNAs that encode secretory proteins and membrane proteins. SNORA73 also interacts with 7SL RNA, part of the signal recognition particle (SRP) required for protein secretion. The mRNA-SNORA73-7SL RNA interactions enhance the association of the SNORA73-target mRNAs with SRP, thereby facilitating the secretion of encoded proteins.

Targeted Covalent Modification Strategies for Drugging the Undruggable Targets

The term "undruggable" refers to proteins or other biological targets that have been historically challenging to target with conventional drugs or therapeutic strategies because of their structural, functional, or dynamic properties. Drugging such undruggable targets is essential to develop new therapies for diseases where current treatment options are limited or nonexistent. Thus, investigating methods to achieve such drugging is an important challenge in medicinal chemistry. Among the numerous methodologies for drug discovery, covalent modification of therapeutic targets has emerged as a transformative strategy. The covalent attachment of diverse functional molecules to targets provides a powerful platform for creating highly potent drugs and chemical tools as well the ability to provide valuable information on the structures and dynamics of undruggable targets. In this review, we summarize recent examples of chemical methods for the covalent modification of proteins and other biomolecules for the development of new therapeutics and to overcome drug discovery challenges and highlight how such methods contribute toward the drugging of undruggable targets. In particular, we focus on the use of covalent chemistry methods for the development of covalent drugs, target identification, drug screening, artificial modulation of post-translational modifications, cancer specific chemotherapies, and nucleic acid-based therapeutics.

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

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

Structural basis for the participation of the SARS-CoV-2 nucleocapsid protein in the template switch mechanism and genomic RNA reorganization

The COVID-19 pandemic has resulted in a significant toll of deaths worldwide, exceeding seven million individuals, prompting intensive research efforts aimed at elucidating the molecular mechanisms underlying the pathogenesis of SARS-CoV-2 infection. Despite the rapid development of effective vaccines and therapeutic interventions, COVID-19 remains a threat to humans due to the emergence of novel variants and largely unknown long-term consequences. Among the viral proteins, the nucleocapsid protein (N) stands out as the most conserved and abundant, playing the primary role in nucleocapsid assembly and genome packaging. The N protein is promiscuous for the recognition of RNA, yet it can perform specific functions. Here, we discuss the structural basis of specificity, which is directly linked to its regulatory role. Notably, the RNA chaperone activity of N is central to its multiple roles throughout the viral life cycle. This activity encompasses double-stranded RNA (dsRNA) annealing and melting and facilitates template switching, enabling discontinuous transcription. N also promotes the formation of membrane-less compartments through liquid-liquid phase separation, thereby facilitating the congregation of the replication and transcription complex. Considering the information available regarding the catalytic activities and binding signatures of the N protein-RNA interaction, this review focuses on the regulatory role of the SARS-CoV-2 N protein. We emphasize the participation of the N protein in discontinuous transcription, template switching, and RNA chaperone activity, including double-stranded RNA melting and annealing activities.

The role of structure in regulatory RNA elements

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

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.

Technological advancements in deciphering RNA-RNA interactions

RNA-RNA interactions (RRIs) can dictate RNA molecules to form intricate higher-order structures and bind their RNA substrates in diverse biological processes. To elucidate the function, binding specificity, and regulatory mechanisms of various RNA molecules, especially the vast repertoire of non-coding RNAs, advanced technologies and methods that globally map RRIs are extremely valuable. In the past decades, many state-of-the-art technologies have been developed for this purpose. This review focuses on those high-throughput technologies for the global mapping of RRIs. We summarize the key concepts and the pros and cons of different technologies. In addition, we highlight the novel biological insights uncovered by these RRI mapping methods and discuss the future challenges for appreciating the crucial roles of RRIs in gene regulation across bacteria, viruses, archaea, and mammals.

Identification of RNA structures and their roles in RNA functions

The development of high-throughput RNA structure profiling methods in the past decade has greatly facilitated our ability to map and characterize different aspects of RNA structures transcriptome-wide in cell populations, single cells and single molecules. The resulting high-resolution data have provided insights into the static and dynamic nature of RNA structures, revealing their complexity as they perform their respective functions in the cell. In this Review, we discuss recent technical advances in the determination of RNA structures, and the roles of RNA structures in RNA biogenesis and functions, including in transcription, processing, translation, degradation, localization and RNA structure-dependent condensates. We also discuss the current understanding of how RNA structures could guide drug design for treating genetic diseases and battling pathogenic viruses, and highlight existing challenges and future directions in RNA structure research.

Structural and biophysical dissection of RNA conformational ensembles

RNA's ability to form and interconvert between multiple secondary and tertiary structures is critical to its functional versatility and the traditional view of RNA structures as static entities has shifted towards understanding them as dynamic conformational ensembles. In this review we discuss RNA structural ensembles and their dynamics, highlighting the concept of conformational energy landscapes as a unifying framework for understanding RNA processes such as folding, misfolding, conformational changes, and complex formation. Ongoing advancements in cryo-electron microscopy and chemical probing techniques are significantly enhancing our ability to investigate multiple structures adopted by conformationally dynamic RNAs, while traditional methods such as nuclear magnetic resonance spectroscopy continue to play a crucial role in providing high-resolution, quantitative spatial and temporal information. We discuss how these methods, when used synergistically, can provide a comprehensive understanding of RNA conformational ensembles, offering new insights into their regulatory functions.

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

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

Chemical and Enzyme-Mediated Chemical Reactions for Studying Nucleic Acids and Their Modifications

Nucleic acids are genetic information-carrying molecules inside cells. Apart from basic nucleotide building blocks, there exist various naturally occurring chemical modifications on nucleobase and ribose moieties, which greatly increase the encoding complexity of nuclei acids, contribute to the alteration of nucleic acid structures, and play versatile regulation roles in gene expression. To study the functions of certain nucleic acids in various biological contexts, robust tools to specifically label and identify these macromolecules and their modifications, and to illuminate their structures are highly necessary. In this review, we summarize recent technique advances of using chemical and enzyme-mediated chemical reactions to study nucleic acids and their modifications and structures. By highlighting the chemical principles of these techniques, we aim to present a perspective on the advancement of the field as well as to offer insights into developing specific chemical reactions and precise enzyme catalysis utilized for nucleic acids and their modifications.

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

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

Potential Transcriptional Enhancers in Coronaviruses: From Infectious Bronchitis Virus to SARS-CoV-2

Coronaviruses constitute a global threat to human and animal health. It is essential to investigate the long-distance RNA-RNA interactions that approximate remote regulatory elements in strategies, including genome circularization, discontinuous transcription, and transcriptional enhancers, aimed at the rapid replication of their large genomes, pathogenicity, and immune evasion. Based on the primary sequences and modeled RNA-RNA interactions of two experimentally defined coronaviral enhancers, we detected via an in silico primary and secondary structural analysis potential enhancers in various coronaviruses, from the phylogenetically ancient avian infectious bronchitis virus (IBV) to the recently emerged SARS-CoV-2. These potential enhancers possess a core duplex-forming region that could transition between closed and open states, as molecular switches directed by viral or host factors. The duplex open state would pair with remote sequences in the viral genome and modulate the expression of downstream crucial genes involved in viral replication and host immune evasion. Consistently, variations in the predicted IBV enhancer region or its distant targets coincide with cases of viral attenuation, possibly driven by decreased open reading frame (ORF)3a immune evasion protein expression. If validated experimentally, the annotated enhancer sequences could inform structural prediction tools and antiviral interventions.

Defective viral genomes: advances in understanding their generation, function, and impact on infection outcomes

Defective viral genomes (DVGs) are truncated derivatives of their parental viral genomes generated during an aberrant round of viral genomic replication. Distinct classes of DVGs have been identified in most families of both positive- and negative-sense RNA viruses. Importantly, DVGs have been detected in clinical samples from virally infected individuals and an emerging body of association studies implicates DVGs in shaping the severity of disease caused by viral infections in humans. Consequently, there is growing interest in understanding the molecular mechanisms of de novo DVG generation, how DVGs interact with the innate immune system, and harnessing DVGs as novel therapeutics and vaccine adjuvants to attenuate viral pathogenesis. This minireview focuses on single-stranded RNA viruses (excluding retroviridae), and summarizes the current knowledge of DVG generation, the functions and diversity of DVG species, the roles DVGs play in influencing disease progression, and their application as antivirals and vaccine adjuvants.

rnaCrosslinkOO: an object-oriented R package for the analysis of RNA structural data generated by RNA crosslinking experiments

SUMMARY

RNA (ribonucleic acid) molecules have secondary and tertiary structures in vivo which play a crucial role in cellular processes such as the regulation of gene expression, RNA processing and localization. The ability to investigate these structures will enhance our understanding of their function and contribute to the diagnosis and treatment of diseases caused by RNA dysregulation. However, there are no mature pipelines or packages for processing and analyzing complex in vivo RNA structural data. Here, we present rnaCrosslinkOO (RNA Crosslink Object-Oriented), a novel software package for the comprehensive analysis of data derived from the COMRADES (Crosslinking of Matched RNA and Deep Sequencing) method. rnaCrosslinkOO offers a comprehensive pipeline from raw sequencing reads to the identification and comparison of RNA structural features. It includes read processing and alignment, clustering of duplexes, data exploration, folding and comparisons of RNA structures. rnaCrosslinkOO also enables comparisons between conditions, the identification of inter-RNA interactions, and the incorporation of reactivity data to improve structure prediction.

AVAILABILITY AND IMPLEMENTATION

rnaCrosslinkOO is freely available to noncommercial users and implemented in R, with the source code and documentation accessible at https://CRAN.R-project.org/package=rnaCrosslinkOO. The software is supported on Linux, macOS, and Windows platforms.

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.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability

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

Bioinformatics Insights on Viral Gene Expression Transactivation: From HIV-1 to SARS-CoV-2

Viruses provide vital insights into gene expression control. Viral transactivators, with other viral and cellular proteins, regulate expression of self, other viruses, and host genes with profound effects on infected cells, underlying inflammation, control of immune responses, and pathogenesis. The multifunctional Tat proteins of lentiviruses (HIV-1, HIV-2, and SIV) transactivate gene expression by recruiting host proteins and binding to transacting responsive regions (TARs) in viral and host RNAs. SARS-CoV-2 nucleocapsid participates in early viral transcription, recruits similar cellular proteins, and shares intracellular, surface, and extracellular distribution with Tat. SARS-CoV-2 nucleocapsid interacting with the replication-transcription complex might, therefore, transactivate viral and cellular RNAs in the transcription and reactivation of self and other viruses, acute and chronic pathogenesis, immune evasion, and viral evolution. Here, we show, by using primary and secondary structural comparisons, that the leaders of SARS-CoV-2 and other coronaviruses contain TAR-like sequences in stem-loops 2 and 3. The coronaviral nucleocapsid C-terminal domains harbor a region of similarity to TAR-binding regions of lentiviral Tat proteins, and coronaviral nonstructural protein 12 has a cysteine-rich metal binding, dimerization domain, as do lentiviral Tat proteins. Although SARS-CoV-1 nucleocapsid transactivated gene expression in a replicon-based study, further experimental evidence for coronaviral transactivation and its possible implications is warranted.

RNA structure: implications in viral infections and neurodegenerative diseases

RNA is an intermediary between DNA and protein, a catalyzer of biochemical reactions, and a regulator of genes and transcripts. RNA structures are essential for complicated functions. Recent years have witnessed rapid advancements in RNA secondary structure probing techniques. These technological strides provided comprehensive insights into RNA structures, which significantly contributed to our understanding of diverse cellular regulatory processes, including gene regulation, epigenetic regulation, and post-transactional regulation. Meanwhile, they have facilitated the creation of therapeutic tools for tackling human diseases. Despite their therapeutic applications, RNA structure probing methods also offer a promising avenue for exploring the mechanisms of human diseases, potentially providing the key to overcoming existing research constraints and obtaining the in-depth information necessary for a deeper understanding of disease mechanisms.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

Decoding the genome of SARS-CoV-2: a pathway to drug development through translation inhibition

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes the coronavirus disease 2019 (COVID-19) pandemic and is continuously spreading globally. The continuous emergence of new SARS-CoV-2 variants keeps posing threats, highlighting the need for fast-acting, mutation-resistant broad-spectrum therapeutics. Protein translation is vital for SARS-CoV-2 replication, producing early non-structural proteins for RNA replication and transcription, and late structural proteins for virion assembly. Targeted blocking of viral protein translation is thus a potential approach to developing effective anti-SARS-CoV-2 drugs. SARS-CoV-2, as an obligate parasite, utilizes the host's translation machinery. Translation-blocking strategies that target the SARS-CoV-2 mRNA, especially those that target its conserved elements are generally preferred. In this review, we discuss the current understanding of SARS-CoV-2 translation, highlighting the important conserved motifs and structures involved in its regulation. We also discuss the current strategies for blocking SARS-CoV-2 translation through viral RNA degradation or RNA element dysfunction.

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.

Dengue and Zika RNA-RNA interactomes reveal pro- and anti-viral RNA in human cells

BACKGROUND

Identifying host factors is key to understanding RNA virus pathogenicity. Besides proteins, RNAs can interact with virus genomes to impact replication.

RESULTS

Here, we use proximity ligation sequencing to identify virus-host RNA interactions for four strains of Zika virus (ZIKV) and one strain of dengue virus (DENV-1) in human cells. We find hundreds of coding and non-coding RNAs that bind to DENV and ZIKV viruses. Host RNAs tend to bind to single-stranded regions along the virus genomes according to hybridization energetics. Compared to SARS-CoV-2 interactors, ZIKV-interacting host RNAs tend to be downregulated upon virus infection. Knockdown of several short non-coding RNAs, including miR19a-3p, and 7SK RNA results in a decrease in viral replication, suggesting that they act as virus-permissive factors. In addition, the 3'UTR of DYNLT1 mRNA acts as a virus-restrictive factor by binding to the conserved dumbbell region on DENV and ZIKV 3'UTR to decrease virus replication. We also identify a conserved set of host RNAs that interacts with DENV, ZIKV, and SARS-CoV-2, suggesting that these RNAs are broadly important for RNA virus infection.

CONCLUSIONS

This study demonstrates that host RNAs can impact virus replication in permissive and restrictive ways, expanding our understanding of host factors and RNA-based gene regulation during viral pathogenesis.

Long-range RNA structures in the human transcriptome beyond evolutionarily conserved regions

RNA structure has been increasingly recognized as a critical player in the biogenesis and turnover of many transcripts classes. In eukaryotes, the prediction of RNA structure by thermodynamic modeling meets fundamental limitations due to the large sizes and complex, discontinuous organization of eukaryotic genes. Signatures of functional RNA structures can be found by detecting compensatory substitutions in homologous sequences, but a comparative approach is applicable only within conserved sequence blocks. Here, we developed a computational pipeline called PHRIC, which is not limited to conserved regions and relies on RNA contacts derived from RNA in situ conformation sequencing (RIC-seq) experiments. It extracts pairs of short RNA fragments surrounded by nested clusters of RNA contacts and predicts long, nearly perfect complementary base pairings formed between these fragments. In application to a panel of RIC-seq experiments in seven human cell lines, PHRIC predicted ~12,000 stable long-range RNA structures with equilibrium free energy below -15 kcal/mol, the vast majority of which fall outside of regions annotated as conserved among vertebrates. These structures, nevertheless, show some level of sequence conservation and remarkable compensatory substitution patterns in other clades. Furthermore, we found that introns have a higher propensity to form stable long-range RNA structures between each other, and moreover that RNA structures tend to concentrate within the same intron rather than connect adjacent introns. These results for the first time extend the application of proximity ligation assays to RNA structure prediction beyond conserved regions.

Versatile toolkit for highly-efficient and scarless overexpression of circular RNAs

Circular RNAs (circRNAs) are a class of single-stranded, covalently closed RNA that contain a unique back-splice junction (bsj) sequence created by the ligation of their 5' and 3' ends via spliceosome-catalyzed back-splicing. A key step in illuminating the cellular roles of specific circRNAs is via increasing their expression. This is frequently done by transfecting cells with plasmid DNA containing cloned exons from which the circRNA is transcribed, flanked by sequences that promote back-splicing. We observed that commonly used plasmids lead to the production of circRNAs with molecular scars at the circRNA bsj. Stepwise redesign of the cloning vector corrected this problem, ensuring bona fide circRNAs are produced with their natural bsj at high efficiency. The fidelity of circRNAs produced from this new construct was validated by RNA sequencing and also functionally validated. To increase the utility of this modified resource for expressing circRNA, we developed an expanded set of vectors incorporating this design that (i) enables selection with a variety of antibiotics and fluorescent proteins, (ii) employs a range of promoters varying in promoter strength and (iii) generated a complementary set of lentiviral plasmids for difficult-to-transfect cells. These resources provide a novel and versatile toolkit for high-efficiency and scarless overexpression of circular RNAs that fulfill a critical need for the investigation of circRNA function.

tRNA renovatio: Rebirth through fragmentation

tRNA function is based on unique structures that enable mRNA decoding using anticodon trinucleotides. These structures interact with specific aminoacyl-tRNA synthetases and ribosomes using 3D shape and sequence signatures. Beyond translation, tRNAs serve as versatile signaling molecules interacting with other RNAs and proteins. Through evolutionary processes, tRNA fragmentation emerges as not merely random degradation but an act of recreation, generating specific shorter molecules called tRNA-derived small RNAs (tsRNAs). These tsRNAs exploit their linear sequences and newly arranged 3D structures for unexpected biological functions, epitomizing the tRNA "renovatio" (from Latin, meaning renewal, renovation, and rebirth). Emerging methods to uncover full tRNA/tsRNA sequences and modifications, combined with techniques to study RNA structures and to integrate AI-powered predictions, will enable comprehensive investigations of tRNA fragmentation products and new interaction potentials in relation to their biological functions. We anticipate that these directions will herald a new era for understanding biological complexity and advancing pharmaceutical engineering.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

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

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

RNA in situ conformation sequencing reveals novel long-range RNA structures with impact on splicing

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

The RNA interactome in the Hallmarks of Cancer

Ribonucleic acid (RNA) molecules are indispensable for cellular homeostasis in healthy and malignant cells. However, the functions of RNA extend well beyond that of a protein-coding template. Rather, both coding and non-coding RNA molecules function through critical interactions with a plethora of cellular molecules, including other RNAs, DNA, and proteins. Deconvoluting this RNA interactome, including the interacting partners, the nature of the interaction, and dynamic changes of these interactions in malignancies has yielded fundamental advances in knowledge and are emerging as a novel therapeutic strategy in cancer. Here, we present an RNA-centric review of recent advances in the field of RNA-RNA, RNA-protein, and RNA-DNA interactomic network analysis and their impact across the Hallmarks of Cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

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.

Machine learning modeling of RNA structures: methods, challenges and future perspectives

The three-dimensional structure of RNA molecules plays a critical role in a wide range of cellular processes encompassing functions from riboswitches to epigenetic regulation. These RNA structures are incredibly dynamic and can indeed be described aptly as an ensemble of structures that shifts in distribution depending on different cellular conditions. Thus, the computational prediction of RNA structure poses a unique challenge, even as computational protein folding has seen great advances. In this review, we focus on a variety of machine learning-based methods that have been developed to predict RNA molecules' secondary structure, as well as more complex tertiary structures. We survey commonly used modeling strategies, and how many are inspired by or incorporate thermodynamic principles. We discuss the shortcomings that various design decisions entail and propose future directions that could build off these methods to yield more robust, accurate RNA structure predictions.

Sequential disruption of SPLASH-identified vRNA-vRNA interactions challenges their role in influenza A virus genome packaging

A fundamental step in the influenza A virus (IAV) replication cycle is the coordinated packaging of eight distinct genomic RNA segments (i.e. vRNAs) into a viral particle. Although this process is thought to be controlled by specific vRNA-vRNA interactions between the genome segments, few functional interactions have been validated. Recently, a large number of potentially functional vRNA-vRNA interactions have been detected in purified virions using the RNA interactome capture method SPLASH. However, their functional significance in coordinated genome packaging remains largely unclear. Here, we show by systematic mutational analysis that mutant A/SC35M (H7N7) viruses lacking several prominent SPLASH-identified vRNA-vRNA interactions involving the HA segment package the eight genome segments as efficiently as the wild-type virus. We therefore propose that the vRNA-vRNA interactions identified by SPLASH in IAV particles are not necessarily critical for the genome packaging process, leaving the underlying molecular mechanism elusive.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

snoRNAs in hematopoiesis and blood malignancies: A comprehensive review

Small nucleolar RNAs (snoRNAs) are noncoding RNA molecules of highly variable size, usually ranging from 60 to 150 nucleotides. They are classified into H/ACA box snoRNAs, C/D box snoRNAs, and scaRNAs. Their functional profile includes biogenesis of ribosomes, processing of rRNAs, 2'-O-methylation and pseudouridylation of RNAs, alternative splicing and processing of mRNAs and the generation of small RNA molecules like miRNA. The snoRNAs have been observed to have an important role in hematopoiesis and malignant hematopoietic conditions including leukemia, lymphoma, and multiple myeloma. Blood malignancies arise in immune system cells or the bone marrow due to chromosome abnormalities. It has been estimated that annually over 1.25 million cases of blood cancer occur worldwide. The snoRNAs often show a differential expression profile in blood malignancies. Recent reports associate the abnormal expression of snoRNAs with the inhibition of apoptosis, uncontrolled cell proliferation, angiogenesis, and metastasis. This implies that targeting snoRNAs could be a potential way to treat hematologic malignancies. In this review, we describe the various functions of snoRNAs, their role in hematopoiesis, and the consequences of their dysregulation in blood malignancies. We also evaluate the potential of the dysregulated snoRNAs as biomarkers and therapeutic targets for blood malignancies.

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.

Chemical crosslinking and ligation methods for in vivo analysis of RNA structures and interactions

RNA structures and interactions in living cells drive a variety of biological processes and play critical roles in physiology and disease states. However, studies of RNA structures and interactions have been challenging due to limitations in available technologies. Direct determination of structures in vitro has been only possible to a small number of RNAs with limited sizes and conformations. We recently introduced two chemical crosslink-ligation techniques that enabled studies of transcriptome-wide secondary and tertiary structures and their dynamics. In a dramatically improved version of the psoralen analysis of RNA interactions and structures (PARIS2) method, we detailed the synthesis and use of amotosalen, a highly soluble psoralen analogue, and enhanced enzymology for higher efficiency duplex capture. We also introduced spatial 2'-hydroxyl acylation reversible crosslinking (SHARC) with exonuclease (exo) trimming, a method which utilizes a novel crosslinker class that targets the 2'-OH to capture three-dimensional (3D) structures. Both are powerful orthogonal approaches for solving in vivo RNA structure and interactions, integrating crosslinking, exo trimming, proximity ligation, and high throughput sequencing. In this chapter, we present a detailed protocol for the methods and highlight steps that outperform existing crosslink-ligation approaches.

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

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

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.

Intragenomic rearrangements involving 5'-untranslated region segments in SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses

BACKGROUND

Variation of the betacoronavirus SARS-CoV-2 has been the bane of COVID-19 control. Documented variation includes point mutations, deletions, insertions, and recombination among closely or distantly related coronaviruses. Here, we describe yet another aspect of genome variation by beta- and alphacoronaviruses that was first documented in an infectious isolate of the betacoronavirus SARS-CoV-2, obtained from 3 patients in Hong Kong that had a 5'-untranslated region segment at the end of the ORF6 gene that in its new location translated into an ORF6 protein with a predicted modified carboxyl terminus. While comparing the amino acid sequences of translated ORF8 genes in the GenBank database, we found a subsegment of the same 5'-UTR-derived amino acid sequence modifying the distal end of ORF8 of an isolate from the United States and decided to carry out a systematic search.

METHODS

Using the nucleotide and in the case of SARS-CoV-2 also the translated amino acid sequence in three reading frames of the genomic termini of coronaviruses as query sequences, we searched for 5'-UTR sequences in regions other than the 5'-UTR in SARS-CoV-2 and reference strains of alpha-, beta-, gamma-, and delta-coronaviruses.

RESULTS

We here report numerous genomic insertions of 5'-untranslated region sequences into coding regions of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses, but not delta- or gammacoronaviruses. To our knowledge this is the first systematic description of such insertions. In many cases, these insertions would change viral protein sequences and further foster genomic flexibility and viral adaptability through insertion of transcription regulatory sequences in novel positions within the genome. Among human Embecorivus betacoronaviruses, for instance, from 65% to all of the surveyed sequences in publicly available databases contain inserted 5'-UTR sequences.

CONCLUSION

The intragenomic rearrangements involving 5'-untranslated region sequences described here, which in several cases affect highly conserved genes with a low propensity for recombination, may underlie the generation of variants homotypic with those of concern or interest and with potentially differing pathogenic profiles. Intragenomic rearrangements thus add to our appreciation of how variants of SARS-CoV-2 and other beta- and alphacoronaviruses may arise.

RNAcontacts: A Pipeline for Predicting Contacts from RNA Proximity Ligation Assays

High-throughput RNA proximity ligation assays are molecular methods that are used to simultaneously analyze the spatial proximity of many RNAs in living cells. Their principle is based on cross-linking, fragmentation, and subsequent religation of RNAs, followed by high-throughput sequencing. The generated fragments have two different types of splits, one resulting from pre-mRNA splicing and the other formed by the ligation of spatially close RNA strands. Here, we present RNAcontacts, a universal pipeline for detecting RNA-RNA contacts in high-throughput RNA proximity ligation assays. RNAcontacts circumvents the inherent problem of mapping sequences with two distinct types of splits using a two-pass alignment, in which splice junctions are inferred from a control RNA-seq experiment on the first pass and then provided to the aligner as bona fide introns on the second pass. Compared to previously developed methods, our approach allows for a more sensitive detection of RNA contacts and has a higher specificity with respect to splice junctions that are present in the biological sample. RNAcontacts automatically extracts contacts, clusters their ligation points, computes the read support, and generates tracks for visualizing through the UCSC Genome Browser. The pipeline is implemented in Snakemake, a reproducible and scalable workflow management system for rapid and uniform processing of multiple datasets. RNAcontacts is a generic pipeline for the detection of RNA contacts that can be used with any proximity ligation method as long as one of the interacting partners is RNA. RNAcontacts is available via the GitHub repository https://github.com/smargasyuk/ RNAcontacts/.

Expression of miRNA-Targeted and Not-Targeted Reporter Genes Shows Mutual Influence and Intercellular Specificity

The regulation of translation by RNA-induced silencing complexes (RISCs) composed of Argonaute proteins and micro-RNAs is well established; however, the mechanisms underlying specific cellular responses to miRNAs and how specific complexes arise are not completely clear. To explore these questions, we performed experiments with Renilla and firefly luciferase reporter genes transfected in a psiCHECK-2 plasmid into human HCT116 or Me45 cells, where only the Renilla gene contained sequences targeted by microRNAs (miRNAs) in the 3'UTR. The effects of targeting were miRNA-specific; miRNA-21-5p caused strong inhibition of translation, whereas miRNA-24-3p or Let-7 family caused no change or an increase in reporter Renilla luciferase synthesis. The mRNA-protein complexes formed by transcripts regulated by different miRNAs differed from each other and were different in different cell types, as shown by sucrose gradient centrifugation. Unexpectedly, the presence of miRNA targets on Renilla transcripts also affected the expression of the co-transfected but non-targeted firefly luciferase gene in both cell types. Renilla and firefly transcripts were found in the same sucrose gradient fractions and specific anti-miRNA oligoribonucleotides, which influenced the expression of the Renilla gene, and also influenced that of firefly gene. These results suggest that, in addition to targeted transcripts, miRNAs may also modulate the expression of non-targeted transcripts, and using the latter to normalize the results may cause bias. We discuss some hypothetical mechanisms which could explain the observed miRNA-induced effects.

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.