Chromatin-associated RNA sequencing (ChAR-seq) maps genome-wide RNA-to-DNA contacts

Unraveling architectural RNAs: Structural and functional blueprints of membraneless organelles and strategies for genome-scale identification

Architectural RNAs (arcRNAs) are long noncoding RNAs that serve as structural scaffolds for membraneless organelles (MLOs), facilitating cellular organization and dynamic responses to stimuli. Acting as blueprints for MLO assembly, arcRNAs recruit specific proteins and nucleic acids to establish and maintain the internal structure of MLOs while coordinating their spatial relationships with other organelles. This organized framework enables precise spatiotemporal regulation, allowing for targeted control of transcription, RNA processing, and cellular responses to stress. Notably, arcRNAs exhibit the "semi-extractable" feature, a property derived from their stable binding to cellular structures, making them partially resistant to conventional RNA extraction methods. This unique feature serves as a useful criterion for identifying novel arcRNAs, providing an opportunity to accelerate research in long noncoding RNAs and deepen our understanding of their functional roles in cellular processes.

Biological functions and affected signaling pathways by Long Non-Coding RNAs in the immune system

Recently, the various regulative functions of long non-coding RNAs (LncRNAs) have been well determined. Recently, the vital role of LncRNAs as gene regulators has been identified in the immune system, especially in the inflammatory response. All cells of the immune system are governed by a complex and ever-changing gene expression program that is regulated through both transcriptional and post-transcriptional processes. LncRNAs regulate gene expression within the cell nucleus by influencing transcription or through post-transcriptional processes that affect the splicing, stability, or translation of messenger RNAs (mRNAs). Recent studies in immunology have revealed substantial alterations in the expression of lncRNAs during the activation of the innate immune system as well as the development, differentiation, and activation of T cells. These lncRNAs regulate key aspects of immune function, including the manufacturing of inflammatory molecules, cellular distinction, and cell movement. They do this by modulating protein-protein interactions or through base pairing with RNA and DNA. Here we review the current understanding of the mechanism of action of lncRNAs as novel immune-related regulators and their impact on physiological and pathological processes related to the immune system, including autoimmune diseases. We also highlight the emerging pattern of gene expression control in important research areas at the intersection between immunology and lncRNA biology.

RADIP technology comprehensively identifies H3K27me3-associated RNA-chromatin interactions

Many RNAs associate with chromatin, either directly or indirectly. Several technologies for mapping regions where RNAs interact across the genome have been developed to investigate the function of these RNAs. Obtaining information on the proteins involved in these RNA-chromatin interactions is critical for further analysis. Here, we developed RADIP [RNA and DNA interacting complexes ligated and sequenced (RADICL-seq) with immunoprecipitation], a novel technology that combines RADICL-seq technology with chromatin immunoprecipitation to characterize RNA-chromatin interactions mediated by individual proteins. Building upon the foundational principles of RADICL-seq, RADIP extends its advantages by increasing genomic coverage and unique mapping rate efficiency compared to existing methods. To demonstrate its effectiveness, we applied an anti-H3K27me3 antibody to the RADIP technology and generated libraries from mouse embryonic stem cells (mESCs). We identified a multitude of RNAs, including RNAs from protein-coding genes and non-coding RNAs, that are associated with chromatin via H3K27me3 and that likely facilitate the spread of Polycomb repressive complexes over broad regions of the mammalian genome, thereby affecting gene expression, chromatin structures and pluripotency of mESCs. Our study demonstrates the applicability of RADIP to investigations of the functions of chromatin-associated RNAs.

Emerging interactions between RNA methylation and chromatin architecture

Epitranscriptomics, the study of chemical modifications of RNA molecules, is increasingly recognized as an important component of gene expression regulation. While the majority of research has focused on N6-methyladenosine (m6A) RNA methylation on mRNAs, emerging evidence has revealed that the m6A modification extends beyond mRNAs to include chromatin-associated RNAs (caRNAs). CaRNAs constitute an important class of RNAs characterized by their interaction with the genome and epigenome. These features allow caRNAs to be actively involved in shaping genome organization. In this review, we bring into focus recent findings of the dynamic interactions between caRNAs and chromatin architecture and how RNA methylation impacts caRNAs' function in this interplay. We highlight several enabling techniques, which were critical for genome-wide profiling of caRNAs and their modifications. Given the nascent stage of the field, we emphasize on the need to address critical gaps in study of these modifications in more relevant biological systems. Overall, these exciting progress have expanded the scope and reach of epitranscriptomics, unveiling new mechanisms that underpin the control of gene expression and cellular phenotypes, with potential therapeutic implications.

HOXDeRNA activates a cancerous transcription program and super enhancers via genome-wide binding

The role of long non-coding RNAs (lncRNAs) in malignant cell transformation remains elusive. We previously identified an enhancer-associated lncRNA, LINC01116 (named HOXDeRNA), as a transformative factor converting human astrocytes into glioma-like cells. Employing a combination of CRISPR editing, chromatin isolation by RNA purification coupled with sequencing (ChIRP-seq), in situ mapping RNA-genome interactions (iMARGI), chromatin immunoprecipitation sequencing (ChIP-seq), HiC, and RNA/DNA FISH, we found that HOXDeRNA directly binds to CpG islands within the promoters of 35 glioma-specific transcription factors (TFs) distributed throughout the genome, including key stem cell TFs SOX2, OLIG2, POU3F2, and ASCL1, liberating them from PRC2 repression. This process requires a distinct RNA quadruplex structure and other segments of HOXDeRNA, interacting with EZH2 and CpGs, respectively. Subsequent transformation activates multiple oncogenes (e.g., EGFR, miR-21, and WEE1), driven by the SOX2- and OLIG2-dependent glioma-specific super enhancers. These results help reconstruct the sequence of events underlying the process of astrocyte transformation, highlighting HOXDeRNA's central genome-wide activity and suggesting a shared RNA-dependent mechanism in otherwise heterogeneous and multifactorial gliomagenesis.

