Systematic Mapping of RNA-Chromatin Interactions In Vivo

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.

RNA in chromatin organization and nuclear architecture

In the cell nucleus, genomic DNA is surrounded by nonmembranous nuclear bodies. This might result from specific regions of the genome being transcribed into long noncoding RNAs (lncRNAs), which tend to remain at the sites of their own transcription. The lncRNAs seed the nuclear bodies by recruiting and concentrating proteins and RNAs, which undergo liquid-liquid-phase separation, and form molecular condensates, the so-called nuclear bodies. These nuclear bodies may provide appropriate environments for gene activation or repression. Notably, lncRNAs also contribute to three-dimensional genome structure by mediating long-range chromatin interactions. In this review, we discuss the mechanisms by which lncRNAs regulate gene expression through shaping chromatin and nuclear architectures. We also explore lncRNAs' potential as a therapeutic target for cancer, because lncRNAs are often expressed in a disease-specific manner.

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.

Long noncoding RNA study: Genome-wide approaches

Long noncoding RNAs (lncRNAs) have been confirmed to play a crucial role in various biological processes across several species. Though many efforts have been devoted to the expansion of the lncRNAs landscape, much about lncRNAs is still unknown due to their great complexity. The development of high-throughput technologies and the constantly improved bioinformatic methods have resulted in a rapid expansion of lncRNA research and relevant databases. In this review, we introduced genome-wide research of lncRNAs in three parts: (i) novel lncRNA identification by high-throughput sequencing and computational pipelines; (ii) functional characterization of lncRNAs by expression atlas profiling, genome-scale screening, and the research of cancer-related lncRNAs; (iii) mechanism research by large-scale experimental technologies and computational analysis. Besides, primary experimental methods and bioinformatic pipelines related to these three parts are summarized. This review aimed to provide a comprehensive and systemic overview of lncRNA genome-wide research strategies and indicate a genome-wide lncRNA research system.

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.

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.

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.

Crosstalk between RNA m6A and DNA methylation regulates transposable element chromatin activation and cell fate in human pluripotent stem cells

Transposable elements (TEs) are parasitic DNA sequences accounting for over half of the human genome. Tight control of the repression and activation states of TEs is critical for genome integrity, development, immunity and diseases, including cancer. However, precisely how this regulation is achieved remains unclear. Here we develop a targeted proteomic proximity labeling approach to capture TE-associated proteins in human embryonic stem cells (hESCs). We find that the RNA N6-methyladenosine (m6A) reader, YTHDC2, occupies genomic loci of the primate-specific TE, LTR7/HERV-H, specifically through its interaction with m6A-modified HERV-H RNAs. Unexpectedly, YTHDC2 recruits the DNA 5-methylcytosine (5mC)-demethylase, TET1, to remove 5mC from LTR7/HERV-H and prevent epigenetic silencing. Functionally, the YTHDC2/LTR7 axis inhibits neural differentiation of hESCs. Our results reveal both an underappreciated crosstalk between RNA m6A and DNA 5mC, the most abundant regulatory modifications of RNA and DNA in eukaryotes, and the fact that in hESCs this interplay controls TE activity and cell fate.

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

Differential nuclear import regulates nuclear RNA inheritance following mitosis

Mitosis results in a dramatic reorganization of chromatin structure to promote chromosome compaction and segregation to daughter cells. Consequently, mitotic entry is accompanied by transcriptional silencing and removal of most chromatin-bound RNA from chromosomes. As cells exit mitosis, chromatin rapidly decondenses and transcription restarts as waves of differential gene expression. However, little is known about the fate of chromatin-bound RNAs following cell division. Here we explored whether nuclear RNA from the previous cell cycle is present in G1 nuclei following mitosis. We found that half of all nuclear RNA is inherited in a transcription-independent manner following mitosis. Interestingly, the snRNA U2 is efficiently inherited by G1 nuclei, while the lncRNAs NEAT1 and MALAT1 show no inheritance following mitosis. We found that the nuclear protein SAF-A, which is hypothesized to tether RNA to DNA, did not play a prominent role in nuclear RNA inheritance, indicating that the mechanism for RNA inheritance may not involve RNA chaperones that have chromatin-binding activity. Instead, we observe that the timing of RNA inheritance indicates that a select group of nuclear RNAs are reimported into the nucleus after the nuclear envelope has reassembled. Our work demonstrates that there is a fraction of nuclear RNA from the previous cell cycle that is reimported following mitosis and suggests that mitosis may serve as a time to reset the interaction of lncRNAs with chromatin.

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.

