Regulated post-transcriptional RNA cleavage diversifies the eukaryotic transcriptome

Core promoterome of barley embryo

Precise localization and dissection of gene promoters are key to understanding transcriptional gene regulation and to successful bioengineering applications. The core RNA polymerase II initiation machinery is highly conserved among eukaryotes, leading to a general expectation of equivalent underlying mechanisms. Still, less is known about promoters in the plant kingdom. In this study, we employed cap analysis of gene expression (CAGE) at three embryonic developmental stages in barley to accurately map, annotate, and quantify transcription initiation events. Unsupervised discovery of de novo sequence clusters grouped promoters based on characteristic initiator and position-specific core-promoter motifs. This grouping was complemented by the annotation of transcription factor binding site (TFBS) motifs. Integration with genome-wide epigenomic data sets and gene ontology (GO) enrichment analysis further delineated the chromatin environments and functional roles of genes associated with distinct promoter categories. The TATA-box presence governs all features explored, supporting the general model of two separate genomic regulatory environments. We describe the extent and implications of alternative transcription initiation events, including those that are specific to developmental stages, which can affect the protein sequence or the presence of regions that regulate translation. The generated promoterome dataset provides a valuable genomic resource for enhancing the functional annotation of the barley genome. It also offers insights into the transcriptional regulation of individual genes and presents opportunities for the informed manipulation of promoter architecture, with the aim of enhancing traits of agronomic importance.

Herpes simplex virus 1 immediate early transcription initiation, pause-release, elongation, and termination in the presence and absence of ICP4

IMPORTANCE

Infection with herpes simplex virus 1 (HSV-1) leads to lifelong infection due to the virus's remarkable ability to control transcription of its own genome, resulting in two transcriptional programs: lytic (highly active) and latent (restricted). The lytic program requires immediate early (IE) proteins to first repress transcription of late viral genes, which then undergo sequential de-repression, leading to a specific sequence of gene expression. Here, we show that the IE ICP4 functions to regulate the cascade by limiting RNA polymerase initiation at immediate early times. However, late viral genes that initiate too early in the absence of ICP4 do not yield mRNA as transcription stalls within gene bodies. It follows that other regulatory steps intercede to prevent elongation of genes at the incorrect time, demonstrating the precise control HSV-1 exerts over its own transcription.

Challenges in identifying mRNA transcript starts and ends from long-read sequencing data

Long-read sequencing (LRS) technologies have the potential to revolutionize scientific discoveries in RNA biology, especially by enabling the comprehensive identification and quantification of full length mRNA isoforms. However, inherently high error rates make the analysis of long-read sequencing data challenging. While these error rates have been characterized for sequence and splice site identification, it is still unclear how accurately LRS reads represent transcript start and end sites. Here, we systematically assess the variability and accuracy of mRNA terminal ends identified by LRS reads across multiple sequencing platforms. We find substantial inconsistencies in both the start and end coordinates of LRS reads spanning a gene, such that LRS reads often fail to accurately recapitulate annotated or empirically derived terminal ends of mRNA molecules. To address this challenge, we introduce an approach to condition reads based on empirically derived terminal ends and identified a subset of reads that are more likely to represent full-length transcripts. Our approach can improve transcriptome analyses by enhancing the fidelity of transcript terminal end identification, but may result in lower power to quantify genes or discover novel isoforms. Thus, it is necessary to be cautious when selecting sequencing approaches and/or interpreting data from long-read RNA sequencing.

Evidence for Existence of Multiple Functional Human Small RNAs Derived from Transcripts of Protein-Coding Genes

The human genome encodes a multitude of different noncoding transcripts that have been traditionally separated on the basis of their lengths into long (>200 nt) or small (<200 nt) noncoding RNAs. The functions, mechanisms of action, and biological relevance of the vast majority of both long and short noncoding transcripts remain unknown. However, according to the functional understanding of the known classes of long and small noncoding RNAs (sncRNAs) that have been shown to play crucial roles in multiple biological processes, it is generally assumed that many unannotated long and small transcripts participate in important cellular functions as well. Nevertheless, direct evidence of functionality is lacking for most noncoding transcripts, especially for sncRNAs that are often dismissed as stable degradation products of longer RNAs. Here, we developed a high-throughput assay to test the functionality of sncRNAs by overexpressing them in human cells. Surprisingly, we found that a significant fraction (>40%) of unannotated sncRNAs appear to have biological relevance. Furthermore, contrary to the expectation, the potentially functional transcripts are not highly abundant and can be derived from protein-coding mRNAs. These results strongly suggest that the small noncoding transcriptome can harbor multiple functional transcripts that warrant future studies.

Noncanonical metabolite RNA caps: Classification, quantification, (de)capping, and function

The 5' cap of eukaryotic mRNA is a hallmark for cellular functions from mRNA stability to translation. However, the discovery of novel 5'-terminal RNA caps derived from cellular metabolites has challenged this long-standing singularity in both eukaryotes and prokaryotes. Reminiscent of the 7-methylguanosine (m7G) cap structure, these noncanonical caps originate from abundant coenzymes such as NAD, FAD, or CoA and from metabolites like dinucleoside polyphosphates (NpnN). As of now, the significance of noncanonical RNA caps is elusive: they differ for individual transcripts, occur in distinct types of RNA, and change in response to environmental stimuli. A thorough comparison of their prevalence, quantity, and characteristics is indispensable to define the distinct classes of metabolite-capped RNAs. This is achieved by a structured analysis of all present studies covering functional, quantitative, and sequencing data which help to uncover their biological impact. The biosynthetic strategies of noncanonical RNA capping and the elaborate decapping machinery reveal the regulation and turnover of metabolite-capped RNAs. With noncanonical capping being a universal and ancient phenomenon, organisms have developed diverging strategies to adapt metabolite-derived caps to their metabolic needs, but ultimately to establish noncanonical RNA caps as another intriguing layer of RNA regulation. This article is categorized under: RNA Processing &gt; Capping and 5' End Modifications RNA Turnover and Surveillance &gt; Turnover/Surveillance Mechanisms RNA Turnover and Surveillance &gt; Regulation of RNA Stability.

Alternative cleavage and polyadenylation generates downstream uncapped RNA isoforms with translation potential

The use of alternative promoters, splicing, and cleavage and polyadenylation (APA) generates mRNA isoforms that expand the diversity and complexity of the transcriptome. Here, we uncovered thousands of previously undescribed 5' uncapped and polyadenylated transcripts (5' UPTs). We show that these transcripts resist exonucleases due to a highly structured RNA and N6-methyladenosine modification at their 5' termini. 5' UPTs appear downstream of APA sites within their host genes and are induced upon APA activation. Strong enrichment in polysomal RNA fractions indicates 5' UPT translational potential. Indeed, APA promotes downstream translation initiation, non-canonical protein output, and consistent changes to peptide presentation at the cell surface. Lastly, we demonstrate the biological importance of 5' UPTs using Bcl2, a prominent anti-apoptotic gene whose entire coding sequence is a 5' UPT generated from 5' UTR-embedded APA sites. Thus, APA is not only accountable for terminating transcripts, but also for generating downstream uncapped RNAs with translation potential and biological impact.

Sequence features around cleavage sites are highly conserved among different species and a critical determinant for RNA cleavage position across eukaryotes

