RNA-chromatin Interactions Research Articles

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.

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

BackgroundLong 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.ResultsWe 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.ConclusionsWe 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.

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.

P487: SMALL NON-CODING RNAS ASSOCIATE WITH CHROMATIN TO MAINTEIN PROPOGATION OF ACUTE MEYLOID LEUKEMIA

Background: Acute myeloid leukemia (AML) is a hematopoietic malignancy characterized by the clonal expansion of myeloid progenitors without terminal differentiation. The interplay between genome and epigenome alterations drives aberrant gene network linking to leukemia transformation. Altered transcription may occur at three-dimensional chromatin level and be orchestrated by a combination of multiple regulatory factors. Although some RNA molecules, especially non-coding RNAs, were reported to modulate transcription in leukemia cells, the regulatory role of genome-wide RNAs on transcription through chromatin association remains largely unexplored. Aims: To identify novel RNA molecules as chromatin-association factors critical for leukemogenesis through modulating transcription Methods: We performed in situ Mapping of RNA-Genome Interactome (iMARGI) in a human acute leukemia cell line – MV4-11 cells. Chromatin association of candidate RNAs was verified using Chromatin Isolation by RNA Purification (ChIRP) by means of biotinylated oligonucleotides for target capture. A series of ChIP-seq assays were performed in multiple AML cell lines to profile chromatin occupancy by Fibrillarin (FBL), a catalytic component of C/D box ribonucleoproteins (snoRNPs), as well as to observe the remodelling of chromatin states marked by H3K27ac and H3K9me3 upon shRNA-mediated knockdown of FBL. Loss-of-function experiments using antisense oligonucleotides were performed on candidate chromatin-associated RNAs in AML cell lines, AML primary samples, and CD34+ healthy bone marrow cells. Results: In total, about 52.1 million RNA-chromatin interactions were identified out of deep sequencing of four iMARGI libraries. Whereas the vast majority (50.5 million, 97%) of RNA-chromatin interactome were intra-chromosome, a small fraction (1.6 million, 3%) occurred inter chromosome. Though protein-coding mRNAs predominate the chromatin-associated RNA species, other non-coding RNAs were also captured. Intriguingly, different from protein-coding mRNAs and long intergenic non-coding RNAs (lincRNAs) which mainly form intra-chromosome interactions, two small non-coding RNA classes, small nucleolar RNAs (snoRNAs) and small nuclear RNAs (snRNAs), were found with prevalent frequency (62% and 92%, respectively) of inter-chromosome associations. SNORD3A (U3) and SNORD118 (U8) were C/D box snoRNAs ranked top by the number of their associated interactions in MV4-11 cells as well as in other cell types. Chromatin association of U3 and U8 was verified using ChIRP, with the observation of unique chromatin binding for both U3 and U8, and of specific binding motifs of transcription factors at their binding sites. ChIP-seq revealed that the snoRNP FBL occupied only a small set of genes (n=118), and shRNA-mediated knockdown of FBL had a minimal impact on H3K27ac and H3K9me3, suggesting that snoRNA-chromatin interactions may function independently on FBL-associated snoRNPs. ASO-mediated loss of U3 or U8 significantly impairs leukemia cells proliferation and colony forming capacity in multiple AML cell lines and primary AML blast samples, whereas only mildly affects CD34+ healthy cells. Summary/Conclusion: We identified chromatin-associated snoRNAs U3 and U8 which were essential for leukemia maintenance and may serve as therapeutic targets for leukemia eradication.

Enhancer-promoter interaction maps provide insights into skeletal muscle-related traits in pig genome

BackgroundGene expression programs are intimately linked to the interplay of active cis regulatory elements mediated by chromatin contacts and associated RNAs. Genome-wide association studies (GWAS) have identified many variants in these regulatory elements that can contribute to phenotypic diversity. However, the functional interpretation of these variants remains nontrivial due to the lack of chromatin contact information or limited contact resolution. Furthermore, the distribution and role of chromatin-associated RNAs in gene expression and chromatin conformation remain poorly understood. To address this, we first present a comprehensive interaction map of nuclear dynamics of 3D chromatin-chromatin interactions (H3K27ac BL-HiChIP) and RNA-chromatin interactions (GRID-seq) to reveal genomic variants that contribute to complex skeletal muscle traits.ResultsIn a genome-wide scan, we provide systematic fine mapping and gene prioritization from GWAS leading signals that underlie phenotypic variability of growth rate, meat quality, and carcass performance. A set of candidate functional variants and 54 target genes previously not detected were identified, with 71% of these candidate functional variants choosing to skip over their nearest gene to regulate the target gene in a long-range manner. The effects of three functional variants regulating KLF6 (related to days to 100 kg), MXRA8 (related to lean meat percentage), and TAF11 (related to loin muscle depth) were observed in two pig populations. Moreover, we find that this multi-omics interaction map consists of functional communities that are enriched in specific biological functions, and GWAS target genes can serve as core genes for exploring peripheral trait-relevant genes.ConclusionsOur results provide a valuable resource of candidate functional variants for complex skeletal muscle-related traits and establish an integrated approach to complement existing 3D genomics by exploiting RNA-chromatin and chromatin-chromatin interactions for future association studies.

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.

