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

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.

Protocol for profiling virus-to-host RNA-RNA interactions in infected cells by RIC-seq

Virus-to-host RNA-RNA interactions directly regulate host mRNA stability and viral replication. However, globally profiling virus-to-host in situ RNA-RNA interactions remains challenging. Here, we present an RNA in situ conformation sequencing (RIC-seq)-based protocol for mapping high-confidence virus-to-host in situ RNA-RNA interactions in infected cells. We detail steps for formaldehyde crosslinking, pCp-biotin labeling, in situ proximity ligation, chimeric RNA enrichment, strand-specific library construction, and data analysis. This protocol allows unbiased identification of virus-to-host RNA-RNA interactions for various RNA viruses and is potentially applicable to DNA virus-derived transcripts. For complete details on the use and execution of this protocol, please refer to Zhao et al.1.

The selective sponging of miRNAs by OIP5-AS1 regulates metabolic reprogramming of pyruvate in adenoma-carcinoma transition of human colorectal cancer

RNA interactomes and their diversified functionalities have recently benefited from critical methodological advances leading to a paradigm shift from a conventional conception on the regulatory roles of RNA in pathogenesis. However, the dynamic RNA interactomes in adenoma-carcinoma sequence of human CRC remain unexplored. The coexistence of adenoma, cancer, and normal tissues in colorectal cancer (CRC) patients provides an appropriate model to address this issue. Here, we adopted an RNA in situ conformation sequencing technology for mapping RNA-RNA interactions in CRC patients. We observed large-scale paired RNA counts and identified some unique RNA complexes including multiple partners RNAs, single partner RNAs, non-overlapping single partner RNAs. We focused on the antisense RNA OIP5-AS1 and found that OIP5-AS1 could sponge different miRNA to regulate the production of metabolites including pyruvate, alanine and lactic acid. Our findings provide novel perspectives in CRC pathogenesis and suggest metabolic reprogramming of pyruvate for the early diagnosis and treatment of CRC.

Global identification of mRNA-interacting circular RNAs by CLiPPR-Seq

Although the functional role of circular RNA (circRNA) interaction with microRNAs and proteins has been studied extensively, circRNA interactions with the protein-coding mRNAs in intact cells remain largely unknown. Here, by employing AMT-mediated proximity ligation of RNA-RNA duplexes followed by circRNA enrichment and deep sequencing, we report a novel Cross-Linking Poly(A) Pulldown RNase R Sequencing (CLiPPR-seq) technology which identified hundreds of mRNA-interacting circRNAs in three different cell types, including βTC6, C2C12 and HeLa cells. Furthermore, CLiPP-seq without RNase R treatment was also performed to identify the mRNA expression in these cells. BLAST analysis of circRNAs in CLiPPR-seq sample with the mRNAs in CLiPP-seq samples determined their potential complementary sequences for circRNA-mRNA interaction. Pulldown of circRNAs and poly(A) RNAs confirmed the direct interaction of circRNAs with target mRNAs. Silencing of mRNA-interacting circRNAs led to the altered expression of target mRNAs in βTC6 cells, suggesting the role of direct interaction of circRNAs with mRNAs in gene expression regulation. CLiPPR-seq thus represents a novel method for illuminating the myriad of uncharacterized circRNA-mRNA hybrids that may regulate gene expression.

SARS-CoV-2 RNA stabilizes host mRNAs to elicit immunopathogenesis

SARS-CoV-2 RNA interacts with host factors to suppress interferon responses and simultaneously induces cytokine release to drive the development of severe coronavirus disease 2019 (COVID-19). However, how SARS-CoV-2 hijacks host RNAs to elicit such imbalanced immune responses remains elusive. Here, we analyzed SARS-CoV-2 RNA in situ structures and interactions in infected cells and patient lung samples using RIC-seq. We discovered that SARS-CoV-2 RNA forms 2,095 potential duplexes with the 3' UTRs of 205 host mRNAs to increase their stability by recruiting RNA-binding protein YBX3 in A549 cells. Disrupting the SARS-CoV-2-to-host RNA duplex or knocking down YBX3 decreased host mRNA stability and reduced viral replication. Among SARS-CoV-2-stabilized host targets, NFKBIZ was crucial for promoting cytokine production and reducing interferon responses, probably contributing to cytokine storm induction. Our study uncovers the crucial roles of RNA-RNA interactions in the immunopathogenesis of RNA viruses such as SARS-CoV-2 and provides valuable host targets for drug development.

Comprehensive landscape and future perspective of long noncoding RNAs in non-small cell lung cancer: it takes a village

