Long non-coding RNAs: definitions, functions, challenges and recommendations.

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

RNA modifications detection by comparative Nanopore direct RNA sequencing.

RNA molecules undergo a vast array of chemical post-transcriptional modifications (PTMs) that can affect their structure and interaction properties. In recent years, a growing number of PTMs have been successfully mapped to the transcriptome using experimental approaches relying on high-throughput sequencing. Oxford Nanopore direct-RNA sequencing has been shown to be sensitive to RNA modifications. We developed and validated Nanocompore, a robust analytical framework that identifies modifications from these data. Our strategy compares an RNA sample of interest against a non-modified control sample, not requiring a training set and allowing the use of replicates. We show that Nanocompore can detect different RNA modifications with position accuracy in vitro, and we apply it to profile m6A in vivo in yeast and human RNAs, as well as in targeted non-coding RNAs. We confirm our results with orthogonal methods and provide novel insights on the co-occurrence of multiple modified residues on individual RNA molecules.

Long non-coding RNAs associated with infection and vaccine-induced immunity.

The immune system responds to infection or vaccination through a dynamic and complex process that involves several molecular and cellular factors. Among these factors, long non-coding RNAs (lncRNAs) have emerged as significant players in all areas of biology, particularly in immunology. Most of the mammalian genome is transcribed in a highly regulated manner, generating a diversity of lncRNAs that impact the differentiation and activation of immune cells and affect innate and adaptive immunity. Here, we have reviewed the range of functions and mechanisms of lncRNAs in response to infectious disease, including pathogen recognition, interferon (IFN) response, and inflammation. We describe examples of lncRNAs exploited by pathogenic agents during infection, which indicate that lncRNAs are a fundamental part of the arms race between hosts and pathogens. We also discuss lncRNAs potentially implicated in vaccine-induced immunity and present examples of lncRNAs associated with the antibody response of subjects receiving Influenza or Yellow Fever vaccines. Elucidating the widespread involvement of lncRNAs in the immune system will improve our understanding of the factors affecting immune response to different pathogenic agents, to better prevent and treat disease.

RNA, Genome Output and Input.

RNA, the transcriptional output of genomes, not only templates protein synthesis or directly engages in catalytic functions, but can feed back to the genome and serve as regulatory input for gene expression. Transcripts affecting the RNA abundance of other genes act by mechanisms similar to and in concert with protein factors that control transcription. Through recruitment or blocking of activating and silencing complexes to specific genomic loci, RNA and protein factors can favor transcription or lower the local gene expression potential. Most regulatory proteins enter nuclei from all directions to start the search for increased affinity to specific DNA sequences or to other proteins nearby genuine gene targets. In contrast, RNAs emerge from spatial point sources within nuclei, their encoding genes. A transcriptional burst can result in the local appearance of multiple nascent RNA copies at once, in turn increasing local nucleic acid density and RNA motif abundance before diffusion into the nuclear neighborhood. The confined initial localization of regulatory RNAs causing accumulation of protein co-factors raises the intriguing possibility that target specificity of non-coding, and probably coding, RNAs is achieved through gene/RNA positioning and spatial proximity to regulated genomic regions. Here we review examples of positional cis conservation of regulatory RNAs with respect to target genes, spatial proximity of enhancer RNAs to promoters through DNA looping and RNA-mediated formation of membrane-less structures to control chromatin structure and expression. We speculate that linear and spatial proximity between regulatory RNA-encoding genes and gene targets could possibly ease the evolutionary pressure on maintaining regulatory RNA sequence conservation.

Long noncoding RNAs are involved in multiple immunological pathways in response to vaccination.

Understanding the mechanisms of vaccine-elicited protection contributes to the development of new vaccines. The emerging field of systems vaccinology provides detailed information on host responses to vaccination and has been successfully applied to study the molecular mechanisms of several vaccines. Long noncoding RNAs (lncRNAs) are crucially involved in multiple biological processes, but their role in vaccine-induced immunity has not been explored. We performed an analysis of over 2,000 blood transcriptome samples from 17 vaccine cohorts to identify lncRNAs potentially involved with antibody responses to influenza and yellow fever vaccines. We have created an online database where all results from this analysis can be accessed easily. We found that lncRNAs participate in distinct immunological pathways related to vaccine-elicited responses. Among them, we showed that the expression of lncRNA FAM30A was high in B cells and correlates with the expression of immunoglobulin genes located in its genomic vicinity. We also identified altered expression of these lncRNAs in RNA-sequencing (RNA-seq) data from a cohort of children following immunization with intranasal live attenuated influenza vaccine, suggesting a common role across several diverse vaccines. Taken together, these findings provide evidence that lncRNAs have a significant impact on immune responses induced by vaccination.

