A genomic analysis of RNA polymerase II modification and chromatin architecture related to 3' end RNA polyadenylation

Intergenic RNA mainly derives from nascent transcripts of known genes

BACKGROUND

Eukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.

RESULTS

We hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the "fuzzy" transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.

CONCLUSIONS

We provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

Regulation of alternative polyadenylation in the yeast Saccharomyces cerevisiae by histone H3K4 and H3K36 methyltransferases

Adjusting DNA structure via epigenetic modifications, and altering polyadenylation (pA) sites at which precursor mRNA is cleaved and polyadenylated, allows cells to quickly respond to environmental stress. Since polyadenylation occurs co-transcriptionally, and specific patterns of nucleosome positioning and chromatin modifications correlate with pA site usage, epigenetic factors potentially affect alternative polyadenylation (APA). We report that the histone H3K4 methyltransferase Set1, and the histone H3K36 methyltransferase Set2, control choice of pA site in Saccharomyces cerevisiae, a powerful model for studying evolutionarily conserved eukaryotic processes. Deletion of SET1 or SET2 causes an increase in serine-2 phosphorylation within the C-terminal domain of RNA polymerase II (RNAP II) and in the recruitment of the cleavage/polyadenylation complex, both of which could cause the observed switch in pA site usage. Chemical inhibition of TOR signaling, which causes nutritional stress, results in Set1- and Set2-dependent APA. In addition, Set1 and Set2 decrease efficiency of using single pA sites, and control nucleosome occupancy around pA sites. Overall, our study suggests that the methyltransferases Set1 and Set2 regulate APA induced by nutritional stress, affect the RNAP II C-terminal domain phosphorylation at Ser2, and control recruitment of the 3′ end processing machinery to the vicinity of pA sites.

Emerging Roles of RNA 3'-end Cleavage and Polyadenylation in Pathogenesis, Diagnosis and Therapy of Human Disorders

A crucial feature of gene expression involves RNA processing to produce 3′ ends through a process termed 3′ end cleavage and polyadenylation (CPA). This ensures the nascent RNA molecule can exit the nucleus and be translated to ultimately give rise to a protein which can execute a function. Further, alternative polyadenylation (APA) can produce distinct transcript isoforms, profoundly expanding the complexity of the transcriptome. CPA is carried out by multi-component protein complexes interacting with multiple RNA motifs and is tightly coupled to transcription, other steps of RNA processing, and even epigenetic modifications. CPA and APA contribute to the maintenance of a multitude of diverse physiological processes. It is therefore not surprising that disruptions of CPA and APA can lead to devastating disorders. Here, we review potential CPA and APA mechanisms involving both loss and gain of function that can have tremendous impacts on health and disease. Ultimately we highlight the emerging diagnostic and therapeutic potential CPA and APA offer.

Distribution pattern of histone marks potentially determines their roles in transcription and RNA processing in rice

Histone marks including histone modifications and histone variants may affect the processes of gene transcription and co-transcriptional RNA processing depending on their specific deposition patterns within genes. Here, we analyzed distribution patterns of rice histone marks and divided them into seven clusters according to their enrichment in promoter, transcription start site (TSS), and gene body regions. Expression levels of the genes in each cluster were explored to disclose the importance of histone marks in the processes of transcription. We show that: a) H3K4me3 and histone acetylation marks show locally different distributions at TSS, implying that they may play different roles in transcription initiation. b) H3K36me1 enriched at TSS has a negative effect on transcription. c) Genes with high level of expression were marked by H3K36me3 at both the TSS and body regions. In addition, we found that H3K4me2, H3K23ac, H3K4ac, and H2A.Z show exon-biased enrichment, suggesting they may be chromatin marks involved in co-transcriptional splicing. Finally, we identified histone marks that discriminate constitutive expression genes (CEGs) from tissue-specific expressed genes (TSEGs). Taken together, the analysis revealed distribution patterns of different histone marks in rice to infer their potential roles in transcription and RNA processing. The results lay foundation for further understanding the mechanism by which histone marks are involved in the regulation of these processes in plants.

Control of RNA Pol II Speed by PNUTS-PP1 and Spt5 Dephosphorylation Facilitates Termination by a "Sitting Duck Torpedo" Mechanism

Control of transcription speed, which influences many co-transcriptional processes, is poorly understood. We report that PNUTS-PP1 phosphatase is a negative regulator of RNA polymerase II (Pol II) elongation rate. The PNUTS W401A mutation, which disrupts PP1 binding, causes genome-wide acceleration of transcription associated with hyper-phosphorylation of the Spt5 elongation factor. Immediately downstream of poly(A) sites, Pol II decelerates from >2 kb/min to <1 kb/min, which correlates with Spt5 dephosphorylation. Pol II deceleration and Spt5 dephosphorylation require poly(A) site recognition and the PNUTS-PP1 complex, which is in turn necessary for transcription termination. These results lead to a model for termination, the "sitting duck torpedo" mechanism, where poly(A) site-dependent deceleration caused by PNUTS-PP1 and Spt5 dephosphorylation is required to convert Pol II into a viable target for the Xrn2 terminator exonuclease. Spt5 and its bacterial homolog NusG therefore have related functions controlling kinetic competition between RNA polymerases and the termination factors that pursue them.

Interactions between the transcription and replication machineries regulate the RNA and DNA synthesis in the herpesviruses

