Rapid turnover of long noncoding RNAs and the evolution of gene expression

Pseudogene-derived lncRNAs: emerging regulators of gene expression

In the more than one decade since the completion of the Human Genome Project, the prevalence of non-protein-coding functional elements in the human genome has emerged as a key revelation in post-genomic biology. Highlighted by the ENCODE (Encyclopedia of DNA Elements) and FANTOM (Functional Annotation of Mammals) consortia, these elements include tens of thousands of pseudogenes, as well as comparably numerous long non-coding RNA (lncRNA) genes. Pseudogene transcription and function remain insufficiently understood. However, the field is of great importance for human disease due to the high sequence similarity between pseudogenes and their parental protein-coding genes, which generates the potential for sequence-specific regulation. Recent case studies have established essential and coordinated roles of both pseudogenes and lncRNAs in development and disease in metazoan systems, including functional impacts of lncRNA transcription at pseudogene loci on the regulation of the pseudogenes’ parental genes. This review synthesizes the nascent evidence for regulatory modalities jointly exerted by lncRNAs and pseudogenes in human disease, and for recent evolutionary origins of these systems.

Decoupling of evolutionary changes in transcription factor binding and gene expression in mammals

To understand the evolutionary dynamics between transcription factor (TF) binding and gene expression in mammals, we compared transcriptional output and the binding intensities for three tissue-specific TFs in livers from four closely related mouse species. For each transcription factor, TF-dependent genes and the TF binding sites most likely to influence mRNA expression were identified by comparing mRNA expression levels between wild-type and TF knockout mice. Independent evolution was observed genome-wide between the rate of change in TF binding and the rate of change in mRNA expression across taxa, with the exception of a small number of TF-dependent genes. We also found that binding intensities are preferentially conserved near genes whose expression is dependent on the TF, and the conservation is shared among binding peaks in close proximity to each other near the TSS. Expression of TF-dependent genes typically showed an increased sensitivity to changes in binding levels as measured by mRNA abundance. Taken together, these results highlight a significant tolerance to evolutionary changes in TF binding intensity in mammalian transcriptional networks and suggest that some TF-dependent genes may be largely regulated by a single TF across evolution.

High expression of long intervening non-coding RNA OLMALINC in the human cortical white matter is associated with regulation of oligodendrocyte maturation

Background

Long intervening non-coding RNAs (lincRNAs) are a recently discovered subclass of non-coding RNAs. LincRNAs are expressed across the mammalian genome and contribute to the pervasive transcription phenomenon. They display a tissue-specific and species-specific mode of expression and are present abundantly in the brain.

Results

Here, we report the expression patterns of oligodendrocyte maturation-associated long intervening non-coding RNA (OLMALINC), which is highly expressed in the white matter (WM) of the human frontal cortex compared to the grey matter (GM) and peripheral tissues. Moreover, we identified a novel isoform of OLMALINC that was also up-regulated in the WM. RNA-interference (RNAi) knockdown of OLMALINC in oligodendrocytes, which are the major cell type in the WM, caused significant changes in the expression of genes regulating cytostructure, cell activation and membrane signaling. Gene ontology enrichment analysis revealed that over 10% of the top 25 up- and down-regulated genes were involved in oligodendrocyte maturation. RNAi experiments in neuronal cells resulted in the perturbation of genes controlling cell proliferation. Furthermore, we identified a novel cis-natural antisense non-coding RNA, which we named OLMALINC-AS, which maps to the first exon of the dominant isoform of OLMALINC.

Conclusions

Our study has demonstrated for the first time that a primate-specific lincRNA regulates the expression of genes critical to human oligodendrocyte maturation, which in turn might be regulated by an antisense counterpart.

Electronic supplementary material

The online version of this article (doi:10.1186/s13041-014-0091-9) contains supplementary material, which is available to authorized users.

Human long noncoding RNAs are substantially less folded than messenger RNAs

Long noncoding RNAs (lncRNAs) do not code for proteins but function as RNAs. Because the functions of an RNA rely on either its sequence or secondary structure, lncRNAs should be folded at least as strongly as messenger RNAs (mRNAs), which serve as messengers for translation and are generally thought to lack secondary structure-dependent RNA-level functions. Contrary to this prediction, analysis of genome-wide experimental data of human RNA folding reveals that lncRNAs are substantially less folded than mRNAs even after the control of expression level and GC% (percentage of guanines and cytosines), although both lncRNAs and mRNAs are more strongly folded than expected by chance. In contrast to mRNAs, lncRNAs show neither the positive correlation between folding strength and expression level nor the negative correlation between folding strength and evolutionary rate. These and other results support that although RNA folding undoubtedly plays a role in RNA biology it is also important in translation and/or protein biology.

Characteristics of long non-coding RNAs in the Brown Norway rat and alterations in the Dahl salt-sensitive rat

Long non-coding RNAs (lncRNAs) are potentially important mediators of genomic regulation. lncRNAs, however, remain poorly characterized in the rat model organism widely used in biomedical research. Using poly(A)-independent and strand-specific RNA-seq, we identified 1,500 to 1,800 lncRNAs expressed in each of the following tissues of Brown Norway rats: the renal cortex, renal outer medulla, liver, cardiac left ventricle, adrenal gland, and hypothalamus. Expression and the binding of histone H3K4me3 to promoter regions were confirmed for several lncRNAs. Rat lncRNA expression appeared to be more tissue-specific than mRNA. Rat lncRNAs had 4.5 times fewer exons and 29% shorter transcripts than mRNA. The median cumulative abundance of rat lncRNAs was 53% of that of mRNA. Approximately 28% of the lncRNAs identified in the renal outer medulla appeared to lack a poly(A) tail. Differential expression of 74 lncRNAs was detected in the renal outer medulla between Dahl SS rats, a model of salt-sensitive hypertension, and salt-insensitive, congenic SS.13^BN26 rats fed a high-salt diet. Two of the differentially expressed lncRNAs, which were confirmed, were located within the congenic region and contained several sequence variants. The study identified genome-wide characteristics of lncRNAs in the rat model and suggested a role of lncRNAs in hypertension.

Identification and diagnostic performance of a small RNA within the PCA3 and BMCC1 gene locus that potentially targets mRNA

BACKGROUND

PCA3 is a long noncoding RNA (lncRNA) with unknown function, upregulated in prostate cancer. LncRNAs may be processed into smaller active species. We hypothesized this for PCA3.

METHODS

We computed feasible RNA hairpins within the BMCC1 gene (encompassing PCA3) and searched a prostate transcriptome for these. We measured expression using qRT-PCR in three cohorts of prostate cancer tissues (n = 60), exfoliated urinary cells (n = 484 with cancer and n = 166 controls), and in cell lines (n = 22). We used in silico predictions and RNA knockup to identify potential mRNA targets of short transcribed RNAs.

RESULTS

We predicted 13 hairpins, of which PCA3-shRNA2 was most abundant within the prostate transcriptome. PCA3-shRNA2 is located within intron 1 of PCA3 and appears regulated by androgens. Expression of PCA3-shRNA2 was upregulated in malignant prostatic tissues, exfoliated urinary cells from men with prostate cancer (13-273 fold change; t test P < 0.003), and closely correlated to PCA3 expression (r = 0.84-0.93; P < 0.001). Urinary PCA3-shRNA2 (C-index, 0.75-0.81) and PCA3 (C-index, 0.78) could predict the presence of cancer in most men. PCA3-shRNA2 knockup altered the expression of predicted target mRNAs, including COPS2, SOX11, WDR48, TEAD1, and Noggin. PCA3-shRNA2 expression was negatively correlated with COPS2 in patient samples (r = -0.32; P < 0.001).

CONCLUSION

We identified a short RNA within PCA3, whose expression is correlated to PCA3, which may target mRNAs implicated in prostate biology.

IMPACT

This short RNA is stable ex vivo, suggesting a role as a robust biomarker. We identify cytoplasmic enrichment of this RNA and potential targeting of mRNAs implicated in prostate carcinogenesis.

Long non-coding RNAs: critical players in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a complex disease with multiple underlying pathogenic mechanisms caused by a variety of etiologic factors. Emerging evidence showed that long non-coding RNAs (lncRNAs), with size larger than 200 nucleotides (nt), play important roles in various types of cancer development and progression. In recent years, some dysregulated lncRNAs in HCC have been revealed and roles for several of them in HCC have been characterized. All these findings point to the potential of lncRNAs as prospective novel therapeutic targets in HCC. In this review, we summarize known dysregulated lncRNAs in HCC, and review potential biological roles and underlying molecular mechanisms of lncRNAs in HCC. Additionally, we discussed prospects of lncRNAs as potential biomarker and therapeutic target for HCC. In conclusion, this paper will help us gain better understanding of molecular mechanisms by which lncRNAs perform their function in HCC and also provide general strategies and directions for future research.

Pervasive transcription: detecting functional RNAs in bacteria

Pervasive, or genome-wide, transcription has been reported in all domains of life. In bacteria, most pervasive transcription occurs antisense to protein-coding transcripts, although recently a new class of pervasive RNAs was identified that originates from within annotated genes. Initially considered to be non-functional transcriptional noise, pervasive transcription is increasingly being recognized as important in regulating gene expression. The function of pervasive transcription is an extensively debated question in the field of transcriptomics and regulatory RNA biology. Here, we highlight the most recent contributions addressing the purpose of pervasive transcription in bacteria and discuss their implications.

Evolutionary dynamics of coding and non-coding transcriptomes

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

Volatile evolution of long noncoding RNA repertoires: mechanisms and biological implications

Thousands of genes encoding long noncoding RNAs (lncRNAs) have been identified in all vertebrate genomes thus far examined. The list of lncRNAs partaking in arguably important biochemical, cellular, and developmental activities is steadily growing. However, it is increasingly clear that lncRNA repertoires are subject to weak functional constraint and rapid turnover during vertebrate evolution. We discuss here some of the factors that may explain this apparent paradox, including relaxed constraint on sequence to maintain lncRNA structure/function, extensive redundancy in the regulatory circuits in which lncRNAs act, as well as adaptive and non-adaptive forces such as genetic drift. We explore the molecular mechanisms promoting the birth and rapid evolution of lncRNA genes, with an emphasis on the influence of bidirectional transcription and transposable elements, two pervasive features of vertebrate genomes. Together these properties reveal a remarkably dynamic and malleable noncoding transcriptome which may represent an important source of robustness and evolvability.

Regulation of transcription by long noncoding RNAs

Over the past decade there has been a greater understanding of genomic complexity in eukaryotes ushered in by the immense technological advances in high-throughput sequencing of DNA and its corresponding RNA transcripts. This has resulted in the realization that beyond protein-coding genes, there are a large number of transcripts that do not encode for proteins and, therefore, may perform their function through RNA sequences and/or through secondary and tertiary structural determinants. This review is focused on the latest findings on a class of noncoding RNAs that are relatively large (>200 nucleotides), display nuclear localization, and use different strategies to regulate transcription. These are exciting times for discovering the biological scope and the mechanism of action for these RNA molecules, which have roles in dosage compensation, imprinting, enhancer function, and transcriptional regulation, with a great impact on development and disease.

Long non-coding RNAs as a source of new peptides

Deep transcriptome sequencing has revealed the existence of many transcripts that lack long or conserved open reading frames (ORFs) and which have been termed long non-coding RNAs (lncRNAs). The vast majority of lncRNAs are lineage-specific and do not yet have a known function. In this study, we test the hypothesis that they may act as a repository for the synthesis of new peptides. We find that a large fraction of the lncRNAs expressed in cells from six different species is associated with ribosomes. The patterns of ribosome protection are consistent with the translation of short peptides. lncRNAs show similar coding potential and sequence constraints than evolutionary young protein coding sequences, indicating that they play an important role in de novo protein evolution.

DOI:http://dx.doi.org/10.7554/eLife.03523.001

Despite the terms being largely interchangeable in modern language, ‘DNA’ and ‘gene’ do not mean the same thing. A gene is made of DNA and contains the instructions to make a protein, and it is the protein that performs the function of the gene. However, cells in the body also contain DNA that does not form genes. Far from being ‘junk’ DNA with no biological purpose; this DNA has a variety of roles, including affecting how other genes are used.

To produce a protein, the DNA sequence of a gene is transcribed into an intermediate molecule called RNA, which is then translated to produce a protein. So-called long non-coding RNA (lncRNA) molecules are also transcribed from DNA, but whether these are translated to make proteins has been a subject of much debate. Indeed, the function of the vast majority of lncRNA molecules is unknown.

Ruiz-Orera et al. analyzed RNA sequences collected from earlier experiments on six different species—humans, mice, fish, flies, yeast, and a plant—and found nearly 2500 as yet unstudied lncRNAs in addition to those previously identified. Many of the lncRNAs that Ruiz-Orera et al. investigated could be found lodged inside the cellular machinery used to translate RNA into proteins. Furthermore, these lncRNA molecules are oriented in the machinery as if they are primed and ready for translation, suggesting that many lncRNAs do produce proteins. However, it is unclear how many of these proteins have a useful function.