RNA degradation is one of the critical steps for control of gene expression, and endonucleolytic cleavage-dependent RNA degradation is conserved among eukaryotes. Some cleavage sites are secondarily capped in the cytoplasm and identified using the Cap analysis of gene expression (CAGE) method. Although uncapped cleavage sites are widespread in eukaryotes, comparatively little information has been obtained about these sites using CAGE-based degradome analysis. Previously, we developed the truncated RNA-end sequencing (TREseq) method in plant species and used it to acquire comprehensive information about uncapped cleavage sites; we observed G-rich sequences near cleavage sites. However, it remains unclear whether this finding is general to other eukaryotes. In this study, we conducted TREseq analyses in fruit flies (Drosophila melanogaster) and budding yeast (Saccharomyces cerevisiae). The results revealed specific sequence features related to RNA cleavage in D. melanogaster and S. cerevisiae that were similar to sequence patterns in Arabidopsis thaliana. Although previous studies suggest that ribosome movements are important for determining cleavage position, feature selection using a random forest classifier showed that sequences around cleavage sites were major determinant for cleaved or uncleaved sites. Together, our results suggest that sequence features around cleavage sites are critical for determining cleavage position, and that sequence-specific endonucleolytic cleavage-dependent RNA degradation is highly conserved across eukaryotes.

Methods for detecting RNA degradation intermediates in plants

RNA degradation is an important process for controlling gene expression and is mediated by decapping / deadenylation-dependent or endonucleolytic cleavage-dependent RNA degradation mechanisms. High-throughput sequencing of RNA degradation intermediates was initially developed in Arabidopsis thaliana and similar RNA degradome sequencing methods were conducted in other eukaryotes. However, interpreting results obtained by these sequencing methods is fragmented, and an overview is needed. Here we review the findings and limitations of these sequencing methods and discuss the missing experiments needed to understand RNA degradation intermediates accurately. This review provides direction for future research on RNA degradation and is a reference for RNA degradome studies in other species.

Towards unveiling the nature of short SERPINA1 transcripts: Avoiding the main ORF control to translate alpha1-antitrypsin C-terminal peptides

Alternative ORFs in-frame with the known genes are challenging to reveal. Yet they may contribute significantly to proteome diversity. Here we focused on the individual expression of the SERPINA1 gene exon 5 leading to direct translation of alpha1-antitrypsin (AAT) C-terminal peptides. The discovery of alternative ways for their production may expand the current understanding of the serpin gene's functioning. We detected short transcripts expressed primarily in hepatocytes. We identified four variants of hepatocyte-specific SERPINA1 short transcripts and individually probed their potential to be translated in living cells. The long mRNA gave the full-length AAT-eGFP fusion, while in case of short transcripts we deduced four active SERPINA1 in-frame alternative ORFs encoding 10, 21, 153 and 169 amino acids AAT C-terminal oligo- and polypeptides. Unlike secretory AAT-eGFP fusion exhibiting classical AAT behavior, truncated AAT-fusions differ by intracellular retention and nuclear enrichment. Immunofluorescence on the endogenous AAT C-terminal epitope showed its accumulation in the cell nucleoli, indicating that short transcripts may be translated in vivo. FANTOM5 CAGE data on SERPINA1 suggest that short transcripts originate from the post-transcriptional cleavage of the spliced mRNA, initiated mainly from the hepatocyte-specific promoter. CONCLUSION: Short SERPINA1 transcripts may represent a source for the direct synthesis of AAT C-terminal peptides with properties uncommon to AAT.

To cap it all off, again: dynamic capping and recapping of coding and non-coding RNAs to control transcript fate and biological activity

The addition of the methyl-7-guanosine (m⁷G) “cap” on the 5' ends of coding and some non-coding RNAs is essential for their protein coding capacity and biochemical activity, respectively. It was previously considered that capping was a constitutive process that generates a complete cap on all transcripts at steady-state. However, development of new methodologies demonstrated that steady-state capping is a dynamic and regulatable feature of many coding and non-coding RNAs. Indeed, capping status of specific RNAs can flux during differentiation and development, thereby impacting on their protein-coding capacity and activity. Moreover, in some primary cancer specimens, capping can be elevated for transcripts encoding proteins involved in proliferation and survival corresponding to their increased protein levels. Overexpression of one of the capping enzymes (RNMT), the transcription factor MYC or the eukaryotic translation initiation factor eIF4E all led to increased levels of steady-state capping of selected transcripts. Additionally, transcripts can be decapped and recapped, allowing these to be sequestered until needed. This review provides a summary of the major advances in enzymatic and affinity-based approaches to quantify m⁷G capping. Further, we summarize the evidence for regulation of capping. Capping has emerged as a significant regulatory step in RNA metabolism which is poised to impact a myriad of biological processes.

Long non-coding RNA linc00665 inhibits CDKN1C expression by binding to EZH2 and affects cisplatin sensitivity of NSCLC cells

Long non-coding RNAs (lncRNAs) can play significant regulatory roles in cells that affect the development and acquired drug resistance of lung cancer. Herein, we report that lncRNA linc00665 is significantly upregulated in non-small cell lung cancer (NSCLC) tissues compared with adjacent normal tissues. linc00665 affects the sensitivity of NSCLC cells to the chemotherapy drug cisplatin (DDP), making it a potential target for the treatment of NSCLC. Functional experiments showed that linc00665 enhanced the proliferation and migration of NSCLC cells in vivo and in vitro, and knocking down linc00665 could enhance the drug sensitivity of NSCLC cells to DDP. Further work revealed that linc00665 could recruit enhancer of zeste homolog 2 (EZH2) to the promoter region of cyclin-dependent kinase inhibitor 1C (CDKN1C) to inhibit its transcription and thus carry out its tumorigenic role. In conclusion, our study elucidated the carcinogenic role of the linc00665-EZH2-CDKN1C axis in NSCLC tumors and its ability to influence the sensitivity of these tumors to DDP. These results suggest that linc00665 may be a potential diagnostic marker and therapeutic target in NSCLC, and they also provide a new direction for the development of clinical reversal methods for acquired drug resistance in patients with NSCLC.

Analysis of splice variants of the human protein disulfide isomerase (P4HB) gene

BACKGROUND

Protein Disulfide Isomerases are thiol oxidoreductase chaperones from thioredoxin superfamily with crucial roles in endoplasmic reticulum proteostasis, implicated in many diseases. The family prototype PDIA1 is also involved in vascular redox cell signaling. PDIA1 is coded by the P4HB gene. While forced changes in P4HB gene expression promote physiological effects, little is known about endogenous P4HB gene regulation and, in particular, gene modulation by alternative splicing. This study addressed the P4HB splice variant landscape.

RESULTS

Ten protein coding sequences (Ensembl) of the P4HB gene originating from alternative splicing were characterized. Structural features suggest that except for P4HB-021, other splice variants are unlikely to exert thiol isomerase activity at the endoplasmic reticulum. Extensive analyses using FANTOM5, ENCODE Consortium and GTEx project databases as RNA-seq data sources were performed. These indicated widespread expression but significant variability in the degree of isoform expression among distinct tissues and even among distinct locations of the same cell, e.g., vascular smooth muscle cells from different origins. P4HB-02, P4HB-027 and P4HB-021 were relatively more expressed across each database, the latter particularly in vascular smooth muscle. Expression of such variants was validated by qRT-PCR in some cell types. The most consistently expressed splice variant was P4HB-021 in human mammary artery vascular smooth muscle which, together with canonical P4HB gene, had its expression enhanced by serum starvation.

CONCLUSIONS

Our study details the splice variant landscape of the P4HB gene, indicating their potential role to diversify the functional reach of this crucial gene. P4HB-021 splice variant deserves further investigation in vascular smooth muscle cells.

Epigenetic regulation of spurious transcription initiation in Arabidopsis

In plants, epigenetic regulation is critical for silencing transposons and maintaining proper gene expression. However, its impact on the genome-wide transcription initiation landscape remains elusive. By conducting a genome-wide analysis of transcription start sites (TSSs) using cap analysis of gene expression (CAGE) sequencing, we show that thousands of TSSs are exclusively activated in various epigenetic mutants of Arabidopsis thaliana and referred to as cryptic TSSs. Many have not been identified in previous studies, of which up to 65% are contributed by transposons. They possess similar genetic features to regular TSSs and their activation is strongly associated with the ectopic recruitment of RNAPII machinery. The activation of cryptic TSSs significantly alters transcription of nearby TSSs, including those of genes important for development and stress responses. Our study, therefore, sheds light on the role of epigenetic regulation in maintaining proper gene functions in plants by suppressing transcription from cryptic TSSs.

Epigenetic regulation can silence transposons and maintain gene expression. Here the authors survey Arabidopsis mutants defective in epigenetic regulation and show ectopic activation of thousands of cryptic TSSs and altered expression of nearby genes demonstrating the importance of suppressing spurious transcription.

AP-TSS: A New Method for the Analysis of RNA Expression from Particular and Challenging Transcription Start Sites

Alternative promoter usage involved in the regulation of transcription, splicing, and translation contributes to proteome diversity and is involved in a large number of diseases, in particular, cancer. Epigenetic mechanisms and cis regulatory elements are involved in alternative promoter activity. Multiple transcript isoforms can be produced from a gene, due to the initiation of transcription at different transcription start sites (TSS). These transcripts may not have regions that allow discrimination during RT-qPCR, making quantification technically challenging. This study presents a general method for the relative quantification of a transcript synthesized from a particular TSS that we called AP-TSS (analysis of particular TSS). AP-TSS is based on the specific elongation of the cDNA of interest, followed by its quantification by qPCR. As proof of principle, AP-TSS was applied to two non-coding RNA: telomeric repeat-containing RNAs (TERRA) from a particular subtelomeric TSS, and Alu transcripts. The treatment of cells with a DNA methylation inhibitor was associated with a global increase of the total TERRA level, but the TERRA expression from the TSS of interest did not change in HT1080 cells, and only modestly increased in HeLa cells. This result suggests that TERRA upregulation induced by global demethylation of the genome is mainly due to activation from sites other than this particular TSS. For Alu RNA, the signal obtained by AP-TSS is specific for the RNA Polymerase III-dependent Alu transcript. In summary, our method provides a tool to study regulation of gene expression from a given transcription start site, in different conditions that could be applied to many genes. In particular, AP-TSS can be used to investigate the epigenetic regulation of alternative TSS usage that is of importance for the development of epigenetic-targeted therapies.

Networks of mRNA Processing and Alternative Splicing Regulation in Health and Disease

mRNA processing events introduce an intricate layer of complexity into gene expression processes, supporting a tremendous level of diversification of the genome's coding and regulatory potential, particularly in vertebrate species. The recent development of massive parallel sequencing methods and their adaptation to the identification and quantification of different RNA species and the dynamics of mRNA metabolism and processing has generated an unprecedented view over the regulatory networks that are established at this level, which contribute to sustain developmental, tissue specific or disease specific gene expression programs. In this chapter, we provide an overview of the recent evolution of transcriptome profiling methods and the surprising insights that have emerged in recent years regarding distinct mRNA processing events - from the 5' end to the 3' end of the molecule.

RNA-binding proteins and heat-shock protein 90 are constituents of the cytoplasmic capping enzyme interactome

The N7-methylguanosine cap is added in the nucleus early in gene transcription and is a defining feature of eukaryotic mRNAs. Mammalian cells also possess cytoplasmic machinery for restoring the cap at uncapped or partially degraded RNA 5' ends. Central to both pathways is capping enzyme (CE) (RNA guanylyltransferase and 5'-phosphatase (RNGTT)), a bifunctional, nuclear and cytoplasmic enzyme. CE is recruited to the cytoplasmic capping complex by binding of a C-terminal proline-rich sequence to the third Src homology 3 (SH3) domain of NCK adapter protein 1 (NCK1). To gain broader insight into the cellular context of cytoplasmic recapping, here we identified the protein interactome of cytoplasmic CE in human U2OS cells through two complementary approaches: chemical cross-linking and recovery with cytoplasmic CE and protein screening with proximity-dependent biotin identification (BioID). This strategy unexpectedly identified 66 proteins, 52 of which are RNA-binding proteins. We found that CE interacts with several of these proteins independently of RNA, mediated by sequences within its N-terminal triphosphatase domain, and we present a model describing how CE-binding proteins may function in defining recapping targets. This analysis also revealed that CE is a client protein of heat shock protein 90 (HSP90). Nuclear and cytoplasmic CEs were exquisitely sensitive to inhibition of HSP90, with both forms declining significantly following treatment with each of several HSP90 inhibitors. Importantly, steady-state levels of capped mRNAs decreased in cells treated with the HSP90 inhibitor geldanamycin, raising the possibility that the cytotoxic effect of these drugs may partially be due to a general reduction in translatable mRNAs.

Diverse Long RNAs Are Differentially Sorted into Extracellular Vesicles Secreted by Colorectal Cancer Cells

The regulation and functional roles of secreted coding and long noncoding RNAs (lncRNAs; >200 nt) are largely unknown. We previously showed that mutant KRAS colorectal cancer (CRC) cells release extracellular vesicles (EVs) containing distinct proteomes, microRNAs (miRNAs), and circular RNAs. Here, we comprehensively identify diverse classes of CRC extracellular long RNAs secreted in EVs and demonstrate differential export of specific RNAs. Distinct noncoding RNAs, including antisense transcripts and transcripts derived from pseudogenes, are enriched in EVs compared to cellular profiles. We detected strong enrichment of Rab13 in mutant KRAS EVs and demonstrate functional delivery of Rab13 mRNA to recipient cells. To assay functional transfer of lncRNAs, we implemented a CRISPR/Cas9-based RNA-tracking system to monitor delivery to recipient cells. We show that gRNAs containing export signals from secreted RNAs can be transferred from donor to recipient cells. Our data support the existence of cellular mechanisms to selectively export diverse classes of RNA.

A recap of RNA recapping

The N7-methylguanosine cap is a hallmark of the 5' end of eukaryotic mRNAs and is required for gene expression. Loss of the cap was believed to lead irreversibly to decay. However, nearly a decade ago, it was discovered that mammalian cells contain enzymes in the cytoplasm that are capable of restoring caps onto uncapped RNAs. In this review, we summarize recent advances in our understanding of cytoplasmic RNA recapping and discuss the biochemistry of this process and its impact on regulating and diversifying the transcriptome. Although most studies focus on mammalian RNA recapping, we also highlight new observations for recapping in disparate eukaryotic organisms, with the trypanosome recapping system appearing to be a fascinating example of convergent evolution. We conclude with emerging insights into the biological significance of RNA recapping and prospects for the future of this evolving area of study. This article is categorized under: RNA Processing > RNA Editing and Modification Translation > Translation Regulation RNA Processing > Capping and 5' End Modifications RNA Turnover and Surveillance > Regulation of RNA Stability.

The State of Long Non-Coding RNA Biology

Transcriptomic studies have demonstrated that the vast majority of the genomes of mammals and other complex organisms is expressed in highly dynamic and cell-specific patterns to produce large numbers of intergenic, antisense and intronic long non-protein-coding RNAs (lncRNAs). Despite well characterized examples, their scaling with developmental complexity, and many demonstrations of their association with cellular processes, development and diseases, lncRNAs are still to be widely accepted as major players in gene regulation. This may reflect an underappreciation of the extent and precision of the epigenetic control of differentiation and development, where lncRNAs appear to have a central role, likely as organizational and guide molecules: most lncRNAs are nuclear-localized and chromatin-associated, with some involved in the formation of specialized subcellular domains. I suggest that a reassessment of the conceptual framework of genetic information and gene expression in the 4-dimensional ontogeny of spatially organized multicellular organisms is required. Together with this and further studies on their biology, the key challenges now are to determine the structure–function relationships of lncRNAs, which may be aided by emerging evidence of their modular structure, the role of RNA editing and modification in enabling epigenetic plasticity, and the role of RNA signaling in transgenerational inheritance of experience.

GALNT5 uaRNA promotes gastric cancer progression through its interaction with HSP90

Recently, long noncoding RNAs (lncRNAs) have been reported to play a pivotal role in the occurrence and progression of cancer because of their unique characteristics and have therefore become an active area of cancer research. The object of this study was to screen lncRNAs that are dysregulated in gastric cancer and to investigate their potential functions. Global expression of lncRNAs in gastric cancer and adjacent normal tissues of patients was profiled using a microarray assay. We identified an lncRNA (GALNT5 uaRNA, UTR-associated RNA) that is derived from the 3'-UTR of GALNT5. This lncRNA was transcribed independently of the coding region of GALNT5 and was determined to be markedly upregulated in human gastric carcinoma relative to their corresponding normal gastric tissues by quantitative RT-PCR (qRT-PCR) analysis of tissues from 122 gastric carcinoma patients. The expression of GALNT5 uaRNA was significantly correlated with the TNM stage and with lymph node metastasis. Further results demonstrated that GALNT5 uaRNA facilitated the proliferation and migration of gastric cancer cells in vitro and promoted tumor growth in a mouse model of human gastric cancer. Our results also indicated that GALNT5 uaRNA might function in gastric cancer by binding with HSP90. Further studies indicated that the 5'-end stem-loop motifs of GALNT5 uaRNA promoted the binding of HSP90 and its client proteins, and thus inhibited ubiquitination of the clients. These results expanded our understanding of GALNT5 uaRNA as a new avenue for therapeutic intervention against gastric cancer progression.

The "dark matter" of DNA and the regulation of bone metabolism: The role of non-coding RNAs

Epigenetics, present a new discipline that attempts to explain significant differences in phenotypes among patients with the same disease. In contrast to the other epigenetic mechanisms that modulate gene transcription, non-coding RNAs act at the post-transcriptional level. They directly modulate the gene expression of mRNA genes leading to mRNA target cleavage and degradation and translation repression. Bioinformatic predictions indicate that non coding RNAs may be involved in the regulation of 60% of the coding genes and each non-coding RNA can have multiple target genes, and each gene may be regulated by more than one non-coding RNAs. In the last decade several studies have shown a significant role of non-coding RNAs in the regulation of bone metabolism and function of bone cells opening a new era in the understanding of bone biology in health and disease.

Non-coding RNA networks in cancer

Thousands of unique non-coding RNA (ncRNA) sequences exist within cells. Work from the past decade has altered our perception of ncRNAs from 'junk' transcriptional products to functional regulatory molecules that mediate cellular processes including chromatin remodelling, transcription, post-transcriptional modifications and signal transduction. The networks in which ncRNAs engage can influence numerous molecular targets to drive specific cell biological responses and fates. Consequently, ncRNAs act as key regulators of physiological programmes in developmental and disease contexts. Particularly relevant in cancer, ncRNAs have been identified as oncogenic drivers and tumour suppressors in every major cancer type. Thus, a deeper understanding of the complex networks of interactions that ncRNAs coordinate would provide a unique opportunity to design better therapeutic interventions.

Alternative mRNA processing sites decrease genetic variability while increasing functional diversity

Recent large-scale RNA sequencing efforts have revealed the extensive diversity of mRNA molecules produced from most eukaryotic coding genes, which arises from the usage of alternative, cryptic or non-canonical splicing and intronic polyadenylation sites. The prevailing view regarding the tremendous diversity of coding gene transcripts is that mRNA processing is a flexible and more-or-less noisy process leading to a diversity of proteins on which natural selection can act depending on protein-mediated cellular functions. However, this concept raises two main questions. First, do alternative mRNA processing pathways have a role other than generating mRNA and protein diversity? Second, is the cellular function of mRNA variants restricted to the biogenesis of functional protein isoforms? Here, I propose that the co-transcriptional use of alternative mRNA processing sites allows first, the resolution of co-transcriptional biophysical constraints that may otherwise result in DNA instability, and second, increases the diversity of cellular functions of mRNAs in a manner that is not restricted to protein synthesis.

Short intron-derived ncRNAs

Introns represent almost half of the human genome, although they are eliminated from transcripts through RNA splicing. Yet, different classes of non-canonical miRNAs have been proposed to originate directly from intron splicing. Here, we considered the alternative splicing of introns as an interesting source of miRNAs, compatible with a developmental switch. We report computational prediction of new Short Intron-Derived ncRNAs (SID), defined as precursors of smaller ncRNAs like miRNAs and snoRNAs produced directly by splicing, and tested their dependence on each key factor in canonical or alternative miRNAs biogenesis (Drosha, DGCR8, DBR1, snRNP70, U2AF65, PRP8, Dicer, Ago2). We found that about half of predicted SID rely on debranching of the excised intron-lariat by the enzyme DBR1, as proposed for mirtrons. However, we identified new classes of SID for which miRNAs biogenesis may rely on intermingling between canonical and alternative pathways. We validated selected SID as putative miRNAs precursors and identified new endogenous miRNAs produced by non-canonical pathways, including one hosted in the first intron of SRA (Steroid Receptor RNA activator). Consistent with increased SRA intron retention during myogenic differentiation, release of SRA intron and its associated mature miRNA decreased in cells from healthy subjects but not from myotonic dystrophy patients with splicing defects.

The 5'UTR in human adenoviruses: leader diversity in late gene expression

Human adenoviruses (HAdVs) shut down host cellular cap-dependent mRNA translation while initiating the translation of viral late mRNAs in a cap-independent manner. HAdV 5′ untranslated regions (5′UTRs) are crucial for cap-independent initiation, and influence mRNA localization and stability. However, HAdV translational regulation remains relatively uncharacterized. The HAdV tripartite leader (TPL), composed of three introns (TPL 1–3), is critical to the translation of HAdV late mRNA. Herein, we annotated and analyzed 72 HAdV genotypes for the HAdV TPL and another previously described leader, the i-leader. Using HAdV species D, type 37 (HAdV-D37), we show by reverse transcription PCR and Sanger sequencing that mRNAs of the HAdV-D37 E3 transcription unit are spliced to the TPL. We also identified a polycistronic mRNA for RID-α and RID-β. Analysis of the i-leader revealed a potential open reading frame within the leader sequence and the termination of this potential protein in TPL3. A potential new leader embedded within the E3 region was also detected and tentatively named the j-leader. These results suggest an underappreciated complexity of post-transcriptional regulation, and the importance of HAdV 5′UTRs for precisely coordinated viral protein expression along the path from genotype to phenotype.

Dynamic landscape of alternative polyadenylation during retinal development

The development of the central nervous system (CNS) is a complex process that must be exquisitely controlled at multiple levels to ensure the production of appropriate types and quantity of neurons. RNA alternative polyadenylation (APA) contributes to transcriptome diversity and gene regulation, and has recently been shown to be widespread in the CNS. However, the previous studies have been primarily focused on the tissue specificity of APA and developmental APA change of whole model organisms; a systematic survey of APA usage is lacking during CNS development. Here, we conducted global analysis of APA during mouse retinal development, and identified stage-specific polyadenylation (pA) sites that are enriched for genes critical for retinal development and visual perception. Moreover, we demonstrated 3'UTR (untranslated region) lengthening and increased usage of intronic pA sites over development that would result in gaining many different RBP (RNA-binding protein) and miRNA target sites. Furthermore, we showed that a considerable number of polyadenylated lncRNAs are co-expressed with protein-coding genes involved in retinal development and functions. Together, our data indicate that APA is highly and dynamically regulated during retinal development and maturation, suggesting that APA may serve as a crucial mechanism of gene regulation underlying the delicate process of CNS development.

Putative RNA-directed adaptive mutations in cancer evolution

Understanding the molecular mechanisms behind the capacity of cancer cells to adapt to the tumor microenvironment and to anticancer therapies is a major challenge. In this context, cancer is believed to be an evolutionary process where random mutations and the selection process shape the mutational pattern and phenotype of cancer cells. This article challenges the notion of randomness of some cancer-associated mutations by describing molecular mechanisms involving stress-mediated biogenesis of mRNA-derived small RNAs able to target and increase the local mutation rate of the genomic loci they originate from. It is proposed that the probability of some mutations at specific loci could be increased in a stress-specific and RNA-depending manner. This would increase the probability of generating mutations that could alleviate stress situations, such as those triggered by anticancer drugs. Such a mechanism is made possible because tumor- and anticancer drug-associated stress situations trigger both cellular reprogramming and inflammation, which leads cancer cells to express molecular tools allowing them to “attack” and mutate their own genome in an RNA-directed manner.

Global analysis of biogenesis, stability and sub-cellular localization of lncRNAs mapping to intragenic regions of the human genome

Long noncoding RNAs (lncRNAs) that map to intragenic regions of the human genome with the same (intronic lncRNAs) or opposite orientation (antisense lncRNAs) relative to protein-coding mRNAs have been largely dismissed from biochemical and functional characterization due to the belief that they are mRNA precursors, byproducts of RNA splicing or simply transcriptional noise. In this work, we used a custom microarray to investigate aspects of the biogenesis, processing, stability, evolutionary conservation, and cellular localization of ∼ 6,000 intronic lncRNAs and ∼ 10,000 antisense lncRNAs. Most intronic (2,903 of 3,427, 85%) and antisense lncRNAs (4,945 of 5,214, 95%) expressed in HeLa cells showed evidence of 5' cap modification, compatible with their transcription by RNAP II. Antisense lncRNAs (median t1/2 = 3.9 h) were significantly (p < 0.0001) more stable than mRNAs (median t1/2 = 3.2 h), whereas intronic lncRNAs (median t1/2 = 2.1 h) comprised a more heterogeneous class that included both stable (t1/2 > 3 h) and unstable (t1/2 < 1 h) transcripts. Intragenic lncRNAs display evidence of evolutionary conservation, have little/no coding potential and were ubiquitously detected in the cytoplasm. Notably, a fraction of the intronic and antisense lncRNAs (13 and 15%, respectively) were expressed from loci at which the corresponding host mRNA was not detected. The abundances of a subset of intronic/antisense lncRNAs were correlated (r ≥ |0.8|) with those of genes encoding proteins involved in cell division and DNA replication. Taken together, the findings of this study contribute novel biochemical and genomic information regarding intronic and antisense lncRNAs, supporting the notion that these classes include independently transcribed RNAs with potentials for exerting regulatory functions in the cell.

Functional annotation of the vlinc class of non-coding RNAs using systems biology approach

Functionality of the non-coding transcripts encoded by the human genome is the coveted goal of the modern genomics research. While commonly relied on the classical methods of forward genetics, integration of different genomics datasets in a global Systems Biology fashion presents a more productive avenue of achieving this very complex aim. Here we report application of a Systems Biology-based approach to dissect functionality of a newly identified vast class of very long intergenic non-coding (vlinc) RNAs. Using highly quantitative FANTOM5 CAGE dataset, we show that these RNAs could be grouped into 1542 novel human genes based on analysis of insulators that we show here indeed function as genomic barrier elements. We show that vlinc RNAs genes likely function in cisto activate nearby genes. This effect while most pronounced in closely spaced vlinc RNA-gene pairs can be detected over relatively large genomic distances. Furthermore, we identified 101 vlinc RNA genes likely involved in early embryogenesis based on patterns of their expression and regulation. We also found another 109 such genes potentially involved in cellular functions also happening at early stages of development such as proliferation, migration and apoptosis. Overall, we show that Systems Biology-based methods have great promise for functional annotation of non-coding RNAs.

Evolution of the unspliced transcriptome

BACKGROUND

Despite their abundance, unspliced EST data have received little attention as a source of information on non-coding RNAs. Very little is know, therefore, about the genomic distribution of unspliced non-coding transcripts and their relationship with the much better studied regularly spliced products. In particular, their evolution has remained virtually unstudied.

RESULTS

We systematically study the evidence on unspliced transcripts available in EST annotation tracks for human and mouse, comprising 104,980 and 66,109 unspliced EST clusters, respectively. Roughly one third of these are located totally inside introns of known genes (TINs) and another third overlaps exonic regions (PINs). Eleven percent are "intergenic", far away from any annotated gene. Direct evidence for the independent transcription of many PINs and TINs is obtained from CAGE tag and chromatin data. We predict more than 2000 3'UTR-associated RNA candidates for each human and mouse. Fifteen to twenty percent of the unspliced EST cluster are conserved between human and mouse. With the exception of TINs, the sequences of unspliced EST clusters evolve significantly slower than genomic background. Furthermore, like spliced lincRNAs, they show highly tissue-specific expression patterns.

CONCLUSIONS