The temporal coordination of viral gene expression is imperative for the regulation of the herpesvirus replication cycle. While the main factors of this transcriptional coordination are known, the subtler control mechanisms of gene expression remain elusive. Recent long read sequencing-based approached have revealed an intricate meshwork of overlaps between the herpesvirus transcripts and the overlap of the replication origins with noncoding RNAs. It has been shown that the transcriptional apparatuses can physically interfere with one another while transcribing overlapping regions. We hypothesize that transcriptional interference regulates the global gene expression across the herpesvirus genome. Additionally, an overall decrease in transcriptional activity in individual viral genes has been observed following the onset of DNA replication. An overlap of the replication origins with specific transcripts has also been described in several herpesviruses. The genome-wide interactions between the transcriptional apparatuses and between the replication and transcriptional machineries suggest the existence of novel layers of genetic regulation.

Nuclear speckles: molecular organization, biological function and role in disease

The nucleoplasm is not homogenous; it consists of many types of nuclear bodies, also known as nuclear domains or nuclear subcompartments. These self-organizing structures gather machinery involved in various nuclear activities. Nuclear speckles (NSs) or splicing speckles, also called interchromatin granule clusters, were discovered as sites for splicing factor storage and modification. Further studies on transcription and mRNA maturation and export revealed a more general role for splicing speckles in RNA metabolism. Here, we discuss the functional implications of the localization of numerous proteins crucial for epigenetic regulation, chromatin organization, DNA repair and RNA modification to nuclear speckles. We highlight recent advances suggesting that NSs facilitate integrated regulation of gene expression. In addition, we consider the influence of abundant regulatory and signaling proteins, i.e. protein kinases and proteins involved in protein ubiquitination, phosphoinositide signaling and nucleoskeletal organization, on pre-mRNA synthesis and maturation. While many of these regulatory proteins act within NSs, direct evidence for mRNA metabolism events occurring in NSs is still lacking. NSs contribute to numerous human diseases, including cancers and viral infections. In addition, recent data have demonstrated close relationships between these structures and the development of neurological disorders.

Posttranscriptional Regulation of HLA-A Protein Expression by Alternative Polyadenylation Signals Involving the RNA-Binding Protein Syncrip

Genomic variation in the untranslated region (UTR) has been shown to influence HLA class I expression level and associate with disease outcomes. Sequencing of the 3'UTR of common HLA-A alleles indicated the presence of two polyadenylation signals (PAS). The proximal PAS is conserved, whereas the distal PAS is disrupted within certain alleles by sequence variants. Using 3'RACE, we confirmed expression of two distinct forms of the HLA-A 3'UTR based on use of either the proximal or the distal PAS, which differ in length by 100 bp. Specific HLA-A alleles varied in the usage of the proximal versus distal PAS, with some alleles using only the proximal PAS, and others using both the proximal and distal PAS to differing degrees. We show that the short and the long 3'UTR produced similar mRNA expression levels. However, the long 3'UTR conferred lower luciferase activity as compared with the short form, indicating translation inhibition of the long 3'UTR. RNA affinity pull-down followed by mass spectrometry analysis as well as RNA coimmunoprecipitation indicated differential binding of Syncrip to the long versus short 3'UTR. Depletion of Syncrip by small interfering RNA increased surface expression of an HLA-A allotype that uses primarily the long 3'UTR, whereas an allotype expressing only the short form was unaffected. Furthermore, specific blocking of the proximal 3'UTR reduced surface expression without decreasing mRNA expression. These data demonstrate HLA-A allele-specific variation in PAS usage, which modulates their cell surface expression posttranscriptionally.

Processing and transcriptome expansion at the mRNA 3' end in health and disease: finding the right end

The human transcriptome is highly dynamic, with each cell type, tissue, and organ system expressing an ensemble of transcript isoforms that give rise to considerable diversity. Apart from alternative splicing affecting the "body" of the transcripts, extensive transcriptome diversification occurs at the 3' end. Transcripts differing at the 3' end can have profound physiological effects by encoding proteins with distinct functions or regulatory properties or by affecting the mRNA fate via the inclusion or exclusion of regulatory elements (such as miRNA or protein binding sites). Importantly, the dynamic regulation at the 3' end is associated with various (patho)physiological processes, including the immune regulation but also tumorigenesis. Here, we recapitulate the mechanisms of constitutive mRNA 3' end processing and review the current understanding of the dynamically regulated diversity at the transcriptome 3' end. We illustrate the medical importance by presenting examples that are associated with perturbations of this process and indicate resulting implications for molecular diagnostics as well as potentially arising novel therapeutic strategies.

Effects of Transcription Elongation Rate and Xrn2 Exonuclease Activity on RNA Polymerase II Termination Suggest Widespread Kinetic Competition

The torpedo model of transcription termination asserts that the exonuclease Xrn2 attacks the 5'PO4-end exposed by nascent RNA cleavage and chases down the RNA polymerase. We tested this mechanism using a dominant-negative human Xrn2 mutant and found that it delayed termination genome-wide. Xrn2 nuclease inactivation caused strong termination defects downstream of most poly(A) sites and modest delays at some histone and U snRNA genes, suggesting that the torpedo mechanism is not limited to poly(A) site-dependent termination. A central untested feature of the torpedo model is that there is kinetic competition between the exonuclease and the pol II elongation complex. Using pol II rate mutants, we found that slow transcription robustly shifts termination upstream, and fast elongation extends the zone of termination further downstream. These results suggest that kinetic competition between elongating pol II and the Xrn2 exonuclease is integral to termination of transcription on most human genes.

Intergenic Alu exonisation facilitates the evolution of tissue-specific transcript ends