Very few lncRNAs were found in more than one species, suggesting that they have evolved recently. The properties of lncRNA molecules also show many similarities with the properties of ‘young’—recently evolved—genes that are known to produce proteins. The combined findings of Ruiz-Orera et al. therefore suggest that lncRNAs are important for developing new proteins. The emergence of proteins with new functions has been an important driving force in evolution, and this work provides important clues into the first steps of this process.

DOI:http://dx.doi.org/10.7554/eLife.03523.002

High-resolution mapping of transcriptional dynamics across tissue development reveals a stable mRNA-tRNA interface

The genetic code is an abstraction of how mRNA codons and tRNA anticodons molecularly interact during protein synthesis; the stability and regulation of this interaction remains largely unexplored. Here, we characterized the expression of mRNA and tRNA genes quantitatively at multiple time points in two developing mouse tissues. We discovered that mRNA codon pools are highly stable over development and simply reflect the genomic background; in contrast, precise regulation of tRNA gene families is required to create the corresponding tRNA transcriptomes. The dynamic regulation of tRNA genes during development is controlled in order to generate an anticodon pool that closely corresponds to messenger RNAs. Thus, across development, the pools of mRNA codons and tRNA anticodons are invariant and highly correlated, revealing a stable molecular interaction interlocking transcription and translation.

Negative regulation of the interferon response by an interferon-induced long non-coding RNA

Long non-coding RNAs (lncRNAs) play critical roles in diverse cellular processes; however, their involvement in many critical aspects of the immune response including the interferon (IFN) response remains poorly understood. To address this gap, we compared the global gene expression pattern of primary human hepatocytes before and at three time points after treatment with IFN-α. Among ∼200 IFN-induced lncRNAs, one transcript showed ∼100-fold induction. This RNA, which we named lncRNA-CMPK2, was a spliced, polyadenylated nuclear transcript that was induced by IFN in diverse cell types from human and mouse. Similar to protein-coding IFN-stimulated genes (ISGs), its induction was dependent on JAK-STAT signaling. Intriguingly, knockdown of lncRNA-CMPK2 resulted in a marked reduction in HCV replication in IFN-stimulated hepatocytes, suggesting that it could affect the antiviral role of IFN. We could show that lncRNA-CMPK2 knockdown resulted in upregulation of several protein-coding antiviral ISGs. The observed upregulation was caused by an increase in both basal and IFN-stimulated transcription, consistent with loss of transcriptional inhibition in knockdown cells. These results indicate that the IFN response involves a lncRNA-mediated negative regulatory mechanism. lncRNA-CMPK2 was strongly upregulated in a subset of HCV-infected human livers, suggesting a role in modulation of the IFN response in vivo.

Long non-coding RNAs in the regulation of the immune response

•

Widespread changes in lncRNA expresssion are associated with the immune response.

•

lncRNAs regulate the inflammatory response following activation of innate immunity.

•

lncRNAs regulate T cell differentiation and migration.

•

The action of long non-coding RNAs is mediated via diverse mechanisms.

It is increasingly clear that long non-coding RNAs (lncRNAs) regulate a variety biological responses, and that they do so by a diverse range of mechanisms. In the field of immunology, recent publications have shown widespread changes in the expression of lncRNAs during the activation of the innate immune response and T cell development, differentiation, and activation. These lncRNAs control important aspects of immunity such as production of inflammatory mediators, differentiation, and cell migration through regulating protein–protein interactions or via their ability to basepair with RNA and DNA. We review the current understanding of the mechanism of action of these immune-related lncRNAs, discuss their impact on physiological and pathological processes, and highlight important areas of inquiry at the intersection between immunology and lncRNA biology.

The transcriptional landscape of mouse beta cells compared to human beta cells reveals notable species differences in long non-coding RNA and protein-coding gene expression

BACKGROUND

Insulin producing beta cell and glucagon producing alpha cells are colocalized in pancreatic islets in an arrangement that facilitates the coordinated release of the two principal hormones that regulate glucose homeostasis and prevent both hypoglycemia and diabetes. However, this intricate organization has also complicated the determination of the cellular source(s) of the expression of genes that are detected in the islet. This reflects a significant gap in our understanding of mouse islet physiology, which reduces the effectiveness by which mice model human islet disease.

RESULTS

To overcome this challenge, we generated a bitransgenic reporter mouse that faithfully labels all beta and alpha cells in mouse islets to enable FACS-based purification and the generation of comprehensive transcriptomes of both populations. This facilitates systematic comparison across thousands of genes between the two major endocrine cell types of the islets of Langerhans whose principal hormones are of cardinal importance for glucose homeostasis. Our data leveraged against similar data for human beta cells reveal a core common beta cell transcriptome of 9900+ genes. Against the backdrop of overall similar beta cell transcriptomes, we describe marked differences in the repertoire of receptors and long non-coding RNAs between mouse and human beta cells.

CONCLUSIONS

The comprehensive mouse alpha and beta cell transcriptomes complemented by the comparison of the global (dis)similarities between mouse and human beta cells represent invaluable resources to boost the accuracy by which rodent models offer guidance in finding cures for human diabetes.

Evolutionary conservation of cold-induced antisense RNAs of FLOWERING LOCUS C in Arabidopsis thaliana perennial relatives

Antisense RNA (asRNA) COOLAIR is expressed at A. thaliana FLOWERING LOCUS C (FLC) in response to winter temperatures. Its contribution to cold-induced silencing of FLC was proposed but its functional and evolutionary significance remain unclear. Here we identify a highly conserved block containing the COOLAIR first exon and core promoter at the 3′ end of several FLC orthologues. Furthermore, asRNAs related to COOLAIR are expressed at FLC loci in the perennials A. alpina and A. lyrata, although some splicing variants differ from A. thaliana. Study of the A. alpina orthologue, PERPETUAL FLOWERING 1 (PEP1), demonstrates that AaCOOLAIR is induced each winter of the perennial life cycle. Introduction of PEP1 into A. thaliana reveals that AaCOOLAIR cis-elements confer cold-inducibility in this heterologous species while the difference between PEP1 and FLC mRNA patterns depends on both cis-elements and species-specific trans-acting factors. Thus, expression of COOLAIR is highly conserved, supporting its importance in FLC regulation.

FLOWERING LOCUS C (FLC) is thought to control the flowering time of A. thaliana in response to winter temperatures, in a process known as vernalization. Here, the authors suggest that the COOLAIR antisense RNA, which is conserved across plant species, acts to repress the expression of FLC during vernalization.

Conserved expression of lincRNA during human and macaque prefrontal cortex development and maturation