MAYEX is an old long noncoding RNA recruited for X chromosome dosage compensation in a reptile

Long noncoding RNAs (lncRNAs) are essential regulatory elements of sex chromosomes that act to equalize gene expression levels between males and females. XIST, RSX, and roX2 regulate X chromosomes in placental mammals, marsupials, and Drosophila, respectively. Because the green anole (Anolis carolinensis) shows complete dosage compensation of its X chromosome, we tested whether a lncRNA was involved. We found an ancient lncRNA, MAYEX, that gained male-specific expression more than 89 million years ago. MAYEX evolved a notable association with the acetylated histone 4 lysine 16 (H4K16ac) epigenetic mark and the ability to loop its locus to the totality of the X chromosome to increase expression levels. MAYEX is the first lncRNA in reptiles linked to a dosage compensation mechanism that balances the expression of sex chromosomes.

Nuclear RNA: a transcription-dependent regulator of chromatin structure

Although the majority of RNAs are retained in the nucleus, their significance is often overlooked. However, it is now becoming clear that nuclear RNA forms a dynamic structure through interacting with various proteins that can influence the three-dimensional structure of chromatin. We review the emerging evidence for a nuclear RNA mesh or gel, highlighting the interplay between DNA, RNA and RNA-binding proteins (RBPs), and assessing the critical role of protein and RNA in governing chromatin architecture. We also discuss a proposed role for the formation and regulation of the nuclear gel in transcriptional control. We suggest that it may concentrate the transcriptional machinery either by direct binding or inducing RBPs to form microphase condensates, nanometre sized membraneless structures with distinct properties to the surrounding medium and an enrichment of particular macromolecules.

The landscape of RNA-chromatin interaction reveals small non-coding RNAs as essential mediators of leukemia maintenance

RNA constitutes a large fraction of chromatin. Spatial distribution and functional relevance of most of RNA-chromatin interactions remain unknown. We established a landscape analysis of RNA-chromatin interactions in human acute myeloid leukemia (AML). In total more than 50 million interactions were captured in an AML cell line. Protein-coding mRNAs and long non-coding RNAs exhibited a substantial number of interactions with chromatin in cis suggesting transcriptional activity. In contrast, small nucleolar RNAs (snoRNAs) and small nuclear RNAs (snRNAs) associated with chromatin predominantly in trans suggesting chromatin specific functions. Of note, snoRNA-chromatin interaction was associated with chromatin modifications and occurred independently of the classical snoRNA-RNP complex. Two C/D box snoRNAs, namely SNORD118 and SNORD3A, displayed high frequency of trans-association with chromatin. The transcription of SNORD118 and SNORD3A was increased upon leukemia transformation and enriched in leukemia stem cells, but decreased during myeloid differentiation. Suppression of SNORD118 and SNORD3A impaired leukemia cell proliferation and colony forming capacity in AML cell lines and primary patient samples. Notably, this effect was leukemia specific with less impact on healthy CD34+ hematopoietic stem and progenitor cells. These findings highlight the functional importance of chromatin-associated RNAs overall and in particular of SNORD118 and SNORD3A in maintaining leukemia propagation.

Exploring the roles of RNAs in chromatin architecture using deep learning

Recent studies have highlighted the impact of both transcription and transcripts on 3D genome organization, particularly its dynamics. Here, we propose a deep learning framework, called AkitaR, that leverages both genome sequences and genome-wide RNA-DNA interactions to investigate the roles of chromatin-associated RNAs (caRNAs) on genome folding in HFFc6 cells. In order to disentangle the cis- and trans-regulatory roles of caRNAs, we have compared models with nascent transcripts, trans-located caRNAs, open chromatin data, or DNA sequence alone. Both nascent transcripts and trans-located caRNAs improve the models' predictions, especially at cell-type-specific genomic regions. Analyses of feature importance scores reveal the contribution of caRNAs at TAD boundaries, chromatin loops and nuclear sub-structures such as nuclear speckles and nucleoli to the models' predictions. Furthermore, we identify non-coding RNAs (ncRNAs) known to regulate chromatin structures, such as MALAT1 and NEAT1, as well as several new RNAs, RNY5, RPPH1, POLG-DT and THBS1-IT1, that might modulate chromatin architecture through trans-interactions in HFFc6. Our modeling also suggests that transcripts from Alus and other repetitive elements may facilitate chromatin interactions through trans R-loop formation. Our findings provide insights and generate testable hypotheses about the roles of caRNAs in shaping chromatin organization.

Molecular insights into regulatory RNAs in the cellular machinery

It is apparent that various functional units within the cellular machinery are derived from RNAs. The evolution of sequencing techniques has resulted in significant insights into approaches for transcriptome studies. Organisms utilize RNA to govern cellular systems, and a heterogeneous class of RNAs is involved in regulatory functions. In particular, regulatory RNAs are increasingly recognized to participate in intricately functioning machinery across almost all levels of biological systems. These systems include those mediating chromatin arrangement, transcription, suborganelle stabilization, and posttranscriptional modifications. Any class of RNA exhibiting regulatory activity can be termed a class of regulatory RNA and is typically represented by noncoding RNAs, which constitute a substantial portion of the genome. These RNAs function based on the principle of structural changes through cis and/or trans regulation to facilitate mutual RNA‒RNA, RNA‒DNA, and RNA‒protein interactions. It has not been clearly elucidated whether regulatory RNAs identified through deep sequencing actually function in the anticipated mechanisms. This review addresses the dominant properties of regulatory RNAs at various layers of the cellular machinery and covers regulatory activities, structural dynamics, modifications, associated molecules, and further challenges related to therapeutics and deep learning.

RNA-DNA triplexes: molecular mechanisms and functional relevance

