Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data

<p>Abstract</p> <p>Background</p> <p>Initiation and elongation of RNA polymerase II (RNAPII) transcription is regulated by both DNA sequence and chromatin signals. Recent breakthroughs make it possible to measure the chromatin state and activity of core promoters genome-wide, but dedicated computational strategies are needed to progress from descriptive annotation of data to quantitative, predictive models.</p> <p>Results</p> <p>Here, we describe a computational framework which with high accuracy can predict the locations of core promoters, the amount of recruited RNAPII at the promoter, the amount of elongating RNAPII in the gene body, the mRNA production originating from the promoter and finally also the stalling characteristics of RNAPII by considering both quantitative and spatial features of histone modifications around the transcription start site (TSS).</p> <p>As the model framework can also pinpoint the signals that are the most influential for prediction, it can be used to infer underlying regulatory biology. For example, we show that the H3K4 di- and tri- methylation signals are strongly predictive for promoter location while the acetylation marks H3K9 and H3K27 are highly important in estimating the promoter usage. All of these four marks are found to be necessary for recruitment of RNAPII but not sufficient for the elongation. We also show that the spatial distributions of histone marks are almost as predictive as the signal strength and that a set of histone marks immediately downstream of the TSS is highly predictive of RNAPII stalling.</p> <p>Conclusions</p> <p>In this study we introduce a general framework to accurately predict the level of RNAPII recruitment, elongation, stalling and mRNA expression from chromatin signals. The versatility of the method also makes it ideally suited to investigate other genomic data.</p

Transcriptional Enhancers in Protein-Coding Exons of Vertebrate Developmental Genes

Many conserved noncoding sequences function as transcriptional enhancers that regulate gene expression. Here, we report that protein-coding DNA also frequently contains enhancers functioning at the transcriptional level. We tested the enhancer activity of 31 protein-coding exons, which we chose based on strong sequence conservation between zebrafish and human, and occurrence in developmental genes, using a Tol2 transposable GFP reporter assay in zebrafish. For each exon we measured GFP expression in hundreds of embryos in 10 anatomies via a novel system that implements the voice-recognition capabilities of a cellular phone. We find that 24/31 (77%) exons drive GFP expression compared to a minimal promoter control, and 14/24 are anatomy-specific (expression in four anatomies or less). GFP expression driven by these coding enhancers frequently overlaps the anatomies where the host gene is expressed (60%), suggesting self-regulation. Highly conserved coding sequences and highly conserved noncoding sequences do not significantly differ in enhancer activity (coding: 24/31 vs. noncoding: 105/147) or tissue-specificity (coding: 14/24 vs. noncoding: 50/105). Furthermore, coding and noncoding enhancers display similar levels of the enhancer-related histone modification H3K4me1 (coding: 9/24 vs noncoding: 34/81). Meanwhile, coding enhancers are over three times as likely to contain an H3K4me1 mark as other exons of the host gene. Our work suggests that developmental transcriptional enhancers do not discriminate between coding and noncoding DNA and reveals widespread dual functions in protein-coding DNA

Regulation of Cyclooxygenase-2 Expression by Heat: A Novel Aspect of Heat Shock Factor 1 Function in Human Cells

The heat-shock response, a fundamental defense mechanism against proteotoxic stress, is regulated by a family of heat-shock transcription factors (HSF). In humans HSF1 is considered the central regulator of heat-induced transcriptional responses. The main targets for HSF1 are specific promoter elements (HSE) located upstream of heat-shock genes encoding cytoprotective heat-shock proteins (HSP) with chaperone function. In addition to its cytoprotective function, HSF1 was recently hypothesized to play a more complex role, regulating the expression of non-HSP genes; however, the non-canonical role of HSF1 is still poorly understood. Herein we report that heat-stress promotes the expression of cyclooxygenase-2 (COX-2), a key regulator of inflammation controlling prostanoid and thromboxane synthesis, resulting in the production of high levels of prostaglandin-E2 in human cells. We show that heat-induced COX-2 expression is regulated at the transcriptional level via HSF1-mediated signaling and identify, by in-vitro reporter gene activity assay and deletion-mutant constructs analysis, the COX-2 heat-responsive promoter region and a new distal cis-acting HSE located at position −2495 from the transcription start site. As shown by ChIP analysis, HSF1 is recruited to the COX-2 promoter rapidly after heat treatment; by using shRNA-mediated HSF1 suppression and HSE-deletion from the COX-2 promoter, we demonstrate that HSF1 plays a central role in the transcriptional control of COX-2 by heat. Finally, COX-2 transcription is also induced at febrile temperatures in endothelial cells, suggesting that HSF1-dependent COX-2 expression could contribute to increasing blood prostaglandin levels during fever. The results identify COX-2 as a human non-classical heat-responsive gene, unveiling a new aspect of HSF1 function

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

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

Chromosomal organization at the level of gene complexes