The 3' untranslated regions (3' UTRs) of transcripts serve as important hubs for posttranscriptional gene expression regulation. Here, we find that the exonisation of intergenic Alu elements introduced new terminal exons and polyadenylation sites during human genome evolution. While Alu exonisation from introns has been described previously, we shed light on a novel mechanism to create alternative 3' UTRs, thereby opening opportunities for differential posttranscriptional regulation. On the mechanistic level, we show that intergenic Alu exonisation can compete both with alternative splicing and polyadenylation in the upstream gene. Notably, the Alu-derived isoforms are often expressed in a tissue-specific manner, and the Alu-derived 3' UTRs can alter mRNA stability. In summary, we demonstrate that intergenic elements can affect processing of preceding genes, and elucidate how intergenic Alu exonisation can contribute to tissue-specific posttranscriptional regulation by expanding the repertoire of 3' UTRs.

Coordination of RNA Polymerase II Pausing and 3' End Processing Factor Recruitment with Alternative Polyadenylation

Most mammalian genes produce transcripts whose 3' ends are processed at multiple alternative positions by cleavage/polyadenylation (CPA). Poly(A) site cleavage frequently occurs cotranscriptionally and is facilitated by CPA factor binding to the RNA polymerase II (Pol II) C-terminal domain (CTD) phosphorylated on Ser2 residues of its heptad repeats (YS2PTSPS). The function of cotranscriptional events in the selection of alternative poly(A) sites is poorly understood. We investigated Pol II pausing, CTD Ser2 phosphorylation, and processing factor CstF recruitment at wild-type and mutant IgM transgenes that use alternative poly(A) sites to produce mRNAs encoding the secreted and membrane-bound forms of the immunoglobulin (Ig) heavy chain. The results show that the sites of Pol II pausing and processing factor recruitment change depending on which poly(A) site is utilized. In contrast, the extent of Pol II CTD Ser2 phosphorylation does not closely correlate with poly(A) site selection. We conclude that changes in properties of the transcription elongation complex closely correlate with utilization of different poly(A) sites, suggesting that cotranscriptional events may influence the decision between alternative modes of pre-mRNA 3' end processing.

In vitro application of ribonucleases: comparison of the effects on mRNA and miRNA stability

Background

MicroRNA has become important in a wide range of research interests. Due to the increasing number of known microRNAs, these molecules are likely to be increasingly seen as a new class of biomarkers. This is driven by the fact that microRNAs are relatively stable when circulating in the plasma. Despite extensive analysis of mechanisms involved in microRNA processing, relatively little is known about the in vitro decay of microRNAs under defined conditions or about the relative stabilities of mRNAs and microRNAs.

Methods

In this in vitro study, equal amounts of total RNA of identical RNA pools were treated with different ribonucleases under defined conditions. Degradation of total RNA was assessed using microfluidic analysis mainly based on ribosomal RNA. To evaluate the influence of the specific RNases on the different classes of RNA (ribosomal RNA, mRNA, miRNA) ribosomal RNA as well as a pattern of specific mRNAs and miRNAs was quantified using RT-qPCR assays. By comparison to the untreated control sample the ribonuclease-specific degradation grade depending on the RNA class was determined.

Results

In the present in vitro study we have investigated the stabilities of mRNA and microRNA with respect to the influence of ribonucleases used in laboratory practice. Total RNA was treated with specific ribonucleases and the decay of different kinds of RNA was analysed by RT-qPCR and miniaturized gel electrophoresis. In addition, we have examined whether the integrity observed for ribosomal RNA is applicable to microRNA and mRNA. Depending on the kind of ribonuclease used, our results demonstrated a higher stability of microRNA relative to mRNA and a limitation of the relevance of ribosomal RNA integrity to the integrity of other RNA groups.

Conclusion

Our results suggest that the degradation status of ribosomal RNA is not always applicable to mRNA and microRNA. In fact, the stabilities of these RNA classes to exposure to ribonucleases are independent from each other, with microRNA being more stable than mRNA. The relative stability of microRNAs supports their potential and further development as biomarkers in a range of applications.

The nucleosome regulates the usage of polyadenylation sites in the human genome

Background

It has been reported that 3' end processing is coupled to transcription and nucleosome depletion near the polyadenylation sites in many species. However, the association between nucleosome occupancy and polyadenylation site usage is still unclear.

Results

By systematic analysis of high-throughput sequencing datasets from the human genome, we found that nucleosome occupancy patterns are different around the polyadenylation sites, and that the patterns associate with both transcription termination and recognition of polyadenylation sites. Upstream of proximal polyadenylation sites, RNA polymerase II accumulated and nucleosomes were better positioned compared with downstream of the sites. Highly used proximal polyadenylation sites had higher upstream nucleosome levels and RNA polymerase II accumulation than lowly used sites. This suggests that nucleosomes positioned upstream of proximal sites function in the recognition of proximal polyadenylation sites and in the preparation for 3' end processing by slowing down transcription speed. Both conserved distal polyadenylation sites and constitutive sites showed stronger nucleosome depletion near polyadenylation sites and had intrinsically better positioned downstream nucleosomes. Finally, there was a higher accumulation of RNA polymerase II downstream of the polyadenylation sites, to guarantee gene transcription termination and recognition of the last polyadenylation sites, if previous sites were missed.

Conclusions

Our study indicates that nucleosome arrays play different roles in the regulation of the usage of polyadenylation sites and transcription termination of protein-coding genes, and form a dual pausing model of RNA polymerase II in the alternative polyadenylation sites’ region, to ensure effective 3' end processing.

Inhibition of cdk9 during herpes simplex virus 1 infection impedes viral transcription

