Search

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

The accurate classification of non-coding RNA (ncRNA) sequences is pivotal for advanced non-coding genome annotation and analysis, a fundamental aspect of genomics that facilitates understanding of ncRNA functions and regulatory mechanisms in various biological processes. While traditional machine learning approaches have been employed for distinguishing ncRNA, these often necessitate extensive feature engineering. Recently, deep learning algorithms have provided advancements in ncRNA classification. This study presents BioDeepFuse, a hybrid deep learning framework integrating convolutional neural networks (CNN) or bidirectional long short-term memory (BiLSTM) networks with handcrafted features for enhanced accuracy. This framework employs a combination of k- mer one-hot, k- mer dictionary, and feature extraction techniques for input representation. Extracted features, when embedded into the deep network, enable optimal utilization of spatial and sequential nuances of ncRNA sequences. Using benchmark datasets and real-world RNA samples from bacterial organisms, we evaluated the performance of BioDeepFuse. Results exhibited high accuracy in ncRNA classification, underscoring the robustness of our tool in addressing complex ncRNA sequence data challenges. The effective melding of CNN or BiLSTM with external features heralds promising directions for future research, particularly in refining ncRNA classifiers and deepening insights into ncRNAs in cellular processes and disease manifestations. In addition to its original application in the context of bacterial organisms, the methodologies and techniques integrated into our framework can potentially render BioDeepFuse effective in various and broader domains.

The Keystone Symposium ' Small Regulatory RNAs: From Bench to Bedside ' was held in Santa Fe, New Mexico from May 1-4, 2022. The symposium was organized by Frank J. Slack, Jörg Vogel, Ivan Martinez and Karyn Schmidt, and brought together scientists working in noncoding RNA biology, therapeutics, and technologies to address mechanistic questions about small regulatory RNAs and facilitate translation of these findings into clinical applications. The conference addressed four specific aims: Aim 1 . Focus on the exciting biology of small regulatory RNAs, highlighting the best current research into the role that small RNAs play in fundamental biological processes; Aim 2 . Focus on the latest efforts to harness the power of these RNAs as agents in the fight against disease and provide the basic understanding that will drive the invention of powerful clinical tools; Aim 3 . Attract leaders from both academia and industry working in small RNAs to one place for critical discussions that will advance the field and accelerate the bench to bedside use of this technology; Aim 4 . Provide a stimulating environment where students, postdoctoral researchers and junior investigators, along with scientists from Biotechnology and Pharmaceutical companies specializing in small regulatory RNAs, can present and discuss their research with the best minds in the field.

Y RNAs (84-112 nt) are non-coding RNAs transcribed by RNA polymerase III and are characterized by a distinctive secondary structure. Human Y RNAs interact with the autoimmune proteins SSB and RO60 that together form a ribonucleoprotein (RNP) complex termed RoRNP and Y RNAs also perform regulatory roles in DNA and RNA replication and stability, which has major implications for diseases including cancer. During cellular stress and apoptosis, Y RNAs are cleaved into 3' and 5' end fragments termed Y RNA-derived small RNAs (ysRNAs). Although some ysRNA functions in stress, apoptosis and cancer have been reported, their fundamental biogenesis has not been described. Here we report that 3' end RNY5 cleavage is structure dependent. In high throughput mutagenesis experiments, cleavage occurred between the 2 nd and 3 rd nt above a double stranded stem comprising high GC content. We demonstrate that an internal loop above stem S3 is critical for producing 3' end ysRNAs (31 nt) with mutants resulting in longer or no ysRNAs. We show a UGGGU sequence motif at position 22 of RNY5 is critical for producing 5' end ysRNAs (22-25 nt). We show that intact RO60 is critical for ysRNA biogenesis. We conclude that ribonuclease L (RNASEL) contributes to Y RNA cleavage in mouse embryonic fibroblasts but is not the only endoribonuclease important in human cells.

Ribosome-associated non-coding RNAs (rancRNAs) have been recognized as an emerging class of regulatory molecules capable of fine-tuning translation in all domains of life. RancRNAs are ideally suited for allowing a swift response to changing environments and are therefore considered pivotal during the first wave of stress adaptation. Previously, we identified an mRNA-derived 18 nucleotides long rancRNA (rancRNA_18) in Saccharomyces cerevisiae that rapidly downregulates protein synthesis during hyperosmotic stress. However, the molecular mechanism of action remained enigmatic. Here, we combine biochemical, genetic, transcriptome-wide and structural evidence, thus revealing rancRNA_18 as global translation inhibitor by targeting the E-site region of the large ribosomal subunit. Ribosomes carrying rancRNA_18 possess decreased affinity for A-site tRNA and impaired structural dynamics. Cumulatively, these discoveries reveal the mode of action of a rancRNA involved in modulating protein biosynthesis at a thus far unequalled precision.

Comparative sequence analysis methods are highly effective for uncovering novel classes of structured noncoding RNAs (ncRNAs) from bacterial genomic DNA sequence datasets. Previously, we developed a computational pipeline to more comprehensively identify structured ncRNA representatives from individual bacterial genomes. This search process exploits the fact that genomic regions serving as templates for the transcription of structured RNAs tend to be present in longer than average noncoding 'intergenic regions' (IGRs) that are enriched in G and C nucleotides compared to the remainder of the genome. In the present study, we apply this computational pipeline to identify structured ncRNA candidates from 26 diverse bacterial species. Numerous novel structured ncRNA motifs were discovered, including several riboswitch candidates, one whose ligand has been identified and others that have yet to be experimentally validated. Our findings support recent predictions that hundreds of novel ribo-switch classes and other ncRNAs remain undiscovered among the limited number of bacterial species whose genomes have been completely sequenced.

Mitochondrial noncoding RNAs (mt-ncRNAs) include noncoding RNAs inside the mitochondria that are transcribed from the mitochondrial genome or nuclear genome, and noncoding RNAs transcribed from the mitochondrial genome that are transported to the cytosol or nucleus. Recent findings have revealed that mt-ncRNAs play important roles in not only mitochondrial functions, but also other cellular activities. This review proposes a classification of mt-ncRNAs and outlines the emerging understanding of mitochondrial circular RNAs (mt-circRNAs), mitochondrial microRNAs (mitomiRs), and mitochondrial long noncoding RNAs (mt-lncRNAs), with an emphasis on their identification and functions.