Interactions of RNA with DNA are principles of gene expression control that have recently gained considerable attention. Among RNA-DNA interactions are R-loops and RNA-DNA hybrid G-quadruplexes, as well as RNA-DNA triplexes. It is proposed that RNA-DNA triplexes guide RNA-associated regulatory proteins to specific genomic locations, influencing transcription and epigenetic decision making. Although triplex formation initially was considered solely an in vitro event, recent progress in computational, biochemical, and biophysical methods support in vivo functionality with relevance for gene expression control. Here, we review the central methodology and biology of triplexes, outline paradigms required for triplex function, and provide examples of physiologically important triplex-forming long non-coding RNAs.

BaRDIC: robust peak calling for RNA-DNA interaction data

Chromatin-associated non-coding RNAs play important roles in various cellular processes by targeting genomic loci. Two types of genome-wide NGS experiments exist to detect such targets: 'one-to-al', which focuses on targets of a single RNA, and 'all-to-al', which captures targets of all RNAs in a sample. As with many NGS experiments, they are prone to biases and noise, so it becomes essential to detect 'peaks'-specific interactions of an RNA with genomic targets. Here, we present BaRDIC-Binomial RNA-DNA Interaction Caller-a tailored method to detect peaks in both types of RNA-DNA interaction data. BaRDIC is the first tool to simultaneously take into account the two most prominent biases in the data: chromatin heterogeneity and distance-dependent decay of interaction frequency. Since RNAs differ in their interaction preferences, BaRDIC adapts peak sizes according to the abundances and contact patterns of individual RNAs. These features enable BaRDIC to make more robust predictions than currently applied peak-calling algorithms and better handle the characteristic sparsity of all-to-all data. The BaRDIC package is freely available at https://github.com/dmitrymyl/BaRDIC.

Advances in the multimodal analysis of the 3D chromatin structure and gene regulation

Recent studies have demonstrated that the three-dimensional conformation of the chromatin plays a crucial role in gene regulation, with aberrations potentially leading to various diseases. Advanced methodologies have revealed a link between the chromatin conformation and biological function. This review divides these methodologies into sequencing-based and imaging-based methodologies, tracing their development over time. We particularly highlight innovative techniques that facilitate the simultaneous mapping of RNAs, histone modifications, and proteins within the context of the 3D architecture of chromatin. This multimodal integration substantially improves our ability to establish a robust connection between the spatial arrangement of molecular components in the nucleus and their functional roles. Achieving a comprehensive understanding of gene regulation requires capturing diverse data modalities within individual cells, enabling the direct inference of functional relationships between these components. In this context, imaging-based technologies have emerged as an especially promising approach for gathering spatial information across multiple components in the same cell.

Single-cell multiplex chromatin and RNA interactions in ageing human brain

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an 'older' transcriptomic signature and Alzheimer's disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

The apparent loss of PRC2 chromatin occupancy as an artifact of RNA depletion

RNA has been implicated in the recruitment of chromatin modifiers, and previous studies have provided evidence in favor and against this idea. RNase treatment of chromatin is commonly used to study RNA-mediated regulation of chromatin modifiers, but the limitations of this approach remain unclear. RNase A treatment during chromatin immunoprecipitation (ChIP) reduces chromatin occupancy of the H3K27me3 methyltransferase Polycomb repressive complex 2 (PRC2). This led to suggestions of an "RNA bridge" between PRC2 and chromatin. Here, we show that RNase A treatment during ChIP causes the apparent loss of all facultative heterochromatin, including both PRC2 and H3K27me3 genome-wide. We track this observation to a gain of DNA from non-targeted chromatin, sequenced at the expense of DNA from facultative heterochromatin, which reduces ChIP signals. Our results emphasize substantial limitations in using RNase A treatment for mapping RNA-dependent chromatin occupancy and invalidate conclusions that were previously established for PRC2 based on this assay.

Single-cell multiplex chromatin and RNA interactions in aging human brain

The dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA (caRNA) 1-3. Chromatin complex compositions change during cellular differentiation and aging, and are expected to be highly heterogeneous among terminally differentiated single cells 4-7. Here we introduce the Multi-Nucleic Acid Interaction Mapping in Single Cell (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression, and RNA-chromatin associations within individual nuclei. Applied to 14 human frontal cortex samples from elderly donors, MUSIC delineates diverse cortical cell types and states. We observed the nuclei exhibiting fewer short-range chromatin interactions are correlated with an "older" transcriptomic signature and with Alzheimer's pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci (cis eQTLs) and a promoter tends to be the cell type where these cis eQTLs specifically affect their target gene's expression. Additionally, the female cortical cells exhibit highly heterogeneous interactions between the XIST non-coding RNA and Chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploring chromatin architecture and transcription at cellular resolution in complex tissues.

The regulatory roles of small nucleolar RNAs within their host locus

Small nucleolar RNAs (snoRNAs) are a class of conserved noncoding RNAs forming complexes with proteins to catalyse site-specific modifications on ribosomal RNA. Besides this canonical role, several snoRNAs are now known to regulate diverse levels of gene expression. While these functions are carried out in trans by mature snoRNAs, evidence has also been emerging of regulatory roles of snoRNAs in cis, either within their genomic locus or as longer transcription intermediates during their maturation. Herein, we review recent findings that snoRNAs can interact in cis with their intron to regulate the expression of their host gene. We also explore the ever-growing diversity of longer host-derived snoRNA extensions and their functional impact across the transcriptome. Finally, we discuss the role of snoRNA duplications into forging these new layers of snoRNA-mediated regulation, as well as their involvement in the genomic imprinting of their host locus.

Integrated Analysis of lncRNA-miRNA-mRNA Regulatory Network in Rapamycin-Induced Cardioprotection against Ischemia/Reperfusion Injury in Diabetic Rabbits

The inhibition of mammalian target of rapamycin (mTOR) with rapamycin (RAPA) provides protection against myocardial ischemia/reperfusion (I/R) injury in diabetes. Since interactions between transcripts, including long non-coding RNA (lncRNA), microRNA(miRNA) and mRNA, regulate the pathophysiology of disease, we performed unbiased miRarray profiling in the heart of diabetic rabbits following I/R injury with/without RAPA treatment to identify differentially expressed (DE) miRNAs and their predicted targets of lncRNAs/mRNAs. Results showed that among the total of 806 unique miRNAs targets, 194 miRNAs were DE after I/R in diabetic rabbits. Specifically, eight miRNAs, including miR-199a-5p, miR-154-5p, miR-543-3p, miR-379-3p, miR-379-5p, miR-299-5p, miR-140-3p, and miR-497-5p, were upregulated and 10 miRNAs, including miR-1-3p, miR-1b, miR-29b-3p, miR-29c-3p, miR-30e-3p, miR-133c, miR-196c-3p, miR-322-5p, miR-499-5p, and miR-672-5p, were significantly downregulated after I/R injury. Interestingly, RAPA treatment significantly reversed these changes in miRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated the participation of miRNAs in the regulation of several signaling pathways related to I/R injury, including MAPK signaling and apoptosis. Furthermore, in diabetic hearts, the expression of lncRNAs, HOTAIR, and GAS5 were induced after I/R injury, but RAPA suppressed these lncRNAs. In contrast, MALAT1 was significantly reduced following I/R injury, with the increased expression of miR-199a-5p and suppression of its target, the anti-apoptotic protein Bcl-2. RAPA recovered MALAT1 expression with its sponging effect on miR-199-5p and restoration of Bcl-2 expression. The identification of novel targets from the transcriptome analysis in RAPA-treated diabetic hearts could potentially lead to the development of new therapeutic strategies for diabetic patients with myocardial infarction.

The complex nature of lncRNA-mediated chromatin dynamics in multiple myeloma

Extensive genome-wide sequencing efforts have unveiled the intricate regulatory potential of long non-protein coding RNAs (lncRNAs) within the domain of haematological malignancies. Notably, lncRNAs have been found to directly modulate chromatin architecture, thereby impacting gene expression and disease progression by interacting with DNA, RNA, and proteins in a tissue- or condition-specific manner. Furthermore, recent studies have highlighted the intricate epigenetic control of lncRNAs in cancer. Consequently, this provides a rationale to explore the possibility of therapeutically targeting lncRNAs themselves or the epigenetic mechanisms that govern their activity. Within the scope of this review, we will assess the current state of knowledge regarding the epigenetic regulation of lncRNAs and how, in turn, lncRNAs contribute to chromatin remodelling in the context of multiple myeloma.

The chromatin - triple helix connection

Mammalian genomes are extensively transcribed, producing a large number of coding and non-coding transcripts. A large fraction of the nuclear RNAs is physically associated with chromatin, functioning in gene activation and silencing, shaping higher-order genome organisation, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions. Different mechanisms allow the tethering of these chromatin-associated RNAs (caRNA) to chromosomes, including RNA binding proteins, the RNA polymerases and R-loops. In this review, we focus on the sequence-specific targeting of RNA to DNA by forming triple helical structures and describe its interplay with chromatin. It turns out that nucleosome positioning at triple helix target sites and the nucleosome itself are essential factors in determining the formation and stability of triple helices. The histone H3-tail plays a critical role in triple helix stabilisation, and the role of its epigenetic modifications in this process is discussed.

Exploring the Roles of RNAs in Chromatin Architecture Using Deep Learning

Recent studies have highlighted the impact of both transcription and transcripts on 3D genome organization, particularly its dynamics. Here, we propose a deep learning framework, called AkitaR, that leverages both genome sequences and genome-wide RNA-DNA interactions to investigate the roles of chromatin-associated RNAs (caRNAs) on genome folding in HFFc6 cells. In order to disentangle the cis- and trans-regulatory roles of caRNAs, we compared models with nascent transcripts, trans-located caRNAs, open chromatin data, or DNA sequence alone. Both nascent transcripts and trans-located caRNAs improved the models' predictions, especially at cell-type-specific genomic regions. Analyses of feature importance scores revealed the contribution of caRNAs at TAD boundaries, chromatin loops and nuclear sub-structures such as nuclear speckles and nucleoli to the models' predictions. Furthermore, we identified non-coding RNAs (ncRNAs) known to regulate chromatin structures, such as MALAT1 and NEAT1, as well as several novel RNAs, RNY5, RPPH1, POLG-DT and THBS1-IT, that might modulate chromatin architecture through trans-interactions in HFFc6. Our modeling also suggests that transcripts from Alus and other repetitive elements may facilitate chromatin interactions through trans R-loop formation. Our findings provide new insights and generate testable hypotheses about the roles of caRNAs in shaping chromatin organization.

Genome-wide analysis of the interplay between chromatin-associated RNA and 3D genome organization in human cells

The interphase genome is dynamically organized in the nucleus and decorated with chromatin-associated RNA (caRNA). It remains unclear whether the genome architecture modulates the spatial distribution of caRNA and vice versa. Here, we generate a resource of genome-wide RNA-DNA and DNA-DNA contact maps in human cells. These maps reveal the chromosomal domains demarcated by locally transcribed RNA, hereafter termed RNA-defined chromosomal domains. Further, the spreading of caRNA is constrained by the boundaries of topologically associating domains (TADs), demonstrating the role of the 3D genome structure in modulating the spatial distribution of RNA. Conversely, stopping transcription or acute depletion of RNA induces thousands of chromatin loops genome-wide. Activation or suppression of the transcription of specific genes suppresses or creates chromatin loops straddling these genes. Deletion of a specific caRNA-producing genomic sequence promotes chromatin loops that straddle the interchromosomal target sequences of this caRNA. These data suggest a feedback loop where the 3D genome modulates the spatial distribution of RNA, which in turn affects the dynamic 3D genome organization.

Global mapping of RNA-chromatin contacts reveals a proximity-dominated connectivity model for ncRNA-gene interactions

Non-coding RNAs (ncRNAs) are transcribed throughout the genome and provide regulatory inputs to gene expression through their interaction with chromatin. Yet, the genomic targets and functions of most ncRNAs are unknown. Here we use chromatin-associated RNA sequencing (ChAR-seq) to map the global network of ncRNA interactions with chromatin in human embryonic stem cells and the dynamic changes in interactions during differentiation into definitive endoderm. We uncover general principles governing the organization of the RNA-chromatin interactome, demonstrating that nearly all ncRNAs exclusively interact with genes in close three-dimensional proximity to their locus and provide a model predicting the interactome. We uncover RNAs that interact with many loci across the genome and unveil thousands of unannotated RNAs that dynamically interact with chromatin. By relating the dynamics of the interactome to changes in gene expression, we demonstrate that activation or repression of individual genes is unlikely to be controlled by a single ncRNA.

Noncoding RNA-chromatin association: Functions and mechanisms

Pervasive transcription of the mammalian genome produces hundreds of thousands of noncoding RNAs (ncRNAs). Numerous studies have suggested that some of these ncRNAs regulate multiple cellular processes and play important roles in physiological and pathological processes. Notably, a large subset of ncRNAs is enriched on chromatin and participates in regulating gene expression and the dynamics of chromatin structure and status. In this review, we summarize recent advances in the functional study of chromatin-associated ncRNAs and mechanistic insights into how these ncRNAs associate with chromatin. We also discuss the potential future challenges which still need to be overcome in this field.

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.

In-depth understanding of higher-order genome architecture in orphan cancer

The human genome is intertwined, folded, condensed, and gradually constitutes the 3D architecture, thereby affecting transcription and widely involving in tumorigenesis. Incidence and mortality rates for orphan cancers increase due to poor early diagnosis and lack of effective medical treatments, which are now getting attention. In-depth understanding in tumorigenesis has fast-tracked over the last decade, however, the further role and mechanism of 3D genome organization in variant orphan tumorigenesis remains to be fully understood. We summarize for the first time that higher-order genome organization can provide novel insights into the occurrence mechanisms of orphan cancers, and discuss probable future research directions for drug development and anti-tumor therapies.

Genome-wide analysis of RNA-chromatin interactions in lizards as a mean for functional lncRNA identification

BACKGROUND

Long non-coding RNAs (lncRNAs) are defined as transcribed molecules longer than 200 nucleotides with little to no protein-coding potential. LncRNAs can regulate gene expression of nearby genes (cis-acting) or genes located on other chromosomes (trans-acting). Several methodologies have been developed to capture lncRNAs associated with chromatin at a genome-wide level. Analysis of RNA-DNA contacts can be combined with epigenetic and RNA-seq data to define potential lncRNAs involved in the regulation of gene expression.

RESULTS

We performed Chromatin Associated RNA sequencing (ChAR-seq) in Anolis carolinensis to obtain the genome-wide map of the associations that RNA molecules have with chromatin. We analyzed the frequency of DNA contacts for different classes of RNAs and were able to define cis- and trans-acting lncRNAs. We integrated the ChAR-seq map of RNA-DNA contacts with epigenetic data for the acetylation of lysine 16 on histone H4 (H4K16ac), a mark connected to actively transcribed chromatin in lizards. We successfully identified three trans-acting lncRNAs significantly associated with the H4K16ac signal, which are likely involved in the regulation of gene expression in A. carolinensis.

CONCLUSIONS

We show that the ChAR-seq method is a powerful tool to explore the RNA-DNA map of interactions. Moreover, in combination with epigenetic data, ChAR-seq can be applied in non-model species to establish potential roles for predicted lncRNAs that lack functional annotations.

The Talented LncRNAs: Meshing into Transcriptional Regulatory Networks in Cancer

As a group of diseases characterized by uncontrollable cell growth, cancer is highly multifaceted in how it overrides checkpoints controlling proliferation. Amongst the regulators of these checkpoints, long non-coding RNAs (lncRNAs) can have key roles in why natural biological processes go haywire. LncRNAs represent a large class of regulatory transcripts that can localize anywhere in cells. They were found to affect gene expression on many levels from transcription to mRNA translation and even protein stability. LncRNA participation in such control mechanisms can depend on cell context, with given transcripts sometimes acting as oncogenes or tumor suppressors. Importantly, the tissue-specificity and low expression levels of lncRNAs make them attractive therapeutic targets or biomarkers. Here, we review the various cellular processes affected by lncRNAs and outline molecular strategies they use to control gene expression, particularly in cancer and in relation to transcription factors.

3plex enables deep computational investigation of triplex forming lncRNAs

Long non-coding RNAs (lncRNAs) regulate gene expression through different molecular mechanisms, including DNA binding via the formation of RNA:DNA:DNA triple helices (TPXs). Despite the increasing amount of experimental evidence, TPXs investigation remains challenging. Here we present 3plex, a software able to predict TPX interactions in silico. Given an RNA sequence and a set of DNA sequences, 3plex integrates 1) Hoogsteen pairing rules that describe the biochemical interactions between RNA and DNA nucleotides, 2) RNA secondary structure prediction and 3) determination of the TPX thermal stability derived from a collection of TPX experimental evidences. We systematically collected and uniformly re-analysed published experimental lncRNA binding sites on human and mouse genomes. We used these data to evaluate 3plex performance and showed that its specific features allow a reliable identification of TPX interactions. We compared 3plex with the other available software and obtained comparable or even better accuracy at a fraction of the computation time. Interestingly, by inspecting collected data with 3plex we found that TPXs tend to be shorter and more degenerated than previously expected and that the majority of analysed lncRNAs can directly bind to the genome by TPX formation. Those results suggest that an important fraction of lncRNAs can exert its biological function through this mechanism. The software is available at https://github.com/molinerisLab/3plex.

