Eukaryotic genomes are pervasively transcribed. Besides protein-coding RNAs, there are different types of non-coding RNAs that modulate complex molecular and cellular processes. RNA sequencing technologies and bioinformatics methods greatly promoted the study of ncRNAs, which revealed ncRNAs’ essential roles in diverse aspects of biological functions. As important key players in gene regulatory networks, ncRNAs work with other biomolecules, including coding and non-coding RNAs, DNAs and proteins. In this review, we discuss the distinct types of ncRNAs, including housekeeping ncRNAs and regulatory ncRNAs, their versatile functions and interactions, transcription, translation, and modification. Moreover, we summarize the integrated networks of ncRNA interactions, providing a comprehensive landscape of ncRNAs regulatory roles.
As a template for protein synthesis, messenger RNAs (mRNAs) have become the major research focus for a long time, while non-coding RNAs (ncRNAs) were considered as by-products of massive transcription with less biological meaning. Since the discovery of ribosomal RNA (rRNA) and transfer RNA (tRNA) in the late 1950s, varieties of RNA species have gradually surfaced, which revealed an unsuspected non-coding world. The combination of large-scale sequencing and computational analysis greatly helps to understand the RNA world. At the beginning of the 21st century, initial sequencing and analysis of human  and mouse genome  revealed that 98% of the “junk” DNAs can be transcribed. Besides mRNAs that have already annotated, most transcripts do not seem to encode proteins. Therefore, those transcripts are commonly referred to as ncRNAs . Soon afterwards, the Human Genome Project (HGP) was achieved in 2005  and abundant long non-coding RNAs (lncRNAs) were detected in mammals , . Later, the widespread application of high throughput sequencing further allowed for more accurate profiling of ncRNAs , . The ENCODE (Encyclopedia of DNA Elements) project that launched in 2005 and its recent reports revealed that up to 80% of the human genome has the capacity to transcribe into ncRNAs , . The large scale datasets generated from those projects have facilitated the establishment of many public databases, including Rfam , NONCODE , miRbase , and circBase .
In general, ncRNAs are found to participate in multiple biological processes, regulate physiological and developmental processes or even disease. They have been identified as tumor suppressors and oncogenic drivers in various cancer types . In plants, ncRNAs have emerged as key regulatory molecules in stress responses . Accumulating evidence indicates that ncRNA analysis has become a frontier research trend, current signs of progress have yielded novel insights into them. With the development of experimental techniques and big data analysis, diverse interactions of ncRNAs have been characterized gradually and form interconnected complex networks. The overall ncRNA network consists of various types of molecules from multi-omics datasets, involved in different subcellular localizations. Considering the rapid growth of ncRNA field, a comprehensive review of the current advances in ncRNAs that involves varieties of ncRNA species and their interrelationship is urgent. In this review, we briefly describe ncRNA transcription and the non-coding members, and then summarize three layers of ncRNA related interactions, which is a help to illustrate those complex ncRNA integrated networks.
2 Various ncRNA Species in Eukaryotes
Eukaryotic transcription from different genomic regions and RNA processing produce various ncRNA species. Several parts of DNAs could be transcribed into ncRNAs, such as protein-coding genes, enhancer regions or transposon elements (Figure 1). Considering their regulatory roles, ncRNA could be divided into two categories (Table 1). Housekeeping ncRNAs are abundantly and ubiquitously expressed in cells, primarily regulate generic cellular functions. While regulatory ncRNAs are usually considered as key regulatory RNA molecules, function as regulators of gene expression at epigenetic, transcriptional, and post-transcriptional levels , , .
|Housekeeping ncRNAs||rRNA||ribosomal RNA||120–4,500 nt|
|tRNA||tansfer RNA||76–90 nt|
|snRNA||small nuclear RNA||100–300 nt|
|snoRNA||small nucleolar RNA||60–400|
|tRF||tRNA-Derived Fragments||16–28 nt|
|tiRNA||tRNA halves||29–50 nt|
|Regulatory ncRNAs||miRNA||microRNA||21–23 nt|
|siRNA||small interfering RNA||20–25 nt|
|piRNA||piwi-interacting RNA||26–32 nt|
|eRNA||enhancer RNA||50–2,000 nt|
|lncRNA||long non-coding RNAs||>200 nt|
|circRNA||circular RNA||100–10,000 nt|
|Y RNA||Y RNA||/|
2.1 Housekeeping ncRNAs
As early discovered ncRNA species, housekeeping ncRNAs have been intensively studied in the past decades, including rRNAs, tRNAs, small nuclear RNAs (snRNAs), small nucleolar RNAs (snoRNAs) and telomerase RNAs , . These ncRNAs are usually small, ranging from 50 nucleotides (nt) to 500 nt, constitutively expressed in all cell types and necessary for cell viability . Besides essential roles like rRNAs and tRNAs in protein synthesis, snRNAs in RNA splicing and snoRNAs in RNA modifications, some housekeeping ncRNAs could perform regulatory roles through cleavage. tRNA-derived RNA fragments (tRFs) and translation interfering tRNAs (tiRNAs) are new classes of small regulatory ncRNAs that are derived from tRNA or pre-tRNA . Studies have suggested that tiRNAs could inhibit translation through recruitment of innovative packed aggregates of proteins and RNAs under stress situation , . Additionally, some small RNAs derived from snoRNAs have been discovered by deep sequencing combing with bioinformatics analysis, such as sno-derived RNAs , sno-miRNAs , and sno-piRNAs .
2.2 Regulatory ncRNAs
Based on their average size, regulatory ncRNAs could be further classified as small non-coding RNAs (sncRNAs) with comprising transcripts fewer than 200 nt, and lncRNAs with greater than 200 nt. The main classes of small ncRNAs are microRNA (miRNA), small interfering RNAs (siRNAs), piwi-interacting RNAs (piRNAs). However, some ncRNAs with variable length might belong to two classifications at the same time, such as promoter-associated transcripts (PATs), enhancer RNAs (eRNAs), and circular RNAs (circRNAs).
miRNAs are the most abundant class of small ncRNAs that are generated from transcribed hairpin loop structures , , . They regulate gene expression in both cytoplasm and nucleus through different mechanisms and mediate gene silencing at the post-transcriptional level . As key players in the complex interplay among diverse RNA species, miRNAs have been considered research hotspot for several years. The latest miRBase release (v22) adds 48 new species, now contains 38,589 hairpin precursors and 48,860 mature miRNAs from 271 organisms, including animals, plants, unicellular algae, and viruses .
siRNAs are a class of double-stranded RNA molecules and play essential roles in RNA interference (RNAi) pathway. Based on their origins, siRNAs can be divided into two subcategories, exogenous siRNA (exo-siRNA), products of exogenous nucleic acid resulting from artificial insertion or virus infections, and endogenous siRNAs (endo-siRNA) that are arising from endogenous genomic locus and mostly transcribed from transposon elements (TEs) , . In plants, the group of siRNAs is very diverse and contains subcategories, including trans-acting siRNAs (tas-siRNAs), heterochromatic siRNAs (hc-siRNAs) , repeat-associated siRNAs (ra-siRNAs), long siRNAs (lsiRNAs) , and natural antisense siRNAs (nat-siRNAs) that are also identified in animals .
Named after PIWI proteins, piRNAs are a class of animal-specific small ncRNAs. Mature piRNAs are derived from two major pathways, precursors are transcribed from piRNA clusters and cleaved by PIWI proteins in the primary pathway, and then modified through “ping-pong” cycle in the second , . Distinct from miRNAs and siRNAs, piRNAs are processed in a Dicer-independent manner from single-stranded RNA precursors  and they only function through binding with PIWI proteins . piRNAs are associated with PIWI proteins and form piRNA-induced silencing complexes, silence their targets at transcriptional and post-transcriptional levels , which acts in a similar manner that endo-siRNAs silence TEs in plants. Endo-siRNAs also silence TEs in animals, but the piRNA pathway is at the forefront of defense against transposons in germ cells .
LncRNAs are defined as transcripts that longer than 200 nt and lack of protein-coding ability. Based on their location with respect to protein-coding genes, lncRNA can be divided into five subcategories . Long intergenic ncRNAs (lincRNAs) are transcribed from both strands in intergenic regions. Long intronic ncRNAs are completely transcribed from introns of protein-coding genes. Sense lncRNAs are transcribed from the sense strand and containing exons of protein-coding genes. In contrast to sense lncRNAs, antisense lncRNAs are transcribed from the antisense strand, overlapping with exonic or intronic regions, or covering the protein-coding sequence through an intron, while bidirectional lncRNAs, localized in the proximity of a coding transcript on the opposite strand . According to their regulatory effects on DNA sequences, lncRNA could be classified as cis-lncRNAs (cis-acting lncRNAs) that regulate the expression of genes in close proximity, and trans-lncRNAs (trans-acting lncRNAs) that regulate distant genes . Some lncRNAs could be further processed to generate small ncRNAs, such as miRNAs, piRNAs, and snoRNAs .
eRNAs are emerging as a group of ncRNAs transcribed from transcriptional enhancers. Besides similar length and protein-coding potentiality, eRNAs share some transcriptional features with lncRNAs, like having both sense and antisense transcripts by bidirectional transcription , , . Bidirectional transcription of enhancer regions generates comparatively shorter (0.5–2kb) and non-polyadenylated eRNAs (polyA-eRNAs) , while long (>4kb) and polyadenylated eRNAs (polyA + eRNAs) are usually unidirectional transcribed from enhancers  (Figure 1). As for eRNAs that are generated by transcription initiation at enhancers in intragenic regions, those are multiexonic polyA + eRNAs (meRNAs) , . Although having similar properties, eRNAs are 5′ capped, unstable, exhibiting short half-lives and easily degraded by the exosome complex, which was quite different from lncRNAs , , , , , .
CircRNAs are a unique class of endogenous ncRNAs that form covalently closed loop structures. Although discovered for decades, advanced deep sequencing and bioinformatics analysis made circRNA gain more attention recently . CircRNAs have a wide range of sizes, ranging from 100 nt to over 10,000 nt, while most of them in mammals and plants are hundreds of nucleotides , . CircRNAs are the circularization product from splicing events, and generate from exons, introns, intergenic regions, untranslated regions (UTRs) or even tRNAs . Similar to other regulatory ncRNAs, circRNAs play important roles in diverse biological processes, such as regulating alternative RNA splicing or transcription , acting as competing endogenous RNAs (ceRNAs)  or miRNA sponges .
3 ncRNA Interactions
3.1 ncRNAs Interact with mRNAs
To define biological functions for ncRNAs, some experimental approaches are employed for the investigation of ncRNA-based targeting, like PAR-CLIP , HITS-CLIP , and CLASH . HITS-CLIP is widely used for functional protein-RNA interaction sites identifying and miRNA-mRNA regulatory interactions analysis , . Compared to the HITS-CLIP, PAR-CLIP improves RNA recovery efficiency and indicates the site of targeting more precisely through identifying the location of the crosslink , . LIGR-seq (LIGation of interacting RNA followed by high-throughput sequencing) is a new developed method that enables the global-scale mapping of RNA-RNA duplexes in vivo , which enlarges the landscape of ncRNA-mRNA interactions. Those contain unexpected interactions between snoRNAs and mRNAs, such as the orphan C/D box snoRNA SNORD83B controls steady-state levels of target mRNAs. These detections may be additional examples for ncRNA-mediated gene regulation. Besides, the ANXA2 mRNA and the ANXA2 Pseudogene ANXA2P2 may have direct or indirect interaction to promote tumor cell invasion through remodeling the hepatocellular carcinoma (HCC) cells movement that could be a novel predictive biomarker for the HCC patients risk evaluation of its recurrence or metastasis .
Furthermore, miRNA-mRNA interactions would result in the silencing of mRNA expression. In neoplasms, upregulated miRNAs could act as oncogenes through silencing the mRNAs encoding tumor suppressor proteins . ncRNAs can also promote the expression of mRNA. The circular RNAs hsa_circ_0032462, hsa_circ_0028173, hsa_circ_0005909 regulate the Cell adhesion molecule 1 (CADM1) gene which is a co-expression mRNA functioning as a miRNA sponge in human osteosarcoma . Thus, these osteosarcoma-related miRNAs could provide a solution in the development of new therapies for human osteosarcoma .
3.2 ncRNA-ncRNA Interactions
ncRNA-ncRNA crosstalk indicates the intersecting non-coding world, and these interactions could impact numerous biological processes, such as epigenetic modifications, transcription, and translation; hence gain a new layer of genomic regulation. Those miRNA response elements (MREs) such as circRNAs, lncRNAs, and eRNA act as competitive endogenous RNAs (ceRNAs), which has abundant implications for gene regulation during multiple physiological and pathophysiological processes at the post-transcriptional level. For example, functional experiments suggested that the circular RNA hsa_circ_0001368 decelerates the growth of gastric cancer by regulating miR-6506-5p/FOXO3 axis by serving as a ceRNAs to sponge miR-6506-5p. Thus hsa_circ_0001368 plays a tumor-suppression role in gastric cancer and may be adapted to the gastric cancer therapy .
In addition to the interaction between different ncRNAs, such as ceRNAs mentioned above, there are interactions between the same kind of ncRNAs, such as miRNA-miRNA and lncRNA-lncRNA. The interactions among multiple miRNAs have enriched the complex mechanisms of post-transcriptional regulations . In mammals, Krek et al. experimentally validated a combination of miR-375, miR-124, and let-7b which lead to synergistic target inhibition . In addition, functional cooperation of miR-125a, miR-125b, and miR-205 can act together to suppress erbB2/erbB3 expression in breast cancer cells, which manifest a novel breast cancer treatment via miRNA-dependent or -independent mechanisms . What’s more, evidence indicated that lncRNA-lncRNA synergistic networks play a vital role in oncogenesis and tumor suppressor pathways. For instance, lncRNAs contribute to the p53 biological network, act as p53-regulator lncRNAs and p53-effector lncRNAs, thus may promote the development of diagnostic methods and therapies focused on these non-coding molecules . We could identify synergistic lncRNA pairs via combining analysis of genome-wide expression datasets and functional information .
3.3 ncRNAs Interact with DNAs
Many ncRNAs serve as direct regulators on chromatin, as exemplified by lncRNAs and tiRNAs. Recent progress in high-throughput RNA sequencing facilitates the detection of promoter-proximal small RNAs, including tiRNAs, promoter-associated small RNAs (PASRs) , and transcription start site-associated RNAs (TSSa RNAs) . They are possibly associated with nucleosome positioning, chromatin marking, and transcriptional regulation, such as tiRNAs are abundant at transcription-initiation chromatin marks.
An RNA sequencing technology called global RNA interactions with DNA by deep sequencing (GRID-seq) has been reported to detect chromatin-RNA interactions, comprehensively . Two highly expressed human lncRNA NEAT1 (nuclear-enriched abundant transcript 1) and MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) were identified by GRID-seq which localize to hundreds of genomic sites in mammals (mostly overactive genes) . NEAT1 and MALAT1 exhibit colocalization, although they display distinct gene body binding patterns at active chromatin. Moreover, the function of enhancer-associated miRNAs has been demonstrated to be able to up-regulate transcription through miRNA-promoter interaction. miR-24-1 has been associated with a function that activates eRNA expression, alters histone modification which results in p300 and RNA polymerase II’s increased enrichment at the enhancer locus . Besides, the lncRNA KHPS1 would induce the synthesis of eRNA Sphk1 towards a triplex-driven recruitment mechanism, which enhances the expression of specific target genes . In addition, the conserved Y RNAs were found to be key factors for the initiation of DNA replication, preferentially with replicated chromatin . Non-coding human Y RNAs can act not only as new cancer biomarkers and also molecular targets for anti-proliferative prevention . hY1 and hY3 RNAs have been found to result in a remarkable cytostatic suppression of cell proliferation.
3.4 ncRNAs Interact with Proteins
ncRNA-protein interactions have critical roles in all aspects of gene expression, therefore, plenty of approaches have been created for comprehensively analysis. Development of deep-sequencing methods together with immunoprecipitation of RNA binding proteins (RBPs) has manifested a great diversity of ncRNA-associated proteins .
Housekeeping ncRNAs interacting with proteins can shape into varied ribonucleoprotein (RNP) complexes that perform various functions. For instance, snRNAs together with multiple proteins form spliceosomes (snRNPs) that participate in canonical splicing and alternative splicing. Many nucleotides in the pre-rRNA, pre-snRNA, and pre-tRNA undergo post-transcriptional modifications via small nucleolar RNP particles . What’s more, interactions between regulatory RNAs and proteins are vital in interceding fundamental cellular processes. Small ncRNAs, like miRNAs, siRNAs, and piRNAs, are well known to affect RNA stability and translation during RNA interference pathway towards interacting with the Argonaute family proteins during RNA interference pathway. LncRNAs achieve their functions as RNPs in several ways like recruitment, inhibition, and act indirectly through genome organization and transcription . There is a well-studied lncRNA Xist which is required during development in female mammals. Three proteins: SHARP, SAF-A, and LBR -interact with Xist are required for transcriptional silencing . Moreover, circRNAs can serve as protein sponges translocate proteins to a specific subcellular location. The expression of circRNA-Foxo3 has significantly increased during cancer cell apoptosis, which is experimentally explored the possibility of inhibiting tumor growth with the delivery of circ-Foxo3 plasmid . It has been clarified that protein recruitment by lncRNAs and circRNAs modifies the cellular concentration and localization of proteins.
4 Integrated ncRNA Networks
More and more evidence proves that ncRNAs provide diverse regulatory functions in the cell. Those RNA mediated interactions are often interconnected, thus, constructing complex regulatory RNA networks. Most ncRNAs, especially housekeeping ncRNAs, have been characterized through their functions in the nucleus and their interactions with various types of nuclear machinery may contribute to their nuclear retention . However, many ncRNAs are also detected in the cytoplasm and play regulatory functions there. The ceRNA network is one of the most representative networks of regulatory ncRNA networks in transcriptomics, while those networks also involve crucial interphases between genomics and proteomics.
4.1 ceRNA Network
Based on the ceRNA hypothesis, the cross-talk between RNAs through MREs could form a large-scale regulatory network across the transcriptome, including both coding and non-coding RNAs . Like other ceRNAs namely pseudogenes, lncRNAs and circRNAs, mRNAs also contain MREs, which are normally located within the 3′UTR. All of these RNAs potentially compete for miRNA binding, thereby regulating miRNA repression. Thus, mRNAs, different types of ceRNAs and miRNAs would form a network of interactions, as called ceRNA network (Figure 2). ceRNAs usually contain carry MREs for multiple miRNAs, each miRNA can regulate multiple ceRNAs and most ceRNAs are regulated by more than one miRNA.
According to ceRNA theory, the existence of this network has three important conditions . First, the relative concentration of ceRNAs and their miRNAs would impact the expression of competing ceRNAs. Considering the density of ceRNA network connections, only the expression level of ceRNAs is large enough could overcome or relieve miRNAs’ repression on downstream mRNA targets. Second, the numbers of MRE normally influence the effectiveness of a ceRNA. Furthermore, in a given cell type or at a specific moment, the identity, concentration, and subcellular distribution of ceRNA and miRNA species impact ceRNA networks. Third, miRNA binding is impact by MREs, not all the MREs on ceRNAs are equal. The specific nucleotide composition of MREs might be partially different even if they are predicted to bind the same miRNA. The binding effectiveness of each MRE is critical for overall ceRNA function. Additionally, if a miRNA is sequestered by ceRNAs, its primary targets would be preferentially affected. Although hundreds of RNAs are predicted to be targeted by a miRNA, their degrees of repression are not same; only a few of them are primary targets while the rest are finely tuned , .
In ceRNA networks, the majority of validated ceRNAs are mRNAs. They compete for miRNA binding to sequester miRNAs from potential alternative targets, this biological role depends on their ability to compete for miRNA binding and are independent of protein-coding function . Zinc finger E-box binding homeobox 2 (ZEB2) mRNA, has been validated as a ceRNA and regulates PTEN expression levels by sequestering several miRNAs, including miR-181, miR-200b, miR-25, and miR-92a. Attenuation of ZEB2 expression leads to inhibition of PTEN in human melanoma . Additionally, 3′UTRs that contain MREs for multiple miRNAs, are critical for the function of mRNAs as ceRNAs. It has been reported that the versican (VCAN) 3′UTR binds to miR-144 and miR-136 to modulate PTEN level, and cell cycle regulator retinoblastoma 1 (Rb1) that targeted by miR-199a-3p and miR-144 could function as a ceRNA for VCAN . Thus, VCAN 3′UTR binds and modulates miRNA activities as a miRNA sponge, and Rb1 and PTEN mRNAs are released for translation. CD34 and FN1 are additional validated VCAN ceRNAs, competing for binding to miR-133a, miR-144, miR-199a-3p, and miR-431 , . As one of the most studied lncRNA, H19 is highly expressed in undifferentiated muscle cells and has canonical and non-canonical binding sites for the let-7 miRNA family. When the expression of H19 is reduced in differentiated cells, the expression of let-7 increases at the same time. H19 has been demonstrated competing with DICER and HMGA2 to modulate let-7 availability, as a molecular sponge . Besides, H19 together with circRNA MYLK and CTDP1 compete for miRNA-29a-3p binding to regulated the expression of its target genes DNMT3B, HAS3, VEGFA, and ITGB1, which would result in cancer growth and metastasis . Pbcas4, the pseudogene of breast carcinoma amplified sequence 4 (BCAS4), is a conserved ceRNA in mouse and human, whose transcripts compete with BCAS4 mRNAs for binding to the common miR-185 .
4.2 ncRNA-Mediated Networks in Genomics and Proteomics
Besides ncRNA-ncRNA interactions, ncRNAs could regulate biological processes through diverse mechanisms. ncRNAs and their mediated networks are an integral part of genomics and proteomics (Figure 3).
ncRNAs could contribute to local and long-range genomic interactions (Figure 3). Functional studies of eRNAs and some lncRNAs have demonstrated that ncRNAs mediate promoter-enhancer interactions to regulate the expression of various protein-coding genes , , . In the nucleus, ncRNAs play critical roles in the organization of the genome to coordinate the expression of gene clusters. Recent studies have shown that the Xist complex explores some larger genomic domains to help spread transcriptional repressor complexes during X-chromosome inactivation , . This strategy could also be used to establish active and inhibitory domains that involve genomic segments separated by a long linear distance on the same chromosome or even from a different chromosome, which in turn may contribute to the organization of the genome in 3D space , .
At the transcriptional and the post-transcriptional levels, regulated gene expression determines the cell type-specific proteome, and ncRNAs are widely involved in protein-mediated interaction networks (Figure 3). There are numerous RNA-dependent protein-protein and protein-DNA interactions in the cell, ncRNAs play pivotal roles in those regulations in which RBPs, lncRNAs, miRNAs, or even circRNAs are involved , . The positive feedforward of RBP HuR is a good example of protein-ncRNA interaction networks. HuR may associate with many mRNAs, which influence cell proliferation, survival, immunity, carcinogenesis, and stress responses , . HuR and lincMD1 are both involved in muscle differentiation, and these two are under the inhibitory control of miR-133. Interestingly, HuR was identified as another component of the lincMD1-regulated circuitry. In that circuitry, HuR binds lincMD1 and protects it from Drosha cleavage at the expense of miR-133b biogenesis .
Benefited from the rapid advance in sequencing technologies and bioinformatics approaches, different classes of ncRNAs have been detected and discovered in the past decades, leading to a complex ncRNA world. Sufficient evidence shows that ncRNA plays significant biological roles in regulation of cellular mechanisms, and is associated to higher eukaryotes complexity. ncRNAs can interact with different species of molecules, and regulate them through various molecular mechanisms generating a complex network. ncRNA-mediated networks not only contained interplay among RNAs, but also involved in genomics and proteomics, in both nucleus and cytoplasm.
The ceRNA hypothesis attributes potentially predictable function to the coding and noncoding transcriptome, which helps to systemically realize the ncRNA network. Although ceRNA networks are one of the most extensively studied ncRNA networks, novel ncRNA species will increase the complexity of these networks. Multiple types of molecules and biology pathways make it more difficult to comprehensively understand ncRNA functions in multi-level networks. Even though the accumulation of information about ncRNA interactions and their mediated networks is relatively slow, the development of tools in silico predictors and network analysis module may accelerate this process. ncRNA-mediated interactions could be captured by experiment detections, and those associations have been systematically analyzed and extended to more connections in various computation method. Based on those intricate interactions, network model could be constructed and predict potential ncRNA related associations.
Funding source: National Key Research and Development Program of China
Award Identifier / Grant number: 2018YFC0310602 and 2016YFA0501704
Funding source: National Natural Science Foundation of China
Award Identifier / Grant number: 31771477 and 31571366
Funding statement: This work has been supported by the National Key Research and Development Program of China (Grant Nos. 2018YFC0310602 and 2016YFA0501704), National Natural Science Foundation of China (Funder Id: Funder Id: http://dx.doi.org/10.13039/501100001809, Grant Nos. 31771477 and 31571366), the Fundamental Research Funds for the Central Universities, Jiangsu Collaborative Innovation Center for Modern Crop Production and 2018 Zhejiang University Academic Award for Outstanding Doctoral Candidates.
Conflict of interest statement: Authors state no conflict of interest. All authors have read the journal’s publication ethics and publication malpractice statement available at the journal’s website and hereby confirm that they comply with all its parts applicable to the present scientific work.
 Carninci P, Kasukawa T, Katayama S, Gough J, Frith M, Maeda N, et al. The transcriptional landscape of the mammalian genome. Science 2005;309:1559–63.10.1126/science.1112014Search in Google Scholar PubMed
 Liu C, Bai B, Skogerbo G, Cai L, Deng W, Zhang Y, et al. NONCODE: an integrated knowledge database of non-coding RNAs. Nucleic Acids Res 2005;33(Database issue):D112–5.10.1093/nar/gki041Search in Google Scholar PubMed
 Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res 2006;34(Database issue):D140–4.10.1093/nar/gkj112Search in Google Scholar PubMed
 Wang J, Meng X, Dobrovolskaya OB, Orlov YL, Chen M. Non-coding RNAs and their roles in stress response in plants. Genomics Proteomics Bioinformatics 2017;15:301–12.10.1016/j.gpb.2017.01.007Search in Google Scholar PubMed
 Ponjavic J, Ponting CP, Lunter G. Functionality or transcriptional noise? Evidence for selection within long noncoding RNAs. Genome Res 2007;17:556–65.10.1101/gr.6036807Search in Google Scholar PubMed
 Cerutti P, Holt JW, Miller N. Detection and determination of 5, 6-dihydrouridine and 4-thiouridine in transfer ribonucleic acid from different sources. J Mol Biol 1968;34:505–18.10.1016/0022-2836(68)90176-9Search in Google Scholar PubMed
 Schaefer M, Pollex T, Hanna K, Tuorto F, Meusburger M, Helm M, et al. RNA methylation by Dnmt2 protects transfer RNAs against stress-induced cleavage. Genes Dev 2010;24:1590–5.10.1101/gad.586710Search in Google Scholar PubMed
 Ono M, Scott MS, Yamada K, Avolio F, Barton GJ, Lamond AI. Identification of human miRNA precursors that resemble box C/D snoRNAs. Nucleic Acids Res 2011;39:3879–91.10.1093/nar/gkq1355Search in Google Scholar
 He X, Chen X, Zhang X, Duan X, Pan T, Hu Q, et al. An Lnc RNA (GAS5)/SnoRNA-derived piRNA induces activation of TRAIL gene by site-specifically recruiting MLL/COMPASS-like complexes. Nucleic Acids Res 2015;43:3712–25.10.1093/nar/gkv214Search in Google Scholar PubMed
 Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993;75:843–54.10.1016/0092-8674(93)90529-YSearch in Google Scholar
 Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, et al. The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature 2000;403:901–6.10.1038/35002607Search in Google Scholar PubMed
 Waterhouse PM, Graham MW, Wang MB. Virus resistance and gene silencing in plants can be induced by simultaneous expression of sense and antisense RNA. Proc Natl Acad Sci USA 1998;95:13959–64.10.1073/pnas.95.23.13959Search in Google Scholar
 Kim JK, Gabel HW, Kamath RS, Tewari M, Pasquinelli A, Rual J-F, et al. Functional genomic analysis of RNA interference in C. elegans. Science 2005;308:1164–7.10.1126/science.1109267Search in Google Scholar
 Vagin VV, Sigova A, Li C, Seitz H, Gvozdev V, Zamore PD. A distinct small RNA pathway silences selfish genetic elements in the germline. Science 2006;313:320–4.10.1126/science.1129333Search in Google Scholar PubMed
 Martinez VD, Vucic EA, Thu KL, Hubaux R, Enfield KS, Pikor LA, et al. Unique somatic and malignant expression patterns implicate PIWI-interacting RNAs in cancer-type specific biology. Sci Rep 2015;5:10423.10.1038/srep10423Search in Google Scholar PubMed
 Kim T-K, Hemberg M, Gray JM, Costa AM, Bear DM, Wu J, et al. Widespread transcription at neuronal activity-regulated enhancers. Nature 2010;465:182.10.1038/nature09033Search in Google Scholar PubMed
 Moghe GD, Lehti-Shiu MD, Seddon AE, Yin S, Chen Y, Juntawong P, et al. Characteristics and significance of intergenic polyadenylated RNA transcription in Arabidopsis. Plant Physiol 2013;161:210–24.10.1104/pp.112.205245Search in Google Scholar PubMed
 Legnini I, Morlando M, Mangiavacchi A, Fatica A, Bozzoni I. A feedforward regulatory loop between HuR and the long noncoding RNA linc-MD1 controls early phases of myogenesis. Mol Cell 2014;53:506–14.10.1016/j.molcel.2013.12.012Search in Google Scholar
 Loraine AE, McCormick S, Estrada A, Patel K, Qin P. RNA-seq of Arabidopsis pollen uncovers novel transcription and alternative splicing. Plant Physiol 2013;162:1092–109.10.1104/pp.112.211441Search in Google Scholar PubMed
 Liu J, Jung C, Xu J, Wang H, Deng S, Bernad L, et al. Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis. Plant Cell 2012;24:4333–45.10.1105/tpc.112.102855Search in Google Scholar PubMed
 Wyers F, Rougemaille M, Badis G, Rousselle J-C, Dufour M-E, Boulay J, et al. Cryptic pol II transcripts are degraded by a nuclear quality control pathway involving a new poly (A) polymerase. Cell 2005;121:725–37.10.1016/j.cell.2005.04.030Search in Google Scholar PubMed
 Natoli G, Andrau J-C. Noncoding transcription at enhancers: general principles and functional models. Annu Rev Genet 2012;46:1–19.10.1146/annurev-genet-110711-155459Search in Google Scholar PubMed
 Lam MT, Cho H, Lesch HP, Gosselin D, Heinz S, Tanaka-Oishi Y, et al. Rev-Erbs repress macrophage gene expression by inhibiting enhancer-directed transcription. Nature 2013;498:511–5.10.1038/nature12209Search in Google Scholar PubMed
 Andersson R, Gebhard C, Miguel-Escalada I, Hoof I, Bornholdt J, Boyd M, et al. An atlas of active enhancers across human cell types and tissues. Nature 2014;507:455–61.10.1038/nature12787Search in Google Scholar PubMed
 Li W, Notani D, Ma Q, Tanasa B, Nunez E, Chen AY, et al. Functional roles of enhancer RNAs for oestrogen-dependent transcriptional activation. Nature 2013;498:516–20.10.1038/nature12210Search in Google Scholar PubMed
 Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PloS One 2012;7:e30733.10.1371/journal.pone.0030733Search in Google Scholar PubMed
 Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013;19:141–57.10.1261/rna.035667.112Search in Google Scholar PubMed
 Huang M, Zhong Z, Lv M, Shu J, Tian Q, Chen J. Comprehensive analysis of differentially expressed profiles of lncRNAs and circRNAs with associated co-expression and ceRNA networks in bladder carcinoma. Oncotarget 2016;7:47186–200.10.18632/oncotarget.9706Search in Google Scholar PubMed
 Hafner M, Landthaler M, Burger L, Khorshid M, Hausser J, Berninger P, et al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell 2010;141:129–41.10.1016/j.cell.2010.03.009Search in Google Scholar PubMed
 Helwak A, Kudla G, Dudnakova T, Tollervey D. Mapping the human miRNA interactome by CLASH reveals frequent noncanonical binding. Cell 2013;153:654–65.10.1016/j.cell.2013.03.043Search in Google Scholar PubMed
 Koo J, Zhang J, Chaterji S. Tiresias: Context-sensitive approach to decipher the presence and strength of MicroRNA regulatory interactions. Theranostics 2018;8:277–91.10.7150/thno.22065Search in Google Scholar PubMed
 Sacar Demirci MD, Yousef M, Allmer J. Computational Prediction of Functional MicroRNA-mRNA Interactions. Methods in molecular biology 2019;1912:175–96.10.1007/978-1-4939-8982-9_7Search in Google Scholar
 Wang QS, Shi LL, Sun F, Zhang YF, Chen RW, Yang SL, et al. High Expression of ANXA2 pseudogene ANXA2P2 promotes an aggressive phenotype in hepatocellular carcinoma. Dis Markers 2019;2019:9267046.10.1155/2019/9267046Search in Google Scholar PubMed
 Chen G, Wang Q, Yang Q, Li Z, Du Z, Ren M, et al. Circular RNAs hsa_circ_0032462, hsa_circ_0028173, hsa_circ_0005909 are predicted to promote CADM1 expression by functioning as miRNAs sponge in human osteosarcoma. PLoS One 2018;13:e0202896.10.1371/journal.pone.0202896Search in Google Scholar PubMed
 Lu J, Zhang P-Y, Li P, Xie J-W, Wang J-B, Lin J-X, et al. Circular RNA hsa_circ_0001368 suppresses the progression of gastric cancer by regulating miR-6506–5p/FOXO3 axis. Biochem Biophys Res Commun 2019;512:29–33.10.1016/j.bbrc.2019.02.111Search in Google Scholar PubMed
 Wang S, Huang J, Lyu H, Lee C, Tan J, Wang J, et al. Functional cooperation of miR-125a, miR-125b, and miR-205 in entinostat-induced downregulation of erbB2/erbB3 and apoptosis in breast cancer cells. Cell Death Dis 2013;4:e556.10.1038/cddis.2013.79Search in Google Scholar PubMed
 Grossi E, Sanchez Y, Huarte M. Expanding the p53 regulatory network: LncRNAs take up the challenge. Biochim Biophys Acta Gene Regul Mech 2016;1859:200–8.10.1016/j.bbagrm.2015.07.011Search in Google Scholar
 Li Y, Chen J, Zhang J, Wang Z, Shao T, Jiang C, et al. Construction and analysis of lncRNA-lncRNA synergistic networks to reveal clinically relevant lncRNAs in cancer. Oncotarget 2015;6:25003–16.10.18632/oncotarget.4660Search in Google Scholar PubMed
 Fejes-Toth K, Sotirova V, Sachidanandam R, Assaf G, Hannon GJ, Kapranov P, et al. Post-transcriptional processing generates a diversity of 5′-modified long and short RNAs: Affymetrix/Cold Spring Harbor Laboratory ENCODE Transcriptome Project. Nature 2009;457:1028–32.10.1038/nature07759Search in Google Scholar PubMed
 West JA, Davis CP, Sunwoo H, Simon MD, Sadreyev RI, Wang PI, et al. The long noncoding RNAs NEAT1 and MALAT1 bind active chromatin sites. Mol Cell 2014;55:791–802.10.1016/j.molcel.2014.07.012Search in Google Scholar
 Blank-Giwojna A, Postepska-Igielska A, Grummt I. lncRNA KHPS1 activates a poised enhancer by triplex-dependent recruitment of epigenomic regulators. Cell Rep 2019;26:2904–15. e4.10.1016/j.celrep.2019.02.059Search in Google Scholar PubMed
 Licatalosi DD, Mele A, Fak JJ, Ule J, Kayikci M, Chi SW, et al. HITS-CLIP yields genome-wide insights into brain alternative RNA processing. Nature 2008;456:464–9.10.1038/nature07488Search in Google Scholar PubMed
 McHugh CA, Chen C-K, Chow A, Surka CF, Tran C, McDonel P, et al. The Xist lncRNA interacts directly with SHARP to silence transcription through HDAC3. Nature 2015;521:232–6.10.1038/nature14443Search in Google Scholar PubMed
 Du WW, Fang L, Yang W, Wu N, Awan FM, Yang Z, et al. Induction of tumor apoptosis through a circular RNA enhancing Foxo3 activity. Cell Death Differ 2017;24:357–70.10.1038/cdd.2016.133Search in Google Scholar PubMed
 Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature 2010;465:1033–8.10.1038/nature09144Search in Google Scholar PubMed
 Karreth FA, Tay Y, Perna D, Ala U, Tan SM, Rust AG, et al. In vivo identification of tumor-suppressive PTEN ceRNAs in an oncogenic BRAF-induced mouse model of melanoma. Cell 2011;147:382–95.10.1016/j.cell.2011.09.032Search in Google Scholar
 Lee DY, Jeyapalan Z, Fang L, Yang J, Zhang Y, Yee AY, et al. Expression of versican 3′-untranslated region modulates endogenous microRNA functions. PloS One 2010;5:e13599.10.1371/journal.pone.0013599Search in Google Scholar PubMed
 Lee DY, Shatseva T, Jeyapalan Z, Du WW, Deng Z, Yang BB. A 3′-untranslated region (3′ UTR) induces organ adhesion by regulating miR-199a* functions. PloS One 2009;4:e4527.10.1371/journal.pone.0004527Search in Google Scholar PubMed
 Fang L, Du WW, Yang X, Chen K, Ghanekar A, Levy G, et al. Versican 3′-untranslated region (3′-UTR) functions as a ceRNA in inducing the development of hepatocellular carcinoma by regulating miRNA activity. FASEB J 2013;27:907–19.10.1096/fj.12-220905Search in Google Scholar PubMed
 Marques AC, Tan J, Lee S, Kong L, Heger A, Ponting CP. Evidence for conserved post-transcriptional roles of unitary pseudogenes and for frequent bifunctionality of mRNAs. Genome Biol 2012;13:R102.10.1186/gb-2012-13-11-r102Search in Google Scholar PubMed
 Lai F, Orom UA, Cesaroni M, Beringer M, Taatjes DJ, Blobel GA, et al. Activating RNAs associate with Mediator to enhance chromatin architecture and transcription. Nature 2013;494:497–501.10.1038/nature11884Search in Google Scholar PubMed
 Halley P, Kadakkuzha BM, Faghihi MA, Magistri M, Zeier Z, Khorkova O, et al. Regulation of the apolipoprotein gene cluster by a long noncoding RNA. Cell Rep 2014;6:222–30.10.1016/j.celrep.2013.12.015Search in Google Scholar PubMed
 Simon MD, Pinter SF, Fang R, Sarma K, Rutenberg-Schoenberg M, Bowman SK, et al. High-resolution Xist binding maps reveal two-step spreading during X-chromosome inactivation. Nature 2013;504:465–9.10.1038/nature12719Search in Google Scholar PubMed
 Engreitz JM, Pandya-Jones A, McDonel P, Shishkin A, Sirokman K, Surka C, et al. The Xist lncRNA exploits three-dimensional genome architecture to spread across the X chromosome. Science 2013;341:1237973.10.1126/science.1237973Search in Google Scholar PubMed
 Hu Q, Kwon YS, Nunez E, Cardamone MD, Hutt KR, Ohgi KA, et al. Enhancing nuclear receptor-induced transcription requires nuclear motor and LSD1-dependent gene networking in interchromatin granules. Proc Natl Acad Sci USA 2008;105:19199–204.10.1073/pnas.0810634105Search in Google Scholar
 Nunez E, Fu XD, Rosenfeld MG. Nuclear organization in the 3D space of the nucleus – cause or consequence? Curr Opin Genet Dev 2009;19:424–36.10.1016/j.gde.2009.07.005Search in Google Scholar PubMed
 Yoon JH, Abdelmohsen K, Srikantan S, Yang X, Martindale JL, De S, et al. LincRNA-p21 suppresses target mRNA translation. Mol Cell 2012;47:648–55.10.1016/j.molcel.2012.06.027Search in Google Scholar PubMed
 You X, Vlatkovic I, Babic A, Will T, Epstein I, Tushev G, et al. Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat Neurosci 2015;18:603–10.10.1038/nn.3975Search in Google Scholar PubMed
© 2019, Peijing Zhang et al., published by Walter de Gruyter GmbH, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 Public License.