The subcellular localization of RNAs correlates with their function and how they are regulated. Most protein-coding mRNAs are exported into the cytoplasm for protein synthesis, while some mRNA species, long noncoding RNAs, and some regulatory element-associated unstable transcripts tend to be retained in the nucleus, where they function as a regulatory unit and/or are regulated by nuclear surveillance pathways. While the mechanisms regulating mRNA export and localization have been well summarized, the mechanisms governing nuclear retention of RNAs, especially of noncoding RNAs, are seldomly reviewed. In this review, we summarize recent advances in the mechanistic study of RNA nuclear retention, especially for noncoding RNAs, from the angle of cis -acting elements embedded in RNA transcripts and their interaction with trans -acting factors. We also try to illustrate the general principles of RNA nuclear retention and we discuss potential areas for future investigation.

As one of the most common forms of RNA modification, N6-methyladenosine (m6A) RNA modification has attracted increasing research interest in recent years. This reversible RNA modification added a new dimension to the post-transcriptional regulation of gene expression. In colorectal cancer (CRC), the role of m6A modification has been extensively studied, not only on mRNAs but also on non-coding RNAs (ncRNAs). In the present review, we depicted the role of m6A modification in CRC, systematically elaborate the interaction between m6A modification and regulatory ncRNAs in function and mechanism. Moreover, we discussed the potential applications in clinical.

Regulation of translation is essential for the diverse biological processes involved in development. Particularly, mammalian oocyte development requires the precisely controlled translation of maternal transcripts to coordinate meiotic and early embryo progression while transcription is silent. It has been recently reported that key components of mRNA translation control are short and long noncoding RNAs (ncRNAs). We found that the ncRNA Brain cytoplasmic 1 (BC1)  has a role in the fully grown germinal vesicle (GV) mouse oocyte, where is highly expressed in the cytoplasm associated with polysomes. Overexpression of  BC1  in GV oocyte leads to a minute decrease in global translation with a significant reduction of specific mRNA translation via interaction with the Fragile X Mental Retardation Protein (FMRP).  BC1  performs a repressive role in translation only in the GV stage oocyte without forming FMRP or Poly(A) granules. In conclusion,  BC1  acts as the translational repressor of specific mRNAs in the GV stage via its binding to a subset of mRNAs and physical interaction with FMRP. The results reported herein contribute to the understanding of the molecular mechanisms of developmental events connected with maternal mRNA translation.

Non-coding RNAs (ncRNAs) are powerful regulators of gene expression but medium-sized (50-300 nts in length) ncRNAs (msRNAs) are barely picked-up precisely by RNA-sequencing. Here we describe msRNA-sequencing (msRNAseq), a modified protocol that associated with a computational analyses pipeline identified about ~1800 msRNA loci, including over 300 putatively novel msRNAs, in human and murine cells. We focused on the identification and initial characterization of three POLIII-derived transcripts. The validation of these uncharacterized msRNAs identified an ncRNA in antisense orientation from the POLR3E locus transcribed by POLIII. This msRNA, termed POLAR (POLR3E Antisense R NA), has a strikingly short half-life, localizes to paraspeckles (PSPs) and associates with PSP-associated proteins indicating that msRNAseq identifies functional msRNAs. Thus, our analyses will pave the way for analysing the roles of msRNAs in cells, development and diseases.

The blood-brain barrier (BBB), which controls permeability into and out of the nervous system, is a tightly connected, structural, and functional separation between the central nervous system (CNS) and circulating blood. CNS diseases, such as Alzheimer's disease, multiple sclerosis, traumatic brain injury, stroke, meningitis, and brain cancers, often develop with the increased BBB permeability and further leads to irreversible CNS injury. Non-coding RNAs (ncRNAs) are functional RNA molecules that generally lack the coding abilities but can actively regulate the mRNA expression and function through different mechanisms. Various types of ncRNAs, including microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs), are highly expressed in brain microvascular endothelial cells and are potential mediators of BBB permeability. Here, we summarized the recent research progress on miRNA, lncRNA, and circRNA roles regulating the BBB permeability in different CNS diseases. Understanding how these ncRNAs affect the BBB permeability shall provide important therapeutic insights into the prevention and control of the BBB dysfunction.

Epstein-Barr virus (EBV) was the first human cancer-causing virus to be discovered over fifty years ago. Given its relatively large genome size for a virus and hence the capacity to store more than mere protein-coding information, EBV also harbours genetic material for producing an array of distinct noncoding (nc)RNAs. The double-stranded nature of its DNA genome allows the utilization of the whole gamut of ncRNA types, including both RNA polymerase II and III transcripts, in devising a sophisticated strategy to ensure its replication upon infection in host cells and evasion of host immune responses. Owing to the development of sensitive technologies in recent years, mostly entailing next-generation sequencing, the list of ncRNA types generated by EBV has expanded now to include two RNAs (EBER1 and EBER2) best categorized as long ncRNAs, dozens of microRNAs, one small nucleolar RNA, stable intronic sequence RNAs, and the most recently discovered circular RNAs. With the application of cutting-edge technology, the molecular mechanisms of some of these noncoding transcripts are beginning to emerge, while others remain yet to be elucidated. As viruses often take advantage of existing molecular pathways established by the host, it is likely that further novel concepts of the greatly unexplored noncoding world can be learned from studying the many EBV ncRNAs.

La and La-related proteins (LARPs) are characterized by a common RNA interaction platform termed the La module. This structural hallmark allows LARPs to pervade various aspects of RNA biology. The metazoan LARP7 protein binds to the 7SK RNA as part of a 7SK small nuclear ribonucleoprotein (7SK snRNP), which inhibits the transcriptional activity of RNA polymerase II (Pol II). Additionally, recent findings revealed unanticipated roles of LARP7 in the assembly of other RNPs, as well as in the modification, processing and cellular transport of RNA molecules. Reduced levels of functional LARP7 have been linked to cancer and Alazami syndrome, two seemingly unrelated human diseases characterized either by hyperproliferation or growth retardation. Here, we review the intricate regulatory networks centered on LARP7 and assess how malfunction of these networks may relate to the etiology of LARP7-linked diseases.

