Systematic Mapping of RNA-Chromatin Interactions In Vivo

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

Telomere-specific chromatin capture using a pyrrole-imidazole polyamide probe for the identification of proteins and non-coding RNAs

Background

Knowing chromatin components at a DNA regulatory element at any given time is essential for understanding how the element works during cellular proliferation, differentiation and development. A region-specific chromatin purification is an invaluable approach to dissecting the comprehensive chromatin composition at a particular region. Several methods (e.g., PICh, enChIP, CAPTURE and CLASP) have been developed for isolating and analyzing chromatin components. However, all of them have some shortcomings in identifying non-coding RNA associated with DNA regulatory elements.

Results

We have developed a new approach for affinity purification of specific chromatin segments employing an N-methyl pyrrole (P)-N-methylimidazole (I) (PI) polyamide probe, which binds to a specific sequence in double-stranded DNA via Watson–Crick base pairing as a minor groove binder. This new technique is called proteomics and RNA-omics of isolated chromatin segments (PI-PRICh). Using PI-PRICh to isolate mouse and human telomeric components, we found enrichments of shelterin proteins, the well-known telomerase RNA component (TERC) and telomeric repeat-containing RNA (TERRA). When PI-PRICh was performed for alternative lengthening of telomere (ALT) cells with highly recombinogenic telomeres, in addition to the conventional telomeric chromatin, we obtained chromatin regions containing telomeric repeat insertions scattered in the genome and their associated RNAs.

Conclusion

PI-PRICh reproducibly identified both the protein and RNA components of telomeric chromatin when targeting telomere repeats. PI polyamide is a promising alternative to simultaneously isolate associated proteins and RNAs of sequence-specific chromatin regions under native conditions, allowing better understanding of chromatin organization and functions within the cell.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00421-8.

Global profiling of RNA-chromatin interactions reveals co-regulatory gene expression networks in Arabidopsis

It is increasingly evident that various RNAs can bind chromatin to regulate gene expression and genome organization. Here we adapted a sequencing-based technique to profile RNA-chromatin interactions at a genome-wide scale in Arabidopsis seedlings. We identified more than 10,000 RNA-chromatin interactions mediated by protein-coding RNAs and non-coding RNAs. Cis and intra-chromosomal interactions are mainly mediated by protein-coding RNAs, whereas inter-chromosomal interactions are primarily mediated by non-coding RNAs. Many RNA-chromatin interactions tend to positively correlate with DNA-DNA interactions, suggesting their mutual influence and reinforcement. We further show that some RNA-chromatin interactions undergo alterations in response to biotic and abiotic stresses and that altered RNA-chromatin interactions form co-regulatory networks. Our study provides a global view on RNA-chromatin interactions in Arabidopsis and a rich resource for future investigations of regulatory roles of RNAs in gene expression and genome organization.

Genome-wide regulation of CpG methylation by ecCEBPα in acute myeloid leukemia

Background: Acute myeloid leukemia (AML) is a hematopoietic malignancy characterized by genetic and epigenetic aberrations that alter the differentiation capacity of myeloid progenitor cells. The transcription factor CEBPα is frequently mutated in AML patients leading to an increase in DNA methylation in many genomic locations. Previously, it has been shown that ecCEBPα (extra coding CEBP α) - a lncRNA transcribed in the same direction as CEBPα gene - regulates DNA methylation of CEBPα promoter in cis. Here, we hypothesize that ecCEBPα could participate in the regulation of DNA methylation in trans.

Method: First, we retrieved the methylation profile of AML patients with mutated CEBPα locus from The Cancer Genome Atlas (TCGA). We then predicted the ecCEBPα secondary structure in order to check the potential of ecCEBPα to form triplexes around CpG loci and checked if triplex formation influenced CpG methylation, genome-wide.

Results: Using DNA methylation profiles of AML patients with a mutated CEBPα locus, we show that ecCEBPα could interact with DNA by forming DNA:RNA triple helices and protect regions near its binding sites from global DNA methylation. Further analysis revealed that triplex-forming oligonucleotides in ecCEBPα are structurally unpaired supporting the DNA-binding potential of these regions. ecCEBPα triplexes supported with the RNA-chromatin co-localization data are located in the promoters of leukemia-linked transcriptional factors such as MLF2.

Discussion: Overall, these results suggest a novel regulatory mechanism for ecCEBPα as a genome-wide epigenetic modulator through triple-helix formation which may provide a foundation for sequence-specific engineering of RNA for regulating methylation of specific genes.

Clinicopathological value of long non-coding RNA profiles in gastrointestinal stromal tumor

Background

Long non-coding RNAs (lncRNAs) have been implicated in diagnosis and prognosis in various cancers. However, few lncRNA signatures have been established for prediction of gastrointestinal stromal tumors (GIST). We aimed to explore a lncRNA signature profile that associated with clinical relevance by mining data from Gene Expression Ominus (GEO) and Surveillance, Epidemiology, and End Results (SEER) Program.

Methods

Using a lncRNA-mining approach, we performed non-negative matrix factorization (NMF) consensus algorithm in Gastrointestinal stromal tumors (GISTs) cohorts (61 patients from GSE8167 and GSE17743) to cluster LncRNA expression profiles. Comparative markers selection, and Gene Set Enrichment Analysis (GSEA) algorithm were performed between distinct molecular subtypes of GIST. The survival rate of GIST patients from SEER stratified by gender were compared by Kaplan-Meier method and log-rank analysis. lncRNA-mRNA co-expression analysis was performed by Pearson correlation coefficients (PCC) using R package LINC. Somatic copy number alterations of GIST patients (GSE40966) were analyzed via web server GenePattern GISTIC2 algorithm.

Results

