Next Article in Journal
Expression of Piwi Genes during the Regeneration of Lineus sanguineus (Nemertea, Pilidiophora, Heteronemertea)
Next Article in Special Issue
RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow
Previous Article in Journal
Approaching Shared Pathophysiology in Immune-Mediated Diseases through Functional Genomics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Research Center of Biotechnology, Institute of Bioengineering, Russian Academy of Science, 119071 Moscow, Russia
2
Department of Biological and Medical Physics, Moscow Institute of Physics and Technology, Dolgoprudny, 141701 Moscow Region, Russia
*
Author to whom correspondence should be addressed.
Genes 2020, 11(12), 1483; https://doi.org/10.3390/genes11121483
Submission received: 24 October 2020 / Revised: 3 December 2020 / Accepted: 4 December 2020 / Published: 10 December 2020
(This article belongs to the Special Issue Algorithms and Workflows in RNA Bioinformatics)

Abstract

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

1. Introduction

Long non-coding RNAs (lncRNAs)—transcripts of at least 200 nucleotides long lacking protein-coding potential—are the most populated class of transcripts in the mammalian genomes. They are expressed from more than 50,000 loci in the human genome [1], which is almost twice the number of protein-coding loci. LncRNAs are usually less conserved [2], lower expressed, and more cell-type-specific [1] compared to protein-coding genes. Still, the transcription of lncRNAs is regulated suggesting their functionality [1,3]. Several hundred human lncRNAs have been extensively characterized so far [4,5,6], revealing their roles in regulating transcription [7] and chromatin state [8,9,10,11]. However, the majority of the non-coding transcriptome remains non-annotated.
Up to date, several genomic screening approaches to functional annotation of lncRNA have been developed. Based on lncRNA knockdown (knockout) followed by cellular or molecular phenotyping, these approaches characterize a fair share of lncRNAs in one series of experiments [12,13,14,15]. In the FANTOM6 pilot project [12] more than 300 RNAs expressed in human primary dermal fibroblasts (HDFs) were suppressed using antisense LNA-modified GapmeR antisense oligonucleotide (ASO) technology. The effect of each knockdown was estimated by Cap Analysis Gene Expression (CAGE) deep sequencing to reveal molecular pathways associated with each lncRNA.
Complementary approaches to lncRNA annotation aim to search for their targets via detecting interactions between lncRNAs and chromatin or other RNA molecules. LncRNAs use different molecular mechanisms to bind to the chromatin, including RNA-DNA binding mediated by protein complexes, direct RNA-DNA hybridization via triplexes, RNA binding to single-stranded DNA regions (known as R-loops) and co-transcriptional RNA-RNA interactions via duplexes with nascent RNA transcripts (reviewed in [16]). Experimental techniques to detect RNA:chromatin interactions (ChIRP-seq [17], MARGI [18], ChAR-seq [19], GRID-seq [20], RADICL-seq [21], Red-C [22]) usually cannot differentiate between various mechanisms of RNA:chromatin interactions and focus solely on the location of such interaction. Direct RNA:RNA interactions can be detected with SPLASH [23], PARIS [24], LIGR-seq [25], MARIO [26] and other methods. Yet, the majority of the contacts obtained by such methods reflect RNA structure leading to very limited information on the inter-molecular lncRNA:RNA interactions available. Due to these limitations computational approaches for detecting RNA:RNA interactions provide valuable information for lncRNA functional annotation.
In this work, we focus on RNA:RNA interactions. We investigated two possible scenarios— lncRNA bound co-transcriptionally to the nascent RNAs via duplexes and lncRNA bound to mature transcripts—by applying the ASSA RNA-RNA interaction prediction tool [27] to the data on expression changes after lncRNA knockdown available from the FANTOM6 pilot project. Therefore, the present study pursued two goals: (i) to identify lncRNAs that are likely to function via direct RNA:RNA interactions with its targets and (ii) for the identified lncRNAs to predict the most likely interaction mode (co-transctiptional or post-transcriptional).

2. Materials and Methods

2.1. Computational Prediction of Inter-Molecular RNA:RNA Interactions