The current annotation of the human genome includes more than 12,000 long intergenic noncoding RNAs (lincRNA). While a handful of lincRNA have been shown to play important regulatory roles, the functionality of most remains unclear. Here, we examined the expression conservation and putative functionality of lincRNA in human and macaque prefrontal cortex (PFC) development and maturation. We analyzed transcriptome sequence (RNA-seq) data from 38 human and 40 macaque individuals covering the entire postnatal development interval. Using the human data set, we detected the expression of 5835 lincRNA annotated in GENCODE and further identified 1888 novel lincRNA. Most of these lincRNA show low DNA sequence conservation, as well as low expression levels. Remarkably, developmental expression patterns of these lincRNA were as conserved between humans and macaques as those of protein-coding genes. Transfection of development-associated lincRNA into human SH-SY5Y cells affected gene expression, indicating their regulatory potential. In brain, expression of these putative target genes correlated with the expression of the corresponding lincRNA during human and macaque PFC development. These results support the potential functionality of lincRNA in primate PFC development.

Long non-coding RNAs and enhancer RNAs regulate the lipopolysaccharide-induced inflammatory response in human monocytes

Early reports indicate that long non-coding RNAs (lncRNAs) are novel regulators of biological responses. However, their role in the human innate immune response, which provides the initial defence against infection, is largely unexplored. To address this issue, here we characterize the long non-coding RNA transcriptome in primary human monocytes using RNA sequencing. We identify 76 enhancer RNAs (eRNAs), 40 canonical lncRNAs, 65 antisense lncRNAs and 35 regions of bidirectional transcription (RBT) that are differentially expressed in response to bacterial lipopolysaccharide (LPS). Crucially, we demonstrate that knockdown of nuclear-localized, NF-κB-regulated, eRNAs (IL1β-eRNA) and RBT (IL1β-RBT46) surrounding the IL1β locus, attenuates LPS-induced messenger RNA transcription and release of the proinflammatory mediators, IL1β and CXCL8. We predict that lncRNAs can be important regulators of the human innate immune response.

Intergenic lncRNAs and the evolution of gene expression

Eukaryote genomes encode a surprisingly large number of noncoding transcripts. Around two-thirds of human transcribed loci do not encode protein, and many are intergenic and produce long (>200 nucleotides) noncoding RNAs (lncRNAs). Extensive analyses using comparative genomics and transcriptomics approaches have established that lncRNA sequence and transcription tend to turn over rapidly during evolution. Our appreciation of the biological roles of lncRNAs, based only on a handful of transcripts with well-characterized functions, is that lncRNAs have diverse roles in regulating gene expression. These proposed roles together with their rapid rates of evolution suggest that lncRNAs could contribute to the divergent expression patterns observed among species and potentially to the origin of new traits.

Characterization of the human cumulus cell transcriptome during final follicular maturation and ovulation

Cumulus expansion and oocyte maturation are central processes in ovulation. Knowledge gained from rodent and other mammalian models has revealed some of the molecular pathways associated with these processes. However, the equivalent pathways in humans have not been thoroughly studied and remain unidentified. Compact cumulus cells (CCs) from germinal vesicle cumulus oocyte complexes (COCs) were obtained from patients undergoing in vitro maturation (IVM) procedures. Expanded CCs from metaphase 2 COC were obtained from patients undergoing IVF/ICSI. Global transcriptome profiles of the samples were obtained using state-of-the-art RNA sequencing techniques. We identified 1746 differentially expressed (DE) genes between compact and expanded CCs. Most of these genes were involved in cellular growth and proliferation, cellular movement, cell cycle, cell-to-cell signaling and interaction, extracellular matrix and steroidogenesis. Out of the DE genes, we found 89 long noncoding RNAs, of which 12 are encoded within introns of genes known to be involved in granulosa cell processes. This suggests that unique noncoding RNA transcripts may contribute to the regulation of cumulus expansion and oocyte maturation. Using global transcriptome sequencing, we were able to generate a library of genes regulated during cumulus expansion and oocyte maturation processes. Analysis of these genes allowed us to identify important new genes and noncoding RNAs potentially involved in COC maturation and cumulus expansion. These results may increase our understanding of the process of oocyte maturation and could ultimately improve the efficacy of IVM treatment.

No gene in the genome makes sense except in the light of evolution

Evolutionary conservation has been an accurate predictor of functional elements across the first decade of metazoan genomics. More recently, there has been a move to define functional elements instead from biochemical annotations. Evolutionary methods are, however, more comprehensive than biochemical approaches can be and can assess quantitatively, especially for subtle effects, how biologically important--how injurious after mutation--different types of elements are. Evolutionary methods are thus critical for understanding the large fraction (up to 10%) of the human genome that does not encode proteins and yet might convey function. These methods can also capture the ephemeral nature of much noncoding functional sequence, with large numbers of functional elements having been gained and lost rapidly along each mammalian lineage. Here, we review how different strengths of purifying selection have impacted on protein-coding and non-protein-coding loci and on transcription factor binding sites in mammalian and fruit fly genomes.

The dark matter rises: the expanding world of regulatory RNAs

The ability to sequence genomes and characterize their products has begun to reveal the central role for regulatory RNAs in biology, especially in complex organisms. It is now evident that the human genome contains not only protein-coding genes, but also tens of thousands of non-protein coding genes that express small and long ncRNAs (non-coding RNAs). Rapid progress in characterizing these ncRNAs has identified a diverse range of subclasses, which vary widely in size, sequence and mechanism-of-action, but share a common functional theme of regulating gene expression. ncRNAs play a crucial role in many cellular pathways, including the differentiation and development of cells and organs and, when mis-regulated, in a number of diseases. Increasing evidence suggests that these RNAs are a major area of evolutionary innovation and play an important role in determining phenotypic diversity in animals.

Identification of novel non-coding RNA-based negative feedback regulating the expression of the oncogenic transcription factor GLI1

Non‐coding RNAs are a complex class of nucleic acids, with growing evidence supporting regulatory roles in gene expression. Here we identify a non‐coding RNA located head‐to‐head with the gene encoding the Glioma‐associated oncogene 1 (GLI1), a transcriptional effector of multiple cancer‐associated signaling pathways. The expression of this three‐exon GLI1 antisense (GLI1AS) RNA in cancer cells was concordant with GLI1 levels. siRNAs knockdown of GLI1AS up‐regulated GLI1 and increased cellular proliferation and tumor growth in a xenograft model system. Conversely, GLI1AS overexpression decreased the levels of GLI1, its target genes PTCH1 and PTCH2, and cellular proliferation. Additionally, we demonstrate that GLI1 knockdown reduced GLI1AS, while GLI1 overexpression increased GLI1AS, supporting the role of GLI1AS as a target gene of the GLI1 transcription factor. Activation of TGFβ and Hedgehog signaling, two known regulators of GLI1 expression, conferred a concordant up‐regulation of GLI1 and GLI1AS in cancer cells. Finally, analysis of the mechanism underlying the interplay between GLI1 and GLI1AS indicates that the non‐coding RNA elicits a local alteration of chromatin structure by increasing the silencing mark H3K27me3 and decreasing the recruitment of RNA polymerase II to this locus. Taken together, the data demonstrate the existence of a novel non‐coding RNA‐based negative feedback loop controlling GLI1 levels, thus expanding the repertoire of mechanisms regulating the expression of this oncogenic transcription factor.