Interaction of Sox2 with RNA binding proteins in mouse embryonic stem cells.

Sox2 is a master transcriptional regulator of embryonic development. In this study, we determined the protein interactome of Sox2 in the chromatin and nucleoplasm of mouse embryonic stem (mES) cells. Apart from canonical interactions with pluripotency-regulating transcription factors, we identified interactions with several chromatin modulators, including members of the heterochromatin protein 1 (HP1) family, suggesting a role for Sox2 in chromatin-mediated transcriptional repression. Sox2 was also found to interact with RNA binding proteins (RBPs), including proteins involved in RNA processing. RNA immunoprecipitation followed by sequencing revealed that Sox2 associates with different messenger RNAs, as well as small nucleolar RNA Snord34 and the non-coding RNA 7SK. 7SK has been shown to regulate transcription at gene regulatory regions, which could suggest a functional interaction with Sox2 for chromatin recruitment. Nevertheless, we found no evidence of Sox2 modulating recruitment of 7SK to chromatin when examining 7SK chromatin occupancy by Chromatin Isolation by RNA Purification (ChIRP) in Sox2 depleted mES cells. In addition, knockdown of 7SK in mES cells did not lead to any change in Sox2 occupancy at 7SK-regulated genes. Thus, our results show that Sox2 extensively interacts with RBPs, and suggest that Sox2 and 7SK co-exist in a ribonucleoprotein complex whose function is not to regulate chromatin recruitment, but could rather regulate other processes in the nucleoplasm.

Genomic positional conservation identifies topological anchor point RNAs linked to developmental loci.

BACKGROUND: The mammalian genome is transcribed into large numbers of long noncoding RNAs (lncRNAs), but the definition of functional lncRNA groups has proven difficult, partly due to their low sequence conservation and lack of identified shared properties. Here we consider promoter conservation and positional conservation as indicators of functional commonality. RESULTS: We identify 665 conserved lncRNA promoters in mouse and human that are preserved in genomic position relative to orthologous coding genes. These positionally conserved lncRNA genes are primarily associated with developmental transcription factor loci with which they are coexpressed in a tissue-specific manner. Over half of positionally conserved RNAs in this set are linked to chromatin organization structures, overlapping binding sites for the CTCF chromatin organiser and located at chromatin loop anchor points and borders of topologically associating domains (TADs). We define these RNAs as topological anchor point RNAs (tapRNAs). Characterization of these noncoding RNAs and their associated coding genes shows that they are functionally connected: they regulate each other's expression and influence the metastatic phenotype of cancer cells in vitro in a similar fashion. Furthermore, we find that tapRNAs contain conserved sequence domains that are enriched in motifs for zinc finger domain-containing RNA-binding proteins and transcription factors, whose binding sites are found mutated in cancers. CONCLUSIONS: This work leverages positional conservation to identify lncRNAs with potential importance in genome organization, development and disease. The evidence that many developmental transcription factors are physically and functionally connected to lncRNAs represents an exciting stepping-stone to further our understanding of genome regulation.

A lncRNA fine tunes the dynamics of a cell state transition involving Lin28, let-7 and de novo DNA methylation.

Execution of pluripotency requires progression from the naïve status represented by mouse embryonic stem cells (ESCs) to a state capacitated for lineage specification. This transition is coordinated at multiple levels. Non-coding RNAs may contribute to this regulatory orchestra. We identified a rodent-specific long non-coding RNA (lncRNA) linc1281, hereafter Ephemeron (Eprn), that modulates the dynamics of exit from naïve pluripotency. Eprn deletion delays the extinction of ESC identity, an effect associated with perduring Nanog expression. In the absence of Eprn, Lin28a expression is reduced which results in persistence of let-7 microRNAs, and the up-regulation of de novo methyltransferases Dnmt3a/b is delayed. Dnmt3a/b deletion retards ES cell transition, correlating with delayed Nanog promoter methylation and phenocopying loss of Eprn or Lin28a. The connection from lncRNA to miRNA and DNA methylation facilitates the acute extinction of naïve pluripotency, a pre-requisite for rapid progression from preimplantation epiblast to gastrulation in rodents. Eprn illustrates how lncRNAs may introduce species-specific network modulations.