Metazoan genomes primarily consist of non-coding DNA in comparison to coding regions. Non-coding fraction of the genome contains cis-regulatory elements, which ensure that the genetic code is read properly at the right time and space during development. Regulatory elements and their target genes define functional landscapes within the genome, and some developmentally important genes evolve by keeping the genes involved in specification of common organs/tissues in clusters and are termed gene complex. The clustering of genes involved in a common function may help in robust spatio-temporal gene expression. Gene complexes are often found to be evolutionarily conserved, and the classic example is the hox complex. The evolutionary constraints seen among gene complexes provide an ideal model system to understand cis and trans-regulation of gene function. This review will discuss the various characteristics of gene regulatory modules found within gene complexes and how they can be characterized

Identification of active transcriptional regulatory elements from GRO-seq data

Modifications to the global run-on and sequencing (GRO-seq) protocol that enrich for 5'-capped RNAs can be used to reveal active transcriptional regulatory elements (TREs) with high accuracy. Here, we introduce discriminative regulatory-element detection from GRO-seq (dREG), a sensitive machine learning method that uses support vector regression to identify active TREs from GRO-seq data without requiring cap-based enrichment (https://github.com/Danko-Lab/dREG/). This approach allows TREs to be assayed together with gene expression levels and other transcriptional features in a single experiment. Predicted TREs are more enriched for several marks of transcriptional activation-including expression quantitative trait loci, disease-associated polymorphisms, acetylated histone 3 lysine 27 (H3K27ac) and transcription factor binding-than those identified by alternative functional assays. Using dREG, we surveyed TREs in eight human cell types and provide new insights into global patterns of TRE function

Dynamics of promoter bivalency and RNAP II pausing in mouse stem and differentiated cells

Mammalian embryonic stem cells display a unique epigenetic and transcriptional state to facilitate pluripotency by maintaining lineage-specification genes in a poised state. Two epigenetic and transcription processes involved in maintaining poised state are bivalent chromatin, characterized by the simultaneous presence of activating and repressive histone methylation marks, and RNA polymerase II (RNAPII) promoter proximal pausing. However, the dynamics of histone modifications and RNAPII at promoters in diverse cellular contexts remains underexplored. We collected genome wide data for bivalent chromatin marks H3K4me3 and H3K27me3, and RNAPII (8WG16) occupancy together with expression profiling in eight different cell types, including ESCs, in mouse. The epigenetic and transcription profiles at promoters grouped in over thirty clusters with distinct functional identities and transcription control. The clustering analysis identified distinct bivalent clusters where genes in one cluster retained bivalency across cell types while in the other were mostly cell type specific, but neither showed a high RNAPII pausing. We noted that RNAPII pausing is more associated with active genes than bivalent genes in a cell type, and was globally reduced in differentiated cell types compared to multipotent

Venous identity requires BMP signalling through ALK3

Venous endothelial cells are molecularly and functionally distinct from their arterial counterparts. Although veins are often considered the default endothelial state, genetic manipulations can modulate both acquisition and loss of venous fate, suggesting that venous identity is the result of active transcriptional regulation. However, little is known about this process. Here we show that BMP signalling controls venous identity via the ALK3/BMPR1A receptor and SMAD1/SMAD5. Perturbations to TGF-β and BMP signalling in mice and zebrafish result in aberrant vein formation and loss of expression of the venous-specific gene Ephb4, with no effect on arterial identity. Analysis of a venous endothelium-specific enhancer for Ephb4 shows enriched binding of SMAD1/5 and a requirement for SMAD binding motifs. Further, our results demonstrate that BMP/SMAD-mediated Ephb4 expression requires the venous-enriched BMP type I receptor ALK3/BMPR1A. Together, our analysis demonstrates a requirement for BMP signalling in the establishment of Ephb4 expression and the venous vasculature

Histone modification profiles are predictive for tissue/cell-type specific expression of both protein-coding and microRNA genes

<p>Abstract</p> <p>Background</p> <p>Gene expression is regulated at both the DNA sequence level and through modification of chromatin. However, the effect of chromatin on tissue/cell-type specific gene regulation (TCSR) is largely unknown. In this paper, we present a method to elucidate the relationship between histone modification/variation (HMV) and TCSR.</p> <p>Results</p> <p>A classifier for differentiating CD4+ T cell-specific genes from housekeeping genes using HMV data was built. We found HMV in both promoter and gene body regions to be predictive of genes which are targets of TCSR. For example, the histone modification types H3K4me3 and H3K27ac were identified as the most predictive for CpG-related promoters, whereas H3K4me3 and H3K79me3 were the most predictive for nonCpG-related promoters. However, genes targeted by TCSR can be predicted using other type of HMVs as well. Such redundancy implies that multiple type of underlying regulatory elements, such as enhancers or intragenic alternative promoters, which can regulate gene expression in a tissue/cell-type specific fashion, may be marked by the HMVs. Finally, we show that the predictive power of HMV for TCSR is not limited to protein-coding genes in CD4+ T cells, as we successfully predicted TCSR targeted genes in muscle cells, as well as microRNA genes with expression specific to CD4+ T cells, by the same classifier which was trained on HMV data of protein-coding genes in CD4+ T cells.</p> <p>Conclusion</p> <p>We have begun to understand the HMV patterns that guide gene expression in both tissue/cell-type specific and ubiquitous manner.</p