A novel negative feedback loop on Hedgehog signaling is demonstrated.

The mechanism involves a non‐coding RNA antisense to the GLI1 gene, GLI1AS.

GLI1AS is shown to be a target gene of the GLI1 transcription factor.

GLI1AS represses gene expression at the GLI1/GLI1AS locus.

GLI1AS acts as an epigenetic modifier eliciting repressive chromatin marks.

Functional and topological characteristics of mammalian regulatory domains

Long-range regulatory interactions play an important role in shaping gene-expression programs. However, the genomic features that organize these activities are still poorly characterized. We conducted a large operational analysis to chart the distribution of gene regulatory activities along the mouse genome, using hundreds of insertions of a regulatory sensor. We found that enhancers distribute their activities along broad regions and not in a gene-centric manner, defining large regulatory domains. Remarkably, these domains correlate strongly with the recently described TADs, which partition the genome into distinct self-interacting blocks. Different features, including specific repeats and CTCF-binding sites, correlate with the transition zones separating regulatory domains, and may help to further organize promiscuously distributed regulatory influences within large domains. These findings support a model of genomic organization where TADs confine regulatory activities to specific but large regulatory domains, contributing to the establishment of specific gene expression profiles.

Evolutionary conservation of long non-coding RNAs; sequence, structure, function

BACKGROUND

Recent advances in genomewide studies have revealed the abundance of long non-coding RNAs (lncRNAs) in mammalian transcriptomes. The ENCODE Consortium has elucidated the prevalence of human lncRNA genes, which are as numerous as protein-coding genes. Surprisingly, many lncRNAs do not show the same pattern of high interspecies conservation as protein-coding genes. The absence of functional studies and the frequent lack of sequence conservation therefore make functional interpretation of these newly discovered transcripts challenging. Many investigators have suggested the presence and importance of secondary structural elements within lncRNAs, but mammalian lncRNA secondary structure remains poorly understood. It is intriguing to speculate that in this group of genes, RNA secondary structures might be preserved throughout evolution and that this might explain the lack of sequence conservation among many lncRNAs.

SCOPE OF REVIEW

Here, we review the extent of interspecies conservation among different lncRNAs, with a focus on a subset of lncRNAs that have been functionally investigated. The function of lncRNAs is widespread and we investigate whether different forms of functionalities may be conserved.

MAJOR CONCLUSIONS

Lack of conservation does not imbue a lack of function. We highlight several examples of lncRNAs where RNA structure appears to be the main functional unit and evolutionary constraint. We survey existing genomewide studies of mammalian lncRNA conservation and summarize their limitations. We further review specific human lncRNAs which lack evolutionary conservation beyond primates but have proven to be both functional and therapeutically relevant.

GENERAL SIGNIFICANCE

Pioneering studies highlight a role in lncRNAs for secondary structures, and possibly the presence of functional "modules", which are interspersed with longer and less conserved stretches of nucleotide sequences. Taken together, high-throughput analysis of conservation and functional composition of the still-mysterious lncRNA genes is only now becoming feasible.

Lineage and species-specific long noncoding RNAs during erythro-megakaryocytic development

Mammals express thousands of long noncoding (lnc) RNAs, a few of which are known to function in tissue development. However, the entire repertoire of lncRNAs in most tissues and species is not defined. Indeed, most lncRNAs are not conserved, raising questions about function. We used RNA sequencing to identify 1109 polyadenylated lncRNAs expressed in erythroblasts, megakaryocytes, and megakaryocyte-erythroid precursors of mice, and 594 in erythroblasts of humans. More than half of these lncRNAs were unannotated, emphasizing the opportunity for new discovery through studies of specialized cell types. Analysis of the mouse erythro-megakaryocytic polyadenylated lncRNA transcriptome indicates that ~75% arise from promoters and 25% from enhancers, many of which are regulated by key transcription factors including GATA1 and TAL1. Erythroid lncRNA expression is largely conserved among 8 different mouse strains, yet only 15% of mouse lncRNAs are expressed in humans and vice versa, reflecting dramatic species-specificity. RNA interference assays of 21 abundant erythroid-specific murine lncRNAs in primary mouse erythroid precursors identified 7 whose knockdown inhibited terminal erythroid maturation. At least 6 of these 7 functional lncRNAs have no detectable expression in human erythroblasts, suggesting that lack of conservation between mammalian species does not predict lack of function.

The evolution of lncRNA repertoires and expression patterns in tetrapods

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

Evolutionary dynamics and tissue specificity of human long noncoding RNAs in six mammals

Long intergenic noncoding RNAs (lincRNAs) play diverse regulatory roles in human development and disease, but little is known about their evolutionary history and constraint. Here, we characterize human lincRNA expression patterns in nine tissues across six mammalian species and multiple individuals. Of the 1898 human lincRNAs expressed in these tissues, we find orthologous transcripts for 80% in chimpanzee, 63% in rhesus, 39% in cow, 38% in mouse, and 35% in rat. Mammalian-expressed lincRNAs show remarkably strong conservation of tissue specificity, suggesting that it is selectively maintained. In contrast, abundant splice-site turnover suggests that exact splice sites are not critical. Relative to evolutionarily young lincRNAs, mammalian-expressed lincRNAs show higher primary sequence conservation in their promoters and exons, increased proximity to protein-coding genes enriched for tissue-specific functions, fewer repeat elements, and more frequent single-exon transcripts. Remarkably, we find that ∼20% of human lincRNAs are not expressed beyond chimpanzee and are undetectable even in rhesus. These hominid-specific lincRNAs are more tissue specific, enriched for testis, and faster evolving within the human lineage.

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.

Chromatin signatures at transcriptional start sites separate two equally populated yet distinct classes of intergenic long noncoding RNAs

BACKGROUND

Mammalian transcriptomes contain thousands of long noncoding RNAs (lncRNAs). Some lncRNAs originate from intragenic enhancers which, when active, behave as alternative promoters producing transcripts that are processed using the canonical signals of their host gene. We have followed up this observation by analyzing intergenic lncRNAs to determine the extent to which they might also originate from intergenic enhancers.

