COMRADES determines in vivo RNA structures and interactions

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.

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.

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.

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.

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.

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.

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.

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.

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/.

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.

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.

Generation and functional analysis of defective viral genomes during SARS-CoV-2 infection.. [PMID: 36172120 PMCID: PMC9516852 DOI: 10.1101/2022.09.22.509123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

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 RNA-seq datasets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hotspots 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 IFN stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the NGS dataset from a published cohort study and observed significantly higher DVG amount and frequency in symptomatic patients than that in asymptomatic patients. Finally, we observed unusually high DVG frequency in one immunosuppressive patient up to 140 days after admitted to hospital due to COVID-19, first-time suggesting 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 how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection.

Detecting RNA-RNA interactome

The last decade has seen a robust increase in various types of novel RNA molecules and their complexity in gene regulation. RNA molecules play a critical role in cellular events by interacting with other biomolecules, including protein, DNA, and RNA. It has been established that RNA-RNA interactions play a critical role in several biological processes by regulating the biogenesis and function of RNA molecules. Interestingly, RNA-RNA interactions regulate the biogenesis of diverse RNA molecules, including mRNAs, microRNAs, tRNAs, and circRNAs, through splicing or backsplicing. Structured RNAs like rRNA, tRNA, and snRNAs achieve their functional conformation by intramolecular RNA-RNA interactions. In addition, functional consequences of many intermolecular RNA-RNA interactions have been extensively studied in the regulation of gene expression. Hence, it is essential to understand the mechanism and functions of RNA-RNA interactions in eukaryotes. Conventionally, RNA-RNA interactions have been identified through diverse biochemical methods for decades. The advent of high-throughput RNA-sequencing technologies has revolutionized the identification of global RNA-RNA interactome in cells and their importance in RNA structure and function in gene expression regulation. Although these technologies revealed tens of thousands of intramolecular and intermolecular RNA-RNA interactions, we further look forward to future unbiased and quantitative high-throughput technologies for detecting transcriptome-wide RNA-RNA interactions. With the ability to detect RNA-RNA interactome, we expect that future studies will reveal the higher-order structures of RNA molecules and multi-RNA hybrids impacting human health and diseases. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

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.

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

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

Discovering functional motifs in long noncoding RNAs

Long noncoding RNAs (lncRNAs) are products of pervasive transcription that closely resemble messenger RNAs on the molecular level, yet function through largely unknown modes of action. The current model is that the function of lncRNAs often relies on specific, typically short, conserved elements, connected by linkers in which specific sequences and/or structures are less important. This notion has fueled the development of both computational and experimental methods focused on the discovery of functional elements within lncRNA genes, based on diverse signals such as evolutionary conservation, predicted structural elements, or the ability to rescue loss-of-function phenotypes. In this review, we outline the main challenges that the different methods need to overcome, describe the recently developed approaches, and discuss their respective limitations. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

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.

Computational approaches for RNA structure ensemble deconvolution from structure probing data

RNA structure probing experiments have emerged over the last decade as a straightforward way to determine the structure of RNA molecules in a number of different contexts. Although powerful, the ability of RNA to dynamically interconvert between, and to simultaneously populate, alternative structural configurations, poses a nontrivial challenge to the interpretation of data derived from these experiments. Recent efforts aimed at developing computational methods for the reconstruction of coexisting alternative RNA conformations from structure probing data are paving the way to the study of RNA structure ensembles, even in the context of living cells. In this review, we critically discuss these methods, their limitations and possible future improvements.

Classification and clustering of RNA crosslink-ligation data reveal complex structures and homodimers

The recent development and application of methods based on the general principle of "crosslinking and proximity ligation" (crosslink-ligation) are revolutionizing RNA structure studies in living cells. However, extracting structure information from such data presents unique challenges. Here, we introduce a set of computational tools for the systematic analysis of data from a wide variety of crosslink-ligation methods, specifically focusing on read mapping, alignment classification, and clustering. We design a new strategy to map short reads with irregular gaps at high sensitivity and specificity. Analysis of previously published data reveals distinct properties and bias caused by the crosslinking reactions. We perform rigorous and exhaustive classification of alignments and discover eight types of arrangements that provide distinct information on RNA structures and interactions. To deconvolve the dense and intertwined gapped alignments, we develop a network/graph-based tool Crosslinked RNA Secondary Structure Analysis using Network Techniques (CRSSANT), which enables clustering of gapped alignments and discovery of new alternative and dynamic conformations. We discover that multiple crosslinking and ligation events can occur on the same RNA, generating multisegment alignments to report complex high-level RNA structures and multi-RNA interactions. We find that alignments with overlapped segments are produced from potential homodimers and develop a new method for their de novo identification. Analysis of overlapping alignments revealed potential new homodimers in cellular noncoding RNAs and RNA virus genomes in the Picornaviridae family. Together, this suite of computational tools enables rapid and efficient analysis of RNA structure and interaction data in living cells.

Global mapping of RNA homodimers in living cells

RNA homodimerization is important for various physiological processes, including the assembly of membraneless organelles, RNA subcellular localization, and packaging of viral genomes. However, understanding RNA dimerization has been hampered by the lack of systematic in vivo detection methods. Here, we show that CLASH, PARIS, and other RNA proximity ligation methods detect RNA homodimers transcriptome-wide as "overlapping" chimeric reads that contain more than one copy of the same sequence. Analyzing published proximity ligation data sets, we show that RNA:RNA homodimers mediated by direct base-pairing are rare across the human transcriptome, but highly enriched in specific transcripts, including U8 snoRNA, U2 snRNA, and a subset of tRNAs. Mutations in the homodimerization domain of U8 snoRNA impede dimerization in vitro and disrupt zebrafish development in vivo, suggesting an evolutionarily conserved role of this domain. Analysis of virus-infected cells reveals homodimerization of SARS-CoV-2 and Zika genomes, mediated by specific palindromic sequences located within protein-coding regions of N gene in SARS-CoV-2 and NS2A gene in Zika. We speculate that regions of viral genomes involved in homodimerization may constitute effective targets for antiviral therapies.

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.

Methods to Analyze the Non-Coding RNA Interactome—Recent Advances and Challenges

Most of the human genome is transcribed to generate a multitude of non-coding RNAs. However, while these transcripts have generated an immense amount of scientific interest, their biological function remains a subject of an intense debate. Understanding mechanisms of action of non-coding RNAs is a key to addressing the issue of biological relevance of these transcripts. Based on some well-understood non-coding RNAs that function inside the cell by interacting with other molecules, it is generally believed many other non-coding transcripts could also function in a similar fashion. Therefore, development of methods that can map RNA interactome is the key to understanding functionality of the extensive cellular non-coding transcriptome. Here, we review the vast progress that has been made in the past decade in technologies that can map RNA interactions with different sites in DNA, proteins or other RNA molecules; the general approaches used to validate the existence of novel interactions; and the challenges posed by interpreting the data obtained using the interactome mapping methods.

Chemical reversible crosslinking enables measurement of RNA 3D distances and alternative conformations in cells

Three-dimensional (3D) structures dictate the functions of RNA molecules in a wide variety of biological processes. However, direct determination of RNA 3D structures in vivo is difficult due to their large sizes, conformational heterogeneity, and dynamics. Here we present a method, Spatial 2′-Hydroxyl Acylation Reversible Crosslinking (SHARC), which uses chemical crosslinkers of defined lengths to measure distances between nucleotides in cellular RNA. Integrating crosslinking, exonuclease (exo) trimming, proximity ligation, and high throughput sequencing, SHARC enables transcriptome-wide tertiary structure contact maps at high accuracy and precision, revealing heterogeneous RNA structures and interactions. SHARC data provide constraints that improves Rosetta-based RNA 3D structure modeling at near-nanometer resolution. Integrating SHARC-exo with other crosslinking-based methods, we discover compact folding of the 7SK RNA, a critical regulator of transcriptional elongation. These results establish a strategy for measuring RNA 3D distances and alternative conformations in their native cellular context.

Determination of RNA 3D structures in vivo is a challenging problem. Here, the authors describe a chemical crosslinking method (SHARC) that they use in combination with RNA sequencing to measure distances between nucleotides in RNA 3D structures.

Promising Assays for Examining a Putative Role of Ribosomal Heterogeneity in COVID-19 Susceptibility and Severity

The heterogeneity of ribosomes, characterized by structural variations, arises from differences in types, numbers, and/or post-translational modifications of participating ribosomal proteins (RPs), ribosomal RNAs (rRNAs) sequence variants plus post-transcriptional modifications, and additional molecules essential for forming a translational machinery. The ribosomal heterogeneity within an individual organism or a single cell leads to preferential translations of selected messenger RNA (mRNA) transcripts over others, especially in response to environmental cues. The role of ribosomal heterogeneity in SARS-CoV-2 coronavirus infection, propagation, related symptoms, or vaccine responses is not known, and a technique to examine these has not yet been developed. Tools to detect ribosomal heterogeneity or to profile translating mRNAs independently cannot identify unique or specialized ribosome(s) along with corresponding mRNA substrate(s). Concurrent characterizations of RPs and/or rRNAs with mRNA substrate from a single ribosome would be critical to decipher the putative role of ribosomal heterogeneity in the COVID-19 disease, caused by the SARS-CoV-2, which hijacks the host ribosome to preferentially translate its RNA genome. Such a protocol should be able to provide a high-throughput screening of clinical samples in a large population that would reach a statistical power for determining the impact of a specialized ribosome to specific characteristics of the disease. These characteristics may include host susceptibility, viral infectivity and transmissibility, severity of symptoms, antiviral treatment responses, and vaccine immunogenicity including its side effect and efficacy. In this study, several state-of-the-art techniques, in particular, chemical probing of ribosomal components or rRNA structures, proximity ligation to generate rRNA-mRNA chimeras for sequencing, nanopore gating of individual ribosomes, nanopore RNA sequencing and/or structural analyses, single-ribosome mass spectrometry, and microfluidic droplets for separating ribosomes or indexing rRNAs/mRNAs, are discussed. The key elements for further improvement and proper integration of the above techniques to potentially arrive at a high-throughput protocol for examining individual ribosomes and their mRNA substrates in a clinical setting are also presented.

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.