Bacterial 6S RNA regulates transcription via binding to the active site of RNA polymerase holoenzymes. 6S RNA has been identified in the majority of bacteria, in most cases encoded by a single gene. Firmicutes including Bacillus subtilis encode two 6S RNA paralogs, 6S-1 and 6S-2 RNA. Hypothesizing that the regulatory role of 6S RNAs may be particularly important under natural, constantly changing environmental conditions, we constructed 6S RNA deletion mutants of the undomesticated B. subtilis wild-type strain NCIB 3610. We observed a strong phenotype for the ∆6S-2 RNA strain that showed increased biofilm formation on solid media and the ability to form surface-attached biofilms in liquid culture. This phenotype remained undetected in derived laboratory strains (168, PY79) that are defective in biofilm formation. Quantitative RT-PCR data revealed transcriptional upregulation of biofilm marker genes such as tasA, epsA and bslA in the ∆6S-2 RNA strain, particularly during transition from exponential to stationary growth phase. Salt stress, which blocks sporulation at a very early stage, was found to override the derepressed biofilm phenotype of the ∆6S-2 RNA strain. Furthermore, the ∆6S-2 RNA strain showed retarded swarming activity and earlier spore formation. Finally, the ∆6S-1&2 RNA double deletion strain showed a prolonged lag phase of growth under oxidative, high salt and alkaline stress conditions, suggesting that the interplay of both 6S RNAs in B. subtilis optimizes and fine-tunes transcriptomic adaptations, thereby contributing to the fitness of B. subtilis under the unsteady and temporarily harsh conditions encountered in natural habitats.

The findings that an RNA is not necessarily either coding or non-coding, or that a precursor RNA can produce different types of mature RNAs, whether coding or non-coding, long or short, have challenged the dichotomous view of the RNA world almost 15 years ago. Since then, and despite an increasing number of studies, the diversity of information that can be conveyed by RNAs is rarely searched for, and when it is known, it remains largely overlooked in further functional studies. Here, we provide an update with prominent examples of multiple functions that are carried by the same RNA or are produced by the same precursor RNA, to emphasize their biological relevance in most living organisms. An important consequence is that the overall function of their locus of origin results from the balance between various RNA species with distinct functions and fates. The consideration of the molecular basis of this multiplicity of information is obviously crucial for downstream functional studies when the targeted functional molecule is often not the one that is believed.

Despite advanced clinical treatments, mortality in patients with metastatic colorectal cancer (CRC) remains high. Three critical determinants in CRC progression include the epithelial proliferation checkpoints, epithelial-to-mesenchymal transition (EMT) and inflammatory cytokines in the tumour microenvironment. Genes involved in these three processes are regulated at the transcriptional and post-transcriptional level. Recent studies revealed previously unappreciated roles of non-coding ribonucleic acids (ncRNAs) in modulating the proliferation checkpoints, EMT, and inflammatory gene expression in CRC. In this review, we will discuss the mechanisms underlying the roles of ncRNAs in CRC as well as examine future perspectives in this field. Better understanding of ncRNA biology will provide novel targets for future therapeutic development.

Enhancers are distal genomic elements critical for gene regulation and cell identify control during development and diseases. Many human cancers were found to associate with enhancer malfunction, due to genetic and epigenetic alterations, which in some cases directly drive tumour growth. Conventionally, enhancers are known to provide DNA binding motifs to recruit transcription factors (TFs) and to control target genes. However, recent progress found that most, if not all, active enhancers pervasively transcribe noncoding RNAs that are referred to as enhancer RNAs (eRNAs). Increasing evidence points to functional roles of at least a subset of eRNAs in gene regulation in both normal and cancer cells, adding new insights into the action mechanisms of enhancers. eRNA expression was observed to be widespread but also specific to tumour types and individual patients, serving as opportunities to exploit them as potential diagnosis markers or therapeutic targets. In this review, we discuss the brief history of eRNA research, their functional mechanisms and importance in cancer gene regulation, as well as their therapeutic and diagnostic values in cancer. We propose that further studies of eRNAs in cancer will offer a promising 'eRNA targeted therapy' for human cancer intervention.

Non-coding RNAs occupy a significant fraction of the human genome. Their biological significance is backed up by a plethora of emerging evidence. One of the most robust approaches to demonstrate non-coding RNA's biological relevance is through their prognostic value. Using the rich gene expression data from The Cancer Genome Altas (TCGA), we designed Advanced Expression Survival Analysis (AESA), a web tool which provides several novel survival analysis approaches not offered by previous tools. In addition to the common single-gene approach, AESA computes the gene expression composite score of a set of genes for survival analysis and utilizes permutation test or cross-validation to assess the significance of log-rank statistic and the degree of over-fitting. AESA offers survival feature selection with post-selection inference and utilizes expanded TCGA clinical data including overall, disease-specific, disease-free, and progression-free survival information. Users can analyse either protein-coding or non-coding regions of the transcriptome. We demonstrated the effectiveness of AESA using several empirical examples. Our analyses showed that non-coding RNAs perform as well as messenger RNAs in predicting survival of cancer patients. These results reinforce the potential prognostic value of non-coding RNAs. AESA is developed as a module in the freely accessible analysis suite MutEx. Abbreviation: ACC: Adrenocortical Carcinoma (n = 92); BLCA: Bladder Urothelial Carcinoma (n = 412); BRCA: Breast Invasive Carcinoma (n = 1098); CESC: Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (n = 307); CHOL: Cholangiocarcinoma (n = 51); COAD: Colon Adenocarcinoma (n = 461); DLBC: Lymphoid Neoplasm Diffuse Large B-cell Lymphoma (n = 58); ESCA: Oesophageal Carcinoma (n = 185); GBM: Glioblastoma Multiforme (n = 617); HNSC: Head and Neck Squamous Cell Carcinoma (n = 528); KICH: Kidney Chromophobe (n = 113); KIRC: Kidney Renal Clear Cell Carcinoma (n = 537); KIRP: Kidney Renal Papillary Cell Carcinoma (n = 291); LAML: Acute Myeloid Leukaemia (n = 200); LGG: Brain Lower Grade Glioma (n = 516); LIHC: Liver Hepatocellular Carcinoma (n = 377); LUAD: Lung Adenocarcinoma (n = 585); LUSC: Lung Squamous Cell Carcinoma (n = 504); MESO: Mesothelioma (n = 87); OV: Ovarian Serous Cystadenocarcinoma (n = 608) PAAD: Pancreatic Adenocarcinoma (n = 185); PCPG: Pheochromocytoma and Paraganglioma (n = 179); PRAD: Prostate Adenocarcinoma (n = 500); READ: Rectum Adenocarcinoma (n = 172); SARC: Sarcoma (n = 261); SKCM: Skin Cutaneous Melanoma (n = 470); STAD: Stomach Adenocarcinoma (n = 443); TGCT: Testicular Germ Cell Tumours (n = 150); THCA: Thyroid Carcinoma (n = 507) THYM: Thymoma (n = 124); UCEC: Uterine Corpus Endometrial Carcinoma (n = 560); UCS: Uterine Carcinosarcoma (n = 57); UVM: Uveal Melanoma (n = 80).