RESULTS

We integrated high-resolution maps of transcriptional initiation and transcription to annotate a conservative set of intergenic lncRNAs expressed in mouse erythroblasts. We subclassified intergenic lncRNAs according to chromatin status at transcriptional initiation regions, defined by relative levels of histone H3K4 mono- and trimethylation. These transcripts are almost evenly divided between those arising from enhancer-associated (elncRNA) or promoter-associated (plncRNA) elements. These two classes of 5' capped and polyadenylated RNA transcripts are indistinguishable with regard to their length, number of exons or transcriptional orientation relative to their closest neighboring gene. Nevertheless, elncRNAs are more tissue-restricted, less highly expressed and less well conserved during evolution. Of considerable interest, we found that expression of elncRNAs, but not plncRNAs, is associated with enhanced expression of neighboring protein-coding genes during erythropoiesis.

CONCLUSIONS

We have determined globally the sites of initiation of intergenic lncRNAs in erythroid cells, allowing us to distinguish two similarly abundant classes of transcripts. Different correlations between the levels of elncRNAs, plncRNAs and expression of neighboring genes suggest that functional lncRNAs from the two classes may play contrasting roles in regulating the transcript abundance of local or distal loci.

Identification of novel transcripts and noncoding RNAs in bovine skin by deep next generation sequencing

BACKGROUND

Deep RNA sequencing (RNAseq) has opened a new horizon for understanding global gene expression. The functional annotation of non-model mammalian genomes including bovines is still poor compared to that of human and mouse. This particularly applies to tissues without direct significance for milk and meat production, like skin, in spite of its multifunctional relevance for the individual. Thus, applying an RNAseq approach, we performed a whole transcriptome analysis of pigmented and nonpigmented bovine skin to describe the comprehensive transcript catalogue of this tissue.

RESULTS

A total of 39,577 unique primary skin transcripts were mapped to the bovine reference genome assembly. The majority of the transcripts were mapped to known transcriptional units (65%). In addition to the reannotation of known genes, a substantial number (10,884) of unknown transcripts (UTs) were discovered, which had not previously been annotated. The classification of UTs was based on the prediction of their coding potential and comparative sequence analysis, subsequently followed by meticulous manual curation. The classification analysis and experimental validation of selected UTs confirmed that RNAseq data can be used to amend the annotation of known genes by providing evidence for additional exons, untranslated regions or splice variants, by approving genes predicted in silico and by identifying novel bovine loci. A large group of UTs (4,848) was predicted to potentially represent long noncoding RNA (lncRNA). Predominantly, potential lncRNAs mapped in intergenic chromosome regions (4,365) and therefore, were classified as potential intergenic lncRNA. Our analysis revealed that only about 6% of all UTs displayed interspecies conservation and discovered a variety of unknown transcripts without interspecies homology but specific expression in bovine skin.

CONCLUSIONS

The results of our study demonstrate a complex transcript pattern for bovine skin and suggest a possible functional relevance of novel transcripts, including lncRNA, in the modulation of pigmentation processes. The results also indicate that the comprehensive identification and annotation of unknown transcripts from whole transcriptome analysis using RNAseq data remains a tremendous future challenge.

Gene regulation by antisense transcription

Antisense transcription, which was initially considered by many as transcriptional noise, is increasingly being recognized as an important regulator of gene expression. It is widespread among all kingdoms of life and has been shown to influence - either through the act of transcription or through the non-coding RNA that is produced - almost all stages of gene expression, from transcription and translation to RNA degradation. Antisense transcription can function as a fast evolving regulatory switch and a modular scaffold for protein complexes, and it can 'rewire' regulatory networks. The genomic arrangement of antisense RNAs opposite sense genes indicates that they might be part of self-regulatory circuits that allow genes to regulate their own expression.

AHT-ChIP-seq: a completely automated robotic protocol for high-throughput chromatin immunoprecipitation

ChIP-seq is an established manually-performed method for identifying DNA-protein interactions genome-wide. Here, we describe a protocol for automated high-throughput (AHT) ChIP-seq. To demonstrate the quality of data obtained using AHT-ChIP-seq, we applied it to five proteins in mouse livers using a single 96-well plate, demonstrating an extremely high degree of qualitative and quantitative reproducibility among biological and technical replicates. We estimated the optimum and minimum recommended cell numbers required to perform AHT-ChIP-seq by running an additional plate using HepG2 and MCF7 cells. With this protocol, commercially available robotics can perform four hundred experiments in five days.

Role of antisense RNAs in evolution of yeast regulatory complexity

Antisense RNAs (asRNAs) are known to regulate gene expression. However, a genome-wide mechanism of asRNA regulation is unclear, and there is no good explanation why partial asRNAs are not functional. To explore its regulatory role, we investigated asRNAs using an evolutionary approach, as genome-wide experimental data are limited. We found that the percentage of genes coupling with asRNAs in Saccharomyces cerevisiae is negatively associated with regulatory complexity and evolutionary age. Nevertheless, asRNAs evolve more slowly when their sense genes are under more complex regulation. Older genes coupling with asRNAs are more likely to demonstrate inverse expression, reflecting the role of these asRNAs as repressors. Our analyses provide novel evidence, suggesting a minor contribution of asRNAs in developing regulatory complexity. Although our results support the leaky hypothesis for asRNA transcription, our evidence also suggests that partial asRNAs may have evolved as repressors. Our study deepens the understanding of asRNA regulatory evolution.

Functional transcriptomics in the post-ENCODE era

The last decade has seen tremendous effort committed to the annotation of the human genome sequence, most notably perhaps in the form of the ENCODE project. One of the major findings of ENCODE, and other genome analysis projects, is that the human transcriptome is far larger and more complex than previously thought. This complexity manifests, for example, as alternative splicing within protein-coding genes, as well as in the discovery of thousands of long noncoding RNAs. It is also possible that significant numbers of human transcripts have not yet been described by annotation projects, while existing transcript models are frequently incomplete. The question as to what proportion of this complexity is truly functional remains open, however, and this ambiguity presents a serious challenge to genome scientists. In this article, we will discuss the current state of human transcriptome annotation, drawing on our experience gained in generating the GENCODE gene annotation set. We highlight the gaps in our knowledge of transcript functionality that remain, and consider the potential computational and experimental strategies that can be used to help close them. We propose that an understanding of the true overlap between transcriptional complexity and functionality will not be gained in the short term. However, significant steps toward obtaining this knowledge can now be taken by using an integrated strategy, combining all of the experimental resources at our disposal.

Integrated analysis of dysregulated lncRNA expression in fetal cardiac tissues with ventricular septal defect