RNA-mediated heterochromatin formation at repetitive elements in mammals

The failure to repress transcription of repetitive genomic elements can lead to catastrophic genome instability and is associated with various human diseases. As such, multiple parallel mechanisms cooperate to ensure repression and heterochromatinization of these elements, especially during germline development and early embryogenesis. A vital question in the field is how specificity in establishing heterochromatin at repetitive elements is achieved. Apart from trans-acting protein factors, recent evidence points to a role of different RNA species in targeting repressive histone marks and DNA methylation to these sites in mammals. Here, we review recent discoveries on this topic and predominantly focus on the role of RNA methylation, piRNAs, and other localized satellite RNAs.

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.

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.

MALAT1 modulates alternative splicing by cooperating with the splicing factors PTBP1 and PSF

Understanding how long noncoding RNAs (lncRNAs) cooperate with splicing factors (SFs) in alternative splicing (AS) control is fundamental to human biology and disease. We show that metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), a well-documented AS-implicated lncRNA, regulates AS via two SFs, polypyrimidine tract-binding protein 1 (PTBP1) and PTB-associated SF (PSF). MALAT1 stabilizes the interaction between PTBP1 and PSF, thereby forming a functional module that affects a network of AS events. The MALAT1-stabilized PTBP1/PSF interaction occurs in multiple cellular contexts; however, the functional module, relative to MALAT1 only, has more dominant pathological significance in hepatocellular carcinoma. MALAT1 also stabilizes the PSF interaction with several heterogeneous nuclear ribonucleoparticle proteins other than PTBP1, hinting a broad role in AS control. We present a model in which MALAT1 cooperates with distinct SFs for AS regulation and pose that, relative to analyses exclusively performed for lncRNAs, a comprehensive consideration of lncRNAs and their binding partners may provide more information about their biological functions.

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.

Triblock copolymer micelle model of spherical paraspeckles

Paraspeckles are nuclear bodies scaffolded by RNP complexes of NEAT1_2 RNA transcripts and multiple RNA-binding proteins. The assembly of paraspeckles is coupled with the transcription of NEAT1_2. Paraspeckles form the core-shell structure, where the two terminal regions of NEAT1_2 RNP complexes compose the shell of the paraspeckle and the middle regions of these complexes compose the core. We here construct a theoretical model of paraspeckles by taking into account the transcription of NEAT1_2 in an extension of the theory of block copolymer micelles. This theory predicts that the core-shell structure of a paraspeckle is assembled by the association of the middle region of NEAT1_2 RNP complexes due to the multivalent interactions between RBPs bound to these regions and by the relative affinity of the terminal regions of the complexes to the nucleoplasm. The latter affinity results in the effective repulsive interactions between terminal regions of the RNA complexes and limits the number of complexes composing the paraspeckle. In the wild type, the repulsive interaction between the middle and terminal block dominates the thermal fluctuation. However, the thermal fluctuation can be significant in a mutant, where a part of the terminal regions of NEAT1_2 is deleted, and distributes the shortened terminal regions randomly between the shell and the core, consistent with our recent experiments. With the upregulated transcription, the shortened terminal regions of NEAT1_2 in a deletion mutant is localized to the core to decrease the repulsive interaction between the terminal regions, while the structure does not change with the upregulation in the wild type. The robustness of the structure of paraspeckles in the wild type results from the polymeric nature of NEAT1_2 complexes.

Epithelial cells-enriched lncRNA SNHG8 regulates chromatin condensation by binding to Histone H1s