The AntiSense Search Approach (ASSA, https://github.com/vanya-antonov/assa) tool was used in this study to predict inter-molecular RNA:RNA interactions [27]. We decided to use ASSA because previously its performance has been evaluated on several well studied lncRNA data sets (RIA-seq for TINCR lncRNA [28], SPLASH for 50 different RNAs [23] as well as 17 lncRNA:RNA interactions from selected publications). We compared the results obtained by ASSA with 17 different RNA:RNA prediction tools. Although accurate computational prediction remains tricky, ASSA outperformed all other tested tools (see [16] for more details).

2.2. Distinguishing Co-Transcriptional and Post-Transcriptional Interactions

To distinguish co-transcriptional and post-transcriptional modes of lncRNA action we prepared two sets of target RNA sequences. We assumed that in the co-transcriptional mode the lncRNA can directly base pair with the nascent transcripts. Thus, the co-transcriptional lncRNA-RNA interactions may regulate gene expression at the transcriptional level. Since the regulation of transcription frequently occurs via gene promoters, we considered the 5’ end of the nascent transcript only. Specifically, we predicted possible interactions between query mature lncRNA (without introns) with the first 1kb of the target gene (including the promoter and all the introns, if any). We included the introns assuming that lncRNA-RNA binding may precede splicing.
By contrast, post-transcriptional lncRNA:RNA hybridizations occur with mature spliced target transcripts in the cytoplasm or the nucleus (depending on cellular localization of the lncRNA). We hypothesized that this type of interaction results in post-transcriptional regulation of gene expression. It should be noted that CAGE-seq technology which has been used in the FANTOM6 project measures expression changes at the promoter level. Thus, to predict post-transcriptional interactions we considered the longest isoform that can be expressed from the corresponding promoter only.
In the present work ASSA was run with default settings in both the co- and post-transcriptional modes.

2.3. Statistical Analysis

To estimate the correspondence between the iMARGI peaks and ASSA predictions for each lncRNA we only considered differentially expressed genes (DEGs) that had predicted antisense interactions with the corresponding lncRNA (ASSA p-values < 0.01 ). It should be noted that the fraction of the ASSA hits among all the DEGs for a particular lncRNA varied depending on the knockdown experiments.
For each lncRNA correspondence between the DEGs with the predicted ASSA hits and the experimentally identified lncRNA:DNA interactions (iMARGI peaks [18]) were computed using the GenometriCorr R package [29]. The tool computes the statistical significance of the observed co-localization between two sets of genomic intervals. The p-values were calculated based on the 100 permutations.
The statistical significance of the gene overlaps was estimated using the hypergeometric distribution (the phyper R function). We used 45480 as the total number of the human genes that was computed as the sum of the protein-coding (19954), long (17957) and small (7569) non-coding RNA genes as annotated in the GENCODE version 35. The human orthologs of the DEGs Chaserr-/- mEFs [30] were identified using the getLDS function of the biomaRt R package. The expression of the CHD2 gene was obtained from the FANTOM5 SSTAR online tool (https://fantom.gsc.riken.jp/5/sstar/EntrezGene:1106) for the “Fibroblast-skin normal, donor2 (nuclear fraction).CNhs12582.14302-155B9” sample.
The code and the corresponding data files that were used to generate the images in the manuscript are located in Supplementary Data S1 and also available at https://github.com/vanya-antonov/article_assa_and_f6.

3. Results

3.1. lncRNA Knockdown Experiments Suggest a Widespread Regulation Via RNA:RNA Interactions

We used the AntiSense Search Approach (ASSA) tool [27] to identify lncRNAs that were likely to regulate gene expression via direct hybridization with the target transcripts. In our previous benchmarking study, we have shown that ASSA outperformed other computational tools both in accuracy and in execution time [16]. In addition to computing the free energy of inter-molecular interaction between two RNAs, ASSA estimates its statistical significance (adjusted p-values) taking into account the lengths and the GC contents of both transcripts.
We analyzed all the lncRNAs knockdown experiments available from the recent pilot FANTOM6 project [12]. In brief, for each lncRNA query we investigated the enrichment of strong ASSA hits among the differentially expressed genes (DEGs) after the knockdown of this lncRNA. We considered two possible scenarios of lncRNA antisense interactions—co- and post-transcriptional modes. To do so we prepared two sets of the target RNA sequences: (i) either the 1 kb downstream regions (including introns, if any) starting from the target promoters (co-transcriptional mode) or (ii) the mature (i.e., without introns) RNA transcripts produced from the corresponding promoters (post-transcriptional mode). This allowed us to investigate possible lncRNA functioning in different cellular compartments–nucleus and cytoplasm (see Figure 1 and Methods for more details).
Among the 337 analyzed knockdown experiments statistically significant enrichment in the number of ASSA hits in the DEGs was identified for 98 and 47 ASOs in the co- and post-transcriptional mode, respectively, suggesting possible regulation of these DEGs via lncRNA:RNA interactions (see Figure 2A,B). The number of the cases of co-transcriptional interactions was more than twice higher than that of the post-transcriptional interactions. Interestingly, the majority (41 out of the 47) of the ASO knockdowns that showed enrichment in the post-transcriptional mode had enrichment in the co-transcriptional mode as well (Figure 2C). These 41 knockdown experiments corresponded to 31 different lncRNAs. It seemed unlikely that all these lncRNAs were able to function in both modes at the same time. To determine the molecular mechanism that was more likely for those 31 lncRNAs we used independent experimental data on the global RNA-chromatin interactions obtained by the iMARGI technology [18,31]. We investigated the co-localization between the iMARGI peaks for CHASERR and the promoters that were differentially expressed upon its knockdown and detected strong ASSA hits in the co-transcriptional mode [29]. Statistically significant co-localization (p-value < 0.01) was observed for 21 of those ambiguous 41 ASO knockdowns (Figure 2D) indicating that the corresponding lncRNAs were likely to function co-transcriptionally. Altogether, 35% and 14% of the FANTOM6 knockdowns showed enrichment of the ASSA hits among the DEGs in the co- and post-transcriptional modes, respectively. This indicated that direct hybridization with the target transcripts may be a common molecular mechanism of the lncRNA action.

3.2. Evolutionary Conserved CHASERR (AC013394.2) lncRNA May Direct the CHD2 Helicase to the Specific Genomic Loci by Interacting with Nascent Transcripts

Two knockdown experiments for the AC013394.2 lncRNA (also known as CHASERR, LINC01578 or LOC100507217) demonstrate coordinated and strong enrichment of the ASSA predicted RNA:RNA duplexes in DEG (Figure 2D). The predicted RNA:RNA contacts are supported by iMARGI co-localization data (Supplementary Figure S4). It should be noted that in the FANTOM6 project there were three independent ASO knockdowns of the CHASERR lncRNA. However, only two of them (ASO_G0272888_AD_07 and ASO_G0272888_AD_10) can be considered successful as the third experiment (ASO_G0272888_04) did not result in a significant reduction of CHASERR levels (Supplementary Figure S1). Importantly, the two successful knockdowns resulted in a change of expression for similar sets of genes. Specifically, more than 50% of the DEGs identified in the ASO_G0272888_AD_10 experiment were also detected in the ASO_G0272888_AD_07 (Supplementary Figure S2, hypergeometric test p-value = 0). This indicated the CHASERR functional significance and reproducibility of the results.
It has recently been shown that the mouse ortholog of CHASERR is mainly localized in the nucleus where it is bound to the chromodomain helicase DNA binding protein 2 (Chd2) [30]. The chromatin-remodeling enzyme CHD2 facilitates disassembly, eviction, sliding, and spacing of nucleosomes [32]. Thus, CHD2 localization to gene promoters opens the chromatin, increases DNA accessibility and expression [33]. Importantly, the primary sequences of human and mouse CHASERR orthologs are highly similar to each other (Supplementary Figure S5). To check the conservation of the CHASERR functionality we obtained human orthologs of all the DEGs identified in the Chaserr-/- mouse embryonic fibroblases (mEFs) (see Methods) and compared them to the lists of the DEGs from the two FANTOM6 knockdowns. We observed a statistically significant overlap between the three sets of genes (hypergeometric test p-values = 0) supporting conservation of the CHASERR target genes (Figure 3A).
Taking into account the sequence and functional similarities between human and mouse CHASERR orthologs as well as the fact that the CHD2 helicase is also expressed in the human fibroblasts (FANTOM5 expression = 113.25), we assumed that CHASERR lncRNA and CHD2 helicase interact in the nucleus of human fibroblasts as well. Our computational analysis suggested that CHASERR may directly interact with the beginning (the 5 part) of the target nascent transcripts and thus localize the CHD2 protein to the promoters of the corresponding genes. The subsequent opening of the chromatin would allow better access for the transcription factors to the promoter regions and activation of the gene expression. Indeed, both CHASERR knockdown experiments in the human fibroblasts resulted mainly in down-regulation of the target genes (Supplementary Figure S3). Moreover, 126 (57.5%) of the 219 DEGs that were common between the two knockdowns and the Chaserr-/- mEFs showed expression reduction in all the three experiments (Figure 3B). Finally, it has been experimentally shown that CHASERR interacts with transcripts of 7 different genes in Hela and HEK293T human cells [24,34]. Two of these 7 genes (TNPO2 and NUP153) were present among the 3599 differentially expressed genes identified in the two successful CHASERR knockdowns in the FANTOM6 project (hypergeometric test p-values = 0.0136). It should also be noted that ASSA predicted statistically significant interactions between CHASERR and transcripts of these two genes (with p-values 2.45 × 10 8 and 1.66 × 10 12 respectively). Altogether, these results indicated that CHASERR mainly enhanced the expression of the target genes possibly by interacting with the nascent transcripts and directing the CHD2 helicase to corresponding promoter regions.

4. Discussion and Conclusions

The main goal of the recent FANTOM6 pilot project was to perform a functional annotation of hundreds of human long non-coding RNAs expressed in human fibroblasts [12]. Although a tremendous amount of work has been done in performing thousands of lncRNA knockdown experiments, subsequent interpretation of the obtained results turned out to be a challenging task. There were several important aspects that complicated the bioinformatics analysis, including (i) relatively low expression of the majority of the selected lncRNAs, (ii) inability of some of the ASOs to efficiently reduce the levels of the corresponding lncRNAs, (iii) non-specific cell response to the performed perturbations, (iv) wide variety of different molecular mechanisms that can be used by lncRNAs to regulate target genes, (v) lack of computationally efficient bioinformatics tools that can produce reliable predictions for lncRNAs.
In the present study we investigated possible regulation via direct lncRNA:RNA interactions. To do so we used ASSA [27]–a tool that not only allows analysis of thousands of possible interactions in a reasonable execution time but also provides an unbiased estimate of the statistical significance of each prediction [16]. The design of the current analysis allowed us to distinguish putative co-transcriptional and post-transcriptional lncRNA interactions. Around 30% of all the FANTOM6 knockdown experiments showed a statistically significant enrichment of lncRNA:RNA duplexes among the DEGs indicating that corresponding lncRNAs may regulate their target genes via direct base-pairing with the transcripts. Surprisingly, the majority of such lncRNAs were likely to function in the nucleus where they can interact with the nascent RNAs co-transcriptionally. This observation was also supported with the co-localization of the predicted lncRNA-RNA interactions and the lncRNA-chromatin peaks experimentally identified by iMARGI technology.
Strong enrichment of predicted lncRNA:RNA duplexes in the DEGs were observed in both successful knockdowns of the CHASERR lncRNA. Moreover, ASSA predictions for these two experiments correlated well with the iMARGI peaks suggesting that CHASERR may regulate its target genes via RNA:RNA interactions co-transcriptionally. Finally, it has previously been shown experimentally that CHD2 had distinct binding profiles across active TSS regions [35]. The authors have concluded that transcription-coupled recruitment of CHD2 occurs at transcribed gene TSSs. Our bioinformatics analysis led us to the same conclusion and highlighted the role of CHASERR lncRNA in this process.
In an attempt to additionally verify the co-transcriptional RNA:RNA interactions we analyzed 9 experimental interactions of CHASERR detected by PARIS method. Unfortunately, these interactions were observed in Hela and HEK293T cells, while FANTOM6 pilot data was obtained in fibroblasts. To the best of our knowledge cell-type specificity of RNA:RNA contacts has never been thoroughly investigated. Yet, lncRNA expression is highly tissue specific, suggesting that contacts obtained in one cell type may not be preserved in another one. On top of that, not all the contacts may lead to functional consequences such as change of gene expression. Given that, the presence of experimental support for two genes TNPO2 and NUP153 is remarkable.
Importantly, a recent study has shown that the Chaserr mouse ortholog is indeed localized in the nucleus and regulates gene expression via interactions with chromatin [30] supporting the validity of our approach. Moreover, in the original study the authors mainly focused on the cis-action of Chaserr. Our results suggested that CHASERR may be able to enhance gene expression in trans by directing the CHD2 helicase to specific genomic loci. The model of CHASERR-guided gene regulation is provided in Figure 4. Since CHASERR enhances gene expression, its knockdown would lead to decreased expression of the regulated genes - the result we obtained by analyzing FANTOM6 data.
To conclude, our global analysis of the FANTOM6 experimental data showed that lncRNA:RNA interactions may be a common targeting mechanism to regulate gene expression (mainly co-transcriptionally). A more detailed analysis of the CHASERR lncRNA suggests the role of lncRNA:nascent RNA co-transcriptional interactions in gene regulation in trans. Therefore, future studies of the nuclear lncRNAs and their interactions with the chromatin may shed a new light on their regulatory mechanisms and further improve their functional annotation.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/12/1483/s1. Supplementary Figure S1: Efficiency of all the ASO knockdowns in the FANTOM6 project; Supplementary Figure S2: The overlap between differentially expressed genes identified in two independent knockdown experiments of the CHASSER lncRNA; Supplementary Figure S3: Types of differentially expressed genes in two independent CHASSER knockdown experiments and in Chasser-/- mouse embryonic fibroblases (mEFs); Supplementary Figure S4: Output of the GenometriCorr tool for the two ASOs corresponding to the lncRNA AC013394.2 (CHASERR); Supplementary Figure S5: Pairwise sequence alignment of the human and mouse CHASSER orthologs; Supplementary Data S1: The code and the corresponding data files that were used to generate the images in the manuscript.

Author Contributions

Y.M. and I.A. formulated and evaluated the ideas. I.A. conducted the research and wrote the original draft. Y.M. supervised the research, reviewed and edited the original draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Russian Science Foundation grant number 18-14-00240 to YAM.

Acknowledgments

We would like to thank Igor Ulitsky (Weizmann Institute of Science, Israel) for providing the details regarding the Chaserr lncRNA analysis [30].

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hon, C.C.; Ramilowski, J.A.; Harshbarger, J.; Bertin, N.; Rackham, O.J.; Gough, J.; Denisenko, E.; Schmeier, S.; Poulsen, T.M.; Severin, J.; et al. An atlas of human long non-coding RNAs with accurate 5’ ends. Nature 2017, 543, 199–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ulitsky, I. Evolution to the rescue: Using comparative genomics to understand long non-coding RNAs. Nat. Rev. Genet. 2016, 17, 601–614. [Google Scholar] [CrossRef] [PubMed]
  3. Alam, T.; Medvedeva, Y.A.; Jia, H.; Brown, J.B.; Lipovich, L.; Bajic, V.B. Promoter analysis reveals globally differential regulation of human long non-coding RNA and protein-coding genes. PLoS ONE 2014, 9, e109443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Quek, X.C.; Thomson, D.W.; Maag, J.L.; Bartonicek, N.; Signal, B.; Clark, M.B.; Gloss, B.S.; Dinger, M.E. lncRNAdb v2. 0: Expanding the reference database for functional long noncoding RNAs. Nucleic Acids Res. 2015, 43, D168–D173. [Google Scholar] [CrossRef] [PubMed]
  5. Volders, P.J.; Verheggen, K.; Menschaert, G.; Vandepoele, K.; Martens, L.; Vandesompele, J.; Mestdagh, P. An update on LNCipedia: A database for annotated human lncRNA sequences. Nucleic Acids Res. 2015, 43, D174–D180. [Google Scholar] [CrossRef] [Green Version]
  6. Ma, L.; Cao, J.; Liu, L.; Du, Q.; Li, Z.; Zou, D.; Bajic, V.B.; Zhang, Z. LncBook: A curated knowledgebase of human long non-coding RNAs. Nucleic Acids Res. 2019, 47, D128–D134. [Google Scholar] [CrossRef] [Green Version]
  7. Engreitz, J.M.; Ollikainen, N.; Guttman, M. Long non-coding RNAs: Spatial amplifiers that control nuclear structure and gene expression. Nat. Rev. Mol. Cell Biol. 2016, 17, 756–770. [Google Scholar] [CrossRef] [Green Version]
  8. Guttman, M.; Donaghey, J.; Carey, B.W.; Garber, M.; Grenier, J.K.; Munson, G.; Young, G.; Lucas, A.B.; Ach, R.; Bruhn, L.; et al. lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature 2011, 477, 295–300. [Google Scholar] [CrossRef] [Green Version]
  9. Guttman, M.; Rinn, J.L. Modular regulatory principles of large non-coding RNAs. Nature 2012, 482, 339–346. [Google Scholar] [CrossRef] [Green Version]
  10. Quinn, J.J.; Chang, H.Y. Unique features of long non-coding RNA biogenesis and function. Nat. Rev. Genet. 2016, 17, 47. [Google Scholar] [CrossRef]
  11. Ransohoff, J.D.; Wei, Y.; Khavari, P.A. The functions and unique features of long intergenic non-coding RNA. Nat. Rev. Mol. Cell Biol. 2018, 19, 143. [Google Scholar] [CrossRef] [PubMed]
  12. Ramilowski, J.A.; Yip, C.W.; Agrawal, S.; Chang, J.C.; Ciani, Y.; Kulakovskiy, I.V.; Mendez, M.; Ooi, J.L.C.; Ouyang, J.F.; Parkinson, N.; et al. Functional annotation of human long noncoding RNAs via molecular phenotyping. Genome Res. 2020, 30, 1060–1072. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, S.J.; Horlbeck, M.A.; Cho, S.W.; Birk, H.S.; Malatesta, M.; He, D.; Attenello, F.J.; Villalta, J.E.; Cho, M.Y.; Chen, Y.; et al. CRISPRi-based genome-scale identification of functional long noncoding RNA loci in human cells. Science 2017, 355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Liu, Y.; Cao, Z.; Wang, Y.; Guo, Y.; Xu, P.; Yuan, P.; Liu, Z.; He, Y.; Wei, W. Genome-wide screening for functional long noncoding RNAs in human cells by Cas9 targeting of splice sites. Nat. Biotechnol. 2018, 36, 1203–1210. [Google Scholar] [CrossRef]
  15. Joung, J.; Engreitz, J.M.; Konermann, S.; Abudayyeh, O.O.; Verdine, V.K.; Aguet, F.; Gootenberg, J.S.; Sanjana, N.E.; Wright, J.B.; Fulco, C.P.; et al. Genome-scale activation screen identifies a lncRNA locus regulating a gene neighbourhood. Nature 2017, 548, 343–346. [Google Scholar] [CrossRef]
  16. Antonov, I.V.; Mazurov, E.; Borodovsky, M.; Medvedeva, Y.A. Prediction of lncRNAs and their interactions with nucleic acids: Benchmarking bioinformatics tools. Briefings Bioinform. 2019, 20, 551–564. [Google Scholar] [CrossRef]
  17. Chu, C.; Qu, K.; Zhong, F.L.; Artandi, S.E.; Chang, H.Y. Genomic maps of long noncoding RNA occupancy reveal principles of RNA-chromatin interactions. Mol. Cell 2011, 44, 667–678. [Google Scholar] [CrossRef] [Green Version]
  18. Sridhar, B.; Rivas-Astroza, M.; Nguyen, T.C.; Chen, W.; Yan, Z.; Cao, X.; Hebert, L.; Zhong, S. Systematic Mapping of RNA-Chromatin Interactions In Vivo. Curr. Biol. 2017, 27, 602–609. [Google Scholar] [CrossRef] [Green Version]
  19. Bell, J.C.; Jukam, D.; Teran, N.A.; Risca, V.I.; Smith, O.K.; Johnson, W.L.; Skotheim, J.; Greenleaf, W.J.; Straight, A. Chromatin-associated RNA sequencing (ChAR-seq) maps genome-wide RNA-to-DNA contacts. Elife 2018, 7, e27024. [Google Scholar] [CrossRef]
  20. Li, X.; Zhou, B.; Chen, L.; Gou, L.T.; Li, H.; Fu, X.D. GRID-seq reveals the global RNA-chromatin interactome. Nat. Biotechnol. 2017, 35, 940–950. [Google Scholar] [CrossRef]
  21. Bonetti, A.; Agostini, F.; Suzuki, A.M.; Hashimoto, K.; Pascarella, G.; Gimenez, J.; Roos, L.; Nash, A.J.; Ghilotti, M.; Cameron, C.J.; et al. RADICL-seq identifies general and cell type–specific principles of genome-wide RNA-chromatin interactions. Nat. Commun. 2020, 11, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Gavrilov, A.A.; Zharikova, A.A.; Galitsyna, A.A.; Luzhin, A.V.; Rubanova, N.M.; Golov, A.K.; Petrova, N.V.; Logacheva, M.D.; Kantidze, O.L.; Ulianov, S.V.; et al. Studying RNA–DNA interactome by Red-C identifies noncoding RNAs associated with various chromatin types and reveals transcription dynamics. Nucleic Acids Res. 2020, 48, 6699–6714. [Google Scholar] [CrossRef] [PubMed]
  23. Aw, J.G.A.; Shen, Y.; Wilm, A.; Sun, M.; Lim, X.N.; Boon, K.L.; Tapsin, S.; Chan, Y.S.; Tan, C.P.; Sim, A.Y.L.; et al. In Vivo Mapping of Eukaryotic RNA Interactomes Reveals Principles of Higher-Order Organization and Regulation. Mol. Cell 2016, 62, 603–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lu, Z.; Zhang, Q.C.; Lee, B.; Flynn, R.A.; Smith, M.A.; Robinson, J.T.; Davidovich, C.; Gooding, A.R.; Goodrich, K.J.; Mattick, J.S.; et al. RNA Duplex Map in Living Cells Reveals Higher-Order Transcriptome Structure. Cell 2016, 165, 1267–1279. [Google Scholar] [CrossRef] [Green Version]
  25. Sharma, E.; Sterne-Weiler, T.; O’Hanlon, D.; Blencowe, B.J. Global Mapping of Human RNA-RNA Interactions. Mol. Cell 2016, 62, 618–626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Nguyen, T.C.; Cao, X.; Yu, P.; Xiao, S.; Lu, J.; Biase, F.H.; Sridhar, B.; Huang, N.; Zhang, K.; Zhong, S. Mapping RNA-RNA interactome and RNA structure in vivo by MARIO. Nat. Commun. 2016, 7, 12023. [Google Scholar] [CrossRef] [Green Version]
  27. Antonov, I.; Marakhonov, A.; Zamkova, M.; Medvedeva, Y. ASSA: Fast identification of statistically significant interactions between long RNAs. J. Bioinform. Comput. Biol. 2018, 16, 1840001. [Google Scholar] [CrossRef]
  28. Kretz, M.; Siprashvili, Z.; Chu, C.; Webster, D.E.; Zehnder, A.; Qu, K.; Lee, C.S.; Flockhart, R.J.; Groff, A.F.; Chow, J.; et al. Control of somatic tissue differentiation by the long non-coding RNA TINCR. Nature 2013, 493, 231–235. [Google Scholar] [CrossRef]
  29. Favorov, A.; Mularoni, L.; Cope, L.M.; Medvedeva, Y.; Mironov, A.A.; Makeev, V.J.; Wheelan, S.J. Exploring massive, genome scale datasets with the GenometriCorr package. PLoS Comput. Biol. 2012, 8, e1002529. [Google Scholar] [CrossRef] [Green Version]
  30. Rom, A.; Melamed, L.; Gil, N.; Goldrich, M.J.; Kadir, R.; Golan, M.; Biton, I.; Perry, R.B.T.; Ulitsky, I. Regulation of CHD2 expression by the Chaserr long noncoding RNA gene is essential for viability. Nat. Commun. 2019, 10, 1–15. [Google Scholar] [CrossRef] [Green Version]
  31. Yan, Z.; Huang, N.; Wu, W.; Chen, W.; Jiang, Y.; Chen, J.; Huang, X.; Wen, X.; Xu, J.; Jin, Q.; et al. Genome-wide colocalization of RNA–DNA interactions and fusion RNA pairs. Proc. Natl. Acad. Sci. USA 2019, 116, 3328–3337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Narlikar, G.J.; Sundaramoorthy, R.; Owen-Hughes, T. Mechanisms and functions of ATP-dependent chromatin-remodeling enzymes. Cell 2013, 154, 490–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Semba, Y.; Harada, A.; Maehara, K.; Oki, S.; Meno, C.; Ueda, J.; Yamagata, K.; Suzuki, A.; Onimaru, M.; Nogami, J.; et al. Chd2 regulates chromatin for proper gene expression toward differentiation in mouse embryonic stem cells. Nucleic Acids Res. 2017, 45, 8758–8772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gong, J.; Shao, D.; Xu, K.; Lu, Z.; Lu, Z.J.; Yang, Y.T.; Zhang, Q.C. RISE: A database of RNA interactome from sequencing experiments. Nucleic Acids Res. 2018, 46, D194–D201. [Google Scholar] [CrossRef] [PubMed]
  35. Siggens, L.; Cordeddu, L.; Rönnerblad, M.; Lennartsson, A.; Ekwall, K. Transcription-coupled recruitment of human CHD1 and CHD2 influences chromatin accessibility and histone H3 and H3.3 occupancy at active chromatin regions. Epigenet. Chromatin 2015, 8, 1–14. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The design of the workflow to detect lncRNAs that may regulate target genes via direct base-pairing co- or post-transcriptionally.
Figure 1. The design of the workflow to detect lncRNAs that may regulate target genes via direct base-pairing co- or post-transcriptionally.
Genes 11 01483 g001
Figure 2. Enrichment of the lncRNA hybridization targets predicted by ASSA among differentially expressed genes identified in the FANTOM6 ASO knockdown experiments in two modes. (A,B) ASSA was applied to the differentially expressed and background transcripts assuming co-transcriptional or post-transcriptional lncRNA:RNA interaction mode. The- log 10 of enrichment p-values were computed using the hypergeometric (HGD) and GSEA tests. Knockdowns with either p-value < 0.01 are marked by the red color. There were 98 and 47 such ASOs for the co- and post-transcriptional modes, respectively. (C) The majority (41 out of the 47) of the ASOs with significant p-values in the post-transcriptional mode also had significant p-values in the co-transcriptional mode. (D) The heatmap showing four p-values for the 41 ASOs that produced strong p-values in both the co- and the post-transcriptional modes. The 41 ASO ids are on the right side and the corresponding 31 lncRNA gene names are on the left side of the heatmap. The p-values of the co-occurrence with the iMARGI peaks are shown on the right.
Figure 2. Enrichment of the lncRNA hybridization targets predicted by ASSA among differentially expressed genes identified in the FANTOM6 ASO knockdown experiments in two modes. (A,B) ASSA was applied to the differentially expressed and background transcripts assuming co-transcriptional or post-transcriptional lncRNA:RNA interaction mode. The- log 10 of enrichment p-values were computed using the hypergeometric (HGD) and GSEA tests. Knockdowns with either p-value < 0.01 are marked by the red color. There were 98 and 47 such ASOs for the co- and post-transcriptional modes, respectively. (C) The majority (41 out of the 47) of the ASOs with significant p-values in the post-transcriptional mode also had significant p-values in the co-transcriptional mode. (D) The heatmap showing four p-values for the 41 ASOs that produced strong p-values in both the co- and the post-transcriptional modes. The 41 ASO ids are on the right side and the corresponding 31 lncRNA gene names are on the left side of the heatmap. The p-values of the co-occurrence with the iMARGI peaks are shown on the right.
Genes 11 01483 g002
Figure 3. Differentially expressed genes (DEGs) observed in CHASERR knockdown experiments in human cells and in Chaserr-/- mouse embryonic fibroblasts (mEFs). (A) Overlap between the DEGs from three independent experiments. (B) Expression changes of the 219 genes that were classified as differentially expressed in all the three cases.
Figure 3. Differentially expressed genes (DEGs) observed in CHASERR knockdown experiments in human cells and in Chaserr-/- mouse embryonic fibroblasts (mEFs). (A) Overlap between the DEGs from three independent experiments. (B) Expression changes of the 219 genes that were classified as differentially expressed in all the three cases.
Genes 11 01483 g003
Figure 4. The proposed model of the CHASERR lncRNA regulation of the target genes involving the CHD2 protein and the RNA-RNA interactions with the nascent transcripts. The figure was created with BioRender.com online tool.
Figure 4. The proposed model of the CHASERR lncRNA regulation of the target genes involving the CHD2 protein and the RNA-RNA interactions with the nascent transcripts. The figure was created with BioRender.com online tool.
Genes 11 01483 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Antonov, I.; Medvedeva, Y. Direct Interactions with Nascent Transcripts Is Potentially a Common Targeting Mechanism of Long Non-Coding RNAs. Genes 2020, 11, 1483. https://doi.org/10.3390/genes11121483

AMA Style

Antonov I, Medvedeva Y. Direct Interactions with Nascent Transcripts Is Potentially a Common Targeting Mechanism of Long Non-Coding RNAs. Genes. 2020; 11(12):1483. https://doi.org/10.3390/genes11121483

Chicago/Turabian Style

Antonov, Ivan, and Yulia Medvedeva. 2020. "Direct Interactions with Nascent Transcripts Is Potentially a Common Targeting Mechanism of Long Non-Coding RNAs" Genes 11, no. 12: 1483. https://doi.org/10.3390/genes11121483

APA Style

Antonov, I., & Medvedeva, Y. (2020). Direct Interactions with Nascent Transcripts Is Potentially a Common Targeting Mechanism of Long Non-Coding RNAs. Genes, 11(12), 1483. https://doi.org/10.3390/genes11121483

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop