Post-transcriptional processing generates a diversity of 5'-modified long and short RNAs. Nature. 2009;457:1028-1032. [PMID: 19169241 DOI: 10.1038/nature07759]

You Don't Muck with MYC

MYC homeostasis is critical for major cellular and organismal processes. The physiological and pathologic patterns of c-myc transcription are programmed by a large number of cis-elements and transfactors (RNAs and proteins). These elements and factors receive inputs from a multitude of intracellular and extracellular pathways. Because c-myc regulation has customarily been dissected element by element and factor by factor, it has been difficult to appreciate how the c-myc promoter and regulatory sequences operate as a system. A full accounting of the regulation of c-myc transcription will require an understanding of the dynamic interplay of these factors and elements with one another, with chromatin, and with the changes in DNA structure and topology that are inevitably coupled with gene activity.

Diverse endonucleolytic cleavage sites in the mammalian transcriptome depend upon microRNAs, Drosha, and additional nucleases

The life span of a mammalian mRNA is determined, in part, by the binding of regulatory proteins and small RNA-guided complexes. The conserved endonuclease activity of Argonaute2 requires extensive complementarity between a small RNA and its target and is not used by animal microRNAs, which pair with their targets imperfectly. Here we investigate the endonucleolytic function of Ago2 and other nucleases by transcriptome-wide profiling of mRNA cleavage products retaining 5' phosphate groups in mouse embryonic stem cells (mESCs). We detect a prominent signature of Ago2-dependent cleavage events and validate several such targets. Unexpectedly, a broader class of Ago2-independent cleavage sites is also observed, indicating participation of additional nucleases in site-specific mRNA cleavage. Within this class, we identify a cohort of Drosha-dependent mRNA cleavage events that functionally regulate mRNA levels in mESCs, including one in the Dgcr8 mRNA. Together, these results highlight the underappreciated role of endonucleolytic cleavage in controlling mRNA fates in mammals.

Genome-wide characterization of methylguanosine-capped and polyadenylated small RNAs in the rice blast fungus Magnaporthe oryzae

Small RNAs are well described in higher eukaryotes such as mammals and plants; however, knowledge in simple eukaryotes such as filamentous fungi is limited. In this study, we discovered and characterized methylguanosine-capped and polyadenylated small RNAs (CPA-sRNAs) by using differential RNA selection, full-length cDNA cloning and 454 transcriptome sequencing of the rice blast fungus Magnaporthe oryzae. This fungus causes blast, a devastating disease on rice, the principle food staple for over half the world’s population. CPA-sRNAs mapped primarily to the transcription initiation and termination sites of protein-coding genes and were positively correlated with gene expression, particularly for highly expressed genes including those encoding ribosomal proteins. Numerous CPA-sRNAs also mapped to rRNAs, tRNAs, snRNAs, transposable elements and intergenic regions. Many other 454 sequence reads could not be mapped to the genome; however, inspection revealed evidence for non-template additions and chimeric sequences. CPA-sRNAs were independently confirmed using a high affinity variant of eIF-4E to capture 5′-methylguanosine-capped RNA followed by 3′-RACE sequencing. These results expand the repertoire of small RNAs in filamentous fungi.

Nuclear-localized tiny RNAs are associated with transcription initiation and splice sites in metazoans

We have recently shown that transcription initiation RNAs (tiRNAs) are derived from sequences immediately downstream of transcription start sites. Here, using cytoplasmic and nuclear small RNA high-throughput sequencing datasets, we report the identification of a second class of nuclear-specific approximately 17- to 18-nucleotide small RNAs whose 3' ends map precisely to the splice donor site of internal exons in animals. These splice-site RNAs (spliRNAs) are associated with highly expressed genes and show evidence of developmental stage- and region-specific expression. We also show that tiRNAs are localized to the nucleus, are enriched at chromatin marks associated with transcription initiation and possess a 3'-nucleotide bias. Additionally, we find that microRNA-offset RNAs (moRNAs), the miR-15/16 cluster previously linked to oncosuppression and most small nucleolar RNA (snoRNA)-derived small RNAs (sdRNAs) are enriched in the nucleus, whereas most miRNAs and two H/ACA sdRNAs are cytoplasmically enriched. We propose that nuclear-localized tiny RNAs are involved in the epigenetic regulation of gene expression.

Tissue-specific transcript annotation and expression profiling with complementary next-generation sequencing technologies

Next-generation sequencing is excellently suited to evaluate the abundance of mRNAs to study gene expression. Here we compare two alternative technologies, cap analysis of gene expression (CAGE) and serial analysis of gene expression (SAGE), for the same RNA samples. Along with quantifying gene expression levels, CAGE can be used to identify tissue-specific transcription start sites, while SAGE monitors 3′-end usage. We used both methods to get more insight into the transcriptional control of myogenesis, studying differential gene expression in differentiated and proliferating C2C12 myoblast cells with statistical evaluation of reproducibility and differential gene expression. Both CAGE and SAGE provided highly reproducible data (Pearson's correlations >0.92 among biological triplicates). With both methods we found around 10 000 genes expressed at levels 2 transcripts per million (0.3 copies per cell), with an overlap of 86%. We identified 4304 and 3846 genes differentially expressed between proliferating and differentiated C2C12 cells by CAGE and SAGE, respectively, with an overlap of 2144. We identified 196 novel regulatory regions with preferential use in proliferating or differentiated cells. Next-generation sequencing of CAGE and SAGE libraries provides consistent expression levels and can enrich current genome annotations with tissue-specific promoters and alternative 3′-UTR usage.

Long noncoding RNAs in mammalian cells: what, where, and why?

Not all long, polyadenylated cellular RNAs encode polypeptides. In recent years, it has become apparent that a number of organisms express abundant amounts of transcripts that lack open reading frames or that are retained in the nucleus. Rather than accumulating silently in the cell, we now know that many of these long noncoding RNAs (lncRNAs) play important roles in nuclear architecture or in the regulation of gene expression. Here, we discuss some recent progress in our understanding of the functions of a number of important lncRNAs in mammalian cells.

Next-generation genomics: an integrative approach

Integrating results from diverse experiments is an essential process in our effort to understand the logic of complex systems, such as development, homeostasis and responses to the environment. With the advent of high-throughput methods--including genome-wide association (GWA) studies, chromatin immunoprecipitation followed by sequencing (ChIP-seq) and RNA sequencing (RNA-seq)--acquisition of genome-scale data has never been easier. Epigenomics, transcriptomics, proteomics and genomics each provide an insightful, and yet one-dimensional, view of genome function; integrative analysis promises a unified, global view. However, the large amount of information and diverse technology platforms pose multiple challenges for data access and processing. This Review discusses emerging issues and strategies related to data integration in the era of next-generation genomics.

Short RNAs are transcribed from repressed polycomb target genes and interact with polycomb repressive complex-2

Polycomb proteins maintain cell identity by repressing the expression of developmental regulators specific for other cell types. Polycomb repressive complex-2 (PRC2) catalyzes trimethylation of histone H3 lysine-27 (H3K27me3). Although repressed, PRC2 targets are generally associated with the transcriptional initiation marker H3K4me3, but the significance of this remains unclear. Here, we identify a class of short RNAs, approximately 50-200 nucleotides in length, transcribed from the 5' end of polycomb target genes in primary T cells and embryonic stem cells. Short RNA transcription is associated with RNA polymerase II and H3K4me3, occurs in the absence of mRNA transcription, and is independent of polycomb activity. Short RNAs form stem-loop structures resembling PRC2 binding sites in Xist, interact with PRC2 through SUZ12, cause gene repression in cis, and are depleted from polycomb target genes activated during cell differentiation. We propose that short RNAs play a role in the association of PRC2 with its target genes.

Linking promoters to functional transcripts in small samples with nanoCAGE and CAGEscan

Large-scale sequencing projects have revealed an unexpected complexity in the origins, structures and functions of mammalian transcripts. Many loci are known to produce overlapping coding and noncoding RNAs with capped 5' ends that vary in size. Methods to identify the 5' ends of transcripts will facilitate the discovery of new promoters and 5' ends derived from secondary capping events. Such methods often require high input amounts of RNA not obtainable from highly refined samples such as tissue microdissections and subcellular fractions. Therefore, we developed nano-cap analysis of gene expression (nanoCAGE), a method that captures the 5' ends of transcripts from as little as 10 ng of total RNA, and CAGEscan, a mate-pair adaptation of nanoCAGE that captures the transcript 5' ends linked to a downstream region. Both of these methods allow further annotation-agnostic studies of the complex human transcriptome.

Re-inspection of small RNA sequence datasets reveals several novel human miRNA genes

Background

miRNAs are key players in gene expression regulation. To fully understand the complex nature of cellular differentiation or initiation and progression of disease, it is important to assess the expression patterns of as many miRNAs as possible. Thereby, identifying novel miRNAs is an essential prerequisite to make possible a comprehensive and coherent understanding of cellular biology.

Methodology/Principal Findings

Based on two extensive, but previously published, small RNA sequence datasets from human embryonic stem cells and human embroid bodies, respectively [1], we identified 112 novel miRNA-like structures and were able to validate miRNA processing in 12 out of 17 investigated cases. Several miRNA candidates were furthermore substantiated by including additional available small RNA datasets, thereby demonstrating the power of combining datasets to identify miRNAs that otherwise may be assigned as experimental noise.

Conclusions/Significance

Our analysis highlights that existing datasets are not yet exhaustedly studied and continuous re-analysis of the available data is important to uncover all features of small RNA sequencing.

A paired-end sequencing strategy to map the complex landscape of transcription initiation

Recent high-throughput sequencing protocols have uncovered the complexity of mammalian transcription by RNA polymerase II, helping to define several initiation patterns in which transcription start sites (TSSs) cluster within both narrow and broad genomic windows. Here, we describe a paired-end sequencing strategy, which enables more robust mapping and characterization of capped transcripts. This strategy was applied to explore the transcription initiation landscape in the Drosophila melanogaster embryo. Extending the previous findings in mammals, we found that fly promoters exhibit distinct initiation patterns, which are linked to specific promoter sequence motifs. Furthermore, we identified a large number of 5′ capped transcripts originating from coding exons; analyses support that they are unlikely the result of alternative TSSs, but rather the product of post-transcriptional modifications. Taken together, paired-end TSS analysis is demonstrated to be a powerful method to uncover the transcriptional complexity of eukaryotic genomes.

Most "dark matter" transcripts are associated with known genes

Short-read RNA sequencing in mouse and human tissues shows that most transcripts are encoded within or nearby known genes and that most of the genome is not transcribed.

A series of reports over the last few years have indicated that a much larger portion of the mammalian genome is transcribed than can be accounted for by currently annotated genes, but the quantity and nature of these additional transcripts remains unclear. Here, we have used data from single- and paired-end RNA-Seq and tiling arrays to assess the quantity and composition of transcripts in PolyA+ RNA from human and mouse tissues. Relative to tiling arrays, RNA-Seq identifies many fewer transcribed regions (“seqfrags”) outside known exons and ncRNAs. Most nonexonic seqfrags are in introns, raising the possibility that they are fragments of pre-mRNAs. The chromosomal locations of the majority of intergenic seqfrags in RNA-Seq data are near known genes, consistent with alternative cleavage and polyadenylation site usage, promoter- and terminator-associated transcripts, or new alternative exons; indeed, reads that bridge splice sites identified 4,544 new exons, affecting 3,554 genes. Most of the remaining seqfrags correspond to either single reads that display characteristics of random sampling from a low-level background or several thousand small transcripts (median length = 111 bp) present at higher levels, which also tend to display sequence conservation and originate from regions with open chromatin. We conclude that, while there are bona fide new intergenic transcripts, their number and abundance is generally low in comparison to known exons, and the genome is not as pervasively transcribed as previously reported.

The human genome was sequenced a decade ago, but its exact gene composition remains a subject of debate. The number of protein-coding genes is much lower than initially expected, and the number of distinct transcripts is much larger than the number of protein-coding genes. Moreover, the proportion of the genome that is transcribed in any given cell type remains an open question: results from “tiling” microarray analyses suggest that transcription is pervasive and that most of the genome is transcribed, whereas new deep sequencing-based methods suggest that most transcripts originate from known genes. We have addressed this discrepancy by comparing samples from the same tissues using both technologies. Our analyses indicate that RNA sequencing appears more reliable for transcripts with low expression levels, that most transcripts correspond to known genes or are near known genes, and that many transcripts may represent new exons or aberrant products of the transcription process. We also identify several thousand small transcripts that map outside known genes; their sequences are often conserved and are often encoded in regions of open chromatin. We propose that most of these transcripts may be by-products of the activity of enhancers, which associate with promoters as part of their role as long-range gene regulatory sites. Overall, however, we find that most of the genome is not appreciably transcribed.

Mining mammalian transcript data for functional long non-coding RNAs

Background

The role of long non-coding RNAs (lncRNAs) in controlling gene expression has garnered increased interest in recent years. Sequencing projects, such as Fantom3 for mouse and H-InvDB for human, have generated abundant data on transcribed components of mammalian cells, the majority of which appear not to be protein-coding. However, much of the non-protein-coding transcriptome could merely be a consequence of ‘transcription noise’. It is therefore essential to use bioinformatic approaches to identify the likely functional candidates in a high throughput manner.

Principal Findings

We derived a scheme for classifying and annotating likely functional lncRNAs in mammals. Using the available experimental full-length cDNA data sets for human and mouse, we identified 78 lncRNAs that are either syntenically conserved between human and mouse, or that originate from the same protein-coding genes. Of these, 11 have significant sequence homology. We found that these lncRNAs exhibit: (i) patterns of codon substitution typical of non-coding transcripts; (ii) preservation of sequences in distant mammals such as dog and cow, (iii) significant sequence conservation relative to their corresponding flanking regions (in 50% cases, flanking regions do not have homology at all; and in the remaining, the degree of conservation is significantly less); (iv) existence mostly as single-exon forms (8/11); and, (v) presence of conserved and stable secondary structure motifs within them. We further identified orthologous protein-coding genes that are contributing to the pool of lncRNAs; of which, genes implicated in carcinogenesis are significantly over-represented.

Conclusion

Our comparative mammalian genomics approach coupled with evolutionary analysis identified a small population of conserved long non-protein-coding RNAs (lncRNAs) that are potentially functional across Mammalia. Additionally, our analysis indicates that amongst the orthologous protein-coding genes that produce lncRNAs, those implicated in cancer pathogenesis are significantly over-represented, suggesting that these lncRNAs could play an important role in cancer pathomechanisms.

Transcriptional regulation of the Alström syndrome gene ALMS1 by members of the RFX family and Sp1

Mutations in the human gene ALMS1 cause Alström syndrome, a disorder characterised by neurosensory degeneration, metabolic defects and cardiomyopathy. ALMS1 encodes a centrosomal protein implicated in the assembly and maintenance of primary cilia. Expression of ALMS1 varies between tissues and recent data suggest that its transcription is modulated during adipogenesis and growth arrest. However the ALMS1 promoter has not been defined. This study focused on identifying and characterising the ALMS1 proximal promoter, initially by using 5' RACE to map transcription start sites. Luciferase reporter assay and EMSA data strongly suggest that ALMS1 transcription is regulated by the ubiquitous factor Sp1. In addition, reporter assay, EMSA, chromatin immunoprecipitation and RNA interference data indicate that ALMS1 transcription is regulated by regulatory factor X (RFX) proteins. These transcription factors are cell-type restricted in their expression profile and known to regulate genes of the ciliogenic pathway. We show binding of RFX proteins to an evolutionarily conserved X-box in the ALMS1 proximal promoter and present evidence that these proteins are responsible for ALMS1 transcription during growth arrest induced by low serum conditions. In summary, this work provides the first data on transcription factors regulating general and context-specific transcription of the disease-associated gene ALMS1.

Small RNAs in the animal gonad: guarding genomes and guiding development

Germ cells must safeguard, apportion, package, and deliver their genomes with exquisite precision to ensure proper reproduction and embryonic development. Classical genetic approaches have identified many genes controlling animal germ cell development, but only recently have some of these genes been linked to the RNA interference (RNAi) pathway, a gene silencing mechanism centered on small regulatory RNAs. Germ cells contain microRNAs (miRNAs), endogenous siRNAs (endo-siRNAs), and Piwi-interacting RNAs (piRNAs); these are bound by members of the Piwi/Argonaute protein family. piwi genes were known to specify germ cell development, but we now understand that mutations disrupting germline development can also affect small RNA accumulation. Small RNA studies in germ cells have revealed a surprising diversity of regulatory mechanisms and a unifying function for germline genes in controlling the spread of transposable elements. Future challenges will be to understand the production of germline small RNAs and to identify the full breadth of gene regulation by these RNAs. Progress in this area will likely impact biomedical goals of manipulating stem cells and preventing diseases caused by the transposition of mobile DNA elements.

RNA dust: where are the genes?

Initial gene discovery efforts through analysis of genome sequences and identification and characterization of expressed RNAs have revealed that only a relatively small portion of the genome is transcribed into protein coding mRNAs in vertebrates. However, in contrast with this paucity of protein coding ‘genes’, there is an enormous complexity in transcription and the protein coding mRNAs contribute to a very small fraction of transcripts in comparison with the different varieties of non-coding RNAs (ncRNAs). This transcriptome complexity may be hypothesized to have a regulatory role that is required for the development and function of organisms as complex as vertebrates. At the same time, it raises the fundamental question of the unequivocal definition of a gene. It is intriguing to postulate that many ncRNAs might finely modulate gene activity by acting as regulatory elements. The emerging hypotheses suggest that the gene regulatory machinery may be deeply interconnected with the world of short RNAs. These RNAs may generally act for fine-tuning of the protein-coding transcriptome.

Non-coding RNAs: regulators of disease

For 50 years the term 'gene' has been synonymous with regions of the genome encoding mRNAs that are translated into protein. However, recent genome-wide studies have shown that the human genome is pervasively transcribed and produces many thousands of regulatory non-protein-coding RNAs (ncRNAs), including microRNAs, small interfering RNAs, PIWI-interacting RNAs and various classes of long ncRNAs. It is now clear that these RNAs fulfil critical roles as transcriptional and post-transcriptional regulators and as guides of chromatin-modifying complexes. Here we review the biology of ncRNAs, focusing on the fundamental mechanisms by which ncRNAs facilitate normal development and physiology and, when dysfunctional, underpin disease. We also discuss evidence that intergenic regions associated with complex diseases express ncRNAs, as well as the potential use of ncRNAs as diagnostic markers and therapeutic targets. Taken together, these observations emphasize the need to move beyond the confines of protein-coding genes and highlight the fact that continued investigation of ncRNA biogenesis and function will be necessary for a comprehensive understanding of human disease.

High DNA melting temperature predicts transcription start site location in human and mouse

The accurate computational prediction of transcription start sites (TSS) in vertebrate genomes is a difficult problem. The physicochemical properties of DNA can be computed in various ways and a many combinations of DNA features have been tested in the past for use as predictors of transcription. We looked in detail at melting temperature, which measures the temperature, at which two strands of DNA separate, considering the cooperative nature of this process. We find that peaks in melting temperature correspond closely to experimentally determined transcription start sites in human and mouse chromosomes. Using melting temperature alone, and with simple thresholding, we can predict TSS with accuracy that is competitive with the most accurate state-of-the-art TSS prediction methods. Accuracy is measured using both experimentally and manually determined TSS. The method works especially well with CpG island containing promoters, but also works when CpG islands are absent. This result is clear evidence of the important role of the physical properties of DNA in the process of transcription. It also points to the importance for TSS prediction methods to include melting temperature as prior information.

Antisense transcription in gammaretroviruses as a mechanism of insertional activation of host genes

Transcription of retroviruses is initiated at the U3-R region boundary in the integrated provirus and continues unidirectionally to produce genomic and mRNA products of positive polarity. Several studies have recently demonstrated the existence of naturally occurring protein-encoding transcripts of negative polarity in complex retroviruses. We report here on the identification of transcripts of negative polarity in simple murine leukemia virus (MLV). In T-cell and B-cell lymphomas induced by SL3-3 and Akv MLV, antisense transcripts initiated in the U3 region of the proviral 5' long terminal repeat (LTR) and continued into the cellular proto-oncogenes Jdp2 and Bach2 to create chimeric transcripts consisting of viral and host sequence. The phenomenon was validated in vivo using a knock-in mouse model homozygous for a single LTR at a position known to activate Nras in B-cell lymphomas. A 5' rapid amplification of cDNA ends (RACE) analysis indicated a broad spectrum of initiation sites within the U3 region of the 5' LTR. Our data show for the first time transcriptional activity of negative polarity initiating in the U3 region of simple retroviruses and suggest a novel mechanism of insertional activation of host genes. Elucidation of the nature and potential regulatory role of 5' LTR antisense transcription will be relevant to the design of therapeutic vectors and may contribute to the increasing recognition of pervasive eukaryotic transcription.

Identification and characterization of microRNAs and endogenous siRNAs in Schistosoma japonicum

BACKGROUND

Small endogenous non-coding RNAs (sncRNAs) such as small interfering RNA (siRNA), microRNA and other small RNA transcripts are derived from distinct loci in the genome and play critical roles in RNA-mediated gene silencing mechanisms in plants and metazoa. They are approximately 22 nucleotides long; regulate mRNA stability through perfect or imperfect match to the targets. The biological activities of sncRNAs have been related to many biological events, from resistance to microbe infections to cellular differentiation. The development of the zoonotic parasite Schistosoma japonicum parasite includes multiple steps of morphological alterations and biological differentiations, which provide a unique model for studies on the functions of small RNAs. Characterization of the genome-wide transcription of the sncRNAs will be a major step in understanding of the parasite biology. The objective of this study is to investigate the transcriptional profile and potential function of the small non-coding RNAs in the development of S. japanicum.

RESULTS

The endogenous siRNAs were found mainly derived from transposable elements (TE) or transposons and the natural antisense transcripts (NAT). In contrast to other organisms, the TE-derived siRNAs in S. japonicum were more predominant than other sncRNAs including microRNAs (miRNAs). Further, there were distinct length and 3'end variations in the sncRNAs, which were associated with the developmental differentiation of the parasite. Among the identified miRNA transcripts, there were 38 unique to S. japonicum and 16 that belonged to 13 miRNA families are common to other metazoan lineages. These miRNAs were either ubiquitously expressed, or they exhibited specific expression patterns related to the developmental stages or sex. Genes that encoded miRNAs are mainly located in clusters within the genome of S. japonicum. However, genes within one cluster could be differentially transcribed, which suggested that individual genes might be regulated by distinct mechanisms during parasite development.

CONCLUSIONS

Many miRNA and endogenous siRNA transcripts were identified in S. japonicum and the amount of siRNA was at least 4.4 and 1.6 times more than that of miRNA in both schistosomulum and adult worm stages respectively. SiRNAs are mainly derived from transposable elements (or transposons); while natural antisense transcripts (NAT)-derived siRNAs were much less. A majority of miRNA transcripts identified in the parasite were species-specific and the expression of certain miRNAs was found developmentally regulated. Both miRNA and siRNAs are potentially important regulators in the development of schistosomal parasites.

Novel RNA-based strategies for therapeutic gene silencing

The past decade has seen intense scientific interest in non-coding RNAs. In particular, the discovery and subsequent exploitation of gene silencing via RNA interference (RNAi) has revolutionized the way in which gene expression is now studied and understood. It is now well established that post-transcriptional gene silencing (PTGS) by the microRNA (miRNA) and other RNAi-associated pathways represents an essential layer of complexity to gene regulation. Gene silencing using RNAi additionally demonstrates huge potential as a therapeutic strategy for eliminating pathogenic gene expression. Yet despite the early promise and excitement of gene-specific silencing, several critical hurdles remain to be overcome before widespread clinical adoption. These include off-target effects, toxicity due to saturation of the endogenous RNAi functions, limited duration of silencing, and effective targeted delivery. In recent years, a range of novel strategies for producing RNA-mediated silencing have been developed that can circumvent many of these hurdles, including small internally segmented interfering RNAs, tandem hairpin RNAs, and pri-miRNA cluster mimics. This review discusses RNA-mediated silencing in light of this recent research, and highlights the benefits and limitations conferred by these novel gene-silencing strategies.

Global analysis of short RNAs reveals widespread promoter-proximal stalling and arrest of Pol II in Drosophila

Emerging evidence indicates that gene expression in higher organisms is regulated by RNA polymerase II stalling during early transcription elongation. To probe the mechanisms responsible for this regulation, we developed methods to isolate and characterize short RNAs derived from stalled RNA polymerase II in Drosophila cells. Significant levels of these short RNAs were generated from more than one-third of all genes, indicating that promoter-proximal stalling is a general feature of early polymerase elongation. Nucleotide composition of the initially transcribed sequence played an important role in promoting transcriptional stalling by rendering polymerase elongation complexes highly susceptible to backtracking and arrest. These results indicate that the intrinsic efficiency of early elongation can greatly affect gene expression.

deepBase: a database for deeply annotating and mining deep sequencing data

Advances in high-throughput next-generation sequencing technology have reshaped the transcriptomic research landscape. However, exploration of these massive data remains a daunting challenge. In this study, we describe a novel database, deepBase, which we have developed to facilitate the comprehensive annotation and discovery of small RNAs from transcriptomic data. The current release of deepBase contains deep sequencing data from 185 small RNA libraries from diverse tissues and cell lines of seven organisms: human, mouse, chicken, Ciona intestinalis, Drosophila melanogaster, Caenhorhabditis elegans and Arabidopsis thaliana. By analyzing ∼14.6 million unique reads that perfectly mapped to more than 284 million genomic loci, we annotated and identified ∼380 000 unique ncRNA-associated small RNAs (nasRNAs), ∼1.5 million unique promoter-associated small RNAs (pasRNAs), ∼4.0 million unique exon-associated small RNAs (easRNAs) and ∼6 million unique repeat-associated small RNAs (rasRNAs). Furthermore, 2038 miRNA and 1889 snoRNA candidates were predicted by miRDeep and snoSeeker. All of the mapped reads can be grouped into about 1.2 million RNA clusters. For the purpose of comparative analysis, deepBase provides an integrative, interactive and versatile display. A convenient search option, related publications and other useful information are also provided for further investigation. deepBase is available at: http://deepbase.sysu.edu.cn/.

The complex eukaryotic transcriptome: unexpected pervasive transcription and novel small RNAs

Over the past few years, techniques have been developed that have allowed the study of transcriptomes without bias from previous genome annotations, which has led to the discovery of a plethora of unexpected RNAs that have no obvious coding capacities. There are many different kinds of products that are generated by this pervasive transcription; this Review focuses on small non-coding RNAs (ncRNAs) that have been found to be associated with promoters in eukaryotes from animals to yeast. After comparing the different classes of such ncRNAs described in various studies, the Review discusses how the models proposed for their origins and their possible functions challenge previous views of the basic transcription process and its regulation.

Annotating non-coding transcription using functional genomics strategies

Non-coding RNA (ncRNA) transcripts are RNA molecules that do not code for proteins, but elicit function by other mechanisms. The vast majority of RNA produced in a cell is non-coding ribosomal RNA, produced from relatively few loci, however more recently complementary DNA (cDNA) cloning, tag sequencing, and genome tiling array studies suggest that ncRNAs also account for the majority of RNA species produced by a cell. ncRNA based regulation has been referred to as a 'hidden layer' of signals or 'dark matter' that control gene expression in cellular processes by poorly described mechanisms. These terms have appeared as ncRNAs until recently have been ignored by expression profiling and cDNA annotation projects and their mode of action is diverse (e.g. influencing chromatin structure and epigenetics, translational silencing, transcriptional silencing). Here, we highlight recent functional genomics strategies toward identifying and assigning function to ncRNA transcription.

Establishing legitimacy and function in the new transcriptome

The last decade has seen an explosion of interest in new classes of non-coding RNA. While some are now firmly established as new categories of legitimate functional RNAs, the purpose and even existence of others remain to be solidified. Here, we discuss the challenges associated with discovery and characterization of non-traditional categories of non-coding RNA.

Evolutionarily stable association of intronic snoRNAs and microRNAs with their host genes

Small nucleolar RNAs (snoRNAs) and microRNAs (miRNAs) are integral to a range of processes, including ribosome biogenesis and gene regulation. Some are intron encoded, and this organization may facilitate coordinated coexpression of host gene and RNA. However, snoRNAs and miRNAs are known to be mobile, so intron-RNA associations may not be evolutionarily stable. We have used genome alignments across 11 mammals plus chicken to examine positional orthology of snoRNAs and miRNAs and report that 21% of annotated snoRNAs and 11% of miRNAs are positionally conserved across mammals. Among RNAs traceable to the bird–mammal common ancestor, 98% of snoRNAs and 76% of miRNAs are intronic. Comparison of the most evolutionarily stable mammalian intronic snoRNAs with those positionally conserved among primates reveals that the former are more overrepresented among host genes involved in translation or ribosome biogenesis and are more broadly and highly expressed. This stability is likely attributable to a requirement for overlap between host gene and intronic snoRNA expression profiles, consistent with an ancestral role in ribosome biogenesis. In contrast, whereas miRNA positional conservation is comparable to that observed for snoRNAs, intronic miRNAs show no obvious association with host genes of a particular functional category, and no statistically significant differences in host gene expression are found between those traceable to mammalian or primate ancestors. Our results indicate evolutionarily stable associations of numerous intronic snoRNAs and miRNAs and their host genes, with probable continued diversification of snoRNA function from an ancestral role in ribosome biogenesis.

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.

Enhancing non-coding RNA information content with ADAR editing

The depth and complexity of the non-coding transcriptome in nervous system tissues provides a rich substrate for adenosine de-amination acting on RNA (ADAR). Non-coding RNAs (ncRNAs) serve diverse regulatory and computational functions, coupling signal flow from the environment to evolutionarily coded analog and digital information elements within the transcriptome. We present a perspective of ADARs interaction with the non-coding transcriptome as a computational matrix, enhancing the information processing power of the cell, adding flexibility, rapid response, and fine tuning to critical pathways. Dramatic increases in ADAR activity during stress response and inflammation result in powerful information processing events that change the functional state of the cell. This review examines the pathways and mechanisms of ADAR interaction with the non-coding transcriptome, and their functional consequences for information processing in nervous system tissues.

Genomic analysis suggests that mRNA destabilization by the microprocessor is specialized for the auto-regulation of Dgcr8

BACKGROUND

The Microprocessor, containing the RNA binding protein Dgcr8 and RNase III enzyme Drosha, is responsible for processing primary microRNAs to precursor microRNAs. The Microprocessor regulates its own levels by cleaving hairpins in the 5'UTR and coding region of the Dgcr8 mRNA, thereby destabilizing the mature transcript.

METHODOLOGY/PRINCIPAL FINDINGS

To determine whether the Microprocessor has a broader role in directly regulating other coding mRNA levels, we integrated results from expression profiling and ultra high-throughput deep sequencing of small RNAs. Expression analysis of mRNAs in wild-type, Dgcr8 knockout, and Dicer knockout mouse embryonic stem (ES) cells uncovered mRNAs that were specifically upregulated in the Dgcr8 null background. A number of these transcripts had evolutionarily conserved predicted hairpin targets for the Microprocessor. However, analysis of deep sequencing data of 18 to 200nt small RNAs in mouse ES, HeLa, and HepG2 indicates that exonic sequence reads that map in a pattern consistent with Microprocessor activity are unique to Dgcr8.

CONCLUSION/SIGNIFICANCE

We conclude that the Microprocessor's role in directly destabilizing coding mRNAs is likely specifically targeted to Dgcr8 itself, suggesting a specialized cellular mechanism for gene auto-regulation.

Re-capping the message

The 5'-cap structure that typifies all polymerase II-transcribed RNAs plays important roles in pre-mRNA processing and mRNA export, translation and quality control. Removal of the cap is a regulated process that is considered to be the first irreversible step in mRNA decay. An emerging view challenges this idea: mRNAs have been identified in mammalian cells that lack sequences from their 5' ends but nevertheless appear to be modified with a cap or cap-like structure. Furthermore, a cytoplasmic form of capping enzyme was recently identified that, together with a novel kinase, generates capped ends from cleaved RNAs. These and other findings provide evidence for re-capping and its possible functions.

The plasticity of the mammalian transcriptome

The dogmatic view of RNA as a mere necessity in the transfer of information between DNA and proteins has during recent years come into question. Novel approaches and new technology has revealed an unprecedented level of inherent complexity in the mammalian transcriptome. Here, the majority of nucleotides are expressed, in sharp contrast to the approximately 1.2% of the human genome harboring protein coding information. Also, >50% of genomic loci contain antisense and interleaved transcription, a conservative estimate since non-coding RNA is highly regulated between tissues and developmental stages, which has only been investigated to a limited extent. Subsequent focus on RNA with no coding potential has revealed numerous species with novel functions, and deep sequencing studies imply that many remain to be discovered. This review gives an overview of the plasticity and dynamics of the mammalian transcriptome and the prevailing interpretation of its effect on the complexity of species.

Mapping accessible chromatin regions using Sono-Seq

Disruptions in local chromatin structure often indicate features of biological interest such as regulatory regions. We find that sonication of cross-linked chromatin, when combined with a size-selection step and massively parallel short-read sequencing, can be used as a method (Sono-Seq) to map locations of high chromatin accessibility in promoter regions. Sono-Seq sites frequently correspond to actively transcribed promoter regions, as evidenced by their co-association with RNA Polymerase II ChIP regions, transcription start sites, histone H3 lysine 4 trimethylation (H3K4me3) marks, and CpG islands; signals over other sites, such as those bound by the CTCF insulator, are also observed. The pattern of breakage by Sono-Seq overlaps with, but is distinct from, that observed for FAIRE and DNase I hypersensitive sites. Our results demonstrate that Sono-Seq can be a useful and simple method by which to map many local alterations in chromatin structure. Furthermore, our results provide insights into the mapping of binding sites by using ChIP-Seq experiments and the value of reference samples that should be used in such experiments.

Endo-siRNAs depend on a new isoform of loquacious and target artificially introduced, high-copy sequences

Colonization of genomes by a new selfish genetic element is detrimental to the host species and must lead to an efficient, repressive response. In vertebrates as well as in Drosophila, piRNAs repress transposons in the germ line, whereas endogenous siRNAs take on this role in somatic cells. We show that their biogenesis depends on a new isoform of the Drosophila TRBP homologue loquacious, which arises by alternative polyadenylation and is distinct from the one that functions during the biogenesis of miRNAs. For endo-siRNAs and piRNAs, it is unclear how an efficient response can be initiated de novo. Our experiments establish that the endo-siRNA pathway will target artificially introduced sequences without the need for a pre-existing template in the genome. This response is also triggered in transiently transfected cells, thus genomic integration is not essential. Deep sequencing showed that corresponding endo-siRNAs are generated throughout the sequence, but preferentially from transcribed regions. One strand of the dsRNA precursor can come from spliced mRNA, whereas the opposite strand derives from independent transcripts in antisense orientation.

MIWI-independent small RNAs (MSY-RNAs) bind to the RNA-binding protein, MSY2, in male germ cells

The germ cell-specific DNA/RNA-binding protein MSY2 binds small RNAs (MSY-RNAs) that are approximately 25-31 nt in length, often initiate with a 5' adenine, and are expressed in both germ cells and somatic cells. MSY-RNA levels do not decrease in Miwi or Msy2 null mice. Most MSY-RNAs map within annotated genes, but some are PIWI-interacting RNA (piRNA)-like and map to piRNA clusters. MSY-RNAs are in both nuclei and cytoplasm. In nuclei, MSY-RNAs are enriched in chromatin, and in the cytoplasm they are detected in both ribonucleoproteins and polysomes.

Repertoire of bovine miRNA and miRNA-like small regulatory RNAs expressed upon viral infection

MicroRNA (miRNA) and other types of small regulatory RNAs play a crucial role in the regulation of gene expression in eukaryotes. Several distinct classes of small regulatory RNAs have been discovered in recent years. To extend the repertoire of small RNAs characterized in mammals and to examine relationship between host miRNA expression and viral infection we used Illumina's ultrahigh throughput sequencing approach. We sequenced three small RNA libraries prepared from cell line derived from the adult bovine kidney under normal conditions and upon infection of the cell line with Bovine herpesvirus 1. We used a bioinformatics approach to distinguish authentic mature miRNA sequences from other classes of small RNAs and short RNA fragments represented in the sequencing data. Using this approach we detected 219 out of 356 known bovine miRNAs and 115 respective miRNA* sequences. In addition we identified five new bovine orthologs of known mammalian miRNAs and discovered 268 new cow miRNAs many of which are not identifiable in other mammalian genomes and thus might be specific to the ruminant lineage. In addition we found seven new bovine mirtron candidates. We also discovered 10 small nucleolar RNA (snoRNA) loci that give rise to small RNA with possible miRNA-like function. Results presented in this study extend our knowledge of the biology and evolution of small regulatory RNAs in mammals and illuminate mechanisms of small RNA biogenesis and function. New miRNA sequences and the original sequencing data have been submitted to miRNA repository (miRBase) and NCBI GEO archive respectively. We envisage that these resources will facilitate functional annotation of the bovine genome and promote further functional and comparative genomics studies of small regulatory RNA in mammals.

Trans-acting antisense RNAs mediate transcriptional gene cosuppression in S. cerevisiae

Homology-dependent gene silencing, a phenomenon described as cosuppression in plants, depends on siRNAs. We provide evidence that in Saccharomyces cerevisiae, which is missing the RNAi machinery, protein coding gene cosuppression exists. Indeed, introduction of an additional copy of PHO84 on a plasmid or within the genome results in the cosilencing of both the transgene and the endogenous gene. This repression is transcriptional and position-independent and requires trans-acting antisense RNAs. Antisense RNAs induce transcriptional gene silencing both in cis and in trans, and the two pathways differ by the implication of the Hda1/2/3 complex. We also show that trans-silencing is influenced by the Set1 histone methyltransferase, which promotes antisense RNA production. Finally we show that although antisense-mediated cis-silencing occurs in other genes, trans-silencing so far depends on features specific to PHO84. All together our data highlight the importance of noncoding RNAs in mediating RNAi-independent transcriptional gene silencing.

Long noncoding RNAs: functional surprises from the RNA world

Most of the eukaryotic genome is transcribed, yielding a complex network of transcripts that includes tens of thousands of long noncoding RNAs with little or no protein-coding capacity. Although the vast majority of long noncoding RNAs have yet to be characterized thoroughly, many of these transcripts are unlikely to represent transcriptional "noise" as a significant number have been shown to exhibit cell type-specific expression, localization to subcellular compartments, and association with human diseases. Here, we highlight recent efforts that have identified a myriad of molecular functions for long noncoding RNAs. In some cases, it appears that simply the act of noncoding RNA transcription is sufficient to positively or negatively affect the expression of nearby genes. However, in many cases, the long noncoding RNAs themselves serve key regulatory roles that were assumed previously to be reserved for proteins, such as regulating the activity or localization of proteins and serving as organizational frameworks of subcellular structures. In addition, many long noncoding RNAs are processed to yield small RNAs or, conversely, modulate how other RNAs are processed. It is thus becoming increasingly clear that long noncoding RNAs can function via numerous paradigms and are key regulatory molecules in the cell.

Epigenetics of human T cells during the G0-->G1 transition

We investigated functional epigenetic changes that occur in primary human T lymphocytes during entry into the cell cycle and mapped these at the single-nucleosome level by ChIP-chip on tiling arrays for chromosomes 1 and 6. We show that nucleosome loss and flanking active histone marks define active transcriptional start sites (TSSs). Moreover, these signatures are already set at many inducible genes in quiescent cells prior to cell stimulation. In contrast, there is a dearth of the inactive histone mark H3K9me3 at the TSS, and under-representation of H3K9me2 and H3K9me3 defines the body of active genes. At the DNA level, cytosine methylation (meC) is enriched for nucleosomes that remain at the TSS, whereas in general there is a dearth of meC at TSSs. Furthermore, a drop in meC also marks 3' transcription termination, and a peak of meC occurs at stop codons. This mimics the 3' nucleosomal distribution in yeast, which we show does not occur in human T cells.

Characterization of two alternative Interleukin(IL)-10 5'UTR mRNA sequences, induced by lipopolysaccharide (LPS) stimulation of peripheral blood mononuclear cells

IL-10 production shows a broad-spectrum of individual response, suggesting a genetic component of approximately 75%. Different polymorphisms located close to, or within the IL-10 gene has been demonstrated to influence its transcription rate whereas the post-transcriptional regulation of IL-10 production has not well elucidated. The main responsible elements at this control level are both the 5'- and 3'-untranslated regions (UTR's) of mRNAs, and as the 3'-UTR regions are mainly involved in the stability and decay rate of mRNAs, the 5'-UTR regions mediate the binding rate of the molecule with ribosomal 40S subunit as a cis-acting element. Herein are report data on the identification of two IL10 mRNA that differ by the length of respective 5'UTR regions (160 and 288 nucleotides, respectively; EMBL accession nrs: EU751618 and EU751619) produced after stimulation of human blood samples with bacterial lipopolysaccharide (LPS). The longer 5'UTR is constitutively expressed in unstimulated PBMC cells cultured at 37 degrees C for 24h, while in LPS stimulated cells an additional IL-10 mRNA molecule, containing a shorter 5'UTR, is synthesized. RNADRAW software (http://www.rnadraw.com/) analysis have indicated that the secondary structures of the shorter 5'UTR IL-10 mRNA region is more available for the binding to the 40S ribosomal subunit. In conclusion, our data seem to suggest that LPS could influence the post-transcriptional control of IL-10 production inducing an alternative mRNA immediately available in response to the inflammatory stimulation.

Studying chromosome-wide transcriptional networks: new insights into disease?

A large amount of experimental data collected over the last decade has shown that genomic organization is very complex and has highlighted the fact that the current set of gene annotations does not fully capture this complexity. Much of the RNA detected in a cell is found to originate from outside the exons of annotated genes. Exons of annotated and unannotated transcripts separated by large genomic distances can be joined together in chimeric transcripts. Any given base-pair in a genome could be traversed by many protein-coding and non-coding RNAs. We discuss the implications of these effects for our understanding of disease.

The transcriptome of human CD34+ hematopoietic stem-progenitor cells

Studying gene expression at different hematopoietic stages provides insights for understanding the genetic basis of hematopoiesis. We analyzed gene expression in human CD34(+) hematopoietic cells that represent the stem-progenitor population (CD34(+) cells). We collected >459,000 transcript signatures from CD34(+) cells, including the de novo-generated 3' ESTs and the existing sequences of full-length cDNAs, ESTs, and serial analysis of gene expression (SAGE) tags, and performed an extensive annotation on this large set of CD34(+) transcript sequences. We determined the genes expressed in CD34(+) cells, verified the known genes and identified the new genes of different functional categories involved in hematopoiesis, dissected the alternative gene expression including alternative transcription initiation, splicing, and adenylation, identified the antisense and noncoding transcripts, determined the CD34(+) cell-specific gene expression signature, and developed the CD34(+) cell-transcription map in the human genome. Our study provides a current view on gene expression in human CD34(+) cells and reveals that early hematopoiesis is an orchestrated process with the involvement of over half of the human genes distributed in various functions. The data generated from our study provide a comprehensive and uniform resource for studying hematopoiesis and stem cell biology.

Ubiquitous antisense transcription in eukaryotes: novel regulatory mechanism or byproduct of opportunistic RNA polymerase?

One of the more notable observations made in the last few years in gene regulation is that eukaryotic genomes appear to be pervasively transcribed. Recent transcriptome mapping studies have shown that much of the genome is transcribed, and in some instances transcripts from both strands of specific genomic loci are detectable. While some of these transcripts map to known RNA polymerase II transcription units [that is, protein encoding open reading frames (ORFs)], many are derived from regions of DNA thought to be non-genic. Parallel chromatin immunoprecipitation studies of template-bound RNA polymerase II have shown that it is indeed resident on those regions found to be transcribed, both ORF and non-ORF. However, the strandedness of these pervasive transcripts has never been measured on a genome-wide basis. Four recent reports have addressed this question and, in the process, have made the startling discovery that many loci of mRNA sense gene transcription are associated with very active antisense or divergent transcription that begins at mapped transcription start sites and proceeds in an upstream direction.

From transcription start site to cell biology

The regulation of transcription is a complex process. Recent novel insights concerning the in vivo regulation and expression of protein-coding and non-coding RNAs have added previously unimagined levels of complexity to these processes.

The regulated retrotransposon transcriptome of mammalian cells

Although repetitive elements pervade mammalian genomes, their overall contribution to transcriptional activity is poorly defined. Here, as part of the FANTOM4 project, we report that 6-30% of cap-selected mouse and human RNA transcripts initiate within repetitive elements. Analysis of approximately 250,000 retrotransposon-derived transcription start sites shows that the associated transcripts are generally tissue specific, coincide with gene-dense regions and form pronounced clusters when aligned to full-length retrotransposon sequences. Retrotransposons located immediately 5' of protein-coding loci frequently function as alternative promoters and/or express noncoding RNAs. More than a quarter of RefSeqs possess a retrotransposon in their 3' UTR, with strong evidence for the reduced expression of these transcripts relative to retrotransposon-free transcripts. Finally, a genome-wide screen identifies 23,000 candidate regulatory regions derived from retrotransposons, in addition to more than 2,000 examples of bidirectional transcription. We conclude that retrotransposon transcription has a key influence upon the transcriptional output of the mammalian genome.