A total of four lncRNA molecular subtypes of GIST were identified with distinct biological pathways and clinical characteristics. LncRNA expression profiles well clustered the GIST samples into small size (<5 mm) and large size tumors (>5 mm), which is a fundamental index for GIST malignancy diagnosis. Several lncRNAs with abundant expression (LRRC75A-AS1, HYMAI, NEAT1, XIST and FTX) were closely associated with tumor size, which may suggest to be biomarkers for the GIST malignancy. Particularly, LRRC75A-AS1 was positively associated with tumor diameters and suggested an oncogene in GIST. Co-expression analysis suggested that chromosome region 17p11.2-p12 may contribute to the oncogenic process in malignant GIST. Interestingly, the gender had a strong influence on clustering by lncRNA expression profile. Data from the Surveillance, Epidemiology, and End Results (SEER) Program were further explored and 7983 patients who were diagnosed with GISTs from 1973 to 2014 were enrolled for analysis. The results also showed the favorable prognosis for female patients. The survival rate between male and female with GIST was statistically significant (P < 0.0001). Gene set enrichment analysis (GSEA) indicated distinct pathways between female and male, and malignant GIST was associated with several cancer metabolism and cell cycle associated pathways.

Conclusions

This lncRNAs-based classification for GISTs may provide a molecular classification applicable to individual GIST that has implications to influence lncRNA markers selection and prediction of tumor progression.

Investigation of the Basic Steps in the Chromosome Conformation Capture Procedure

A constellation of chromosome conformation capture methods (С-methods) are an important tool for biochemical analysis of the spatial interactions between DNA regions that are separated in the primary sequence. All these methods are based on the long sequence of basic steps of treating cells, nuclei, chromatin, and finally DNA, thus representing a significant technical challenge. Here, we present an in-depth study of the basic steps in the chromatin conformation capture procedure (3С), which was performed using Drosophila Schneider 2 cells as a model. We investigated the steps of cell lysis, nuclei washing, nucleoplasm extraction, chromatin treatment with SDS/Triton X-100, restriction enzyme digestion, chromatin ligation, reversion of cross-links, DNA extraction, treatment of a 3C library with RNases, and purification of the 3C library. Several options were studied, and optimal conditions were found. Our work contributes to the understanding of the 3C basic steps and provides a useful guide to the 3C procedure.

Epigenetic Regulation of the Vascular Endothelium by Angiogenic LncRNAs

The functional properties of the vascular endothelium are diverse and heterogeneous between vascular beds. This is especially evident when new blood vessels develop from a pre-existing closed cardiovascular system, a process termed angiogenesis. Endothelial cells are key drivers of angiogenesis as they undergo a highly choreographed cascade of events that has both exogenous (e.g., hypoxia and VEGF) and endogenous regulatory inputs. Not surprisingly, angiogenesis is critical in health and disease. Diverse therapeutics target proteins involved in coordinating angiogenesis with varying degrees of efficacy. It is of great interest that recent work on non-coding RNAs, especially long non-coding RNAs (lncRNAs), indicates that they are also important regulators of the gene expression paradigms that underpin this cellular cascade. The protean effects of lncRNAs are dependent, in part, on their subcellular localization. For instance, lncRNAs enriched in the nucleus can act as epigenetic modifiers of gene expression in the vascular endothelium. Of great interest to genetic disease, they are undergoing rapid evolution and show extensive inter- and intra-species heterogeneity. In this review, we describe endothelial-enriched lncRNAs that have robust effects in angiogenesis.

Non-Coding Variants in Cancer: Mechanistic Insights and Clinical Potential for Personalized Medicine

The cancer genome is characterized by extensive variability, in the form of Single Nucleotide Polymorphisms (SNPs) or structural variations such as Copy Number Alterations (CNAs) across wider genomic areas. At the molecular level, most SNPs and/or CNAs reside in non-coding sequences, ultimately affecting the regulation of oncogenes and/or tumor-suppressors in a cancer-specific manner. Notably, inherited non-coding variants can predispose for cancer decades prior to disease onset. Furthermore, accumulation of additional non-coding driver mutations during progression of the disease, gives rise to genomic instability, acting as the driving force of neoplastic development and malignant evolution. Therefore, detection and characterization of such mutations can improve risk assessment for healthy carriers and expand the diagnostic and therapeutic toolbox for the patient. This review focuses on functional variants that reside in transcribed or not transcribed non-coding regions of the cancer genome and presents a collection of appropriate state-of-the-art methodologies to study them.

Non-coding RNAs in chromatin folding and nuclear organization

One of the most intriguing questions facing modern biology concerns how the genome directs the construction of cells, tissues, and whole organisms. It is tempting to suggest that the part of the genome that does not encode proteins contains architectural plans. We are still far from understanding how these plans work at the level of building tissues and the body as a whole. However, the results of recent studies demonstrate that at the cellular level, special non-coding RNAs serve as scaffolds for the construction of various intracellular structures. The term "architectural RNAs" was proposed to designate this subset of non-coding RNAs. In this review, we discuss the role of architectural RNAs in the construction of the cell nucleus and maintenance of the three-dimensional organization of the genome.

Decoding LncRNAs

Non-coding RNAs (ncRNAs) have been considered as unimportant additions to the transcriptome. Yet, in light of numerous studies, it has become clear that ncRNAs play important roles in development, health and disease. Long-ignored, long non-coding RNAs (lncRNAs), ncRNAs made of more than 200 nucleotides have gained attention due to their involvement as drivers or suppressors of a myriad of tumours. The detailed understanding of some of their functions, structures and interactomes has been the result of interdisciplinary efforts, as in many cases, new methods need to be created or adapted to characterise these molecules. Unlike most reviews on lncRNAs, we summarize the achievements on lncRNA studies by taking into consideration the approaches for identification of lncRNA functions, interactomes, and structural arrangements. We also provide information about the recent data on the involvement of lncRNAs in diseases and present applications of these molecules, especially in medicine.

Very long intergenic non-coding (vlinc) RNAs directly regulate multiple genes in cis and trans

BACKGROUND

The majority of the human genome is transcribed in the form of long non-coding (lnc) RNAs. While these transcripts have attracted considerable interest, their molecular mechanisms of function and biological significance remain controversial. One of the main reasons behind this lies in the significant challenges posed by lncRNAs requiring the development of novel methods and concepts to unravel their functionality. Existing methods often lack cross-validation and independent confirmation by different methodologies and therefore leave significant ambiguity as to the authenticity of the outcomes. Nonetheless, despite all the caveats, it appears that lncRNAs may function, at least in part, by regulating other genes via chromatin interactions. Therefore, the function of a lncRNA could be inferred from the function of genes it regulates. In this work, we present a genome-wide functional annotation strategy for lncRNAs based on identification of their regulatory networks via the integration of three distinct types of approaches: co-expression analysis, mapping of lncRNA-chromatin interactions, and assaying molecular effects of lncRNA knockdowns obtained using an inducible and highly specific CRISPR/Cas13 system.