Recent transcriptome-wide studies have identified a diverse pool of transfer RNA (tRNA)-derived RNAs or tRNA-derived fragments (tRFs). Some of these RNAs have been demonstrated to be functional and involved in multiple biological processes ranging from the regulation of gene expression to transgenerational epigenetic inheritance. Post-transcriptional maturation of tRNAs includes various processing events including extensive decoration by various RNA modifications, which are required for correct tRNA folding and stability. Moreover, tRNA modifications determine the pattern and specificity of tRNA cleavage. The major drawbacks of many studies in the field of tRFs are that most of them used synthetic RNAs that closely mimic endogenous tRFs in their sequence, yet lack RNA modification that is found in vivo . Here, we developed a simple method to isolate tRNA-derived stress-induced RNAs (tiRNAs), a specific subset of tRFs. Our approach is scalable, cost-effective and relies on the purification of individual tiRNAs based on a sequence-specific RNA/DNA isolation technique using DNA probes. Our method facilitates functional studies of tiRNAs by addressing how physiological RNA modifications within tRNA fragments affect their biological activities. Here, we report pilot functional studies on selected endogenous tiRNAs, namely tiRNA Ala and tiRNA Gly . We show that natural 5'-tiRNA Ala molecules assemble into G-quadruplex structures, and endogenous 5'-tiRNA Gly possesses translation inhibition activity.