Unspliced long non-coding RNAs are an important, rapidly evolving, component of mammalian transcriptomes. Their analysis is complicated by their preferential association with complex transcribed loci that usually also harbor a plethora of spliced transcripts. Unspliced EST data, although typically disregarded in transcriptome analysis, can be used to gain insights into this rarely investigated transcriptome component. The frequently postulated connection between lack of splicing and nuclear retention and the surprising overlap of chromatin-associated transcripts suggests that this class of transcripts might be involved in chromatin organization and possibly other mechanisms of epigenetic control.

mRIN for direct assessment of genome-wide and gene-specific mRNA integrity from large-scale RNA-sequencing data

The volume of RNA-Seq data sets in public repositories has been expanding exponentially, providing unprecedented opportunities to study gene expression regulation. Because degraded RNA samples, such as those collected from post-mortem tissues, can result in distinct expression profiles with potential biases, a particularly important step in mining these data is quality control. Here we develop a method named mRIN to directly assess mRNA integrity from RNA-Seq data at the sample and individual gene level. We systematically analyse large-scale RNA-Seq data sets of the human brain transcriptome generated by different consortia. Our analysis demonstrates that 3′ bias resulting from partial RNA fragmentation in post-mortem tissues has a marked impact on global expression profiles, and that mRIN effectively identifies samples with different levels of mRNA degradation. Unexpectedly, this process has a reproducible and gene-specific component, and transcripts with different stabilities are associated with distinct functions and structural features reminiscent of mRNA decay in living cells.

With the rapid increase in the volume of publically available RNA-seq data, quality control is an increasingly important consideration. Here Feng et al. develop mRIN, a method to directly assess mRNA integrity, and show that RNA degradation in post-mortem samples has a strong impact on global expression profiles.

Exosomes Secreted from Human Cancer Cell Lines Contain Inhibitors of Apoptosis (IAP)

Exosomes are endosomal-derived nanovesicles released by normal and tumor cells which have been shown to transfer functionally active protein, lipids, mRNAs and miRNAs between cells. Varying in molecular profiles, biological roles, functional roles and protein contents, exosomes have been described as "multi-purpose carriers" playing a role in supporting the survival and growth of tumor cells. The IAP Survivin has been found to be present in tumor exosomes. However, the existence of other IAPs in tumor exosomes is still unknown. Survivin, cIAP1, cIAP2 and XIAP mRNA and protein are differently expressed in a panel of tumor cell lines: DLCL2, HeLa, MCF-7, Panc-1, and PC3. Exosomes were isolated from conditioned media collected from the cells from which RNA and protein were extracted. Our results provide evidence that like Survivin, XIAP, cIAP1 and cIAP2 proteins are found in tumor exosomes. The mRNA expression, however, is differentially expressed across the tumor cell lines. The presence of these bioactive molecules in exosomes may not only serve as warning signals, but also play a role in providing protection to the cancer cells against changes that are constantly occurring in the tumor microenvironment.

MicroRNA-directed cleavage of targets: mechanism and experimental approaches

MicroRNAs (miRNAs) are a large family of post-transcriptional regulators, which are 21-24 nt in length and play a role in a wide variety of biological processes in eukaryotes. The past few years have seen rapid progress in our understanding of miRNA biogenesis and the mechanism of action, which commonly entails a combination of target degradation and translational repression. The target degradation mediated by Argonaute-catalyzed endonucleolytic cleavage exerts a significant repressive effect on target mRNA expression, particularly during rapid developmental transitions. This review outlines the current understanding of the mechanistic aspects of this important process and discusses several different experimental approaches to identify miRNA cleavage targets. [BMB Reports 2014; 47(8): 417-423]

The many faces of Dicer: the complexity of the mechanisms regulating Dicer gene expression and enzyme activities

There is increasing evidence indicating that the production of small regulatory RNAs is not the only process in which ribonuclease Dicer can participate. For example, it has been demonstrated that this enzyme is also involved in chromatin structure remodelling, inflammation and apoptotic DNA degradation. Moreover, it has become increasingly clear that cellular transcript and protein levels of Dicer must be strictly controlled because even small changes in their accumulation can initiate various pathological processes, including carcinogenesis. Accordingly, in recent years, a number of studies have been performed to identify the factors regulating Dicer gene expression and protein activity. As a result, a large amount of complex and often contradictory data has been generated. None of these data have been subjected to an exhaustive review or critical discussion. This review attempts to fill this gap by summarizing the current knowledge of factors that regulate Dicer gene transcription, primary transcript processing, mRNA translation and enzyme activity. Because of the high complexity of this topic, this review mainly concentrates on human Dicer. This review also focuses on an additional regulatory layer of Dicer activity involving the interactions of protein and RNA factors with Dicer substrates.

A 3'UTR-associated RNA, FLJ11812 maintains stemness of human embryonic stem cells by targeting miR-4459

The 3' untranslated region (UTR)-associated RNAs (uaRNAs) have important roles in various biological processes, especially in development. However, since they overlap with protein-coding mRNAs, uaRNAs are difficult to study by RNA interference techniques. We recently identified a chemical molecule, 3-benzyl-5-((2-nitrophenoxy) methyl)-dihydrofuran-2(3H)-one (3BDO), that could efficiently induce human embryonic stem cells (hESCs) differentiation, and meanwhile selectively and efficiently downregulate the uaRNA FLJ11812. By acting as a competing endogenous RNA, downregulated FLJ11812 by 3BDO further increased miR-4459 level in hESCs. miR-4459 could decrease the expression of its targets, CDC20B and ATG13, and thus altered stemness via cell cycle and autophagy. Our results revealed that FLJ11812 played a key role in maintenance of stemness of hESCs for the first time. The findings provide new clues and a powerful tool for investigating the action mechanism of FLJ11812 in early development.

Exosomes in human semen carry a distinctive repertoire of small non-coding RNAs with potential regulatory functions

Semen contains relatively ill-defined regulatory components that likely aid fertilization, but which could also interfere with defense against infection. Each ejaculate contains trillions of exosomes, membrane-enclosed subcellular microvesicles, which have immunosuppressive effects on cells important in the genital mucosa. Exosomes in general are believed to mediate inter-cellular communication, possibly by transferring small RNA molecules. We found that seminal exosome (SE) preparations contain a substantial amount of RNA from 20 to 100 nucleotides (nts) in length. We sequenced 20-40 and 40-100 nt fractions of SE RNA separately from six semen donors. We found various classes of small non-coding RNA, including microRNA (21.7% of the RNA in the 20-40 nt fraction) as well as abundant Y RNAs and tRNAs present in both fractions. Specific RNAs were consistently present in all donors. For example, 10 (of ∼2600 known) microRNAs constituted over 40% of mature microRNA in SE. Additionally, tRNA fragments were strongly enriched for 5'-ends of 18-19 or 30-34 nts in length; such tRNA fragments repress translation. Thus, SE could potentially deliver regulatory signals to the recipient mucosa via transfer of small RNA molecules.

The rise of regulatory RNA

Discoveries over the past decade portend a paradigm shift in molecular biology. Evidence suggests that RNA is not only functional as a messenger between DNA and protein but also involved in the regulation of genome organization and gene expression, which is increasingly elaborate in complex organisms. Regulatory RNA seems to operate at many levels; in particular, it plays an important part in the epigenetic processes that control differentiation and development. These discoveries suggest a central role for RNA in human evolution and ontogeny. Here, we review the emergence of the previously unsuspected world of regulatory RNA from a historical perspective.