The Evx1/Evx1as gene locus regulates anterior-posterior patterning during gastrulation.

Thousands of sense-antisense mRNA-lncRNA gene pairs occur in the mammalian genome. While there is usually little doubt about the function of the coding transcript, the function of the lncRNA partner is mostly untested. Here we examine the function of the homeotic Evx1-Evx1as gene locus. Expression is tightly co-regulated in posterior mesoderm of mouse embryos and in embryoid bodies. Expression of both genes is enhanced by BMP4 and WNT3A, and reduced by Activin. We generated a suite of deletions in the locus by CRISPR-Cas9 editing. We show EVX1 is a critical downstream effector of BMP4 and WNT3A with respect to patterning of posterior mesoderm. The lncRNA, Evx1as arises from alternative promoters and is difficult to fully abrogate by gene editing or siRNA approaches. Nevertheless, we were able to generate a large 2.6 kb deletion encompassing the shared promoter with Evx1 and multiple additional exons of Evx1as. This led to an identical dorsal-ventral patterning defect to that generated by micro-deletion in the DNA-binding domain of EVX1. Thus, Evx1as has no function independent of EVX1, and is therefore unlikely to act in trans. We predict many antisense lncRNAs have no specific trans function, possibly only regulating the linked coding genes in cis.

The non-coding snRNA 7SK controls transcriptional termination, poising, and bidirectionality in embryonic stem cells.

BACKGROUND: Pluripotency is characterized by a unique transcriptional state, in which lineage-specification genes are poised for transcription upon exposure to appropriate stimuli, via a bivalency mechanism involving the simultaneous presence of activating and repressive methylation marks at promoter-associated histones. Recent evidence suggests that other mechanisms, such as RNA polymerase II pausing, might be operational in this process, but their regulation remains poorly understood. RESULTS: Here we identify the non-coding snRNA 7SK as a multifaceted regulator of transcription in embryonic stem cells. We find that 7SK represses a specific cohort of transcriptionally poised genes with bivalent or activating chromatin marks in these cells, suggesting a novel poising mechanism independent of Polycomb activity. Genome-wide analysis shows that 7SK also prevents transcription downstream of polyadenylation sites at several active genes, indicating that 7SK is required for normal transcriptional termination or control of 3'-UTR length. In addition, 7SK suppresses divergent upstream antisense transcription at more than 2,600 loci, including many that encode divergent long non-coding RNAs, a finding that implicates the 7SK snRNA in the control of transcriptional bidirectionality. CONCLUSIONS: Our study indicates that a single non-coding RNA, the snRNA 7SK, is a gatekeeper of transcriptional termination and bidirectional transcription in embryonic stem cells and mediates transcriptional poising through a mechanism independent of chromatin bivalency.

Non-coding RNAs in homeostasis, disease and stress responses: an evolutionary perspective.

Cells and organisms are subject to challenges and perturbations in their environment and physiology in all stages of life. The molecular response to such changes, including insulting conditions such as pathogen infections, involves coordinated modulation of gene expression programmes and has not only homeostatic but also ecological and evolutionary importance. Although attention has been primarily focused on signalling pathways and protein networks, non-coding RNAs (ncRNAs), which comprise a significant output of the genomes of prokaryotes and especially eukaryotes, are increasingly implicated in the molecular mechanisms of these responses. Long and short ncRNAs not only regulate development and cell physiology, they are also involved in disease states, including cancers, in host-pathogen interactions, and in a variety of stress responses. Indeed, regulatory RNAs are part of genetically encoded response networks and also underpin epigenetic processes, which are emerging as key mechanisms of adaptation and transgenerational inheritance. Here we present the growing evidence that ncRNAs are intrinsically involved in cellular and organismal adaptation processes, in both robustness and protection to stresses, as well as in mechanisms generating evolutionary change.

Pinstripe: a suite of programs for integrating transcriptomic and proteomic datasets identifies novel proteins and improves differentiation of protein-coding and non-coding genes.