Long noncoding RNAs (lncRNAs) are a distinct subtype of RNA that lack protein-coding capacity but exert significant influence on various cellular processes. In non-small cell lung cancer (NSCLC), dysregulated lncRNAs act as either oncogenes or tumor suppressors, contributing to tumorigenesis and tumor progression. LncRNAs directly modulate gene expression, act as competitive endogenous RNAs by interacting with microRNAs or proteins, and associate with RNA binding proteins. Moreover, lncRNAs can reshape the tumor immune microenvironment and influence cellular metabolism, cancer cell stemness, and angiogenesis by engaging various signaling pathways. Notably, lncRNAs have shown great potential as diagnostic or prognostic biomarkers in liquid biopsies and therapeutic strategies for NSCLC. This comprehensive review elucidates the significant roles and diverse mechanisms of lncRNAs in NSCLC. Furthermore, we provide insights into the clinical relevance, current research progress, limitations, innovative research approaches, and future perspectives for targeting lncRNAs in NSCLC. By summarizing the existing knowledge and advancements, we aim to enhance the understanding of the pivotal roles played by lncRNAs in NSCLC and stimulate further research in this field. Ultimately, unraveling the complex network of lncRNA-mediated regulatory mechanisms in NSCLC could potentially lead to the development of novel diagnostic tools and therapeutic strategies.

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.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

lhCLIP reveals the in vivo RNA-RNA interactions recognized by hnRNPK

RNA-RNA interactions play a crucial role in regulating gene expression and various biological processes, but identifying these interactions on a transcriptomic scale remains a challenge. To address this, we have developed a new biochemical technique called pCp-biotin labelled RNA hybrid and ultraviolet crosslinking and immunoprecipitation (lhCLIP) that enables the transcriptome-wide identification of intra- and intermolecular RNA-RNA interactions mediated by a specific RNA-binding protein (RBP). Using lhCLIP, we have uncovered a diverse landscape of intermolecular RNA interactions recognized by hnRNPK in human cells, involving all major classes of noncoding RNAs (ncRNAs) and mRNA. Notably, hnRNPK selectively binds with snRNA U4, U11, and U12, and shapes the secondary structure of these snRNAs, which may impact RNA splicing. Our study demonstrates the potential of lhCLIP as a user-friendly and widely applicable method for discovering RNA-RNA interactions mediated by a particular protein of interest and provides a valuable tool for further investigating the role of RBPs in gene expression and biological processes.

CRIC-seq protocol for in situ profiling of proximal RNA-RNA contacts associated with RNA-binding proteins

RNA-binding proteins (RBPs) can bind and mediate RNA-RNA contacts. However, identifying specific RBP-organized RNA-RNA contacts remains challenging. Here, we present a capture RIC-seq (CRIC-seq) technique to map specific RBP-associated RNA-RNA contacts globally. We describe steps for formaldehyde cross-linking to fix RNA in situ conformation, pCp-biotin labeling to mark RNA juncture, and in situ proximity ligation to join proximal RNAs. We then detail immunoprecipitation to isolate specific RBP-associated RNA-RNA contacts, biotin-streptavidin selection to enrich chimeric RNAs, and library construction for paired-end sequencing. For complete information on the generation and use of this protocol, please refer to Ye et al.1.

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.

Complementary Alu sequences mediate enhancer-promoter selectivity

Enhancers determine spatiotemporal gene expression programs by engaging with long-range promoters^1-4. However, it remains unknown how enhancers find their cognate promoters. We recently developed a RNA in situ conformation sequencing technology to identify enhancer-promoter connectivity using pairwise interacting enhancer RNAs and promoter-derived noncoding RNAs^5,6. Here we apply this technology to generate high-confidence enhancer-promoter RNA interaction maps in six additional cell lines. Using these maps, we discover that 37.9% of the enhancer-promoter RNA interaction sites are overlapped with Alu sequences. These pairwise interacting Alu and non-Alu RNA sequences tend to be complementary and potentially form duplexes. Knockout of Alu elements compromises enhancer-promoter looping, whereas Alu insertion or CRISPR-dCasRx-mediated Alu tethering to unregulated promoter RNAs can create new loops to homologous enhancers. Mapping 535,404 noncoding risk variants back to the enhancer-promoter RNA interaction maps enabled us to construct variant-to-function maps for interpreting their molecular functions, including 15,318 deletions or insertions in 11,677 Alu elements that affect 6,497 protein-coding genes. We further demonstrate that polymorphic Alu insertion at the PTK2 enhancer can promote tumorigenesis. Our study uncovers a principle for determining enhancer-promoter pairing specificity and provides a framework to link noncoding risk variants to their molecular functions.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

Capture RIC-seq reveals positional rules of PTBP1-associated RNA loops in splicing regulation