An mRNA-derived noncoding RNA targets and regulates the ribosome

The structural and functional repertoire of small non-protein-coding RNAs (ncRNAs) is central for establishing gene regulation networks in cells and organisms. Here, we show that an mRNA-derived 18-nucleotide-long ncRNA is capable of downregulating translation in Saccharomyces cerevisiae by targeting the ribosome. This 18-mer ncRNA binds to polysomes upon salt stress and is crucial for efficient growth under hyperosmotic conditions. Although the 18-mer RNA originates from the TRM10 locus, which encodes a tRNA methyltransferase, genetic analyses revealed the 18-mer RNA nucleotide sequence, rather than the mRNA-encoded enzyme, as the translation regulator. Our data reveal the ribosome as a target for a small regulatory ncRNA and demonstrate the existence of a yet unkown mechanism of translation regulation. Ribosome-targeted small ncRNAs are found in all domains of life and represent a prevalent but so far largely unexplored class of regulatory molecules.

•

This study reveals the yeast ribosome as direct target for small regulatory ncRNAs

•

An 18-nt-long exon-derived RNA fragment from the TRM10 locus binds to ribosomes

•

This 18-mer ncRNA inhibits global protein biosynthesis in vivo and in vitro

•

This translation attenuation is crucial for adaption under hyperosmotic stress

Cell cycle, oncogenic and tumor suppressor pathways regulate numerous long and macro non-protein-coding RNAs

Background

The genome is pervasively transcribed but most transcripts do not code for proteins, constituting non-protein-coding RNAs. Despite increasing numbers of functional reports of individual long non-coding RNAs (lncRNAs), assessing the extent of functionality among the non-coding transcriptional output of mammalian cells remains intricate. In the protein-coding world, transcripts differentially expressed in the context of processes essential for the survival of multicellular organisms have been instrumental in the discovery of functionally relevant proteins and their deregulation is frequently associated with diseases. We therefore systematically identified lncRNAs expressed differentially in response to oncologically relevant processes and cell-cycle, p53 and STAT3 pathways, using tiling arrays.

Results

We found that up to 80% of the pathway-triggered transcriptional responses are non-coding. Among these we identified very large macroRNAs with pathway-specific expression patterns and demonstrated that these are likely continuous transcripts. MacroRNAs contain elements conserved in mammals and sauropsids, which in part exhibit conserved RNA secondary structure. Comparing evolutionary rates of a macroRNA to adjacent protein-coding genes suggests a local action of the transcript. Finally, in different grades of astrocytoma, a tumor disease unrelated to the initially used cell lines, macroRNAs are differentially expressed.

Conclusions

It has been shown previously that the majority of expressed non-ribosomal transcripts are non-coding. We now conclude that differential expression triggered by signaling pathways gives rise to a similar abundance of non-coding content. It is thus unlikely that the prevalence of non-coding transcripts in the cell is a trivial consequence of leaky or random transcription events.

Slicing tRNAs to boost functional ncRNA diversity

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. Substrates for such RNA fragmentations are coding as well as non-protein-coding RNAs. In particular, fragments derived from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA fragments possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules. Here we review recent reports on tRNA cleavage and attempt to categorize tRNA pieces according to their origin and cellular function. The biological scope of tRNA-derived fragments ranges from translation control, over RNA silencing, to regulating apoptosis, and thus clearly enlarges the functional repertoire of ncRNA biology.

Dynamic regulation of the transcription initiation landscape at single nucleotide resolution during vertebrate embryogenesis

Spatiotemporal control of gene expression is central to animal development. Core promoters represent a previously unanticipated regulatory level by interacting with cis-regulatory elements and transcription initiation in different physiological and developmental contexts. Here, we provide a first and comprehensive description of the core promoter repertoire and its dynamic use during the development of a vertebrate embryo. By using cap analysis of gene expression (CAGE), we mapped transcription initiation events at single nucleotide resolution across 12 stages of zebrafish development. These CAGE-based transcriptome maps reveal genome-wide rules of core promoter usage, structure, and dynamics, key to understanding the control of gene regulation during vertebrate ontogeny. They revealed the existence of multiple classes of pervasive intra- and intergenic post-transcriptionally processed RNA products and their developmental dynamics. Among these RNAs, we report splice donor site-associated intronic RNA (sRNA) to be specific to genes of the splicing machinery. For the identification of conserved features, we compared the zebrafish data sets to the first CAGE promoter map of Tetraodon and the existing human CAGE data. We show that a number of features, such as promoter type, newly discovered promoter properties such as a specialized purine-rich initiator motif, as well as sRNAs and the genes in which they are detected, are conserved in mammalian and Tetraodon CAGE-defined promoter maps. The zebrafish developmental promoterome represents a powerful resource for studying developmental gene regulation and revealing promoter features shared across vertebrates.

A global profiling of uncapped mRNAs under cold stress reveals specific decay patterns and endonucleolytic cleavages in Brachypodium distachyon

Background

mRNA degradation is a critical factor in determining mRNA abundance and enables rapid adjustment of gene expression in response to environmental stress. The involvement of processing bodies in stress response suggests a role for decapping-mediated mRNA degradation. However, little is known about the role of mRNA degradation under stressful environmental conditions.

Results

Here, we perform a global study of uncapped mRNAs, via parallel analysis of RNA ends (PARE), under cold stress in Brachypodium distachyon. Enrichment analysis indicates that degradation products detected by PARE are mainly generated by the decapping pathway. Endonucleolytic cleavages are detected, uncovering another way of modulating gene expression. PARE and RNA-Seq analyses identify four types of mRNA decay patterns. Type II genes, for which light-harvesting processes are over-represented in gene ontology analyses, show unchanged transcript abundance and altered uncapped transcript abundance. Uncapping-mediated transcript stability of light harvesting-related genes changes significantly in response to cold stress, which may allow rapid adjustments in photosynthetic activity in response to cold stress. Transcript abundance and uncapped transcript abundance for type III genes changes in opposite directions in response to cold stress, indicating that uncapping-mediated mRNA degradation plays a role in regulating gene expression.

Conclusion

To our knowledge, this is the first global analysis of mRNA degradation under environmental stress conditions in Brachypodium distachyon. We uncover specific degradation and endonucleolytic cleavage patterns under cold stress, which will deepen our understanding of mRNA degradation under stressful environmental conditions, as well as the cold stress response mechanism in monocots.

Exosomes secreted by human cells transport largely mRNA fragments that are enriched in the 3'-untranslated regions

Small secreted membrane vesicles called exosomes have recently attracted a great interest after the discovery that they transfer mRNA that can be translated into protein in recipient cells. Surprisingly, we found that for the majority of exosomal mRNAs only a fraction of their corresponding probes is detectable on the expression microarrays. Exosomal mRNA fragmentation is characterized with a specific structural pattern. The closer to the 3'-end of the transcript the fragments are localized, the larger fraction among the secreted RNAs they constitute. Since the 3'-ends of transcripts contain elements conferring subcellular localization of mRNA and are rich in miRNA-binding sites, exosomal RNA may act as competing RNA to regulate stability, localization and translation activity of mRNAs in recipient cells.

REVIEWERS

This article was reviewed by Neil Smalheiser and Sandor Pongor.

Widespread and extensive lengthening of 3' UTRs in the mammalian brain