RESULTS

We applied the strategy to annotate 407 very long intergenic non-coding (vlinc) RNAs belonging to a novel widespread subclass of lncRNAs. We show that vlincRNAs indeed appear to regulate multiple genes encoding proteins predominantly involved in RNA- and development-related functions, cell cycle, and cellular adhesion via a mechanism involving proximity between vlincRNAs and their targets in the nucleus. A typical vlincRNAs can be both a positive and negative regulator and regulate multiple genes both in trans and cis. Finally, we show vlincRNAs and their regulatory networks potentially represent novel components of DNA damage response and are functionally important for the ability of cancer cells to survive genotoxic stress.

CONCLUSIONS

This study provides strong evidence for the regulatory role of the vlincRNA class of lncRNAs and a potentially important role played by these transcripts in the hidden layer of RNA-based regulation in complex biological systems.

The Toolbox for Untangling Chromosome Architecture in Immune Cells

The code of life is not only encrypted in the sequence of DNA but also in the way it is organized into chromosomes. Chromosome architecture is gradually being recognized as an important player in regulating cell activities (e.g., controlling spatiotemporal gene expression). In the past decade, the toolbox for elucidating genome structure has been expanding, providing an opportunity to explore this under charted territory. In this review, we will introduce the recent advancements in approaches for mapping spatial organization of the genome, emphasizing applications of these techniques to immune cells, and trying to bridge chromosome structure with immune cell activities.

Battle for supremacy: nucleic acid interactions between viruses and cells

Since the COVID-19 pandemic swept across the globe, researchers have been trying to understand its origin, life cycle, and pathogenesis. There is a striking variability in the phenotypic response to infection with SARS-CoV-2 that may reflect differences in host genetics and/or immune response. It is known that the human epigenome is influenced by ethnicity, age, lifestyle, and environmental factors, including previous viral infections. This Review examines the influence of viruses on the host epigenome. We describe general lessons and methodologies that can be used to understand how the virus evades the host immune response. We consider how variation in the epigenome may contribute to heterogeneity in the response to SARS-CoV-2 and may identify a precision medicine approach to treatment.

Evolution of Genome-Organizing Long Non-coding RNAs in Metazoans

Long non-coding RNAs (lncRNAs) have important regulatory functions across eukarya. It is now clear that many of these functions are related to gene expression regulation through their capacity to recruit epigenetic modifiers and establish chromatin interactions. Several lncRNAs have been recently shown to participate in modulating chromatin within the spatial organization of the genome in the three-dimensional space of the nucleus. The identification of lncRNA candidates is challenging, as it is their functional characterization. Conservation signatures of lncRNAs are different from those of protein-coding genes, making identifying lncRNAs under selection a difficult task, and the homology between lncRNAs may not be readily apparent. Here, we review the evidence for these higher-order genome organization functions of lncRNAs in animals and the evolutionary signatures they display.

The lncRNA Toolkit: Databases and In Silico Tools for lncRNA Analysis

Long non-coding RNAs (lncRNAs) are a rapidly expanding field of research, with many new transcripts identified each year. However, only a small subset of lncRNAs has been characterized functionally thus far. To aid investigating the mechanisms of action by which new lncRNAs act, bioinformatic tools and databases are invaluable. Here, we review a selection of computational tools and databases for the in silico analysis of lncRNAs, including tissue-specific expression, protein coding potential, subcellular localization, structural conformation, and interaction partners. The assembled lncRNA toolkit is aimed primarily at experimental researchers as a useful starting point to guide wet-lab experiments, mainly containing multi-functional, user-friendly interfaces. With more and more new lncRNA analysis tools available, it will be essential to provide continuous updates and maintain the availability of key software in the future.

Direct Interactions with Nascent Transcripts Is Potentially a Common Targeting Mechanism of Long Non-Coding RNAs

Although thousands of mammalian long non-coding RNAs (lncRNAs) have been reported in the last decade, their functional annotation remains limited. A wet-lab approach to detect functions of a novel lncRNA usually includes its knockdown followed by RNA sequencing and identification of the deferentially expressed genes. However, identification of the molecular mechanism(s) used by the lncRNA to regulate its targets frequently becomes a challenge. Previously, we developed the ASSA algorithm that detects statistically significant inter-molecular RNA-RNA interactions. Here we designed a workflow that uses ASSA predictions to estimate the ability of an lncRNA to function via direct base pairing with the target transcripts (co- or post-transcriptionally). The workflow was applied to 300+ lncRNA knockdown experiments from the FANTOM6 pilot project producing statistically significant predictions for 71 unique lncRNAs (104 knockdowns). Surprisingly, the majority of these lncRNAs were likely to function co-transcriptionally, i.e., hybridize with the nascent transcripts of the target genes. Moreover, a number of the obtained predictions were supported by independent iMARGI experimental data on co-localization of lncRNA and chromatin. We detected an evolutionarily conserved lncRNA CHASERR (AC013394.2 or LINC01578) that could regulate target genes co-transcriptionally via interaction with a nascent transcript by directing CHD2 helicase. The obtained results suggested that this nuclear lncRNA may be able to activate expression of the target genes in trans by base-pairing with the nascent transcripts and directing the CHD2 helicase to the regulated promoters leading to open the chromatin and active transcription. Our study highlights the possible importance of base-pairing between nuclear lncRNAs and nascent transcripts for the regulation of gene expression.

Chromatin-Associated Membraneless Organelles in Regulation of Cellular Differentiation

Membrane-free intracellular biocondensates are enclosures of proteins and nucleic acids that form by phase separation. Extensive ensembles of nuclear "membraneless organelles" indicate their involvement in genome regulation. Indeed, nuclear bodies have been linked to regulation of gene expression by formation of condensates made of chromatin and RNA processing factors. Important questions pertain to the involvement of membraneless organelles in determining cell identity through their cell-type-specific composition and function. Paraspeckles provide a prism to these questions because they exhibit striking cell-type-specific patterns and since they are crucial in embryogenesis. Here, we outline known interactions between paraspeckles and chromatin, and postulate how such interactions may be important in regulation of cell fate transitions. Moreover, we propose long non-coding RNAs (lncRNAs) as candidates for similar regulation because many form foci that resemble biocondensates and exhibit dynamic patterns during differentiation. Finally, we outline approaches that could ascertain how chromatin-associated membraneless organelles regulate cellular differentiation.