RedChIP identifies noncoding RNAs associated with genomic sites occupied by Polycomb and CTCF proteins

Nuclear noncoding RNAs (ncRNAs) are key regulators of gene expression and chromatin organization. The progress in studying nuclear ncRNAs depends on the ability to identify the genome-wide spectrum of contacts of ncRNAs with chromatin. To address this question, a panel of RNA–DNA proximity ligation techniques has been developed. However, neither of these techniques examines proteins involved in RNA–chromatin interactions. Here, we introduce RedChIP, a technique combining RNA–DNA proximity ligation and chromatin immunoprecipitation for identifying RNA–chromatin interactions mediated by a particular protein. Using antibodies against architectural protein CTCF and the EZH2 subunit of the Polycomb repressive complex 2, we identify a spectrum of cis- and trans-acting ncRNAs enriched at Polycomb- and CTCF-binding sites in human cells, which may be involved in Polycomb-mediated gene repression and CTCF-dependent chromatin looping. By providing a protein-centric view of RNA–DNA interactions, RedChIP represents an important tool for studies of nuclear ncRNAs.

SAF-A mutants disrupt chromatin structure through dominant negative effects on RNAs associated with chromatin

Here we provide a brief review of relevant background before presenting results of our investigation into the interplay between scaffold attachment factor A (SAF-A), chromatin-associated RNAs, and DNA condensation. SAF-A, also termed heterogenous nuclear protein U (hnRNP U), is a ubiquitous nuclear scaffold protein that was implicated in XIST RNA localization to the inactive X-chromosome (Xi) but also reported to maintain open DNA packaging in euchromatin. Here we use several means to perturb SAF-A and examine potential impacts on the broad association of RNAs on euchromatin, and on chromatin compaction. SAF-A has an N-terminal DNA binding domain and C-terminal RNA binding domain, and a prominent model has been that the protein provides a single-molecule bridge between XIST RNA and chromatin. Here analysis of the impact of SAF-A on broad RNA-chromatin interactions indicate greater biological complexity. We focus on SAF-A’s role with repeat-rich C0T-1 hnRNA (repeat-rich heterogeneous nuclear RNA), shown recently to comprise mostly intronic sequences of pre-mRNAs and diverse long non-coding RNAs (lncRNAs). Our results show that SAF-A mutants cause dramatic changes to cytological chromatin condensation through dominant negative effects on C0T-1 RNA’s association with euchromatin, and likely other nuclear scaffold factors. In contrast, depletion of SAF-A by RNA interference (RNAi) had no discernible impact on C0T-1 RNA, nor did it cause similarly marked chromatin changes as did three different SAF-A mutations. Overall results support the concept that repeat-rich, chromatin-associated RNAs interact with multiple RNA binding proteins (RBPs) in a complex dynamic meshwork that is integral to larger-scale chromatin architecture and collectively influences cytological-scale DNA condensation.

A novel RNA-mediated mechanism causing down-regulation of insulating promoter interactions in human embryonic stem cells

The genome-wide promoter interactome is primarily maintained and regulated by architectural proteins such as CTCF and cohesin. However, some studies suggest a role for non-coding RNAs (ncRNAs) in this process. We aimed to characterise the regulatory role of RNA-mediated promoter interactions in the control of gene expression. We integrated genome-wide datasets of RNA-chromatin and promoter-genome interactions in human embryonic stem cells (hESCs) to identify putative RNA-mediated promoter interactions. We discovered that CTCF sites were enriched in RNA-PIRs (promoter interacting regions co-localising with RNA-chromatin interaction sites) and genes interacting with RNA-PIRs containing CTCF sites showed higher expression levels. One of the long noncoding RNAs (lncRNAs) expressed in hESCs, Syntaxin 18-Antisense 1 (STX18-AS1), appeared to be involved in an insulating promoter interaction with the neighbouring gene, MSX1. By knocking down STX18-AS1, the MSX1 promoter-PIR interaction was intensified and the target gene (MSX1) expression was down-regulated. Conversely, reduced MSX1 promoter-PIR interactions, resulting from CRISPR-Cas9 deletion of the PIR, increased the expression of MSX1. We conclude that STX18-AS1 RNA antagonised local CTCF-mediated insulating promoter interactions to augment gene expression. Such down-regulation of the insulating promoter interactions by this novel mechanism may explain the higher expression of genes interacting with RNA-PIRs linked to CTCF sites.

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.