MOTIVATION: Comparing transcriptomic data with proteomic data to identify protein-coding sequences is a long-standing challenge in molecular biology, one that is exacerbated by the increasing size of high-throughput datasets. To address this challenge, and thereby to improve the quality of genome annotation and understanding of genome biology, we have developed an integrated suite of programs, called Pinstripe. We demonstrate its application, utility and discovery power using transcriptomic and proteomic data from publicly available datasets. RESULTS: To demonstrate the efficacy of Pinstripe for large-scale analysis, we applied Pinstripe's reverse peptide mapping pipeline to a transcript library including de novo assembled transcriptomes from the human Illumina Body Atlas (IBA2) and GENCODE v10 gene annotations, and the EBI Proteomics Identifications Database (PRIDE) peptide database. This analysis identified 736 canonical open reading frames (ORFs) supported by three or more PRIDE peptide fragments that are positioned outside any known coding DNA sequence (CDS). Because of the unfiltered nature of the PRIDE database and high probability of false discovery, we further refined this list using independent evidence for translation, including the presence of a Kozak sequence or functional domains, synonymous/non-synonymous substitution ratios and ORF length. Using this integrative approach, we observed evidence of translation from a previously unknown let7e primary transcript, the archetypical lncRNA H19, and a homolog of RD3. Reciprocally, by exclusion of transcripts with mapped peptides or significant ORFs (>80 codon), we identify 32 187 loci with RNAs longer than 2000 nt that are unlikely to encode proteins. AVAILABILITY AND IMPLEMENTATION: Pinstripe (pinstripe.matticklab.com) is freely available as source code or a Mono binary. Pinstripe is written in C# and runs under the Mono framework on Linux or Mac OS X, and both under Mono and .Net under Windows. CONTACT: m.dinger@garvan.org.au or j.mattick@garvan.org.au SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

lncRNAdb: a reference database for long noncoding RNAs.

Large numbers of long RNAs with little or no protein-coding potential [long noncoding RNAs (lncRNAs)] are being identified in eukaryotes. In parallel, increasing data describing the expression profiles, molecular features and functions of individual lncRNAs in a variety of systems are accumulating. To enable the systematic compilation and updating of this information, we have developed a database (lncRNAdb) containing a comprehensive list of lncRNAs that have been shown to have, or to be associated with, biological functions in eukaryotes, as well as messenger RNAs that have regulatory roles. Each entry contains referenced information about the RNA, including sequences, structural information, genomic context, expression, subcellular localization, conservation, functional evidence and other relevant information. lncRNAdb can be searched by querying published RNA names and aliases, sequences, species and associated protein-coding genes, as well as terms contained in the annotations, such as the tissues in which the transcripts are expressed and associated diseases. In addition, lncRNAdb is linked to the UCSC Genome Browser for visualization and Noncoding RNA Expression Database (NRED) for expression information from a variety of sources. lncRNAdb provides a platform for the ongoing collation of the literature pertaining to lncRNAs and their association with other genomic elements. lncRNAdb can be accessed at: http://www.lncrnadb.org/.

Pervasive transcription of the eukaryotic genome: functional indices and conceptual implications.

Genome-wide analyses of the eukaryotic transcriptome have revealed that the majority of the genome is transcribed, producing large numbers of non-protein-coding RNAs (ncRNAs). This surprising observation challenges many assumptions about the genetic programming of higher organisms and how information is stored and organized within the genome. Moreover, the rapid advances in genomics have given little opportunity for biologists to integrate these emerging findings into their intellectual and experimental frameworks. This problem has been compounded by the perception that genome-wide studies often generate more questions than answers, which in turn has led to confusion and controversy. In this article, we address common questions associated with the phenomenon of pervasive transcription and consider the indices that can be used to evaluate the function (or lack thereof) of the resulting ncRNAs. We suggest that many lines of evidence, including expression profiles, conservation signatures, chromatin modification patterns and examination of increasing numbers of individual cases, argue in favour of the widespread functionality of non-coding transcription. We also discuss how informatic and experimental approaches used to analyse protein-coding genes may not be applicable to ncRNAs and how the general perception that protein-coding genes form the main informational output of the genome has resulted in much of the misunderstanding surrounding pervasive transcription and its potential significance. Finally, we present the conceptual implications of the majority of the eukaryotic genome being functional and describe how appreciating this perspective will provide considerable opportunity to further understand the molecular basis of development and complex diseases.

Complex architecture and regulated expression of the Sox2ot locus during vertebrate development.