Remarkable advances in techniques for gene expression profiling have radically changed our knowledge of the transcriptome. Recently, the mammalian brain was reported to express many long intergenic noncoding (lincRNAs) from loci downstream from protein-coding genes. Our experimental tests failed to validate specific accumulation of lincRNA transcripts, and instead revealed strongly distal 3' UTRs generated by alternative cleavage and polyadenylation (APA). With this perspective in mind, we analyzed deep mammalian RNA-seq data using conservative criteria, and identified 2035 mouse and 1847 human genes that utilize substantially distal novel 3' UTRs. Each of these extends at least 500 bases past the most distal 3' termini available in Ensembl v65, and collectively they add 6.6 Mb and 5.1 Mb to the mRNA space of mouse and human, respectively. Extensive Northern analyses validated stable accumulation of distal APA isoforms, including transcripts bearing exceptionally long 3' UTRs (many >10 kb and some >18 kb in length). The Northern data further illustrate that the extensions we annotated were not due to unprocessed transcriptional run-off events. Global tissue comparisons revealed that APA events yielding these extensions were most prevalent in the mouse and human brain. Finally, these extensions collectively contain thousands of conserved miRNA binding sites, and these are strongly enriched for many well-studied neural miRNAs. Altogether, these new 3' UTR annotations greatly expand the scope of post-transcriptional regulatory networks in mammals, and have particular impact on the central nervous system.

Regional differences in gene expression and promoter usage in aged human brains

To characterize the promoterome of caudate and putamen regions (striatum), frontal and temporal cortices, and hippocampi from aged human brains, we used high-throughput cap analysis of gene expression to profile the transcription start sites and to quantify the differences in gene expression across the 5 brain regions. We also analyzed the extent to which methylation influenced the observed expression profiles. We sequenced more than 71 million cap analysis of gene expression tags corresponding to 70,202 promoter regions and 16,888 genes. More than 7000 transcripts were differentially expressed, mainly because of differential alternative promoter usage. Unexpectedly, 7% of differentially expressed genes were neurodevelopmental transcription factors. Functional pathway analysis on the differentially expressed genes revealed an overrepresentation of several signaling pathways (e.g., fibroblast growth factor and wnt signaling) in hippocampus and striatum. We also found that although 73% of methylation signals mapped within genes, the influence of methylation on the expression profile was small. Our study underscores alternative promoter usage as an important mechanism for determining the regional differences in gene expression at old age.

Identification of cytoplasmic capping targets reveals a role for cap homeostasis in translation and mRNA stability

The notion that decapping leads irreversibly to messenger RNA (mRNA) decay was contradicted by the identification of capped transcripts missing portions of their 5' ends and a cytoplasmic complex that can restore the cap on uncapped mRNAs. In this study, we used accumulation of uncapped transcripts in cells inhibited for cytoplasmic capping to identify the targets of this pathway. Inhibition of cytoplasmic capping results in the destabilization of some transcripts and the redistribution of others from polysomes to nontranslating messenger ribonucleoproteins, where they accumulate in an uncapped state. Only a portion of the mRNA transcriptome is affected by cytoplasmic capping, and its targets encode proteins involved in nucleotide binding, RNA and protein localization, and the mitotic cell cycle. The 3' untranslated regions of recapping targets are enriched for AU-rich elements and microRNA binding sites, both of which function in cap-dependent mRNA silencing. These findings identify a cyclical process of decapping and recapping that we term cap homeostasis.

Computational identification and functional predictions of long noncoding RNA in Zea mays

BACKGROUND

Computational analysis of cDNA sequences from multiple organisms suggests that a large portion of transcribed DNA does not code for a functional protein. In mammals, noncoding transcription is abundant, and often results in functional RNA molecules that do not appear to encode proteins. Many long noncoding RNAs (lncRNAs) appear to have epigenetic regulatory function in humans, including HOTAIR and XIST. While epigenetic gene regulation is clearly an essential mechanism in plants, relatively little is known about the presence or function of lncRNAs in plants.

METHODOLOGY/PRINCIPAL FINDINGS

To explore the connection between lncRNA and epigenetic regulation of gene expression in plants, a computational pipeline using the programming language Python has been developed and applied to maize full length cDNA sequences to identify, classify, and localize potential lncRNAs. The pipeline was used in parallel with an SVM tool for identifying ncRNAs to identify the maximal number of ncRNAs in the dataset. Although the available library of sequences was small and potentially biased toward protein coding transcripts, 15% of the sequences were predicted to be noncoding. Approximately 60% of these sequences appear to act as precursors for small RNA molecules and may function to regulate gene expression via a small RNA dependent mechanism. ncRNAs were predicted to originate from both genic and intergenic loci. Of the lncRNAs that originated from genic loci, ∼20% were antisense to the host gene loci.

CONCLUSIONS/SIGNIFICANCE

Consistent with similar studies in other organisms, noncoding transcription appears to be widespread in the maize genome. Computational predictions indicate that maize lncRNAs may function to regulate expression of other genes through multiple RNA mediated mechanisms.

Identification of the human PMR1 mRNA endonuclease as an alternatively processed product of the gene for peroxidasin-like protein

The PMR1 endonuclease was discovered in Xenopus liver and identified as a member of the large and diverse peroxidase gene family. The peroxidase genes arose from multiple duplication and rearrangement events, and their high degree of sequence similarity confounded attempts to identify human PMR1. The functioning of PMR1 in mRNA decay depends on the phosphorylation of a tyrosine in the C-terminal polysome targeting domain by c-Src. The sequences of regions that are required for c-Src binding and phosphorylation of Xenopus PMR1 were used to inform a bioinformatics search that identified two related genes as potential candidates for human PMR1: peroxidasin homolog (PXDN) and peroxidasin homolog-like (PXDNL) protein. Although each of these genes is predicted to encode a large, multidomain membrane-bound peroxidase, alternative splicing of PXDNL pre-mRNA yields a transcript whose predicted product is a 57-kDa protein with 42% sequence identity to Xenopus PMR1. Results presented here confirm the existence of the predicted 57-kDa protein, show this is the only form of PXDNL detected in any of the human cell lines examined, and confirm its identity as human PMR1. Like the Xenopus protein, human PMR1 binds to c-Src, is tyrosine phosphorylated, sediments on polysomes, and catalyzes the selective decay of a PMR1 substrate mRNA. Importantly, the expression of human PMR1 stimulates cell motility in a manner similar to that of the Xenopus PMR1 expressed in human cells, thus providing definitive evidence linking endonuclease decay to the regulation of cell motility.

Dark matter RNA: an intelligent scaffold for the dynamic regulation of the nuclear information landscape

Perhaps no other topic in contemporary genomics has inspired such diverse viewpoints as the 95+% of the genome, previously known as “junk DNA,” that does not code for proteins. Here, we present a theory in which dark matter RNA plays a role in the generation of a landscape of spatial micro-domains coupled to the information signaling matrix of the nuclear landscape. Within and between these micro-domains, dark matter RNAs additionally function to tether RNA interacting proteins and complexes of many different types, and by doing so, allow for a higher performance of the various processes requiring them at ultra-fast rates. This improves signal to noise characteristics of RNA processing, trafficking, and epigenetic signaling, where competition and differential RNA binding among proteins drives the computational decisions inherent in regulatory events.