Ventricular septal defects (VSD) are the most common form of congenital heart disease, which is the leading non-infectious cause of death in children; nevertheless, the exact cause of VSD is not yet fully understood. Long non-coding RNAs (lncRNAs) have been shown to play key roles in various biological processes, such as imprinting control, circuitry controlling pluripotency and differentiation, immune responses and chromosome dynamics. Notably, a growing number of lncRNAs have been implicated in disease etiology, although an association with VSD has not been reported. In the present study, we conducted an integrated analysis of dysregulated lncRNAs, focusing specifically on the identification and characterization of lncRNAs potentially involving in initiation of VSD. Comparison of the transcriptome profiles of cardiac tissues from VSD-affected and normal hearts was performed using a second-generation lncRNA microarray, which covers the vast majority of expressed RefSeq transcripts (29,241 lncRNAs and 30,215 coding transcripts). In total, 880 lncRNAs were upregulated and 628 were downregulated in VSD. Furthermore, our established filtering pipeline indicated an association of two lncRNAs, ENST00000513542 and RP11-473L15.2, with VSD. This dysregulation of the lncRNA profile provides a novel insight into the etiology of VSD and furthermore, illustrates the intricate relationship between coding and ncRNA transcripts in cardiac development. These data may offer a background/reference resource for future functional studies of lncRNAs related to VSD.

Long non-coding RNA expression profiling of mouse testis during postnatal development

Mammalian testis development and spermatogenesis play critical roles in male fertility and continuation of a species. Previous research into the molecular mechanisms of testis development and spermatogenesis has largely focused on the role of protein-coding genes and small non-coding RNAs, such as microRNAs and piRNAs. Recently, it has become apparent that large numbers of long (>200 nt) non-coding RNAs (lncRNAs) are transcribed from mammalian genomes and that lncRNAs perform important regulatory functions in various developmental processes. However, the expression of lncRNAs and their biological functions in post-natal testis development remain unknown. In this study, we employed microarray technology to examine lncRNA expression profiles of neonatal (6-day-old) and adult (8-week-old) mouse testes. We found that 8,265 lncRNAs were expressed above background levels during post-natal testis development, of which 3,025 were differentially expressed. Candidate lncRNAs were identified for further characterization by an integrated examination of genomic context, gene ontology (GO) enrichment of their associated protein-coding genes, promoter analysis for epigenetic modification, and evolutionary conservation of elements. Many lncRNAs overlapped or were adjacent to key transcription factors and other genes involved in spermatogenesis, such as Ovol1, Ovol2, Lhx1, Sox3, Sox9, Plzf, c-Kit, Wt1, Sycp2, Prm1 and Prm2. Most differentially expressed lncRNAs exhibited epigenetic modification marks similar to protein-coding genes and tend to be expressed in a tissue-specific manner. In addition, the majority of differentially expressed lncRNAs harbored evolutionary conserved elements. Taken together, our findings represent the first systematic investigation of lncRNA expression in the mammalian testis and provide a solid foundation for further research into the molecular mechanisms of lncRNAs function in mammalian testis development and spermatogenesis.

Sense-antisense gene pairs: sequence, transcription, and structure are not conserved between human and mouse

Previous efforts to characterize conservation between the human and mouse genomes focused largely on sequence comparisons. These studies are inherently limited because they don't account for gene structure differences, which may exist despite genomic sequence conservation. Recent high-throughput transcriptome studies have revealed widespread and extensive overlaps between genes, and transcripts, encoded on both strands of the genomic sequence. This overlapping gene organization, which produces sense-antisense (SAS) gene pairs, is capable of effecting regulatory cascades through established mechanisms. We present an evolutionary conservation assessment of SAS pairs, on three levels: genomic, transcriptomic, and structural. From a genome-wide dataset of human SAS pairs, we first identified orthologous loci in the mouse genome, then assessed their transcription in the mouse, and finally compared the genomic structures of SAS pairs expressed in both species. We found that approximately half of human SAS loci have single orthologous locations in the mouse genome; however, only half of those orthologous locations have SAS transcriptional activity in the mouse. This suggests that high human-mouse gene conservation overlooks widespread distinctions in SAS pair incidence and expression. We compared gene structures at orthologous SAS loci, finding frequent differences in gene structure between human and orthologous mouse SAS pair members. Our categorization of human SAS pairs with respect to mouse conservation of expression as well as structure points to limitations of mouse models. Gene structure differences, including at SAS loci, may account for some of the phenotypic distinctions between primates and rodents. Genes in non-conserved SAS pairs may contribute to evolutionary lineage-specific regulatory outcomes.

lincRNAs: genomics, evolution, and mechanisms

Long intervening noncoding RNAs (lincRNAs) are transcribed from thousands of loci in mammalian genomes and might play widespread roles in gene regulation and other cellular processes. This Review outlines the emerging understanding of lincRNAs in vertebrate animals, with emphases on how they are being identified and current conclusions and questions regarding their genomics, evolution and mechanisms of action.

Long non-coding RNAs and complex human diseases

Long non-coding RNAs (lncRNAs) are a heterogeneous class of RNAs that are generally defined as non-protein-coding transcripts longer than 200 nucleotides. Recently, an increasing number of studies have shown that lncRNAs can be involved in various critical biological processes, such as chromatin remodeling, gene transcription, and protein transport and trafficking. Moreover, lncRNAs are dysregulated in a number of complex human diseases, including coronary artery diseases, autoimmune diseases, neurological disorders, and various cancers, which indicates their important roles in these diseases. Here, we reviewed the current understanding of lncRNAs, including their definition and subclassification, regulatory functions, and potential roles in different types of complex human diseases.

Identification and characterization of long non-coding RNAs related to mouse embryonic brain development from available transcriptomic data

Long non-coding RNAs (lncRNAs) as a key group of non-coding RNAs have gained widely attention. Though lncRNAs have been functionally annotated and systematic explored in higher mammals, few are under systematical identification and annotation. Owing to the expression specificity, known lncRNAs expressed in embryonic brain tissues remain still limited. Considering a large number of lncRNAs are only transcribed in brain tissues, studies of lncRNAs in developmental brain are therefore of special interest. Here, publicly available RNA-sequencing (RNA-seq) data in embryonic brain are integrated to identify thousands of embryonic brain lncRNAs by a customized pipeline. A significant proportion of novel transcripts have not been annotated by available genomic resources. The putative embryonic brain lncRNAs are shorter in length, less spliced and show less conservation than known genes. The expression of putative lncRNAs is in one tenth on average of known coding genes, while comparable with known lncRNAs. From chromatin data, putative embryonic brain lncRNAs are associated with active chromatin marks, comparable with known lncRNAs. Embryonic brain expressed lncRNAs are also indicated to have expression though not evident in adult brain. Gene Ontology analysis of putative embryonic brain lncRNAs suggests that they are associated with brain development. The putative lncRNAs are shown to be related to possible cis-regulatory roles in imprinting even themselves are deemed to be imprinted lncRNAs. Re-analysis of one knockdown data suggests that four regulators are associated with lncRNAs. Taken together, the identification and systematic analysis of putative lncRNAs would provide novel insights into uncharacterized mouse non-coding regions and the relationships with mammalian embryonic brain development.