The Sox2 gene is a key regulator of pluripotency embedded within an intron of a long noncoding RNA (ncRNA), termed Sox2 overlapping transcript (Sox2ot), which is transcribed in the same orientation. However, this ncRNA remains uncharacterized. Here we show that Sox2ot has multiple transcription start sites associated with genomic features that indicate regulated expression, including highly conserved elements (HCEs) and chromatin marks characteristic of gene promoters. To identify biological processes in which Sox2ot may be involved, we analyzed its expression in several developmental systems, compared to expression of Sox2. We show that Sox2ot is a stable transcript expressed in mouse embryonic stem cells, which, like Sox2, is down-regulated upon induction of embryoid body (EB) differentiation. However, in contrast to Sox2, Sox2ot is up-regulated during EB mesoderm-lineage differentiation. In adult mouse, Sox2ot isoforms were detected in tissues where Sox2 is expressed, as well as in different tissues, supporting independent regulation of expression of the ncRNA. Sox2dot, an isoform of Sox2ot transcribed from a distal HCE located >500 kb upstream of Sox2, was detected exclusively in the mouse brain, with enrichment in regions of adult neurogenesis. In addition, Sox2ot isoforms are transcribed from HCEs upstream of Sox2 in other vertebrates, including in several regions of the human brain. We also show that Sox2ot is dynamically regulated during chicken and zebrafish embryogenesis, consistently associated with central nervous system structures. These observations provide insight into the structure and regulation of the Sox2ot gene, and suggest conserved roles for Sox2ot orthologs during vertebrate development.

RNA regulation of epigenetic processes.

There is increasing evidence that dynamic changes to chromatin, chromosomes and nuclear architecture are regulated by RNA signalling. Although the precise molecular mechanisms are not well understood, they appear to involve the differential recruitment of a hierarchy of generic chromatin modifying complexes and DNA methyltransferases to specific loci by RNAs during differentiation and development. A significant fraction of the genome-wide transcription of non-protein coding RNAs may be involved in this process, comprising a previously hidden layer of intermediary genetic information that underpins developmental ontogeny and the differences between species, ecotypes and individuals. It is also evident that RNA editing is a primary means by which hardwired genetic information in animals can be altered by environmental signals, especially in the brain, indicating a dynamic RNA-mediated interplay between the transcriptome, the environment and the epigenome. Moreover, RNA-directed regulatory processes may also transfer epigenetic information not only within cells but also between cells and organ systems, as well as across generations.

MEN epsilon/beta nuclear-retained non-coding RNAs are up-regulated upon muscle differentiation and are essential components of paraspeckles.

Studies of the transcriptional output of the human and mouse genomes have revealed that there are many more transcripts produced than can be accounted for by predicted protein-coding genes. Using a custom microarray, we have identified 184 non-coding RNAs that exhibit more than twofold up- or down-regulation upon differentiation of C2C12 myoblasts into myotubes. Here, we focus on the Men epsilon/beta locus, which is up-regulated 3.3-fold during differentiation. Two non-coding RNA isoforms are produced from a single RNA polymerase II promoter, differing in the location of their 3' ends. Men epsilon is a 3.2-kb polyadenylated RNA, whereas Men beta is an approximately 20-kb transcript containing a genomically encoded poly(A)-rich tract at its 3'-end. The 3'-end of Men beta is generated by RNase P cleavage. The Men epsilon/beta transcripts are localized to nuclear paraspeckles and directly interact with NONO. Knockdown of MEN epsilon/beta expression results in the disruption of nuclear paraspeckles. Furthermore, the formation of paraspeckles, after release from transcriptional inhibition by DRB treatment, was suppressed in MEN epsilon/beta-depleted cells. Our findings indicate that the MEN epsilon/beta non-coding RNAs are essential structural/organizational components of paraspeckles.

Noncoding RNA in development.

Non-protein-coding sequences increasingly dominate the genomes of multicellular organisms as their complexity increases, in contrast to protein-coding genes, which remain relatively static. Most of the mammalian genome and indeed that of all eukaryotes is expressed in a cell- and tissue-specific manner, and there is mounting evidence that much of this transcription is involved in the regulation of differentiation and development. Different classes of small and large noncoding RNAs (ncRNAs) have been shown to regulate almost every level of gene expression, including the activation and repression of homeotic genes and the targeting of chromatin-remodeling complexes. ncRNAs are involved in developmental processes in both simple and complex eukaryotes, and we illustrate this in the latter by focusing on the animal germline, brain, and eye. While most have yet to be systematically studied, the emerging evidence suggests that there is a vast hidden layer of regulatory ncRNAs that constitutes the majority of the genomic programming of multicellular organisms and plays a major role in controlling the epigenetic trajectories that underlie their ontogeny.