During herpes simplex virus 1 (HSV-1) infection there is a loss of the serine-2 phosphorylated form of RNA polymerase II (RNAP II) found in elongation complexes. This occurs in part because RNAP II undergoes ubiquitination and proteasomal degradation during times of highly active viral transcription, which may result from stalled elongating complexes. In addition, a viral protein, ICP22, was reported to trigger a loss of serine-2 RNAP II. These findings have led to some speculation that the serine-2 phosphorylated form of RNAP II may not be required for HSV-1 transcription, although this form is required for cellular transcription elongation and RNA processing. Cellular kinase cdk9 phosphorylates serine-2 in the C-terminal domain (CTD) of RNAP II. To determine if serine-2 phosphorylated RNAP II is required for HSV-1 transcription, we inhibited cdk9 during HSV-1 infection and measured viral gene expression. Inhibition was achieved by adding cdk9 inhibitors 5,6-dichlorobenzimidazone-1-β-D-ribofuranoside (DRB) or flavopiridol (FVP) or by expression of a dominant–negative cdk9 or HEXIM1, which in conjunction with 7SK snRNA inhibits cdk9 in complex with cyclin 1. Here we report that inhibition of cdk9 resulted in decreased viral yields and levels of late proteins, poor formation of viral transcription-replication compartments, reduced levels of poly(A)+ mRNA and decreased RNA synthesis as measured by uptake of 5-bromouridine into nascent RNA. Importantly, a global reduction in viral mRNAs was seen as determined by microarray analysis. We conclude that serine-2 phosphorylation of the CTD of RNAP II is required for HSV-1 transcription.

Distinct polyadenylation landscapes of diverse human tissues revealed by a modified PA-seq strategy

Background

Polyadenylation is a key regulatory step in eukaryotic gene expression and one of the major contributors of transcriptome diversity. Aberrant polyadenylation often associates with expression defects and leads to human diseases.

Results

To better understand global polyadenylation regulation, we have developed a polyadenylation sequencing (PA-seq) approach. By profiling polyadenylation events in 13 human tissues, we found that alternative cleavage and polyadenylation (APA) is prevalent in both protein-coding and noncoding genes. In addition, APA usage, similar to gene expression profiling, exhibits tissue-specific signatures and is sufficient for determining tissue origin. A 3′ untranslated region shortening index (USI) was further developed for genes with tandem APA sites. Strikingly, the results showed that different tissues exhibit distinct patterns of shortening and/or lengthening of 3′ untranslated regions, suggesting the intimate involvement of APA in establishing tissue or cell identity.

Conclusions

This study provides a comprehensive resource to uncover regulated polyadenylation events in human tissues and to characterize the underlying regulatory mechanism.

An alternative polyadenylation mechanism coopted to the Arabidopsis RPP7 gene through intronic retrotransposon domestication

Transposable elements (TEs) can drive evolution by creating genetic and epigenetic variation. Although examples of adaptive TE insertions are accumulating, proof that epigenetic information carried by such "domesticated" TEs has been coopted to control host gene function is still limited. We show that COPIA-R7, a TE inserted into the Arabidopsis thaliana disease resistance gene RPP7 recruited the histone mark H3K9me2 to this locus. H3K9me2 levels at COPIA-R7 affect the choice between two alternative RPP7 polyadenylation sites in the pre-mRNA and, thereby, influence the critical balance between RPP7-coding and non-RPP7-coding transcript isoforms. Function of RPP7 is fully dependent on high levels of H3K9me2 at COPIA-R7. We present a direct in vivo demonstration for cooption of a TE-associated histone mark to the epigenetic control of pre-mRNA processing and establish a unique mechanism for regulation of plant immune surveillance gene expression. Our results functionally link a histone mark to alternative polyadenylation and the balance between distinct transcript isoforms from a single gene.

Chromatin loop organization of the junb locus in mouse dendritic cells

The junb gene behaves as an immediate early gene in bacterial lipopolysaccharide (LPS)-stimulated dendritic cells (DCs), where its transient transcriptional activation is necessary for the induction of inflammatory cytokines. junb is a short gene and its transcriptional activation by LPS depends on the binding of NF-κB to an enhancer located just downstream of its 3′ UTR. Here, we have addressed the mechanisms underlying the transcriptional hyper-reactivity of junb. Using transfection and pharmacological assays to complement chromatin immunoprecipitation analyses addressing the localization of histones, polymerase II, negative elongation factor (NELF)-, DRB sensitivity-inducing factor (DSIF)- and Positive Transcription Factor b complexes, we demonstrate that junb is a RNA Pol II-paused gene where Pol II is loaded in the transcription start site domain but poorly active. Moreover, High salt-Recovered Sequence, chromosome conformation capture (3C)- and gene transfer experiments show that (i) junb is organized in a nuclear chromatin loop bringing into close spatial proximity the upstream promoter region and the downstream enhancer and (ii) this configuration permits immediate Pol II release on the junb body on binding of LPS-activated NF-κB to the enhancer. Thus, our work unveils a novel topological framework underlying fast junb transcriptional response in DCs. Moreover, it also points to a novel layer of complexity in the modes of action of NF-κB.

Alterations in polyadenylation and its implications for endocrine disease

INTRODUCTION

Polyadenylation is the process in which the pre-mRNA is cleaved at the poly(A) site and a poly(A) tail is added - a process necessary for normal mRNA formation. Genes with multiple poly(A) sites can undergo alternative polyadenylation (APA), producing distinct mRNA isoforms with different 3' untranslated regions (3' UTRs) and in some cases different coding regions. Two thirds of all human genes undergo APA. The efficiency of the polyadenylation process regulates gene expression and APA plays an important part in post-transcriptional regulation, as the 3' UTR contains various cis-elements associated with post-transcriptional regulation, such as target sites for micro-RNAs and RNA-binding proteins. Implications of alterations in polyadenylation for endocrine disease: Alterations in polyadenylation have been found to be causative of neonatal diabetes and IPEX (immune dysfunction, polyendocrinopathy, enteropathy, X-linked) and to be associated with type I and II diabetes, pre-eclampsia, fragile X-associated premature ovarian insufficiency, ectopic Cushing syndrome, and many cancer diseases, including several types of endocrine tumor diseases.