Linker histone H1 proteins contain many variants in mammalian and can stabilize the condensed state of chromatin by binding to nucleosomes and promoting a more inaccessible structure of DNA. However, it is poorly understood how the binding of histone H1s to chromatin DNA is regulated. Screened as one of a collection of epithelial cells-enriched long non-coding RNAs (lncRNAs), here we found that small nucleolar RNA host gene 8 (SNHG8) is a chromatin-localized lncRNA and presents strong interaction and phase separation with histone H1 variants. Moreover, SNHG8 presents stronger ability to bind H1s than linker DNA, and outcompetes linker DNA for H1 binding. Consequently, loss of SNHG8 increases the amount of H1s that bind to chromatin, promotes chromatin condensation, and induces an epithelial differentiation-associated gene expression pattern. Collectively, our results propose that the highly abundant SNHG8 in epithelial cells keeps histone H1 variants out of nucleosome and its loss contributes to epithelial cell differentiation.

Nuclear speckles - a driving force in gene expression

Nuclear speckles are dynamic membraneless bodies located in the cell nucleus. They harbor RNAs and proteins, many of which are splicing factors, that together display complex biophysical properties dictating nuclear speckle formation and maintenance. Although these nuclear bodies were discovered decades ago, only recently has in-depth genomic analysis begun to unravel their essential functions in modulation of gene activity. Major advancements in genomic mapping techniques combined with microscopy approaches have enabled insights into the roles nuclear speckles may play in enhancing gene expression, and how gene positioning to specific nuclear landmarks can regulate gene expression and RNA processing. Some studies have drawn a link between nuclear speckles and disease. Certain maladies either involve nuclear speckles directly or dictate the localization and reorganization of many nuclear speckle factors. This is most striking during viral infection, as viruses alter the entire nuclear architecture and highjack host machinery. As discussed in this Review, nuclear speckles represent a fascinating target of study not only to reveal the links between gene positioning, genome subcompartments and gene activity, but also as a potential target for therapeutics.

The era of 3D and spatial genomics

Over a decade ago the advent of high-throughput chromosome conformation capture (Hi-C) sparked a new era of 3D genomics. Since then the number of methods for mapping the 3D genome has flourished, enabling an ever-increasing understanding of how DNA is packaged in the nucleus and how the spatiotemporal organization of the genome orchestrates its vital functions. More recently, the next generation of spatial genomics technologies has begun to reveal how genome sequence and 3D genome organization vary between cells in their tissue context. We summarize how the toolkit for charting genome topology has evolved over the past decade and discuss how new technological developments are advancing the field of 3D and spatial genomics.

Ribosomal RNA regulates chromosome clustering during mitosis

Noncoding RNAs are known to associate with mitotic chromosomes, but the identities and functions of chromosome-associated RNAs in mitosis remain elusive. Here, we show that rRNA species associate with condensed chromosomes during mitosis. In particular, pre-rRNAs such as 45S, 32S, and 30S are highly enriched on mitotic chromosomes. Immediately following nucleolus disassembly in mitotic prophase, rRNAs are released and associate with and coat each condensed chromosome at prometaphase. Using unbiased mass spectrometry analysis, we further demonstrate that chromosome-bound rRNAs are associated with Ki-67. Moreover, the FHA domain and the repeat region of Ki-67 recognize and anchor rRNAs to chromosomes. Finally, suppression of chromosome-bound rRNAs by RNA polymerase I inhibition or by using rRNA-binding-deficient Ki-67 mutants impair mitotic chromosome dispersion during prometaphase. Our study thus reveals an important role of rRNAs in preventing chromosome clustering during mitosis.

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.

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.

The 4D Nucleome Data Portal as a resource for searching and visualizing curated nucleomics data