Transfer RNA-derived fragments (tRFs) are a new class of small non-coding RNAs whose biological roles in cancers are not well understood. Emerging evidence suggests that tRFs are involved in gene regulation at multiple levels. In this study, we constructed an integrative database, OncotRF (http://bioinformatics.zju.edu.cn/OncotRF), for in silico exploration of tRF functions, and identification of diagnostic and prognostic biomarkers in cancers. The database contains an analysis pipeline for tRF identification and characterization, analysis results of 11,211 small RNA sequencing samples and 8,776 RNA sequencing samples, and clinicopathologic annotation data from The Cancer Genome Atlas (TCGA). The results include: tRF identification and quantification across 33 cancers, abnormally expressed tRFs and genes, tRF-gene correlations, tRF-gene networks, survival analyses, and tRF-related functional enrichment analyses. Users are also able to identify differentially expressed tRFs, predict their functions, and assess the relevance of the tRF expression levels to the clinical outcome according to user-defined groups. Additionally, an online Kaplan-Meier plotter is available in OncotRF for plotting survival curves according to user-defined groups. OncotRF will be a valuable online database and functional annotation tool for researchers studying the roles, functions, and mechanisms of tRFs in human cancers.

Extracellular vesicles (EVs) are cell-derived nanoparticles that act as natural carriers of nucleic acids between cells. They offer advantages as delivery vehicles for therapeutic nucleic acids such as small RNAs. Loading of desired nucleic acids into EVs can be achieved by electroporation or transfection once purified. An attractive alternative is to transfect cells with the desired small RNAs and harness the cellular machinery for RNA sorting into the EVs. This possibility has been less explored because cells are believed to secrete only specific RNAs. However, we hypothesized that, even in the presence of selective secretion, concentration-driven RNA sorting to EVs would still be feasible. To show this, we transfected cells with glycine 5' tRNA halves, which we have previously shown to better resist RNases. We then measured their levels in EVs and in recipient cells and found that, in contrast to unstable RNAs of random sequence, these tRNA halves were present in vesicles and in recipient cells in amounts proportional to the concentration of RNA used for transfection. Similar efficiencies were obtained with other stable oligonucleotides of random sequence. Our results demonstrate that RNA stability is a key factor needed to maintain high intracellular concentrations, a prerequisite for efficient non-selective RNA sorting to EVs and delivery to cells. Given that glycine 5' tRNA halves belong to the group of stress-induced tRNA fragments frequently detected in extracellular space and biofluids, we propose that upregulation of extracellular tRNA fragments is consequential to cellular stress and might be involved in intercellular signalling.

Transfer RNA-derived fragments (tRFs) exist in all branches of life. They are involved in RNA degradation, regulation of gene expression, ribosome biogenesis. In archaebacteria, kinetoplastid, yeast, and human cells, they were also shown to regulate translation. In Arabidopsis, the tRFs population fluctuates under developmental or environmental conditions but their functions are yet poorly understood. Here, we show that populations of long (30-35 nt) or short (19-25 nt) tRFs produced from Arabidopsis tRNAs can inhibit in vitro translation of a reporter gene. Analysing a series of oligoribonucleotides mimicking natural tRFs, we demonstrate that only a limited set of tRFs possess the ability to affect protein synthesis. Out of a dozen of tRFs, only two deriving from tRNA Ala (AGC) and tRNA Asn (GUU) strongly attenuate translation in vitro . Contrary to human tRF(Ala), the 4 Gs present at the 5' extremity of Arabidopsis tRF(Ala) are not implicated in this inhibition while the G 18 and G 19 residues are essential. Protein synthesis inhibition by tRFs does not require complementarity with the translated mRNA but, having the capability to be associated with polyribosomes, tRFs likely act as general modulation factors of the translation process in plants.

tiRNAs are small non-coding RNAs produced when tRNA is cleaved under stress. tRNA methylation modifications has emerged in recent years as important regulators for tRNA structural stability and sensitivity to cleavage and tiRNA generation during stress, however, the specificity and higher regulation of such a process is not fully understood. Alkbh1 is a m 1 A demethylase that leads to destabilization of tRNA and enhanced tRNA cleavage. We examined the impact of Alkbh1 targeting via gene knockdown or overexpression on B35 rat neuroblastoma cell line fate following stresses and on tRNA cleavage. We show that Alkbh1 impact on cell fate and tRNA cleavage is a stress specific process that is impacted by the demethylating capacity of the cellular stress in question. We also show that not all tRNAs are cleaved equally following Alkbh1 manipulation and stress, and that Alkbh1 KD fails to rescue tRNAs from cleavage following demethylating stresses. These findings shed a light on the specificity and higher regulation of tRNA cleavage and should act as a guide for future work exploring the utility of Alkbh1 as a therapeutic target for cancers or ischaemic insult.

Transfer RNA fragments (tRFs) are an emerging class of small RNA molecules derived from mature or precursor tRNAs. They are found across a wide range of organisms and tissues, in small RNA fraction or loaded to Argonaute in numbers comparable to microRNAs. Their functions and mechanisms of action are largely unknown, and results obtained on individual tRFs are often hard to generalize. Here we predicted binding mechanisms and specific target interaction sites of 26 human Argonaute-loaded tRFs of different types using large-scale meta-analyses of available experimental data. Strikingly, our findings matched all interaction sites detected in a recent experimental screen, confirming the validity of our computational approach. Such sites are primarily located on the 5' end of tRFs and often involve additional binding along the tRF length, similar to microRNAs. Indicative of multiple layers of regulation, diverse regulatory non-coding RNAs comprised a third of the tRF targets, with the rest being protein-coding transcripts. In the latter, coding sequence and 3' UTRs were the likely primary target regions, although we observed interactions of tRFs with 5' UTRs. Another novel phenomenon we report, a large number of putative interactions between tRFs and introns, is compatible with the roles of Argonaute in the nucleus. Further, observed tRF-intron binding modes suggest a mechanism of interaction of tRFs with Argonaute-dependent introns, and we predict here >20 candidate introns of this type. Taken together, these results present tRFs as regulatory molecules with a rich functional spectrum.

Direct single-molecule sequencing of full-length transcripts allows efficient identification of gene isoforms, which is apt to alternative splicing (AS), polyadenylation, and long non-coding RNA analyses. However, the identification of gene isoforms and long non-coding RNAs with novel regulatory functions remains challenging, especially for species without a reference genome. Here, we present a comprehensive analysis of a combined long-read and short-read transcriptome sequencing in Camellia japonica . Through a novel bioinformatic pipeline of reverse-tracing the split-sites, we have uncovered 257,692 AS sites from 61,838 transcripts; and 13,068 AS isoforms have been validated by aligning the short reads. We have identified the tissue-specific AS isoforms along with 6,373 AS events that were found in all tissues. Furthermore, we have analysed the polyadenylation (polyA) patterns of transcripts, and found that the preference for polyA signals was different between the AS and non-AS transcripts. Moreover, we have predicted the phased small interfering RNA (phasiRNA) loci through integrative analyses of transcriptome and small RNA sequencing. We have shown that a newly evolved phasiRNA locus from lipoxygenases generated 12 consecutive 21 bp secondary RNAs, which were responsive to cold and heat stress in Camellia . Our studies of the isoform transcriptome provide insights into gene splicing and functions that may facilitate the mechanistic understanding of plants.

Noncoding RNAs (ncRNAs) play critical roles in many critical biological processes and have become a novel class of potential targets and bio-markers for disease diagnosis, therapy, and prognosis. Annotating and analysing ncRNA-disease association data are essential but challenging. Current computational resources lack comprehensive database platforms to consistently interpret and prioritize ncRNA-disease association data for biomedical investigation and application. Here, we present the ncRPheno database platform (http://lilab2.sysu.edu.cn/ncrpheno), which comprehensively integrates and annotates ncRNA-disease association data and provides novel searches, visualizations, and utilities for association identification and validation. ncRPheno contains 482,751 non-redundant associations between 14,494 ncRNAs and 3,210 disease phenotypes across 11 species with supporting evidence in the literature. A scoring model was refined to prioritize the associations based on evidential metrics. Moreover, ncRPheno provides user-friendly web interfaces, novel visualizations, and programmatic access to enable easy exploration, analysis, and utilization of the association data. A case study through ncRPheno demonstrated a comprehensive landscape of ncRNAs dysregulation associated with 22 cancers and uncovered 821 cancer-associated common ncRNAs. As a unique database platform, ncRPheno outperforms the existing similar databases in terms of data coverage and utilities, and it will assist studies in encoding ncRNAs associated with phenotypes ranging from genetic disorders to complex diseases., Abbreviations: APIs: application programming interfaces; circRNA: circular RNA; ECO: Evidence & Conclusion Ontology; EFO: Experimental Factor Ontology; FDR: false discovery rate; GO: Gene Ontology; GWAS: genome wide association studies; HPO: Human Phenotype Ontology; ICGC: International Cancer Genome Consortium; lncRNA: long noncoding RNA; miRNA: micro RNA; ncRNA: noncoding RNA; NGS: next generation sequencing; OMIM: Online Mendelian Inheritance in Man; piRNA: piwi-interacting RNA; snoRNA: small nucleolar RNA; TCGA: The Cancer Genome Atlas.

The tight gene expression regulation controls the development and pathogenesis of human malaria parasite Plasmodium falciparum throughout the complex life cycle. Recent studies have revealed the pervasive nascent transcripts in the genome of P. falciparum , suggesting the existence of a hidden transcriptome involved in the dynamic gene expression. However, the landscape and related biological functions of nascent non-coding RNAs (ns-ncRNAs) are still poorly explored. Here we profiled the transcription dynamics of nascent RNAs by rRNA-depleted and stranded RNA sequencing over the course of 48-h intraerythrocytic developmental cycle (IDC). We identified the genome-wide sources of a total of 2252 ns-ncRNAs, mostly originating from intergenic and untranslated regions of annotated genes. By integrating the nascent RNA abundances with ATAC-seq and ChIP-seq analysis, we uncovered the euchromatic microenvironment surrounding the ns-ncRNA loci, and revealed a positive correlation between ns-ncRNAs and corresponding mRNA abundances. Finally, by gene knock-down strategy, we showed that the cooperation of RNA exosome catalytic subunit PfDis3 and PfMtr4 cofactor played a major role in ns-ncRNAs degradation. Collectively, this study contributes to understanding of the potential roles of short-lived nascent ncRNAs in regulating gene expression in malaria parasites.

DEAD-box RNA-helicases catalyze the reorganization of structured RNAs and the formation of RNP complexes. The cyanobacterium Synechocystis sp. PCC 6803 encodes a single DEAD-box RNA helicase, CrhR (Slr0083), whose expression is regulated by abiotic stresses that alter the redox potential of the photosynthetic electron transport chain, including temperature downshift. Despite its proposed effect on RNA metabolism and its known relevance in cold-stress adaptation, the reported impact of a CrhR knockout on the cold adaption of the transcriptome only identified eight affected genes. Here, we utilized a custom designed microarray to assess the impact of the absence of CrhR RNA helicase activity on the transcriptome, independent of cold stress. CrhR truncation impacts an RNA subset comprising ~10% of the ncRNA and also ~10% of the mRNA transcripts. While equal numbers of mRNAs showed increased as well as decreased abundance, more than 90% of the ncRNAs showed enhanced expression in the absence of CrhR, indicative of a negative effect on ncRNA transcription or stability. We further tested the effect of CrhR on the stability of strongly responding RNAs that identify examples of post-transcriptional and transcriptional regulation. The data suggest that CrhR impacts multiple aspects of RNA metabolism in Synechocystis .

Leishmaniasis is a worldwide public health problem caused by protozoan parasites of the genus Leishmania. Leishmania braziliensis is the most important species responsible for tegumentary leishmaniases in Brazil. An understanding of the molecular mechanisms underlying the success of this parasite is urgently needed. An in-depth study on the modulation of gene expression across the life cycle stages of L. braziliensis covering coding and noncoding RNAs (ncRNAs) was missing and is presented herein. Analyses of differentially expressed (DE) genes revealed that most prominent differences were observed between the transcriptomes of insect and mammalian proliferative forms (6,576 genes). Gene ontology (GO) analysis indicated stage-specific enriched biological processes. A computational pipeline and 5 ncRNA predictors allowed the identification of 11,372 putative ncRNAs. Most of the DE ncRNAs were found between the transcriptomes of insect and mammalian proliferative stages (38%). Of the DE ncRNAs, 295 were DE in all three stages and displayed a wide range of lengths, chromosomal distributions and locations; many of them had a distinct expression profile compared to that of their protein-coding neighbors. Thirty-five putative ncRNAs were submitted to northern blotting analysis, and one or more hybridization-positive signals were observed in 22 of these ncRNAs. This work presents an overview of the L. braziliensis transcriptome and its adjustments throughout development. In addition to determining the general features of the transcriptome at each life stage and the profile of protein-coding transcripts, we identified and characterized a variety of noncoding transcripts. The novel putative ncRNAs uncovered in L. braziliensis might be regulatory elements to be further investigated.

Non-coding Y RNAs and stem-bulge RNAs are homologous small RNAs in vertebrates and nematodes, respectively. They share a conserved function in the replication of chromosomal DNA in these two groups of organisms. However, functional homologues have not been found in insects, despite their common early evolutionary history. Here, we describe the identification and functional characterization of two sbRNAs in Drosophila melanogaster, termed Dm1 and Dm2. The genes coding for these two RNAs were identified by a computational search in the genome of D. melanogaster for conserved sequence motifs present in nematode sbRNAs. The predicted secondary structures of Dm1 and Dm2 partially resemble nematode sbRNAs and show stability in molecular dynamics simulations. Both RNAs are phylogenetically closer related to nematode sbRNAs than to vertebrate Y RNAs. Dm1, but not Dm2 sbRNA is abundantly expressed in D. melanogaster S2 cells and adult flies. Only Dm1, but not Dm2 sbRNA can functionally replace Y RNAs in a human cell-free DNA replication initiation system. Therefore, Dm1 is the first functional sbRNA described in insects, allowing future investigations into the physiological roles of sbRNAs in the genetically tractable model organism D. melanogaster.

DIS3L2, in which mutations have been linked to Perlman syndrome, is an RNA-binding protein with 3'-5' exoribonuclease activity. It contains two CSD domains and one S1 domain, all of which are RNA-binding domains, and one RNB domain that is responsible for the exoribonuclease activity. The 3' polyuridine of RNA substrates can serve as a degradation signal for DIS3L2. Because DIS3L2 is predominantly localized in the cytoplasm, it can recognize, bind, and mediate the degradation of cytoplasmic uridylated RNA, including pre-microRNA, mature microRNA, mRNA, and some other non-coding RNAs. Therefore, DIS3L2 plays an important role in cytoplasmic RNA surveillance and decay. DIS3L2 is involved in multiple biological and physiological processes such as cell division, proliferation, differentiation, and apoptosis. Nonetheless, the function of DIS3L2, especially its association with cancer, remains largely unknown. We summarize here the RNA substrates degraded by DIS3L2 with its exonucleolytic activity, together with the corresponding biological functions it is implicated in. Furthermore, we discuss whether DIS3L2 can function independently of its 3'-5' exoribonuclease activity, as well as its potential tumor-suppressive or oncogenic roles during cancer progression.

Cancer is characterized by multiple genetic and epigenetic alterations, including a higher prevalence of mutations of oncogenes and/or tumor suppressors. Mounting evidences have shown that noncoding RNAs (ncRNAs) are involved in the epigenetic regulation of cancer genes and their associated pathways. The clustered regularly interspaced short palindromic repeats (CRISPR)-associated nuclease 9 (CRISPR/Cas9) system, a revolutionary genome-editing technology, has shed light on ncRNA-based cancer therapy. Here, we briefly introduce the classifications and mechanisms of CRISPR/Cas9 system. Importantly, we mainly focused on the applications of CRISPR/Cas9 system as a molecular tool for ncRNA (microRNA, long noncoding RNA and circular RNA, etc.) editing in human cancers, and the novel techniques that are based on CRISPR/Cas9 system. Additionally, the off-target effects and the corresponding solutions as well as the challenges toward CRISPR/Cas9 were also evaluated and discussed. Long- and short-ncRNAs have been employed as targets in precision oncology, and CRISPR/Cas9-mediated ncRNA editing may provide an excellent way to cure cancer.

Epstein-Barr virus (EBV) expresses an abundant nuclear noncoding RNA called EBER2, which interacts with and acts as a guide RNA for the host transcription factor PAX5. This ribonucleoprotein complex localizes to the terminal repeat (TR) regions of the EBV genome via RNA-RNA interactions between EBER2 and nascent transcripts originating from these target sites. Given the fact that EBER2 base pairs with a viral RNA, we developed a protocol to identify EBER2-interacting RNAs in a transcriptome-wide manner. Our approach entails psoralen-mediated crosslinking, selection with antisense oligonucleotides targeting EBER2, and RNase V1 digestion coupled to next-generation sequencing. The use of RNase V1 circumvents the need of extensive computational analysis post data acquisition to search for predicted RNA hybrids, as the RNase V1 cleavage site marks the region of RNA duplex formation. As proof of principle, we show that our approach correctly identifies the known EBER2 interaction with TR RNAs. Moreover, we identify the host functional noncoding RNAs MRP, H1, and 7SL RNAs as well as three putative enhancer RNAs as candidate EBER2-interacting RNAs. As all of these gene loci exhibit PAX5 occupancy, we propose that EBER2 is recruited to these sites through its binding partner PAX5 and forms RNA-RNA interactions with nascent transcripts on chromatin. Thus, our novel approach facilitates the identification of targeted RNA-RNA-interactions and minimizes the need of downstream computational analyses to predict RNA duplexes.

The Arabidopsis GUT15 RNA belongs to a class of noncoding RNAs that are expressed from the intergenic regions of protein-coding genes. We show that the RNA polymerase II transcribed GUT15 transcript serves as a precursor for two stable RNA species, a tRNA-like molecule and GUT15-tRF-F5, which are both encoded by the final intron in the GUT15 gene. The GUT15-encoded tRNA-like molecule cannot be autonomously transcribed by RNA polymerase III. However, this molecule contains a CCA motif, suggesting that it may enter the tRNA maturation pathway. The GUT15-encoded tRNA-like sequence has an inhibiting effect on the splicing of its host intron. Moreover, we demonstrate that the canonical tRNA genes nested within introns do not affect the splicing patterns of their host protein-coding transcripts.

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

Viruses masterfully regulate host gene expression during infection. Many do so, in part, by expressing non-coding RNAs. Recent work has shown that HSUR 2, a viral non-coding RNA expressed by the oncogenic Herpesvirus saimiri, regulates mRNA expression through a novel mechanism. HSUR 2 base pairs with both target mRNAs and host miRNAs in infected cells. This results in HSUR 2-dependent recruitment of host miRNAs and associated Ago proteins to target mRNAs, and the subsequent destabilization of target mRNAs. Using this mechanism, this virus regulates key cellular pathways during viral infection. Here I discuss the evolution of our thinking about HSUR function and explore the implications of recent findings in relation to the current views on the functions of interactions between miRNAs and other classes of non-coding RNAs, the potential advantages of this mechanism of regulation of gene expression, and the evolutionary origin of HSUR 2.

Bordetella pertussis is the causative agent of whooping cough, a respiratory disease still considered as a major public health threat and for which recent re-emergence has been observed. Constant reshuffling of Bordetella pertussis genome organization was observed during evolution. These rearrangements are essentially mediated by Insertion Sequences (IS), a mobile genetic elements present in more than 230 copies in the genome, which are supposed to be one of the driving forces enabling the pathogen to escape from vaccine-induced immunity. Here we use high-throughput sequencing approaches (RNA-seq and differential RNA-seq), to decipher Bordetella pertussis transcriptome characteristics and to evaluate the impact of IS elements on transcriptome architecture. Transcriptional organization was determined by identification of transcription start sites and revealed also a large variety of non-coding RNAs including sRNAs, leaderless mRNAs or long 3' and 5'UTR including seven riboswitches. Unusual topological organizations, such as overlapping 5'- or 3'-extremities between oppositely orientated mRNA were also unveiled. The pivotal role of IS elements in the transcriptome architecture and their effect on the transcription of neighboring genes was examined. This effect is mediated by the introduction of IS harbored promoters or by emergence of hybrid promoters. This study revealed that in addition to their impact on genome rearrangements, most of the IS also impact on the expression of their flanking genes. Furthermore, the transcripts produced by IS are strain-specific due to the strain to strain variation in IS copy number and genomic context.

Post-testicular sperm maturation and storage within the epididymis is a key determinant of gamete quality and fertilization competence. Here we demonstrate that mouse spermatozoa possess a complex small non-protein-coding RNA (sRNA) profile, the composition of which is markedly influenced by their epididymal transit. Thus, although microRNAs (miRNAs) are highly represented in the spermatozoa of the proximal epididymis, this sRNA class is largely diminished in mature spermatozoa of the distal epididymis. Coincident with this, a substantial enrichment in Piwi-interacting RNA (piRNA) abundance in cauda spermatozoa was detected. Further, features of cauda piRNAs, including; predominantly 29-31 nts in length; preference for uracil at their 5' terminus; no adenine enrichment at piRNA nt 10, and; predominantly mapping to intergenic regions of the mouse genome, indicate that these piRNAs are generated by the PIWIL1-directed primary piRNA production pathway. Accordingly, PIWIL1 was detected via immunoblotting and mass spectrometry in epididymal spermatozoa. These data provide insight into the complexity and dynamic nature of the sRNA profile of spermatozoa and raise the intriguing prospect that piRNAs are generated in situ in maturing spermatozoa. Such information is of particular interest in view of the potential role for paternal sRNAs in influencing conception, embryo development and intergenerational inheritance.

Adenosine-to-inosine RNA editing is a conserved process, which is performed by ADAR enzymes. By changing nucleotides in coding regions of genes and altering codons, ADARs expand the cell's protein repertoire. This function of the ADAR enzymes is essential for human brain development. However, most of the known editing sites are in non-coding repetitive regions in the transcriptome and the purpose of editing in these regions is unclear. Recent studies, which have shown that editing levels of transcripts vary between tissues and developmental stages in many organisms, suggest that the targeted RNA and ADAR editing are both regulated. We discuss the implications of these findings, and the possible role of RNA editing in innate immunity.

Cis-natural antisense transcripts (cis-NATs) are long noncoding RNAs transcribed from the opposite strand and overlapping coding and noncoding genes on the sense strand. cis-NATs are widely present in the human genome and can be involved in multiple mechanisms of gene regulation. Here, we describe the presence of cis-NATs in the 3' distal region of the c-MYC locus and investigate their impact on transcriptional regulation of this key oncogene in human cancers. We found that cis-NATs are produced as consequence of the activation of cryptic transcription initiation sites in the 3' distal region downstream of the c-MYC 3'UTR. The process is tightly regulated and leads to the formation of two main transcripts, NAT6531 and NAT6558, which differ in their ability to fold into stem-loop secondary structures. NAT6531 acts as a substrate for DICER and as a source of small RNAs capable of modulating c-MYC transcription. This complex system, based on the interplay between cis-NATs and NAT-derived small RNAs, may represent an important layer of epigenetic regulation of the expression of c-MYC and other genes in human cells.

Enhancers are cis- regulatory genetic elements crucial for controlling temporal and cell-type specific patterns of gene expression. Active enhancers generate bi-directional non-coding RNA transcripts called enhancer RNAs (eRNAs). eRNAs are important for stimulating gene expression, but precise mechanisms for this ability remain unclear. Here we highlight recent findings that demonstrate a direct interaction between RNAs and the transcriptional co-activator Creb-binding protein (CBP). Notably, RNA binding could stimulate the core histone acetyltransferase activity of the enzyme, observable in cells as a link between eRNA production, CBP-dependent histone acetylation and expression of genes regulated by specific enhancers. Although RNA binding was independent of RNA sequence, specificity arises in a locus-specific manner at transcribed sites where CBP was bound to chromatin. The results suggest a functional role for eRNAs as regulatory molecules that are able to stimulate the activity of a key epigenetic regulatory enzyme, thereby promoting gene expression. Furthermore, they suggest an intriguing role for eRNAs: by modulating the activity of chromatin modifying enzymes, they could directly impact transcription by altering the chromatin environment.

The function of 6S RNA, a global regulator of transcription, was studied in the photosynthetic α-proteobacterium Rhodobacter sphaeroides. The cellular levels of R. sphaeroides 6S RNA peak toward the transition to stationary phase and strongly decrease during extended stationary phase. The synthesis of so-called product RNA transcripts (mainly 12-16-mers) on 6S RNA as template by RNA polymerase was found to be highest in late exponential phase. Product RNA ≥ 13-mers are expected to trigger the dissociation of 6S RNA:RNA polymerase complexes. A 6S RNA deletion in R. sphaeroides had no impact on growth under various metabolic and oxidative stress conditions (with the possible exception of tert-butyl hydroperoxide stress). However, the 6S RNA knockout resulted in a robust growth defect under high salt stress (0.25 M NaCl). Remarkably, the sspA gene encoding the putative salt stress-induced membrane protein SspA and located immediately downstream of the 6S RNA (ssrS) gene on the antisense strand was expressed at elevated levels in the ΔssrS strain when grown in the presence of 250 mM NaCl.

For a long time, eukaryotic stand-alone pseudouridine synthases (Pus enzymes) were neglected as non-essential enzymes adding seemingly simple modifications to tRNAs and small nuclear RNAs. Most studies were limited to the identification and initial characterization of the yeast Pus enzymes. However, recent transcriptome-wide mapping of pseudouridines in yeast and humans revealed pervasive modification of mRNAs and other non-coding RNAs by Pus enzymes which is dynamically regulated in response to cellular stress. Moreover, mutations in at least 2 genes encoding human Pus enzymes cause inherited diseases affecting muscle and brain function. Together, the recent findings suggest a broader-than-anticipated role of the Pus enzymes which are emerging as potential regulators of gene expression. In this review, we summarize the current knowledge on Pus enzymes, generate hypotheses regarding their cellular function and outline future areas of research of pseudouridine synthases.

Methodological advances that allow deeper characterization of non-coding elements in the genome have started to reveal the full spectrum of deregulation in cancer. We generated an inducible cell model to track transcriptional changes after induction of a well-known leukemia-inducing fusion gene, ETV6-RUNX1. Our data revealed widespread transcriptional alterations outside coding elements in the genome. This adds to the growing list of various alterations in the non-coding genome in cancer and pinpoints their role in diseased cellular state.

Chromosomal fusions are common in normal and cancer cells and can produce aberrant gene products that promote transformation. The mechanisms driving these fusions are poorly understood, but recurrent fusions are widespread. This suggests an underlying mechanism, and some authors have proposed a possible role for RNA in this process. The unicellular eukaryote Oxytricha trifallax displays an exorbitant capacity for natural genome editing, when it rewrites its germline genome to form a somatic epigenome. This developmental process provides a powerful model system to directly test the influence of small noncoding RNAs on chromosome fusion events during somatic differentiation. Here we show that small RNAs are capable of inducing chromosome fusions in 4 distinct cases (out of 4 tested), including one fusion of 3 chromosomes. We further show that these RNA-mediated chromosome fusions are heritable over multiple sexual generations and that transmission of the acquired fusion is associated with endogenous production of novel piRNA molecules that target the fused junction. We also demonstrate the capacity of a long noncoding RNA (lncRNA) to induce chromosome fusion of 2 distal germline loci. These results underscore the ability of short-lived, aberrant RNAs to act as drivers of chromosome fusion events that can be stably transmitted to future generations.

Eukaryotic cells produce a variety of non-coding RNAs (ncRNAs), many of which have been shown to play pivotal roles in biological processes such as differentiation, maintenance of pluripotency of stem cells, and cellular response to various stresses. Genome-wide analyses have revealed that many ncRNAs are transcribed around regulatory DNA elements located proximal or distal to gene promoters, but their biological functions are largely unknown. Recently, it has been demonstrated in yeast and mouse that ncRNA transcription around gene promoters and enhancers facilitates DNA binding of transcription factors to their target sites. These results suggest universal roles of promoter/enhancer-associated ncRNAs in the recruitment of transcription factors to their binding sites.

Cells secrete extracellular RNA (exRNA) to their surrounding environment and exRNA has been found in many body fluids such as blood, breast milk and cerebrospinal fluid. However, there are conflicting results regarding the nature of exRNA. Here, we have separated 2 distinct exRNA profiles released by mast cells, here termed high-density (HD) and low-density (LD) exRNA. The exRNA in both fractions was characterized by microarray and next-generation sequencing. Both exRNA fractions contained mRNA and miRNA, and the mRNAs in the LD exRNA correlated closely with the cellular mRNA, whereas the HD mRNA did not. Furthermore, the HD exRNA was enriched in lincRNA, antisense RNA, vault RNA, snoRNA, and snRNA with little or no evidence of full-length 18S and 28S rRNA. The LD exRNA was enriched in mitochondrial rRNA, mitochondrial tRNA, tRNA, piRNA, Y RNA, and full-length 18S and 28S rRNA. The proteomes of the HD and LD exRNA-containing fractions were determined with LC-MS/MS and analyzed with Gene Ontology term finder, which showed that both proteomes were associated with the term extracellular vesicles and electron microscopy suggests that at least a part of the exRNA is associated with exosome-like extracellular vesicles. Additionally, the proteins in the HD fractions tended to be associated with the nucleus and ribosomes, whereas the LD fraction proteome tended to be associated with the mitochondrion. We show that the 2 exRNA signatures released by a single cell type can be separated by floatation on a density gradient. These results show that cells can release multiple types of exRNA with substantial differences in RNA species content. This is important for any future studies determining the nature and function of exRNA released from different cells under different conditions.