RNA-binding proteins (RBPs) bind at different positions of the pre-mRNA molecules to promote or reduce the usage of a particular exon. Seeking to understand the working principle of these positional effects, we develop a capture RIC-seq (CRIC-seq) method to enrich specific RBP-associated in situ proximal RNA-RNA fragments for deep sequencing. We determine hnRNPA1-, SRSF1-, and PTBP1-associated proximal RNA-RNA contacts and regulatory mechanisms in HeLa cells. Unexpectedly, the 3D RNA map analysis shows that PTBP1-associated loops in individual introns preferentially promote cassette exon splicing by accelerating asymmetric intron removal, whereas the loops spanning across cassette exon primarily repress splicing. These "positional rules" can faithfully predict PTBP1-regulated splicing outcomes. We further demonstrate that cancer-related splicing quantitative trait loci can disrupt RNA loops by reducing PTBP1 binding on pre-mRNAs to cause aberrant splicing in tumors. Our study presents a powerful method for exploring the functions of RBP-associated RNA-RNA proximal contacts in gene regulation and disease.

Approaches for Modes of Action Study of Long Non-Coding RNAs: From Single Verification to Genome-Wide Determination

Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides (nt) that are not translated into known functional proteins. This broad definition covers a large collection of transcripts with diverse genomic origins, biogenesis, and modes of action. Thus, it is very important to choose appropriate research methodologies when investigating lncRNAs with biological significance. Multiple reviews to date have summarized the mechanisms of lncRNA biogenesis, their localization, their functions in gene regulation at multiple levels, and also their potential applications. However, little has been reviewed on the leading strategies for lncRNA research. Here, we generalize a basic and systemic mind map for lncRNA research and discuss the mechanisms and the application scenarios of ‘up-to-date’ techniques as applied to molecular function studies of lncRNAs. Taking advantage of documented lncRNA research paradigms as examples, we aim to provide an overview of the developing techniques for elucidating lncRNA interactions with genomic DNA, proteins, and other RNAs. In the end, we propose the future direction and potential technological challenges of lncRNA studies, focusing on techniques and applications.

From molecular basics to agronomic benefits: Insights into noncoding RNA-mediated gene regulation in plants

The development of plants is largely dependent on their growth environment. To better adapt to a particular habitat, plants have evolved various subtle regulatory mechanisms for altering gene expression. Non coding RNAs (ncRNAs) constitute a major portion of the transcriptomes of eukaryotes. Various ncRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycles of plants. In this review, we summarize the current understanding of the biogenesis and contributions of small nucle olar RNA (snoRNA)- and regulatory long non coding RNA (lncRNA)-mediated gene regulation in plant development and environmental responses. Many regulatory ncRNAs appear to be associated with increased yield, quality and disease resistance of various species and cultivars. These ncRNAs may potentially be used as genetic resources for improving agronomic traits and for molecular breeding. The challenges in understanding plant ncRNA biology and the possibilities to make better use of these valuable gene resources in the future are discussed in this review.

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.

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.

Current Status of Epitranscriptomic Marks Affecting lncRNA Structures and Functions

Long non-coding RNAs (lncRNAs) belong to a class of non-protein-coding RNAs with their lengths longer than 200 nucleotides. Most of the mammalian genome is transcribed as RNA, yet only a small percent of the transcribed RNA corresponds to exons of protein-coding genes. Thus, the number of lncRNAs is predicted to be several times higher than that of protein-coding genes. Because of sheer number of lncRNAs, it is often difficult to elucidate the functions of all lncRNAs, especially those arising from their relationship to their binding partners, such as DNA, RNA, and proteins. Due to their binding to other macromolecules, it has become evident that the structures of lncRNAs influence their functions. In this regard, the recent development of epitranscriptomics (the field of study to investigate RNA modifications) has become important to further elucidate the structures and functions of lncRNAs. In this review, the current status of lncRNA structures and functions influenced by epitranscriptomic marks is discussed.

CLASH Analyst: A Web Server to Identify In Vivo RNA–RNA Interactions from CLASH Data

Non-coding RNAs, such as miRNAs and piRNAs, play critical roles in gene regulation through base-pairing interactions with their target molecules. The recent development of the crosslinking, ligation, and sequencing of hybrids (CLASH) method has allowed scientists to map transcriptome-wide RNA–RNA interactions by identifying chimeric reads consisting of fragments from regulatory RNAs and their targets. However, analyzing CLASH data requires scientists to use advanced bioinformatics, and currently available tools are limited for users with little bioinformatic experience. In addition, many published CLASH studies do not show the full scope of RNA–RNA interactions that were captured, highlighting the importance of reanalyzing published data. Here, we present CLASH Analyst, a web server that can analyze raw CLASH data within a fully customizable and easy-to-use interface. CLASH Analyst accepts raw CLASH data as input and identifies the RNA chimeras containing the regulatory and target RNAs according to the user’s interest. Detailed annotation of the captured RNA–RNA interactions is then presented for the user to visualize within the server or download for further analysis. We demonstrate that CLASH Analyst can identify miRNA- and piRNA-targeting sites reported from published CLASH data and should be applicable to analyze other RNA–RNA interactions. CLASH Analyst is freely available for academic use.

Architecture of RNA-RNA interactions

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