PERSPECTIVES

Recent developments in high-throughput sequencing have made it possible to characterize polyadenylation genome-wide. Antisense elements inhibiting or enhancing specific poly(A) site usage can induce desired alterations in polyadenylation, and thus hold the promise of new therapeutic approaches.

SUMMARY

This review gives a detailed description of alterations in polyadenylation in endocrine disease, an overview of the current literature on polyadenylation and summarizes the clinical implications of the current state of research in this field.

Bromodomain-containing protein 4 (BRD4) regulates RNA polymerase II serine 2 phosphorylation in human CD4+ T cells

Transcriptional elongation by RNA polymerase II (Pol II) is regulated by positive transcription elongation factor b (P-TEFb) in association with bromodomain-containing protein 4 (BRD4). We used genome-wide chromatin immunoprecipitation sequencing in primary human CD4+ T cells to reveal that BRD4 co-localizes with Ser-2-phosphorylated Pol II (Pol II Ser-2) at both enhancers and promoters of active genes. Disruption of bromodomain-histone acetylation interactions by JQ1, a small-molecule bromodomain inhibitor, resulted in decreased BRD4 binding, reduced Pol II Ser-2, and reduced expression of lineage-specific genes in primary human CD4+ T cells. A large number of JQ1-disrupted BRD4 binding regions exhibited diacetylated H4 (lysine 5 and -8) and H3K27 acetylation (H3K27ac), which correlated with the presence of histone acetyltransferases and deacetylases. Genes associated with BRD4/H3K27ac co-occupancy exhibited significantly higher activity than those associated with H3K27ac or BRD4 binding alone. Comparison of BRD4 binding in T cells and in human embryonic stem cells revealed that enhancer BRD4 binding sites were predominantly lineage-specific. Our findings suggest that BRD4-driven Pol II phosphorylation at serine 2 plays an important role in regulating lineage-specific gene transcription in human CD4+ T cells.

Alternative transcription and alternative splicing in cancer

In recent years, the notion of "one gene makes one protein that functions in one signaling pathway" in mammalian cells has been shown to be overly simplistic. Recent genome-wide studies suggest that at least half of the human genes, including many therapeutic target genes, produce multiple protein isoforms through alternative splicing and alternative usage of transcription initiation and/or termination. For example, alternative splicing of the vascular endothelial growth factor gene (VEGFA) produces multiple protein isoforms, which display either pro-angiogenic or anti-angiogenic activities. Similarly, for the majority of human genes, the inclusion or exclusion of exonic sequences enhances the generation of transcript variants and/or protein isoforms that can vary in structure and functional properties. Many of the isoforms produced in this manner are tightly regulated during normal development but are misregulated in cancer cells. Altered expression of transcript variants and protein isoforms for numerous genes is linked with disease and its prognosis, and cancer cells manipulate regulatory mechanisms to express specific isoforms that confer drug resistance and survival advantages. Emerging insights indicate that modulating the expression of transcript and protein isoforms of a gene may hold the key to impeding tumor growth and act as a model for efficient targeting of disease-associated genes at the isoform level. This review highlights the role and regulation of alternative transcription and splicing mechanisms in generating the transcriptome, and the misuse and diagnostic/prognostic potential of alternative transcription and splicing in cancer.

RNA polymerase II pausing downstream of core histone genes is different from genes producing polyadenylated transcripts

Recent genome-wide chromatin immunoprecipitation coupled high throughput sequencing (ChIP-seq) analyses performed in various eukaryotic organisms, analysed RNA Polymerase II (Pol II) pausing around the transcription start sites of genes. In this study we have further investigated genome-wide binding of Pol II downstream of the 3′ end of the annotated genes (EAGs) by ChIP-seq in human cells. At almost all expressed genes we observed Pol II occupancy downstream of the EAGs suggesting that Pol II pausing 3′ from the transcription units is a rather common phenomenon. Downstream of EAGs Pol II transcripts can also be detected by global run-on and sequencing, suggesting the presence of functionally active Pol II. Based on Pol II occupancy downstream of EAGs we could distinguish distinct clusters of Pol II pause patterns. On core histone genes, coding for non-polyadenylated transcripts, Pol II occupancy is quickly dropping after the EAG. In contrast, on genes, whose transcripts undergo polyA tail addition [poly(A)⁺], Pol II occupancy downstream of the EAGs can be detected up to 4–6 kb. Inhibition of polyadenylation significantly increased Pol II occupancy downstream of EAGs at poly(A)⁺ genes, but not at the EAGs of core histone genes. The differential genome-wide Pol II occupancy profiles 3′ of the EAGs have also been confirmed in mouse embryonic stem (mES) cells, indicating that Pol II pauses genome-wide downstream of the EAGs in mammalian cells. Moreover, in mES cells the sharp drop of Pol II signal at the EAG of core histone genes seems to be independent of the phosphorylation status of the C-terminal domain of the large subunit of Pol II. Thus, our study uncovers a potential link between different mRNA 3′ end processing mechanisms and consequent Pol II transcription termination processes.

Transcription goes digital

Transcription is a complex process that integrates the state of the cell and its environment to generate adequate responses for cell fitness and survival. Recent microscopy experiments have been able to monitor transcription from single genes in individual cells. These observations have revealed two striking features: transcriptional activity can vary markedly from one cell to another, and is subject to large changes over time, sometimes within minutes. How the chromatin structure, transcription machinery assembly and signalling networks generate such patterns is still unclear. In this review, we present the techniques used to investigate transcription from single genes, introduce quantitative modelling tools, and discuss transcription mechanisms and their implications for gene expression regulation.