Analysis of Enhancer-Promoter Interactions using CAGE and RADICL-Seq Technologies.

Regulation of gene expression is a key feature for higher eukaryotes and how chromatin topology relates to gene activation is an intense area of research. Enhancer-promoter interactions are believed to mediate activation of target genes. Bidirectional transcription represents one hallmark of active enhancers that can be measured using transcriptome technologies such as Cap analysis of gene expression (CAGE). Recently, we have developed RNA and DNA interacting complexes ligated and sequenced (RADICL-Seq) a novel methodology to map genome-wide RNA-chromatin interactions in intact nuclei. Here, we describe how CAGE and RADICL-Seq data can be used to characterize enhancer elements and identify their target genes.

Stress-induced RNA\u2013chromatin interactions promote endothelial dysfunction

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

Genome-Wide Technologies to Study RNA-Chromatin Interactions.

An increasing number of studies have revealed that long non-coding RNAs (lncRNAs) play important roles in gene regulation and nuclear organization. Although the mechanisms are still largely unknown, many lncRNAs have been shown to interact with chromatin. Thus, one approach to understanding the function of these lncRNAs is to identify their sites of genomic interaction. Hybridization capture methods using oligonucleotide probes have been used for years to study chromatin-associated RNA. Recently, several groups have developed novel methods based on proximity ligation to investigate RNA–chromatin interactions at a genome-wide scale. This review discusses these technologies and highlights their advantages and disadvantages for the consideration of potential users.

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

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

Author Correction: GRID-seq for comprehensive analysis of global RNA-chromatin interactions.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Analyzing RNA-DNA Triplex Formation in Chromatin.

A significant fraction of non-coding RNAs (ncRNAs) is associated with chromatin, shown to regulate gene expression and to organize nuclear architecture. Mechanisms of direct and indirect RNA-chromatin interactions have been described, including the sequence-specific formation of triple helix structures. Triplexes are formed by the sequence-specific binding of RNA to the bases located in the major groove of DNA. We recently showed that triplexes do exist in the context of cellular chromatin and that these structures are stabilized by the histone H3 tail of adjacent nucleosomes. The in vitro characterization of the specificity and binding affinity of triplex sequences next to nucleosomes are essential parameters to identify potential sites of RNA-chromatin interaction in vivo. Here we provide a detailed protocol to determine the influence of nucleosome positioning on triple helix formation. This assay allows the comparative quantification of triplex formation and specificity for triplex targeting sequences relative to the spatial nucleosome position.

Mapping RNA-Chromatin Interactions In Vivo with RNA-DamID.

Long-noncoding RNAs (lncRNAs) are emerging as regulators of development and disease. lncRNAs are expressed in exquisitely precise expression patterns in vivo and many interact with chromatin to regulate gene expression. However, the limited sensitivity of RNA-purification techniques has precluded the identification of genomic targets of cell-type specific lncRNAs. RNA-DamID is a powerful new approach to understand the mechanisms by which lncRNAs act in vivo. RNA-DamID is highly sensitive and accurate, and can resolve cell-type-specific chromatin binding patterns without cell isolation. The determinants of RNA-chromatin interactions can be identified with RNA-DamID by analyzing RNA and protein cofactor mutants. Here we describe how to implement RNA-DamID and the design considerations to take into account to accurately identify lncRNA-chromatin interactions in vivo.

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

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

Mapping RNA-chromatin interactions by sequencing with iMARGI.

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

29 Introduction to functional annotation of animal genomes consortium

Abstract The annotation of the functional components of genomes is required for both basic and applied genomic research, yet farmed animal genomics are deficient in such annotation. This talk will introduce the Functional Annotation of ANimal Genomes (FAANG) consortium, which was initiated in 2014 to address this knowledge deficit. An overarching theme of early FAANG efforts is the collaborative approach required for such comprehensive research, thus another emphasis of this talk will be our community efforts to develop the environment for successful FAANG research. FAANG scientists have created a coordinated data collection and analysis enterprise crucial for success of this global effort, and funding of ~$20 million has been secured for current projects in several species, with an anticipated doubling of this support during 2019. Member FAANG consortium labs are producing genome-wide datasets on RNA expression, DNA methylation, chromatin modification, chromatin accessibility, and chromatin interactions for many agriculturally relevant animal species. The goal of a first-generation global chromatin state map for cattle, chicken, pig, and potentially other species is projected for completion in the next 3–5 years. These data will be used both to better understand animal genome function at the epigenetic level, as well as improve the precision and sensitivity of genomic selection for animal improvement. The functional annotation delivered by the FAANG initiative will add value and utility to the greatly improved genome sequences being established for domesticated animal species.