Unraveling the 3D Genome Architecture in Plants: Present and Future

The eukaryotic genome has a hierarchical three-dimensional (3D) organization with functional implications for DNA replication, DNA repair, and transcriptional regulation. Over the past decade, scientists have endeavored to elucidate the spatial characteristics and functions of plant genome architecture using high-throughput chromatin conformation capturing technologies such as Hi-C, ChIA-PET, and HiChIP. Here, we systematically review current understanding of chromatin organization in plants at multiple scales. We also discuss the emerging opinions and concepts in 3D genome research, focusing on state-of-the-art 3D genome techniques, RNA-chromatin interactions, liquid-liquid phase separation, and dynamic chromatin alterations. We propose the application of single-cell/single-molecule multi-omics, multiway (DNA-DNA, DNA-RNA, and RNA-RNA interactions) chromatin conformation capturing methods, and proximity ligation-independent 3D genome-mapping technologies to explore chromatin organization structure and function in plants. Such methods could reveal the spatial interactions between trait-related SNPs and their target genes at various spatiotemporal resolutions, and elucidate the molecular mechanisms of the interactions among DNA elements, RNA molecules, and protein factors during the formation of key traits in plants.

Stress-induced RNA-chromatin interactions promote endothelial dysfunction

Chromatin-associated RNA (caRNA) has been proposed as a type of epigenomic modifier. Here, we test whether environmental stress can induce cellular dysfunction through modulating RNA-chromatin interactions. We induce endothelial cell (EC) dysfunction with high glucose and TNFα (H + T), that mimic the common stress in diabetes mellitus. We characterize the H + T-induced changes in gene expression by single cell (sc)RNA-seq, DNA interactions by Hi-C, and RNA-chromatin interactions by iMARGI. H + T induce inter-chromosomal RNA-chromatin interactions, particularly among the super enhancers. To test the causal relationship between H + T-induced RNA-chromatin interactions and the expression of EC dysfunction-related genes, we suppress the LINC00607 RNA. This suppression attenuates the expression of SERPINE1, a critical pro-inflammatory and pro-fibrotic gene. Furthermore, the changes of the co-expression gene network between diabetic and healthy donor-derived ECs corroborate the H + T-induced RNA-chromatin interactions. Taken together, caRNA-mediated dysregulation of gene expression modulates EC dysfunction, a crucial mechanism underlying numerous diseases.

Global interaction of chromatin-associated RNAs and DNA can be identified in situ. Here the authors report the genome-wide increase of interchromosomal RNA-DNA interactions and demonstrate the importance of such RNA-DNA contacts exemplified by LINC00607 RNA and SERPINE1 gene’s super enhancer in dysfunctional endothelial cell models.

Decoding the plant genome: From epigenome to 3D organization

The linear genome of eukaryotes is partitioned into diverse chromatin states and packaged into a three-dimensional (3D) structure, which has functional implications in DNA replication, DNA repair, and transcriptional regulation. Over the past decades, research on plant functional genomics and epigenomics has made great progress, with thousands of genes cloned and molecular mechanisms of diverse biological processes elucidated. Recently, 3D genome research has gradually attracted great attention of many plant researchers. Herein, we briefly review the progress in genomic and epigenomic research in plants, with a focus on Arabidopsis and rice, and summarize the currently used technologies and advances in plant 3D genome organization studies. We also discuss the relationships between one-dimensional linear genome sequences, epigenomic states, and the 3D chromatin architecture. This review provides basis for future research on plant 3D genomics.

Weak interactions in higher-order chromatin organization

The detailed principles of the hierarchical folding of eukaryotic chromosomes have been revealed during the last two decades. Along with structures composing three-dimensional (3D) genome organization (chromatin compartments, topologically associating domains, chromatin loops, etc.), the molecular mechanisms that are involved in their establishment and maintenance have been characterized. Generally, protein-protein and protein-DNA interactions underlie the spatial genome organization in eukaryotes. However, it is becoming increasingly evident that weak interactions, which exist in biological systems, also contribute to the 3D genome. Here, we provide a snapshot of our current understanding of the role of the weak interactions in the establishment and maintenance of the 3D genome organization. We discuss how weak biological forces, such as entropic forces operating in crowded solutions, electrostatic interactions of the biomolecules, liquid-liquid phase separation, DNA supercoiling, and RNA environment participate in chromosome segregation into structural and functional units and drive intranuclear functional compartmentalization.

Studying RNA-DNA interactome by Red-C identifies noncoding RNAs associated with various chromatin types and reveals transcription dynamics

Non-coding RNAs (ncRNAs) participate in various biological processes, including regulating transcription and sustaining genome 3D organization. Here, we present a method termed Red-C that exploits proximity ligation to identify contacts with the genome for all RNA molecules present in the nucleus. Using Red-C, we uncovered the RNA-DNA interactome of human K562 cells and identified hundreds of ncRNAs enriched in active or repressed chromatin, including previously undescribed RNAs. Analysis of the RNA-DNA interactome also allowed us to trace the kinetics of messenger RNA production. Our data support the model of co-transcriptional intron splicing, but not the hypothesis of the circularization of actively transcribed genes.

Long non-coding RNAs and transposable elements: A functional relationship

Long non-coding RNAs (lncRNAs) have become increasingly important in the past decade. They are known to regulate gene expression and to interact with chromatin, proteins and other coding and non-coding RNAs. The study of lncRNAs has been challenging due to their low expression and the lack of tools developed to adapt to their particular features. Studies on lncRNAs performed to date have largely focused on cellular functions, whereas details on the mechanism of action has only been thoroughly investigated for a small number of lncRNAs. Nevertheless, some studies have highlighted the potential of these transcripts to contain functional domains, following the same accepted trend as proteins. Interestingly, many of these identified "domains" are attributed to functional units derived from transposable elements. Here, we review several types of functions of lncRNAs and relate these functions to lncRNA-embedded transposable elements.