The 4D Nucleome (4DN) Network aims to elucidate the complex structure and organization of chromosomes in the nucleus and the impact of their disruption in disease biology. We present the 4DN Data Portal (https://data.4dnucleome.org/), a repository for datasets generated in the 4DN network and relevant external datasets. Datasets were generated with a wide range of experiments, including chromosome conformation capture assays such as Hi-C and other innovative sequencing and microscopy-based assays probing chromosome architecture. All together, the 4DN data portal hosts more than 1800 experiment sets and 36000 files. Results of sequencing-based assays from different laboratories are uniformly processed and quality-controlled. The portal interface allows easy browsing, filtering, and bulk downloads, and the integrated HiGlass genome browser allows interactive visualization and comparison of multiple datasets. The 4DN data portal represents a primary resource for chromosome contact and other nuclear architecture data for the scientific community.

This paper describes the ‘4DN Data Portal’ that hosts data generated by the 4D Nucleome network, including Hi-C and other chromatin conformation capture assays, as well as various sequencing-based and imaging-based assays. Raw data have been uniformly processed to increase comparability and the portal is implemented with visualization tools to browse the data without download.

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.

Pairs and Pairix: a file format and a tool for efficient storage and retrieval for Hi-C read pairs

SUMMARY

As the amount of 3D chromosomal interaction data continues to increase, storing and accessing such data efficiently becomes paramount. We introduce Pairs, a block-compressed text file format for storing paired genomic coordinates from Hi-C data, and Pairix, an open-source C application to index and query Pairs files. Pairix (also available in Python and R) extends the functionalities of Tabix to paired coordinates data. We have also developed PairsQC, a collapsible HTML quality control report generator for Pairs files.

AVAILABILITY AND IMPLEMENTATION

The format specification and source code are available at https://github.com/4dn-dcic/pairix, https://github.com/4dn-dcic/Rpairix and https://github.com/4dn-dcic/pairsqc.

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.

Long Non-coding RNA MALAT1: A Key Player in Liver Diseases

Long non-coding RNAs (lncRNAs) exceed 200 nucleotides in length are considered to be involved in both developmental processes and various diseases. Here, we focus on lncRNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1), which was one of the most important lncRNAs in proliferation, apoptosis, and migration. MALAT1 plays a regulatory role in liver diseases, including hepatic fibrosis, liver regeneration, liver cancer, and fatty liver diseases. In the current review, we summarize the latest literature about the function roles of MALAT1 in liver disorders. Probing the regulatory mechanism and cross talk of MALAT1 with other signaling pathways of pathological processes would improve the prognosis, diagnosis of liver diseases, and offer a promising candidate target for therapeutic interventions.

The regulatory activities of MALAT1 in the development of bone and cartilage diseases

Long non-coding RNAs (lncRNAs) have been comprehensively implicated in various cellular functions by mediating transcriptional or post-transcriptional activities. MALAT1 is involved in the differentiation, proliferation, and apoptosis of multiple cell lines, including BMSCs, osteoblasts, osteoclasts, and chondrocytes. Interestingly, MALAT1 may interact with RNAs or proteins, regulating cellular processes. Recently, MALAT1 has been reported to be associated with the development of bone and cartilage diseases by orchestrating the signaling network. The involvement of MALAT1 in the pathological development of bone and cartilage diseases makes it available to be a potential biomarker for clinical diagnosis or prognosis. Although the potential mechanisms of MALAT1 in mediating the cellular processes of bone and cartilage diseases are still needed for further elucidation, MALAT1 shows great promise for drug development.

Hi-C Data Formats

Processing, storing, and visualizing high-resolution Hi-C data required development of efficient data formats. A sparse matrix format saving only nonzero values has become the norm. A "zoomable" matrix style also became popular, storing multiple resolutions in a single file for interactive visualization. This chapter discusses the latest matrix file formats such as .hic and .mcool, and other intermediate formats including SAM/BAM and random-accessible contact lists.

Deciphering the RNA universe in sperm in its role as a vertical information carrier

The inheritance of neurophysiologic and neuropsychologic complex diseases can only partly be explained by the Mendelian concept of genetic inheritance. Previous research showed that both psychological disorders like post-traumatic stress disorder and metabolic diseases are more prevalent in the progeny of affected parents. This could suggest an epigenetic mode of transmission. Human studies give first insight into the scope of intergenerational influence of stressors but are limited in exploring the underlying mechanisms. Animal models have elucidated the mechanistic underpinnings of epigenetic transmission. In this review, we summarize progress on the mechanisms of paternal intergenerational transmission by means of sperm RNA in mouse models. We discuss relevant details for the modelling of RNA-mediated transmission, point towards currently unanswered questions and propose experimental considerations for tackling these questions.