Immunoglobulin heavy chain gene regulation through polyadenylation and splicing competition

The immunoglobulin heavy chain (IgH) genes, which encode one of the two chains of antibody molecules, were the first cellular genes shown to undergo developmentally regulated alternative RNA processing. These genes produce two different mRNAs from a single primary transcript. One mRNA is cleaved and polyadenylated at an upstream poly(A) signal while the other mRNA removes this poly(A) signal by RNA splicing and is cleaved and polyadenylated at a downstream poly(A) site. A broad range of studies have been performed to understand the mechanism of IgH RNA processing regulation during B lymphocyte development. The model that has emerged is much more complex than envisioned by the earliest view of regulation through poly(A) signal choice. Regulation requires that the IgH gene contain competing splice and cleavage-polyadenylation reactions with balanced efficiencies. Because non-IgH genes with these structural features also can be regulated, IgH gene-specific sequence elements are not required for regulation. Changes in cleavage-polyadenylation and RNA splicing, as well as pol II elongation, all contribute to IgH developmental RNA processing regulation. Multiple factors are likely involved in the regulation during B lymphocyte maturation. Additional biologically relevant factors that contribute to IgH regulation remain to be identified and incorporated into a mechanistic model for regulation. Much of the work to date confirms the complex nature of IgH mRNA regulation and suggests that a thorough understanding of this control will remain a challenge. However, it is also likely that such understanding will help elucidate novel mechanisms of RNA processing regulation.

CREPT accelerates tumorigenesis by regulating the transcription of cell-cycle-related genes

Tumorigenesis is caused by an uncontrolled cell cycle and the altered expression of many genes. Here, we report a gene CREPT that is preferentially expressed in diverse human tumors. Overexpression of CREPT accelerates tumor growth, whereas depletion of CREPT demonstrates a reversed effect. CREPT regulates cyclin D1 expression by binding to its promoter, enhancing its transcription both in vivo and in vitro, and interacting with RNA polymerase II (RNAPII). Interestingly, CREPT promotes the formation of a chromatin loop and prevents RNAPII from reading through the 3' end termination site of the gene. Our findings reveal a mechanism where CREPT increases cyclin D1 transcription during tumorigenesis, through enhancing the recruitment of RNAPII to the promoter region, possibly, as well as chromatin looping.

Long conserved fragments upstream of Mammalian polyadenylation sites

Polyadenylation is a cotranscriptional nuclear RNA processing event involving endonucleolytic cleavage of the nascent, emerging pre-messenger RNA (pre-mRNA) from the RNA polymerase, immediately followed by the polymerization of adenine ribonucleotides, called the poly(A) tail, to the cleaved 3' end of the polyadenylation site (PAS). This apparently simple molecular processing step has been discovered to be connected to transcription and splicing therefore increasing its potential for regulation of gene expression. Here, through a bioinformatic analysis of cis-PAS-regulatory elements in mammals that includes taking advantage of multiple evolutionary time scales, we find unexpected selection pressure much further upstream, up to 200 nt, from the PAS than previously thought. Strikingly, close to 3,000 long (30-500 nt) noncoding conserved fragments (CFs) were discovered in the PAS flanking region of three remotely related mammalian species, human, mouse, and cow. When an even more remote transitional mammal, platypus, was included, still over a thousand CFs were found in the proximity of the PAS. Even though the biological function of these CFs remains unknown, their considerable sizes makes them unlikely to serve as protein recognition sites, which are typically ≤15 nt. By harnessing genome wide DNaseI hypersensitivity data, we have discovered that the presence of CFs correlates with chromatin accessibility. Our study is important in highlighting novel experimental targets, which may provide new understanding about the regulatory aspects of polyadenylation.

Epigenomic and RNA structural correlates of polyadenylation

Polyadenylation (poly(A)) of mRNA plays a critical role in regulating gene expression. Identifying the sequence, structural, and epigenomic determinants of poly(A) site usage is an important long term goal. Several cis elements that mediate poly(A) regulation have been identified. Highly used poly(A) sites are also known to have a greater nucleosome occupancy in the immediate downstream. However, a detailed exploration of additional epigenomic and mRNA structural correlates of poly(A) site usage has not been reported. Importantly, functional interaction between sequence, structure, and the epigenome in determining the poly(A) site usage is not known. We show that highly used poly(A) sites are positively associated with an mRNA structure that is energetically more favorable and one that better exposes a critical polyadenylation cis element. In exploring potential interplay between RNA and chromatin structure, we found that a stronger nucleosome occupancy downstream of poly(A) site strongly correlated with (1) a more favorable mRNA structure, and (2) a greater accumulation of RNA Polymerase II (PolII) at the poly(A) site. Further analysis suggested a causal relationship pointing from PolII accumulation to a stable RNA structure. Additionally, we found that distinct patterns of histone modifications characterize poly(A) sites and these epigenetic patterns alone can distinguish true poly(A) sites with ~76% accuracy and also discriminate between high and low usage poly(A) sites with ~74% accuracy. Our results suggest a causative link between chromatin structure and mRNA structure whereby a compacted chromatin downstream of the poly(A) site slows down the elongating transcript, thus facilitating the folding of nascent mRNA in a favorable structure at poly(A) site during transcription. Additionally we report hitherto unknown epigenomic correlates for poly(A) site usage.

Tiling array data analysis: a multiscale approach using wavelets

BACKGROUND

Tiling array data is hard to interpret due to noise. The wavelet transformation is a widely used technique in signal processing for elucidating the true signal from noisy data. Consequently, we attempted to denoise representative tiling array datasets for ChIP-chip experiments using wavelets. In doing this, we used specific wavelet basis functions, Coiflets, since their triangular shape closely resembles the expected profiles of true ChIP-chip peaks.