Functional annotation of human long noncoding RNAs via molecular phenotyping

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

Mapping Transcriptome-Wide and Genome-Wide RNA-DNA Contacts with Chromatin-Associated RNA Sequencing (ChAR-seq)

RNAs play key roles in the cell as molecular intermediates for protein synthesis and as regulators of nuclear processes such as splicing, posttranscriptional regulation, or chromatin remodeling. Various classes of non-coding RNAs, including long non-coding RNAs (lncRNAs), can bind chromatin either directly or via interaction with chromatin binding proteins. It has been proposed that lncRNAs regulate cell-state-specific genes by coordinating the locus-dependent activity of chromatin-modifying complexes. Yet, the vast majority of lncRNAs have unknown functions, and we know little about the specific loci they regulate. A key step toward understanding chromatin regulation by RNAs is to map the genomic loci with which every nuclear RNA interacts and, reciprocally, to identify all RNAs that target a given locus. Our ability to generate such data has been limited, until recently, by the lack of methods to probe the genomic localization of more than a few RNAs at a time. Here, we describe a protocol for ChAR-seq, an RNA-DNA proximity ligation method that maps the binding loci for thousands of RNAs at once and without the need for specific RNA or DNA probe sequences. The ChAR-seq approach generates chimeric RNA-DNA molecules in situ and then converts those chimeras to DNA for next-generation sequencing. Using ChAR-seq we detect many types of chromatin-associated RNA, both coding and non-coding. Understanding the RNA-DNA interactome and its changes during differentiation or disease with ChAR-seq will likely provide key insights into chromatin and RNA biology.

The New Frontier of Functional Genomics: From Chromatin Architecture and Noncoding RNAs to Therapeutic Targets

Common diseases are complex, multifactorial disorders whose pathogenesis is influenced by the interplay of genetic predisposition and environmental factors. Genome-wide association studies have interrogated genetic polymorphisms across genomes of individuals to test associations between genotype and susceptibility to specific disorders, providing insights into the genetic architecture of several complex disorders. However, genetic variants associated with the susceptibility to common diseases are often located in noncoding regions of the genome, such as tissue-specific enhancers or long noncoding RNAs, suggesting that regulatory elements might play a relevant role in human diseases. Enhancers are cis-regulatory genomic sequences that act in concert with promoters to regulate gene expression in a precise spatiotemporal manner. They can be located at a considerable distance from their cognate target promoters, increasing the difficulty of their identification. Genomes are organized in domains of chromatin folding, namely topologically associating domains (TADs). Identification of enhancer-promoter interactions within TADs has revealed principles of cell-type specificity across several organisms and tissues. The vast majority of mammalian genomes are pervasively transcribed, accounting for a previously unappreciated complexity of the noncoding RNA fraction. Particularly, long noncoding RNAs have emerged as key players for the establishment of chromatin architecture and regulation of gene expression. In this perspective, we describe the new advances in the fields of transcriptomics and genome organization, focusing on the role of noncoding genomic variants in the predisposition of common diseases. Finally, we propose a new framework for the identification of the next generation of pharmacological targets for common human diseases.

MALAT1 Long Non-Coding RNA: Functional Implications

The mammalian genome is pervasively transcribed and the functional significance of many long non-coding RNA (lncRNA) transcripts are gradually being elucidated. Metastasis Associated Lung Adenocarcinoma Transcript 1 (MALAT1) is one of the most well-studied lncRNAs. MALAT1 is a highly conserved nuclear retained lncRNA that is abundantly expressed in cells and tissues and has been shown to play a role in regulating genes at both the transcriptional and post-transcriptional levels in a context-dependent manner. However, Malat1 has been shown to be dispensable for normal development and viability in mice. Interestingly, accumulating evidence suggests that MALAT1 plays an important role in numerous diseases including cancer. Here, we discuss the current state-of-knowledge in regard to MALAT1 with respect to its function, role in diseases, and the potential therapeutic opportunities for targeting MALAT1 using antisense oligonucleotides and small molecules.

Functions and properties of nuclear lncRNAs-from systematically mapping the interactomes of lncRNAs

Protein and DNA have been considered as the major components of chromatin. But beyond that, an increasing number of studies show that RNA occupies a large amount of chromatin and acts as a regulator of nuclear architecture. A significant fraction of long non-coding RNAs (lncRNAs) prefers to stay in the nucleus and cooperate with protein complexes to modulate epigenetic regulation, phase separation, compartment formation, and nuclear organization. An RNA strand also can invade into double-stranded DNA to form RNA:DNA hybrids (R-loops) in living cells, contributing to the regulation of gene expression and genomic instability. In this review, we discuss how nuclear lncRNAs orchestrate cellular processes through their interactions with proteins and DNA and summarize the recent genome-wide techniques to study the functions of lncRNAs by revealing their interactomes in vivo.

RADICL-seq identifies general and cell type-specific principles of genome-wide RNA-chromatin interactions

Mammalian genomes encode tens of thousands of noncoding RNAs. Most noncoding transcripts exhibit nuclear localization and several have been shown to play a role in the regulation of gene expression and chromatin remodeling. To investigate the function of such RNAs, methods to massively map the genomic interacting sites of multiple transcripts have been developed; however, these methods have some limitations. Here, we introduce RNA And DNA Interacting Complexes Ligated and sequenced (RADICL-seq), a technology that maps genome-wide RNA-chromatin interactions in intact nuclei. RADICL-seq is a proximity ligation-based methodology that reduces the bias for nascent transcription, while increasing genomic coverage and unique mapping rate efficiency compared with existing methods. RADICL-seq identifies distinct patterns of genome occupancy for different classes of transcripts as well as cell type-specific RNA-chromatin interactions, and highlights the role of transcription in the establishment of chromatin structure.

Short Tandem Repeat-Enriched Architectural RNAs in Nuclear Bodies: Functions and Associated Diseases