Evolutionarily conserved long intergenic non-coding RNAs in the eye

The discovery that the mammalian transcriptome encodes thousands of long intergenic non-coding (linc) RNA transcripts, together with recent evidence that lincRNAs can regulate protein-coding genes, has added a new level of complexity to cellular transcriptional/translational regulation. Indeed several reports now link mutations in lincRNAs to heritable human disorders. Here, we identified a subset of lincRNAs in terminally differentiated adult human retinal neurons based on their sequence conservation across species. RNA sequencing of eye tissue from several mammalian species with varied rod/cone photoreceptor content identified 18 lincRNAs that were highly conserved across these species. Sixteen of the 18 were conserved in human retinal tissue with 14 of these also conserved in the macular region. A subset of lincRNAs exhibited restricted tissue expression profiles in mice, with preferential expression in the retina. Mouse models with different populations of retinal cells as well as in situ hybridization provided evidence that these lincRNAs localized to specific retinal compartments, most notably to the photoreceptor neuronal layer. Computational genomic loci and promoter region analyses provided a basis for regulated expression of these conserved lincRNAs in retinal post-mitotic neurons. This combined approach identified several lincRNAs that could be critical for retinal and visual maintenance in adults.

Detecting and comparing non-coding RNAs in the high-throughput era

In recent years there has been a growing interest in the field of non-coding RNA. This surge is a direct consequence of the discovery of a huge number of new non-coding genes and of the finding that many of these transcripts are involved in key cellular functions. In this context, accurately detecting and comparing RNA sequences has become important. Aligning nucleotide sequences is a key requisite when searching for homologous genes. Accurate alignments reveal evolutionary relationships, conserved regions and more generally any biologically relevant pattern. Comparing RNA molecules is, however, a challenging task. The nucleotide alphabet is simpler and therefore less informative than that of amino-acids. Moreover for many non-coding RNAs, evolution is likely to be mostly constrained at the structural level and not at the sequence level. This results in very poor sequence conservation impeding comparison of these molecules. These difficulties define a context where new methods are urgently needed in order to exploit experimental results to their full potential. This review focuses on the comparative genomics of non-coding RNAs in the context of new sequencing technologies and especially dealing with two extremely important and timely research aspects: the development of new methods to align RNAs and the analysis of high-throughput data.

Differential expression of long noncoding RNAs in the livers of female B6C3F1 mice exposed to the carcinogen furan

The mammalian genome is transcribed into mRNAs that code for protein and a broad spectrum of other noncoding (nc) RNA products. Long ncRNAs (lncRNA), defined as ncRNA species > 200 nucleotides long, are emerging as important epigenetic regulators of gene expression that are involved in a spectrum of biological processes of relevance to toxicology. We conducted a gene expression profiling study in the livers of female B6C3F1 mice exposed to the carcinogen furan at 0.0, 1.0, and 2.0mg/kg (noncarcinogenic doses) and at 4.0 and 8.0mg/kg (carcinogenic doses) for 3 weeks. LncRNA differential expression showed a nonlinear dose response with none differentially expressed at 1.0 or 2.0mg/kg, 2 lncRNAs at 4.0mg/kg furan, and 83 at 8mg/kg, representing 13.3% (83/632) of the total number of differentially expressed transcripts. Among the lncRNAs observed, two lncRNAs examined showed transcriptional clustering with nearby protein-coding genes. LincRNA-p21 is an antisense transcript that is 15kb downstream from Cdkn1a locus and appears to be cotranscribed with the protein coding gene Cdkn1a at 8.0mg/kg furan. In a separate independent study, RNA samples from the livers of mice administered benzo(a)pyrene also demonstrated increased levels of Cdkn1a and the antisense lincRNA-p21 transcript. These data demonstrate that lncRNAs are transcriptional targets of furan exposures associated with levels of furan that are cytotoxic and induce cell proliferation. In addition, certain lncRNA transcripts are associated with the expression of nearby coding protein genes. We hypothesize that lncRNAs have potential as epigenetic biomarkers of carcinogenic exposures.

Mutations within lncRNAs are effectively selected against in fruitfly but not in human

Background

Previous studies in Drosophila and mammals have revealed levels of long non-coding RNAs (lncRNAs) sequence conservation that are intermediate between neutrally evolving and protein-coding sequence. These analyses compared conservation between species that diverged up to 75 million years ago. However, analysis of sequence polymorphisms within a species' population can provide an understanding of essentially contemporaneous selective constraints that are acting on lncRNAs and can quantify the deleterious effect of mutations occurring within these loci.

Results

We took advantage of polymorphisms derived from the genome sequences of 163 Drosophila melanogaster strains and 174 human individuals to calculate the distribution of fitness effects of single nucleotide polymorphisms occurring within intergenic lncRNAs and compared this to distributions for SNPs present within putatively neutral or protein-coding sequences. Our observations show that in D.melanogaster there is a significant excess of rare frequency variants within intergenic lncRNAs relative to neutrally evolving sequences, whereas selection on human intergenic lncRNAs appears to be effectively neutral. Approximately 30% of mutations within these fruitfly lncRNAs are estimated as being weakly deleterious.

Conclusions

These contrasting results can be attributed to the large difference in effective population sizes between the two species. Our results suggest that while the sequences of lncRNAs will be well conserved across insect species, such loci in mammals will accumulate greater proportions of deleterious changes through genetic drift.

New gene expression pipelines gush lncRNAs

Genome-wide techniques provide robust and comprehensive identification of lncRNAs in adult mouse neural stem cells and their derivatives, illuminating the functions of these underappreciated transcripts.

Long noncoding RNAs in biology and hematopoiesis

Genome and transcriptome sequencing have revealed a rich assortment of noncoding RNAs in eukaryote cells, including long noncoding RNAs (lncRNAs), which regulate gene expression independent of protein coding potential. LncRNAs modulate protein coding gene expression in many cell types by regulating multiple processes, including epigenetic control of transcription, mRNA stability, and protein localization. Although little is known about lncRNAs in hematopoiesis, they are likely to exert widespread roles in this process.