RESULTS

In our wavelet-transformed data, we observed that noise tends to be confined to small scales while the useful signal-of-interest spans multiple large scales. We were also able to show that wavelet coefficients due to non-specific cross-hybridization follow a log-normal distribution, and we used this fact in developing a thresholding procedure. In particular, wavelets allow one to set an unambiguous, absolute threshold, which has been hard to define in ChIP-chip experiments. One can set this threshold by requiring a similar confidence level at different length-scales of the transformed signal. We applied our algorithm to a number of representative ChIP-chip data sets, including those of Pol II and histone modifications, which have a diverse distribution of length-scales of biochemical activity, including some broad peaks.

CONCLUSIONS

Finally, we benchmarked our method in comparison to other approaches for scoring ChIP-chip data using spike-ins on the ENCODE Nimblegen tiling array. This comparison demonstrated excellent performance, with wavelets getting the best overall score.

Alternative polyadenylation of antisense RNAs and flowering time control

Flowering time is controlled by precision in gene regulation mediated by different pathways. Two Arabidopsis thaliana components of the autonomous flowering pathway, FCA and FPA, function as genetically independent trans-acting regulators of alternative cleavage and polyadenylation. FCA and FPA directly associate with chromatin at the locus encoding the floral repressor FLC, but appear to control FLC transcription by mediating alternative polyadenylation of embedded non-coding antisense RNAs. These findings prompt the re-examination of how other factors control FLC expression, as it is formally possible that they function primarily to control alternative processing of antisense RNAs. As co-expressed sense and antisense gene pairs are widespread in eukaryotes, alternative processing of antisense RNAs may represent a significant form of gene regulation.

Alternative mRNA polyadenylation in eukaryotes: an effective regulator of gene expression