RNA-DNA interactomes of three prokaryotes uncovered by proximity ligation

Proximity ligation approaches, which are widely used to study the spatial organization of the genome, also make it possible to reveal patterns of RNA-DNA interactions. Here, we use RedC, an RNA-DNA proximity ligation approach, to assess the distribution of major RNA types along the genomes of E. coli, B. subtilis, and thermophilic archaeon T. adornatum. We find that (i) messenger RNAs preferentially interact with their cognate genes and the genes located downstream in the same operon, which is consistent with polycistronic transcription; (ii) ribosomal RNAs preferentially interact with active protein-coding genes in both bacteria and archaea, indicating co-transcriptional translation; and (iii) 6S noncoding RNA, a negative regulator of bacterial transcription, is depleted from active genes in E. coli and B. subtilis. We conclude that the RedC data provide a rich resource for studying both transcription dynamics and the function of noncoding RNAs in microbial organisms.

RNA-Chrom: a manually curated analytical database of RNA-chromatin interactome

Every year there is more and more evidence that non-coding RNAs play an important role in biological processes affecting various levels of organization of living systems: from the cellular (regulation of gene expression, remodeling and maintenance of chromatin structure, co-transcriptional suppression of transposons, splicing, post-transcriptional RNA modifications, etc.) to cell populations and even organismal ones (development, aging, cancer, cardiovascular and many other diseases). The development and creation of mutually complementary databases that will aggregate, unify and structure different types of data can help to reach the system level of studying non-coding RNAs. Here we present the RNA-Chrom manually curated analytical database, which contains the coordinates of billions of contacts of thousands of human and mouse RNAs with chromatin. Through the user-friendly web interface (https://rnachrom2.bioinf.fbb.msu.ru/), two approaches to the analysis of the RNA-chromatin interactome were implemented. Firstly, to find out whether the RNA of interest to a user contacts with chromatin, and if so, with which genes or DNA loci? Secondly, to find out which RNAs are in contact with the DNA locus of interest to a user (and probably participate in its regulation), and if there are such, what is the nature of their interaction? For a more detailed study of contact maps and their comparison with other data, the web interface allows a user to view them in the UCSC Genome Browser. Database URL https://genome.ucsc.edu/.

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.

Two Novel lncRNAs Regulate Primordial Germ Cell Development in Zebrafish

Long noncoding RNAs (lncRNAs) are regulatory transcripts in various biological processes. However, the role of lncRNAs in germline development remains poorly understood, especially for fish primordial germ cell (PGC) development. In this study, the lncRNA profile of zebrafish PGC was revealed by single cell RNA-sequencing and bioinformatic prediction. We established the regulation network of lncRNA-mRNA associated with PGC development, from which we identified three novel lncRNAs-lnc172, lnc196, and lnc304-highly expressing in PGCs and gonads. Fluorescent in situ hybridization indicated germline-specific localization of lnc196 and lnc304 in the cytoplasm and nucleus of spermatogonia, spermatocyte, and occyte, and they were co-localized with vasa in the cytoplasm of the spermatogonia. By contrast, lnc172 was localized in the cytoplasm of male germline, myoid cells and ovarian somatic cells. Loss- and gain-of-function experiments demonstrated that knockdown and PGC-specific overexpression of lnc304 as well as universal overexpression of lnc172 significantly disrupted PGC development. In summary, the present study revealed the lncRNA profile of zebrafish PGC and identified two novel lncRNAs associated with PGC development, providing new insights for understanding the regulatory mechanism of PGC development.

Recent Development in Biomedical Applications of Oligonucleotides with Triplex-Forming Ability

A DNA structure, known as triple-stranded DNA, is made up of three oligonucleotide chains that wind around one another to form a triple helix (TFO). Hoogsteen base pairing describes how triple-stranded DNA may be built at certain conditions by the attachment of the third strand to an RNA, PNA, or DNA, which might all be employed as oligonucleotide chains. In each of these situations, the oligonucleotides can be employed as an anchor, in conjunction with a specific bioactive chemical, or as a messenger that enables switching between transcription and replication through the triplex-forming zone. These data are also considered since various illnesses have been linked to the expansion of triplex-prone sequences. In light of metabolic acidosis and associated symptoms, some consideration is given to the impact of several low-molecular-weight compounds, including pH on triplex production in vivo. The review is focused on the development of biomedical oligonucleotides with triplexes.

The chromatin-associated RNAs in gene regulation and cancer

Eukaryotic genomes are prevalently transcribed into many types of RNAs that translate into proteins or execute gene regulatory functions. Many RNAs associate with chromatin directly or indirectly and are called chromatin-associated RNAs (caRNAs). To date, caRNAs have been found to be involved in gene and transcriptional regulation through multiple mechanisms and have important roles in different types of cancers. In this review, we first present different categories of caRNAs and the modes of interaction between caRNAs and chromatin. We then detail the mechanisms of chromatin-associated nascent RNAs, chromatin-associated noncoding RNAs and emerging m⁶A on caRNAs in transcription and gene regulation. Finally, we discuss the roles of caRNAs in cancer as well as epigenetic and epitranscriptomic mechanisms contributing to cancer, which could provide insights into the relationship between different caRNAs and cancer, as well as tumor treatment and intervention.

Computational Methods to Study DNA:DNA:RNA Triplex Formation by lncRNAs

Long non-coding RNAs (lncRNAs) impact cell function via numerous mechanisms. In the nucleus, interactions between lncRNAs and DNA and the consequent formation of non-canonical nucleic acid structures seems to be particularly relevant. Along with interactions between single-stranded RNA (ssRNA) and single-stranded DNA (ssDNA), such as R-loops, ssRNA can also interact with double-stranded DNA (dsDNA) to form DNA:DNA:RNA triplexes. A major challenge in the study of DNA:DNA:RNA triplexes is the identification of the precise RNA component interacting with specific regions of the dsDNA. As this is a crucial step towards understanding lncRNA function, there exist several computational methods designed to predict these sequences. This review summarises the recent progress in the prediction of triplex formation and highlights important DNA:DNA:RNA triplexes. In particular, different prediction tools (Triplexator, LongTarget, TRIPLEXES, Triplex Domain Finder, TriplexFFP, TriplexAligner and Fasim-LongTarget) will be discussed and their use exemplified by selected lncRNAs, whose DNA:DNA:RNA triplex forming potential was validated experimentally. Collectively, these tools revealed that DNA:DNA:RNA triplexes are likely to be numerous and make important contributions to gene expression regulation.

Non-coding RNAs in human health and disease: potential function as biomarkers and therapeutic targets

Human diseases have been a critical threat from the beginning of human history. Knowing the origin, course of action and treatment of any disease state is essential. A microscopic approach to the molecular field is a more coherent and accurate way to explore the mechanism, progression, and therapy with the introduction and evolution of technology than a macroscopic approach. Non-coding RNAs (ncRNAs) play increasingly important roles in detecting, developing, and treating all abnormalities related to physiology, pathology, genetics, epigenetics, cancer, and developmental diseases. Noncoding RNAs are becoming increasingly crucial as powerful, multipurpose regulators of all biological processes. Parallel to this, a rising amount of scientific information has revealed links between abnormal noncoding RNA expression and human disorders. Numerous non-coding transcripts with unknown functions have been found in addition to advancements in RNA-sequencing methods. Non-coding linear RNAs come in a variety of forms, including circular RNAs with a continuous closed loop (circRNA), long non-coding RNAs (lncRNA), and microRNAs (miRNA). This comprises specific information on their biogenesis, mode of action, physiological function, and significance concerning disease (such as cancer or cardiovascular diseases and others). This study review focuses on non-coding RNA as specific biomarkers and novel therapeutic targets.

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

Methods to Study RNA-Chromatin Interactions

RNA plays a fundamental role in the organization of chromatin as well as the regulation of gene expression. Although the chromatin is pervasively attached by both coding and noncoding RNAs, the impact of these chromatin-associated RNAs (caRNAs) on gene expression and cellular functions and their underlying mechanisms have just begun to be unraveled. One approach to understand the potential mechanism of gene regulation by caRNAs is to identify the caRNA-associated genomic regions. Several groups have developed methods to capture RNA-chromatin interactions in either one RNA vs the whole genome, i.e., "one-to-all" or all RNAs vs the whole genome, i.e., "all-to-all" manner. In this chapter, we discuss several state-of-the-art methods highlighting the principles behind them, the experimental procedures, the advantages and limitations, and their applications. Our goal is to provide an overview and guide to researchers interested in exploring caRNAs using these techniques.

Novel biological insights revealed from the investigation of multiscale genome architecture

Gene expression and cell fate determination require precise and coordinated epigenetic regulation. The complex three-dimensional (3D) genome organization plays a critical role in transcription in myriad biological processes. A wide range of architectural features of the 3D genome, including chromatin loops, topologically associated domains (TADs), chromatin compartments, and phase separation, together regulate the chromatin state and transcriptional activity at multiple levels. With the help of 3D genome informatics, recent biochemistry and imaging approaches based on different strategies have revealed functional interactions among biomacromolecules, even at the single-cell level. Here, we review the occurrence, mechanistic basis, and functional implications of dynamic genome organization, and outline recent experimental and computational approaches for profiling multiscale genome architecture to provide robust tools for studying the 3D genome.

NPInter v5.0: ncRNA interaction database in a new era

Noncoding RNAs (ncRNAs) play key regulatory roles in biological processes by interacting with other biomolecules. With the development of high-throughput sequencing and experimental technologies, extensive ncRNA interactions have been accumulated. Therefore, we updated the NPInter database to a fifth version to document these interactions. ncRNA interaction entries were doubled from 1 100 618 to 2 596 695 by manual literature mining and high-throughput data processing. We integrated global RNA-DNA interactions from iMARGI, ChAR-seq and GRID-seq, greatly expanding the number of RNA-DNA interactions (from 888 915 to 8 329 382). In addition, we collected different types of RNA interaction between SARS-CoV-2 virus and its host from recently published studies. Long noncoding RNA (lncRNA) expression specificity in different cell types from tumor single cell RNA-seq (scRNA-seq) data were also integrated to provide a cell-type level view of interactions. A new module named RBP was built to display the interactions of RNA-binding proteins with annotations of localization, binding domains and functions. In conclusion, NPInter v5.0 (http://bigdata.ibp.ac.cn/npinter5/) provides informative and valuable ncRNA interactions for biological researchers.

Profiling RNA at chromatin targets in situ by antibody-targeted tagmentation

Whereas techniques to map chromatin-bound proteins are well developed, mapping chromatin-associated RNAs remains a challenge. Here, we describe Reverse Transcribe and Tagment (RT&Tag), in which RNAs associated with a chromatin epitope are targeted by an antibody followed by a protein A-Tn5 transposome. Localized reverse transcription generates RNA/cDNA hybrids that are subsequently tagmented by Tn5 transposases for downstream sequencing. We demonstrate the utility of RT&Tag in Drosophila cells for capturing the noncoding RNA roX2 with the dosage compensation complex and maturing transcripts associated with silencing histone modifications. We also show that RT&Tag can detect N6-methyladenosine-modified mRNAs, and show that genes producing methylated transcripts are characterized by extensive promoter pausing of RNA polymerase II. The high efficiency of in situ antibody tethering and tagmentation makes RT&Tag especially suitable for rapid low-cost profiling of chromatin-associated RNAs.

A universal model of RNA.DNA:DNA triplex formation accurately predicts genome-wide RNA-DNA interactions

RNA.DNA:DNA triple helix (triplex) formation is a form of RNA-DNA interaction which regulates gene expression but is difficult to study experimentally in vivo. This makes accurate computational prediction of such interactions highly important in the field of RNA research. Current predictive methods use canonical Hoogsteen base pairing rules, which whilst biophysically valid, may not reflect the plastic nature of cell biology. Here, we present the first optimization approach to learn a probabilistic model describing RNA-DNA interactions directly from motifs derived from triplex sequencing data. We find that there are several stable interaction codes, including Hoogsteen base pairing and novel RNA-DNA base pairings, which agree with in vitro measurements. We implemented these findings in TriplexAligner, a program that uses the determined interaction codes to predict triplex binding. TriplexAligner predicts RNA-DNA interactions identified in all-to-all sequencing data more accurately than all previously published tools in human and mouse and also predicts previously studied triplex interactions with known regulatory functions. We further validated a novel triplex interaction using biophysical experiments. Our work is an important step towards better understanding of triplex formation and allows genome-wide analyses of RNA-DNA interactions.

Detecting TAD-like domains from RNA-associated interactions

Topologically associated domains (TADs) are crucial chromatin structural units. Evidence has illustrated that RNA-chromatin and RNA-RNA spatial interactions, so-called RNA-associated interactions (RAIs), may be associated with TAD-like domains (TLDs). To decode hierarchical TLDs from RAIs, we proposed SuperTLD, a domain detection algorithm incorporating imputation. We applied SuperTLD on four RAI data sets and compared TLDs with the TADs identified from the corresponding Hi-C datasets. The TLDs and TADs share a moderate similarity of hierarchies ≥ 0.5312 and the finest structures ≥ 0.8295. Comparison between boundaries and domains further demonstrated the novelty of TLDs. Enrichment analysis of epigenetic characteristics illustrated that the novel TLDs exhibit an enriched CTCF by 0.6245 fold change and H3 histone marks enriched within domains. GO analysis on the TLD novel boundaries exhibited enriched diverse terms, revealing TLDs' formation mechanism related closely to gene regulation.

Chromatin-Based Transcriptional Reprogramming in Plants under Abiotic Stresses

Plants' stress response machinery is characterized by an intricate network of signaling cascades that receive and transmit environmental cues and ultimately trigger transcriptional reprogramming. The family of epigenetic regulators that are the key players in the stress-induced signaling cascade comprise of chromatin remodelers, histone modifiers, DNA modifiers and regulatory non-coding RNAs. Changes in the histone modification and DNA methylation lead to major alterations in the expression level and pattern of stress-responsive genes to adjust with abiotic stress conditions namely heat, cold, drought and salinity. The spotlight of this review falls primarily on the chromatin restructuring under severe abiotic stresses, crosstalk between epigenetic regulators along with a brief discussion on stress priming in plants.

Enhancer RNAs (eRNAs) in Cancer: The Jacks of All Trades

Simple Summary

This review focuses on eRNAs and the several mechanisms by which they can regulate gene expression. In particular we describe here the most recent examples of eRNAs dysregulated in cancer or involved in the immune escape of tumor cells.

Abstract

Enhancer RNAs (eRNAs) are non-coding RNAs (ncRNAs) transcribed in enhancer regions. They play an important role in transcriptional regulation, mainly during cellular differentiation. eRNAs are tightly tissue- and cell-type specific and are induced by specific stimuli, activating promoters of target genes in turn. eRNAs usually have a very short half-life but in some cases, once activated, they can be stably expressed and acquire additional functions. Due to their critical role, eRNAs are often dysregulated in cancer and growing number of interactions with chromatin modifiers, transcription factors, and splicing machinery have been described. Enhancer activation and eRNA transcription have particular relevance also in inflammatory response, placing the eRNAs at the interplay between cancer and immune cells. Here, we summarize all the possible molecular mechanisms recently reported in association with eRNAs activity.

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.

HiMoRNA: A Comprehensive Database of Human lncRNAs Involved in Genome-Wide Epigenetic Regulation

Long non-coding RNAs (lncRNAs) play an important role in genome regulation. Specifically, many lncRNAs interact with chromatin, recruit epigenetic complexes and in this way affect large-scale gene expression programs. However, the experimental data about lncRNA-chromatin interactions is still limited. The majority of experimental protocols do not provide any insight into the mechanics of lncRNA-based genome-wide epigenetic regulation. Here we present the HiMoRNA (Histone-Modifying RNA) database, a resource containing correlated lncRNA–epigenetic changes in specific genomic locations genome-wide. HiMoRNA integrates a large amount of multi-omics data to characterize the effects of lncRNA on epigenetic modifications and gene expression. The current release of HiMoRNA includes more than five million associations in humans for ten histone modifications in multiple genomic loci and 4145 lncRNAs. HiMoRNA provides a user-friendly interface to facilitate browsing, searching and retrieving of lncRNAs associated with epigenetic profiles of various chromatin loci. Analysis of the HiMoRNA data suggests that several lncRNA including JPX might be involved not only in regulation of XIST locus but also in direct establishment or maintenance of X-chromosome inactivation. We believe that HiMoRNA is a convenient and valuable resource that can provide valuable biological insights and greatly facilitate functional annotation of lncRNAs.

The landscape of promoter-centred RNA-DNA interactions in rice

Chromatin-associated RNAs play key roles in various biological processes. However, both their repository and conjugation genomic loci and potential functions remain largely unclear. Here, we develop an effective method for mapping of chromatin-associated RNA-DNA interactions, followed by paired-end-tag sequencing (ChRD-PET) in rice. We present a comprehensive interaction map between RNAs and H3K4me3-marked regions based on H3K4me3 ChRD-PET data, showing three types of RNA-DNA interactions-local, proximal and distal. We further characterize the origin and composition of the RNA strand in R-loop RNA-DNA hybrids and identify that extensive cis and trans RNAs, including trans-non-coding RNAs, are prevalently involved in the R-loop. Integrative analysis of rice epigenome and three-dimensional genome data suggests that both coding and non-coding RNAs engage extensively in the formation of chromatin loops and chromatin-interacting domains. In summary, ChRD-PET is an efficient method for studying the features of RNA-chromatin interactions, and the resulting datasets constitute a valuable resource for the study of RNAs and their biological functions.