Nuclear bodies are membraneless, phase-separated compartments that concentrate specific proteins and RNAs in the nucleus. They are believed to serve as sites for the modification, sequestration, and storage of specific factors, and to act as organizational hubs of chromatin structure to control gene expression and cellular function. Architectural (arc) RNA, a class of long noncoding RNA (lncRNA), plays essential roles in the formation of nuclear bodies. Herein, we focus on specific arcRNAs containing short tandem repeat-enriched sequences and introduce their biological functions and recently elucidated underlying molecular mechanism. In various neurodegenerative diseases, abnormal nuclear and cytoplasmic bodies are built on disease-causing RNAs or toxic RNAs with aberrantly expanded short tandem repeat-enriched sequences. We discuss the possible analogous functions of natural arcRNAs and toxic RNAs with short tandem repeat-enriched sequences. Finally, we describe the technical utility of short tandem repeat-enriched arcRNAs as a model for exploring the structures and functions of nuclear bodies, as well as the pathogenic mechanisms of neurodegenerative diseases.

Practical Guidance in Genome-Wide RNA:DNA Triple Helix Prediction

Long noncoding RNAs (lncRNAs) play a key role in many cellular processes including chromatin regulation. To modify chromatin, lncRNAs often interact with DNA in a sequence-specific manner forming RNA:DNA triple helices. Computational tools for triple helix search do not always provide genome-wide predictions of sufficient quality. Here, we used four human lncRNAs (MEG3, DACOR1, TERC and HOTAIR) and their experimentally determined binding regions for evaluating triplex parameters that provide the highest prediction accuracy. Additionally, we combined triplex prediction with the lncRNA secondary structure and demonstrated that considering only single-stranded fragments of lncRNA can further improve DNA-RNA triplexes prediction.

RNAInter in 2020: RNA interactome repository with increased coverage and annotation

Research on RNA-associated interactions has exploded in recent years, and increasing numbers of studies are not limited to RNA-RNA and RNA-protein interactions but also include RNA-DNA/compound interactions. To facilitate the development of the interactome and promote understanding of the biological functions and molecular mechanisms of RNA, we updated RAID v2.0 to RNAInter (RNA Interactome Database), a repository for RNA-associated interactions that is freely accessible at http://www.rna-society.org/rnainter/ or http://www.rna-society.org/raid/. Compared to RAID v2.0, new features in RNAInter include (i) 8-fold more interaction data and 94 additional species; (ii) more definite annotations organized, including RNA editing/localization/modification/structure and homology interaction; (iii) advanced functions including fuzzy/batch search, interaction network and RNA dynamic expression and (iv) four embedded RNA interactome tools: RIscoper, IntaRNA, PRIdictor and DeepBind. Consequently, RNAInter contains >41 million RNA-associated interaction entries, involving more than 450 thousand unique molecules, including RNA, protein, DNA and compound. Overall, RNAInter provides a comprehensive RNA interactome resource for researchers and paves the way to investigate the regulatory landscape of cellular RNAs.

PIRCh-seq: functional classification of non-coding RNAs associated with distinct histone modifications

We develop PIRCh-seq, a method which enables a comprehensive survey of chromatin-associated RNAs in a histone modification-specific manner. We identify hundreds of chromatin-associated RNAs in several cell types with substantially less contamination by nascent transcripts. Non-coding RNAs are found enriched on chromatin and are classified into functional groups based on the patterns of their association with specific histone modifications. We find single-stranded RNA bases are more chromatin-associated, and we discover hundreds of allele-specific RNA-chromatin interactions. These results provide a unique resource to globally study the functions of chromatin-associated lncRNAs and elucidate the basic mechanisms of chromatin-RNA interactions.

Understanding Long Noncoding RNA and Chromatin Interactions: What We Know So Far

With the evolution of technologies that deal with global detection of RNAs to probing of lncRNA-chromatin interactions and lncRNA-chromatin structure regulation, we have been updated with a comprehensive repertoire of chromatin interacting lncRNAs, their genome-wide chromatin binding regions and mode of action. Evidence from these new technologies emphasize that chromatin targeting of lncRNAs is a prominent mechanism and that these chromatin targeted lncRNAs exert their functionality by fine tuning chromatin architecture resulting in an altered transcriptional readout. Currently, there are no unifying principles that define chromatin association of lncRNAs, however, evidence from a few chromatin-associated lncRNAs show presence of a short common sequence for chromatin targeting. In this article, we review how technological advancements contributed in characterizing chromatin associated lncRNAs, and discuss the potential mechanisms by which chromatin associated lncRNAs execute their functions.

Multimodal Long Noncoding RNA Interaction Networks: Control Panels for Cell Fate Specification

Lineage specification in early development is the basis for the exquisitely precise body plan of multicellular organisms. It is therefore critical to understand cell fate decisions in early development. Moreover, for regenerative medicine, the accurate specification of cell types to replace damaged/diseased tissue is strongly dependent on identifying determinants of cell identity. Long noncoding RNAs (lncRNAs) have been shown to regulate cellular plasticity, including pluripotency establishment and maintenance, differentiation and development, yet broad phenotypic analysis and the mechanistic basis of their function remains lacking. As components of molecular condensates, lncRNAs interact with almost all classes of cellular biomolecules, including proteins, DNA, mRNAs, and microRNAs. With functions ranging from controlling alternative splicing of mRNAs, to providing scaffolding upon which chromatin modifiers are assembled, it is clear that at least a subset of lncRNAs are far from the transcriptional noise they were once deemed. This review highlights the diversity of lncRNA interactions in the context of cell fate specification, and provides examples of each type of interaction in relevant developmental contexts. Also highlighted are experimental and computational approaches to study lncRNAs.

Principles and Practices of Hybridization Capture Experiments to Study Long Noncoding RNAs That Act on Chromatin

The diverse roles of cellular RNAs can be studied by purifying RNAs of interest together with the biomolecules they bind. Biotinylated antisense oligonucleotides that hybridize specifically to the RNA of interest provide a general approach to develop affinity reagents for these experiments. Such oligonucleotides can be used to enrich endogenous RNAs from cross-linked chromatin extracts to study the genomic binding sites of RNAs. These hybridization capture protocols are evolving modular experiments that are compatible with a range of cross-linkers and conditions. This review discusses the principles of these hybridization capture experiments as well as considerations and controls necessary to interpret the resulting data without being misled by artifactual signals.

Mapping RNA-chromatin interactions by sequencing with iMARGI