Long non-coding RNAs as possible therapeutic targets in protozoa, and in Schistosoma and other helminths

Long non-coding RNAs (lncRNAs) emerged in the past 20 years due to massive amounts of scientific data regarding transcriptomic analyses. They have been implicated in a plethora of cellular processes in higher eukaryotes. However, little is known about lncRNA possible involvement in parasitic diseases, with most studies only detecting their presence in parasites of human medical importance. Here, we review the progress on lncRNA studies and their functions in protozoans and helminths. In addition, we show an example of knockdown of one lncRNA in Schistosoma mansoni, SmLINC156349, which led to in vitro parasite adhesion, motility, and pairing impairment, with a 20% decrease in parasite viability and 33% reduction in female oviposition. Other observed phenotypes were a decrease in the proliferation rate of both male and female worms and their gonads, and reduced female lipid and vitelline droplets that are markers for well-developed vitellaria. Impairment of female worms' vitellaria in SmLINC156349-silenced worms led to egg development deficiency. All those results demonstrate the great potential of the tools and methods to characterize lncRNAs as potential new therapeutic targets. Further, we discuss the challenges and limitations of current methods for studying lncRNAs in parasites and possible solutions to overcome them, and we highlight the future directions of this exciting field.

RNA promotes the formation of spatial compartments in the nucleus

RNA, DNA, and protein molecules are highly organized within three-dimensional (3D) structures in the nucleus. Although RNA has been proposed to play a role in nuclear organization, exploring this has been challenging because existing methods cannot measure higher-order RNA and DNA contacts within 3D structures. To address this, we developed RNA & DNA SPRITE (RD-SPRITE) to comprehensively map the spatial organization of RNA and DNA. These maps reveal higher-order RNA-chromatin structures associated with three major classes of nuclear function: RNA processing, heterochromatin assembly, and gene regulation. These data demonstrate that hundreds of ncRNAs form high-concentration territories throughout the nucleus, that specific RNAs are required to recruit various regulators into these territories, and that these RNAs can shape long-range DNA contacts, heterochromatin assembly, and gene expression. These results demonstrate a mechanism where RNAs form high-concentration territories, bind to diffusible regulators, and guide them into compartments to regulate essential nuclear functions.

Single-cell Hi-C data analysis: safety in numbers

Over the past decade, genome-wide assays for chromatin interactions in single cells have enabled the study of individual nuclei at unprecedented resolution and throughput. Current chromosome conformation capture techniques survey contacts for up to tens of thousands of individual cells, improving our understanding of genome function in 3D. However, these methods recover a small fraction of all contacts in single cells, requiring specialised processing of sparse interactome data. In this review, we highlight recent advances in methods for the interpretation of single-cell genomic contacts. After discussing the strengths and limitations of these methods, we outline frontiers for future development in this rapidly moving field.

Three-dimensional genome organization via triplex-forming RNAs

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF.

The macro and micro of chromosome conformation capture

The 3D organization of the genome facilitates gene regulation, replication, and repair, making it a key feature of genomic function and one that remains to be properly understood. Over the past two decades, a variety of chromosome conformation capture (3C) methods have delineated genome folding from megabase-scale compartments and topologically associating domains (TADs) down to kilobase-scale enhancer-promoter interactions. Understanding the functional role of each layer of genome organization is a gateway to understanding cell state, development, and disease. Here, we discuss the evolution of 3C-based technologies for mapping 3D genome organization. We focus on genomics methods and provide a historical account of the development from 3C to Hi-C. We also discuss ChIP-based techniques that focus on 3D genome organization mediated by specific proteins, capture-based methods that focus on particular regions or regulatory elements, 3C-orthogonal methods that do not rely on restriction digestion and proximity ligation, and methods for mapping the DNA-RNA and RNA-RNA interactomes. We consider the biological discoveries that have come from these methods, examine the mechanistic contributions of CTCF, cohesin, and loop extrusion to genomic folding, and detail the 3D genome field's current understanding of nuclear architecture. Finally, we give special consideration to Micro-C as an emerging frontier in chromosome conformation capture and discuss recent Micro-C findings uncovering fine-scale chromatin organization in unprecedented detail. This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

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.