Alternative RNA processing mechanisms, including alternative splicing and alternative polyadenylation, are increasingly recognized as important regulators of gene expression. This article will focus on what has recently been described about alternative polyadenylation in development, differentiation, and disease in higher eukaryotes. We will also describe how the evolving global methodologies for examining the cellular transcriptome, both experimental and bioinformatic, are revealing new details about the complex nature of alternative 3(') end formation as well as interactions with other RNA-mediated and RNA processing mechanisms.

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.

Alternative polyadenylation blooms

Alternative polyadenylation generates mRNAs with 3' untranscribed regions of different lengths, often affecting transcript stability. Hornyik et al., in this issue of Developmental Cell, and Liu et al. now demonstrate a role for alternative polyadenylation in gene silencing and the regulation of flowering time in Arabidopsis thaliana.

Detection of novel 3' untranslated region extensions with 3' expression microarrays

Background

The 3' untranslated regions (UTRs) of transcripts are not well characterized for many genes and often extend beyond the annotated regions. Since Affymetrix 3' expression arrays were designed based on expressed sequence tags, many probesets map to intergenic regions downstream of genes. We used expression information from these probesets to predict transcript extension beyond currently known boundaries.

Results

Based on our dataset encompassing expression in 22 different murine tissues, we identified 845 genes with predicted 3'UTR extensions. These extensions have a similar conservation as known 3'UTRs, which is distinctly higher than intergenic regions. We verified 8 of the predictions by PCR and found all of the predicted regions to be expressed. The method can be extended to other 3' expression microarray platforms as we demonstrate with human data. Additional confirming evidence was obtained from public paired end read data.

Conclusions

We show that many genes have 3'UTR regions extending beyond currently known gene regions and provide a method to identify such regions based on microarray expression data. Since 3' UTR contain microRNA binding sites and other stability determining regions, identification of the full length 3' UTR is important to elucidate posttranscriptional regulation.

Stability of histone modifications across mammalian genomes: implications for 'epigenetic' marking

The combination of chromatin immunoprecipitation (ChIP) with microarray analysis (ChIP-chip) or high-throughput sequencing (ChIP-seq and ChIP-SAGE) has provided maps of a wide variety of site-specific histone modifications across mammalian genomes in various cell types. Although distinct genomic regions and functional elements have been strongly associated with specific histone modifications, an overwhelming number of combinatorial patterns have also been observed across the genome. While peaks of enrichment in ChIP-chip and ChIP-seq data may suggest stable and predictive 'landmarks' across the genomic landscape, studies from transcribed genes indicate a more dynamic model of how these data may be interpreted. In light of such studies, which show highly dynamic methylation, acetylation and phosphorylation of histone H3 during gene transcription, we consider the extent to which genome-wide maps of chromatin state could be interpreted as 'snapshots' of heterogeneous profiles deriving from dynamic modification processes. Rather than acting as static 'epigenetic' landmarks, histone modifications may function as dynamic and transient operational marks supporting specific steps in diverse processes throughout the mammalian genome.

Poly(A) signal-dependent degradation of unprocessed nascent transcripts accompanies poly(A) signal-dependent transcriptional pausing in vitro

The poly(A) signal has long been known for its role in directing the cleavage and polyadenylation of eukaryotic mRNA. In recent years its additional coordinating role in multiple related aspects of gene expression has also become increasingly clear. Here we use HeLa nuclear extracts to study two of these activities, poly(A) signal-dependent transcriptional pausing, which was originally proposed as a surveillance checkpoint, and poly(A) signal-dependent degradation (PDD) of unprocessed transcripts from weak poly(A) signals. We confirm directly, by measuring the length of RNA within isolated transcription elongation complexes, that a newly transcribed poly(A) signal reduces the rate of elongation by RNA polymerase II and causes the accumulation of elongation complexes downstream from the poly(A) signal. We then show that if the RNA in these elongation complexes contains a functional but unprocessed poly(A) signal, degradation of the transcripts ensues. The degradation depends on the unprocessed poly(A) signal being functional, and does not occur if a mutant poly(A) signal is used. We suggest that during normal 3'-end processing the uncleaved poly(A) signal continuously samples competing reaction pathways for processing and for degradation, and that in the case of weak poly(A) signals, where poly(A) site cleavage is slow, the default pathway to degradation predominates.

Biased chromatin signatures around polyadenylation sites and exons

Core RNA-processing reactions in eukaryotic cells occur cotranscriptionally in a chromatin context, but the relationship between chromatin structure and pre-mRNA processing is poorly understood. We observed strong nucleosome depletion around human polyadenylation sites (PAS) and nucleosome enrichment just downstream of PAS. In genes with multiple alternative PAS, higher downstream nucleosome affinity was associated with higher PAS usage, independently of known PAS motifs that function at the RNA level. Conversely, exons were associated with distinct peaks in nucleosome density. Exons flanked by long introns or weak splice sites exhibited stronger nucleosome enrichment, and incorporation of nucleosome density data improved splicing simulation accuracy. Certain histone modifications, including H3K36me3 and H3K27me2, were specifically enriched on exons, suggesting active marking of exon locations at the chromatin level. Together, these findings provide evidence for extensive functional connections between chromatin structure and RNA processing.

The HBS1L-MYB intergenic interval associated with elevated HbF levels shows characteristics of a distal regulatory region in erythroid cells

HBS1L-MYB intergenic polymorphism (HMIP) on chromosome 6q23 is associated with elevated fetal hemoglobin levels and has pleiotropic effects on several hematologic parameters. To investigate potential regulatory activity in the region, we have measured sensitivity of the sequences to DNase I cleavage that identified 3 tissue-specific DNase I hypersensitive sites in the core intergenic interval. Chromatin immunoprecipitation with microarray (ChIP-chip) analysis showed strong histone acetylation in a defined interval of 65 kb corresponding to the core HBS1L-MYB intergenic region in primary human erythroid cells but not in non–MYB-expressing HeLa cells. ChIP-chip analysis also identified several potential cis-regulatory elements as strong GATA-1 signals that coincided with the DNase I hypersensitive sites present in MYB-expressing erythroid cells. We suggest that HMIP contains regulatory sequences that could be important in hematopoiesis by controlling MYB expression. This study provides the functional link between genetic association of HMIP with control of fetal hemoglobin and other hematologic parameters. We also present a large-scale analysis of histone acetylation as well as RNA polymerase II and GATA-1 interactions on chromosome 6q, and α and β globin gene loci. The data suggest that GATA-1 regulates numerous genes of various functions on chromosome 6q.

Transcription termination by nuclear RNA polymerases

Gene transcription in the cell nucleus is a complex and highly regulated process. Transcription in eukaryotes requires three distinct RNA polymerases, each of which employs its own mechanisms for initiation, elongation, and termination. Termination mechanisms vary considerably, ranging from relatively simple to exceptionally complex. In this review, we describe the present state of knowledge on how each of the three RNA polymerases terminates and how mechanisms are conserved, or vary, from yeast to human.

HBO1 HAT complexes target chromatin throughout gene coding regions via multiple PHD finger interactions with histone H3 tail

The HBO1 HAT protein is the major source of histone H4 acetylation in vivo and has been shown to play critical roles in gene regulation and DNA replication. A distinctive characteristic of HBO1 HAT complexes is the presence of three PHD finger domains in two different subunits: tumor suppressor proteins ING4/5 and JADE1/2/3. Biochemical and functional analyses indicate that these domains interact with histone H3 N-terminal tail region, but with a different specificity toward its methylation status. Their combinatorial action is essential in regulating chromatin binding and substrate specificity of HBO1 complexes, as well as cell growth. Importantly, localization analyses on the human genome indicate that HBO1 complexes are enriched throughout the coding regions of genes, supporting a role in transcription elongation. These results underline the importance and versatility of PHD finger domains in regulating chromatin association and histone modification crosstalk within a single protein complex.

Nascent RNA sequencing reveals widespread pausing and divergent initiation at human promoters

RNA polymerases are highly regulated molecular machines. We present a method (global run-on sequencing, GRO-seq) that maps the position, amount, and orientation of transcriptionally engaged RNA polymerases genome-wide. In this method, nuclear run-on RNA molecules are subjected to large-scale parallel sequencing and mapped to the genome. We show that peaks of promoter-proximal polymerase reside on approximately 30% of human genes, transcription extends beyond pre-messenger RNA 3' cleavage, and antisense transcription is prevalent. Additionally, most promoters have an engaged polymerase upstream and in an orientation opposite to the annotated gene. This divergent polymerase is associated with active genes but does not elongate effectively beyond the promoter. These results imply that the interplay between polymerases and regulators over broad promoter regions dictates the orientation and efficiency of productive transcription.

Alternative polyadenylation: a twist on mRNA 3' end formation

Regulation of gene expression by RNA processing mechanisms is now understood to be an important level of control in mammalian cells. Regulation at the level of RNA transcription, splicing, polyadenylation, nucleo-cytoplasmic transport, and translation into polypeptides has been well-studied. Alternative RNA processing events, such as alternative splicing, also have been recognized as key contributors to the complexity of mammalian gene expression. Pre-messenger RNAs (pre-mRNAs) may be polyadenylated in several different ways due to more than one polyadenylation signal, allowing a single gene to encode multiple mRNA transcripts. However, alternative polyadenylation has only recently taken the field as a major player in gene regulation. This review summarizes what is currently known about alternative polyadenylation. It covers results from bioinformatics, as well as those from investigations of viral and tissue-specific studies and, importantly, will set the stage for what is yet to come.