RNA-chromatin interactions represent an important aspect of the transcriptional regulation of genes and transposable elements. However, analyses of chromatin-associated RNAs (caRNAs) are often limited to one caRNA at a time. Here, we describe the iMARGI (in situ mapping of RNA-genome interactome) technique, which is used to discover caRNAs and reveal their respective genomic interaction loci. iMARGI starts with in situ crosslinking and genome fragmentation, followed by converting each proximal RNA-DNA pair into an RNA-linker-DNA chimeric sequence. These chimeric sequences are subsequently converted into a sequencing library suitable for paired-end sequencing. A standardized bioinformatic software package, iMARGI-Docker, is provided to decode the paired-end sequencing data into caRNA-DNA interactions. Compared to its predecessor MARGI (mapping RNA-genome interactions), the number of input cells for iMARGI is 3-5 million (a 100-fold reduction), experimental time is reduced, and clear checkpoints have been established. It takes a few hours a day and a total of 8 d to complete the construction of an iMARGI sequencing library and 1 d to carry out data processing with iMARGI-Docker.

Isolation and genome-wide characterization of cellular DNA:RNA triplex structures

RNA can directly bind to purine-rich DNA via Hoogsteen base pairing, forming a DNA:RNA triple helical structure that anchors the RNA to specific sequences and allows guiding of transcription regulators to distinct genomic loci. To unravel the prevalence of DNA:RNA triplexes in living cells, we have established a fast and cost-effective method that allows genome-wide mapping of DNA:RNA triplex interactions. In contrast to previous approaches applied for the identification of chromatin-associated RNAs, this method uses protein-free nucleic acids isolated from chromatin. High-throughput sequencing and computational analysis of DNA-associated RNA revealed a large set of RNAs which originate from non-coding and coding loci, including super-enhancers and repeat elements. Combined analysis of DNA-associated RNA and RNA-associated DNA identified genomic DNA:RNA triplex structures. The results suggest that triplex formation is a general mechanism of RNA-mediated target-site recognition, which has major impact on biological functions.

Chromatin-associated RNAs as facilitators of functional genomic interactions

Mammalian genomes are extensively transcribed, which produces a large number of both coding and non-coding transcripts. Various RNAs are physically associated with chromatin, through being either retained in cis at their site of transcription or recruited in trans to other genomic regions. Driven by recent technological innovations for detecting chromatin-associated RNAs, diverse roles are being revealed for these RNAs and associated RNA-binding proteins (RBPs) in gene regulation and genome function. Such functions include locus-specific roles in gene activation and silencing, as well as emerging roles in higher-order genome organization, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions.

A Polymer Physics Framework for the Entropy of Arbitrary Pseudoknots

The accurate prediction of RNA secondary structure from primary sequence has had enormous impact on research from the past 40 years. Although many algorithms are available to make these predictions, the inclusion of non-nested loops, termed pseudoknots, still poses challenges arising from two main factors: 1) no physical model exists to estimate the loop entropies of complex intramolecular pseudoknots, and 2) their NP-complete enumeration has impeded their study. Here, we address both challenges. First, we develop a polymer physics model that can address arbitrarily complex pseudoknots using only two parameters corresponding to concrete physical quantities-over an order of magnitude fewer than the sparsest state-of-the-art phenomenological methods. Second, by coupling this model to exhaustive enumeration of the set of possible structures, we compute the entire free energy landscape of secondary structures resulting from a primary RNA sequence. We demonstrate that for RNA structures of ∼80 nucleotides, with minimal heuristics, the complete enumeration of possible secondary structures can be accomplished quickly despite the NP-complete nature of the problem. We further show that despite our loop entropy model's parametric sparsity, it performs better than or on par with previously published methods in predicting both pseudoknotted and non-pseudoknotted structures on a benchmark data set of RNA structures of ≤80 nucleotides. We suggest ways in which the accuracy of the model can be further improved.

Interaction of Sox2 with RNA binding proteins in mouse embryonic stem cells

Sox2 is a master transcriptional regulator of embryonic development. In this study, we determined the protein interactome of Sox2 in the chromatin and nucleoplasm of mouse embryonic stem (mES) cells. Apart from canonical interactions with pluripotency-regulating transcription factors, we identified interactions with several chromatin modulators, including members of the heterochromatin protein 1 (HP1) family, suggesting a role for Sox2 in chromatin-mediated transcriptional repression. Sox2 was also found to interact with RNA binding proteins (RBPs), including proteins involved in RNA processing. RNA immunoprecipitation followed by sequencing revealed that Sox2 associates with different messenger RNAs, as well as small nucleolar RNA Snord34 and the non-coding RNA 7SK. 7SK has been shown to regulate transcription at gene regulatory regions, which could suggest a functional interaction with Sox2 for chromatin recruitment. Nevertheless, we found no evidence of Sox2 modulating recruitment of 7SK to chromatin when examining 7SK chromatin occupancy by Chromatin Isolation by RNA Purification (ChIRP) in Sox2 depleted mES cells. In addition, knockdown of 7SK in mES cells did not lead to any change in Sox2 occupancy at 7SK-regulated genes. Thus, our results show that Sox2 extensively interacts with RBPs, and suggest that Sox2 and 7SK co-exist in a ribonucleoprotein complex whose function is not to regulate chromatin recruitment, but could rather regulate other processes in the nucleoplasm.

Prediction of lncRNAs and their interactions with nucleic acids: benchmarking bioinformatics tools

The genomes of mammalian species are pervasively transcribed producing as many noncoding as protein-coding RNAs. There is a growing body of evidence supporting their functional role. Long noncoding RNA (lncRNA) can bind both nucleic acids and proteins through several mechanisms. A reliable computational prediction of the most probable mechanism of lncRNA interaction can facilitate experimental validation of its function. In this study, we benchmarked computational tools capable to discriminate lncRNA from mRNA and predict lncRNA interactions with other nucleic acids. We assessed the performance of 9 tools for distinguishing protein-coding from noncoding RNAs, as well as 19 tools for prediction of RNA-RNA and RNA-DNA interactions. Our conclusions about the considered tools were based on their performances on the entire genome/transcriptome level, as it is the most common task nowadays. We found that FEELnc and CPAT distinguish between coding and noncoding mammalian transcripts in the most accurate manner. ASSA, RIBlast and LASTAL, as well as Triplexator, turned out to be the best predictors of RNA-RNA and RNA-DNA interactions, respectively. We showed that the normalization of the predicted interaction strength to the transcript length and GC content may improve the accuracy of inferring RNA interactions. Yet, all the current tools have difficulties to make accurate predictions of short-trans RNA-RNA interactions-stretches of sparse contacts. All over, there is still room for improvement in each category, especially for predictions of RNA interactions.

Profiling the long noncoding RNA interaction network in the regulatory elements of target genes by chromatin in situ reverse transcription sequencing

Long noncoding RNAs (lncRNAs) can regulate the activity of target genes by participating in the organization of chromatin architecture. We have devised a “chromatin-RNA in situ reverse transcription sequencing” (CRIST-seq) approach to profile the lncRNA interaction network in gene regulatory elements by combining the simplicity of RNA biotin labeling with the specificity of the CRISPR/Cas9 system. Using gene-specific gRNAs, we describe a pluripotency-specific lncRNA interacting network in the promoters of Sox2 and Pou5f1, two critical stem cell factors that are required for the maintenance of pluripotency. The promoter-interacting lncRNAs were specifically activated during reprogramming into pluripotency. Knockdown of these lncRNAs caused the stem cells to exit from pluripotency. In contrast, overexpression of the pluripotency-associated lncRNA activated the promoters of core stem cell factor genes and enhanced fibroblast reprogramming into pluripotency. These CRIST-seq data suggest that the Sox2 and Pou5f1 promoters are organized within a unique lncRNA interaction network that determines the fate of pluripotency during reprogramming. This CRIST approach may be broadly used to map lncRNA interaction networks at target loci across the genome.

GRID-seq for comprehensive analysis of global RNA-chromatin interactions

Chromatin in higher eukaryotic nuclei is extensively bound by various RNA species. We recently developed a method for in situ capture of global RNA interactions with DNA by deep sequencing (GRID-seq) of fixed permeabilized nuclei that allows identification of the entire repertoire of chromatin-associated RNAs in an unbiased manner. The experimental design of GRID-seq is related to those of two recently published strategies (MARGI (mapping RNA-genome interactions) and ChAR-seq (chromatin-associated RNA sequencing)), which also use a bivalent linker to ligate RNA and DNA in proximity. Importantly, however, GRID-seq also implements a combined experimental and computational approach to control nonspecific RNA-DNA interactions that are likely to occur during library construction, which is critical for accurate interpretation of detected RNA-DNA interactions. GRID-seq typically finds both coding and non-coding RNAs (ncRNAs) that interact with tissue-specific promoters and enhancers, especially super-enhancers, from which a global promoter-enhancer connectivity map can be deduced. Here, we provide a detailed protocol for GRID-seq that includes nuclei preparation, chromatin fragmentation, RNA and DNA in situ ligation with a bivalent linker, PCR amplification and high-throughput sequencing. To further enhance the utility of GRID-seq, we include a pipeline for data analysis, called GridTools, into which key steps such as background correction and inference of genomic element proximity are integrated. For researchers experienced in molecular biology with minimal bioinformatics skills, the protocol typically takes 4-5 d from cell fixation to library construction and 2-3 d for data processing.

Understanding the 3D genome: Emerging impacts on human disease

Recent burst of new technologies that allow for quantitatively delineating chromatin structure has greatly expanded our understanding of how the genome is organized in the three-dimensional (3D) space of the nucleus. It is now clear that the hierarchical organization of the eukaryotic genome critically impacts nuclear activities such as transcription, replication, as well as cellular and developmental events such as cell cycle, cell fate decision and embryonic development. In this review, we discuss new insights into how the structural features of the 3D genome hierarchy are established and maintained, how this hierarchy undergoes dynamic rearrangement during normal development and how its perturbation will lead to human disease, highlighting the accumulating evidence that links the diverse 3D genome architecture components to a multitude of human diseases and the emerging mechanisms by which 3D genome derangement causes disease phenotypes.

Cellular functions of long noncoding RNAs

A diverse catalog of long noncoding RNAs (lncRNAs), which lack protein-coding potential, are transcribed from the mammalian genome. They are emerging as important regulators in gene expression networks by controlling nuclear architecture and transcription in the nucleus and by modulating mRNA stability, translation and post-translational modifications in the cytoplasm. In this Review, we highlight recent progress in cellular functions of lncRNAs at the molecular level in mammalian cells.

lncRedibly versatile: biochemical and biological functions of long noncoding RNAs

Long noncoding RNAs (lncRNAs) are transcripts that do not code for proteins, but nevertheless exert regulatory effects on various biochemical pathways, in part via interactions with proteins, DNA, and other RNAs. LncRNAs are thought to regulate transcription and other biological processes by acting, for example, as guides that target proteins to chromatin, scaffolds that facilitate protein-protein interactions and complex formation, and orchestrators of phase-separated compartments. The study of lncRNAs has reached an exciting time, as recent advances in experimental and computational methods allow for genome-wide interrogation of biochemical and biological mechanisms of these enigmatic transcripts. A better appreciation for the biochemical versatility of lncRNAs has allowed us to begin closing gaps in our knowledge of how they act in diverse cellular and organismal contexts, including development and disease.

Chromatin-Associated RNA Sequencing (ChAR-seq)

RNA is a fundamental component of chromatin. Noncoding RNAs (ncRNAs) can associate with chromatin to influence gene expression and chromatin state; many also act at long distances from their transcriptional origin. Yet we know almost nothing about the functions or sites of action for most ncRNAs. Current methods to identify sites of RNA interaction with the genome are limited to the study of a single RNA at a time. Here we describe a protocol for ChAR-seq, a strategy to identify all chromatin-associated RNAs and map their DNA contacts genome-wide. In ChAR-seq, proximity ligation of RNA and DNA to a linker molecule is used to construct a chimeric RNA-DNA molecule that is converted to DNA for sequencing. In a single assay, ChAR-seq can discover de novo chromatin interactions of distinct RNAs, including nascent transcripts, splicing RNAs, and long noncoding RNAs (lncRNAs). Resulting "maps" of genome-bound RNAs should provide new insights into RNA biology. © 2019 by John Wiley & Sons, Inc.