Characteristic bimodal profiles of RNA polymerase II at thousands of active mammalian promoters

Regulation of Pol II Pausing during Daily Gene Transcription in Mouse Liver

Cell autonomous circadian oscillation is present in central and various peripheral tissues. The intrinsic tissue clock and various extrinsic cues drive gene expression rhythms. Transcription regulation is thought to be the main driving force for gene rhythms. However, how transcription rhythms arise remains to be fully characterized due to the fact that transcription is regulated at multiple steps. In particular, Pol II recruitment, pause release, and premature transcription termination are critical regulatory steps that determine the status of Pol II pausing and transcription output near the transcription start site (TSS) of the promoter. Recently, we showed that Pol II pausing exhibits genome-wide changes during daily transcription in mouse liver. In this article, we review historical as well as recent findings on the regulation of transcription rhythms by the circadian clock and other transcription factors, and the potential limitations of those results in explaining rhythmic transcription at the TSS. We then discuss our results on the genome-wide characteristics of daily changes in Pol II pausing, the possible regulatory mechanisms involved, and their relevance to future research on circadian transcription regulation.

Premature transcription termination at the expanded GAA repeats and aberrant alternative polyadenylation contributes to the Frataxin transcriptional deficit in Friedreich's ataxia

Frataxin deficiency in Friedreich's ataxia results from transcriptional downregulation of the FXN gene caused by expansion of the intronic trinucleotide guanine-adenine-adenine (GAA) repeats. We used multiple transcriptomic approaches to determine the molecular mechanism of transcription inhibition caused by long GAAs. We uncovered that transcription of FXN in patient cells is prematurely terminated upstream of the expanded repeats leading to the formation of a novel, truncated and stable RNA. This FXN early terminated transcript (FXN-ett) undergoes alternative, non-productive splicing and does not contribute to the synthesis of functional frataxin. The level the FXN-ett RNA directly correlates with the length of the longer of the two expanded GAA tracts. Targeting GAAs with antisense oligonucleotides or excision of the repeats eliminates the transcription impediment, diminishes expression of the aberrant FXN-ett, while increasing levels of FXN mRNA and frataxin. Non-productive transcription may represent a common phenomenon and attractive therapeutic target in diseases caused by repeat-mediated transcription aberrations.

Molecular Insights into the Recognition of Acetylated Histone Modifications by the BRPF2 Bromodomain

HBO1 [HAT bound to the origin recognition complex (ORC)], a member of the MYST family of histone acetyltransferases (HATs), was initially identified as a binding partner of ORC that acetylates free histone H3, H4, and nucleosomal H3. It functions as a quaternary complex with the BRPF (BRPF1/2/3) scaffolding protein and two accessory proteins, ING4/5 and Eaf6. Interaction of BRPF2 with HBO1 has been shown to be important for regulating H3K14 acetylation during embryonic development. However, how BRPF2 directs the HBO1 HAT complex to chromatin to regulate its HAT activity toward nucleosomal substrates remains unclear. Our findings reveal novel interacting partners of the BRPF2 bromodomain that recognizes different acetyllysine residues on the N-terminus of histone H4, H3, and H2A and preferentially binds to H4K5ac, H4K8ac, and H4K5acK12ac modifications. In addition, mutational analysis of the BRPF2 bromodomain coupled with isothermal titration calorimetry binding and pull-down assays on the histone substrates identified critical residues responsible for acetyllysine binding. Moreover, the BRPF2 bromodomain could enrich H4K5ac mark-bearing mononucleosomes compared to other acetylated H4 marks. Consistent with this, ChIP-seq analysis revealed that BRPF2 strongly co-localizes with HBO1 at histone H4K5ac and H4K8ac marks near the transcription start sites in the genome. Our study provides novel insights into how the histone binding function of the BRPF2 bromodomain directs the recruitment of the HBO1 HAT complex to chromatin to regulate gene expression.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

Corticosterone pattern-dependent glucocorticoid receptor binding and transcriptional regulation within the liver

Ultradian glucocorticoid rhythms are highly conserved across mammalian species, however, their functional significance is not yet fully understood. Here we demonstrate that pulsatile corticosterone replacement in adrenalectomised rats induces a dynamic pattern of glucocorticoid receptor (GR) binding at ~3,000 genomic sites in liver at the pulse peak, subsequently not found during the pulse nadir. In contrast, constant corticosterone replacement induced prolonged binding at the majority of these sites. Additionally, each pattern further induced markedly different transcriptional responses. During pulsatile treatment, intragenic occupancy by active RNA polymerase II exhibited pulsatile dynamics with transient changes in enrichment, either decreased or increased depending on the gene, which mostly returned to baseline during the inter-pulse interval. In contrast, constant corticosterone exposure induced prolonged effects on RNA polymerase II occupancy at the majority of gene targets, thus acting as a sustained regulatory signal for both transactivation and repression of glucocorticoid target genes. The nett effect of these differences were consequently seen in the liver transcriptome as RNA-seq analysis indicated that despite the same overall amount of corticosterone infused, twice the number of transcripts were regulated by constant corticosterone infusion, when compared to pulsatile. Target genes that were found to be differentially regulated in a pattern-dependent manner were enriched in functional pathways including carbohydrate, cholesterol, glucose and fat metabolism as well as inflammation, suggesting a functional role for dysregulated glucocorticoid rhythms in the development of metabolic dysfunction.

Adrenal glucocorticoid hormones are released in a characteristic ultradian rhythm that becomes dysregulated during chronic stress, disease, or synthetic corticosteroid treatment. Metabolic dysfunction is a comorbidity associated with all these conditions, but the role that altered glucocorticoid dynamics play is unknown. As the liver is a major site of glucocorticoid action on metabolic homeostasis regulated by the glucocorticoid receptor, we have assessed how different patterns of hormone replacement in adrenalectomised rats differentially regulate gene pathways involved in type II diabetes, cirrhosis, and fatty liver development, via altering the pattern of glucocorticoid receptor binding to regulatory sites. We believe our findings have important implications for therapies that can reproduce the endogenous glucocorticoid rhythm and thus minimize adverse metabolic side-effects in patients.

Characterizing Genetic Regulatory Elements in Ovine Tissues

The Ovine Functional Annotation of Animal Genomes (FAANG) project, part of the broader livestock species FAANG initiative, aims to identify and characterize gene regulatory elements in domestic sheep. Regulatory element annotation is essential for identifying genetic variants that affect health and production traits in this important agricultural species, as greater than 90% of variants underlying genetic effects are estimated to lie outside of transcribed regions. Histone modifications that distinguish active or repressed chromatin states, CTCF binding, and DNA methylation were used to characterize regulatory elements in liver, spleen, and cerebellum tissues from four yearling sheep. Chromatin immunoprecipitation with sequencing (ChIP-seq) was performed for H3K4me3, H3K27ac, H3K4me1, H3K27me3, and CTCF. Nine chromatin states including active promoters, active enhancers, poised enhancers, repressed enhancers, and insulators were characterized in each tissue using ChromHMM. Whole-genome bisulfite sequencing (WGBS) was performed to determine the complement of whole-genome DNA methylation with the ChIP-seq data. Hypermethylated and hypomethylated regions were identified across tissues, and these locations were compared with chromatin states to better distinguish and validate regulatory elements in these tissues. Interestingly, chromatin states with the poised enhancer mark H3K4me1 in the spleen and cerebellum and CTCF in the liver displayed the greatest number of hypermethylated sites. Not surprisingly, active enhancers in the liver and spleen, and promoters in the cerebellum, displayed the greatest number of hypomethylated sites. Overall, chromatin states defined by histone marks and CTCF occupied approximately 22% of the genome in all three tissues. Furthermore, the liver and spleen displayed in common the greatest percent of active promoter (65%) and active enhancer (81%) states, and the liver and cerebellum displayed in common the greatest percent of poised enhancer (53%), repressed enhancer (68%), hypermethylated sites (75%), and hypomethylated sites (73%). In addition, both known and de novo CTCF-binding motifs were identified in all three tissues, with the highest number of unique motifs identified in the cerebellum. In summary, this study has identified the regulatory regions of genes in three tissues that play key roles in defining health and economically important traits and has set the precedent for the characterization of regulatory elements in ovine tissues using the Rambouillet reference genome.

LETR1 is a lymphatic endothelial-specific lncRNA governing cell proliferation and migration through KLF4 and SEMA3C

Recent studies have revealed the importance of long noncoding RNAs (lncRNAs) as tissue-specific regulators of gene expression. There is ample evidence that distinct types of vasculature undergo tight transcriptional control to preserve their structure, identity, and functions. We determine a comprehensive map of lineage-specific lncRNAs in human dermal lymphatic and blood vascular endothelial cells (LECs and BECs), combining RNA-Seq and CAGE-Seq. Subsequent antisense oligonucleotide-knockdown transcriptomic profiling of two LEC- and two BEC-specific lncRNAs identifies LETR1 as a critical gatekeeper of the global LEC transcriptome. Deep RNA-DNA, RNA-protein interaction studies, and phenotype rescue analyses reveal that LETR1 is a nuclear trans-acting lncRNA modulating, via key epigenetic factors, the expression of essential target genes, including KLF4 and SEMA3C, governing the growth and migratory ability of LECs. Together, our study provides several lines of evidence supporting the intriguing concept that every cell type expresses precise lncRNA signatures to control lineage-specific regulatory programs.

Looping of upstream cis-regulatory elements is required for CFTR expression in human airway epithelial cells

The CFTR gene lies within an invariant topologically associated domain (TAD) demarcated by CTCF and cohesin, but shows cell-type specific control mechanisms utilizing different cis-regulatory elements (CRE) within the TAD. Within the respiratory epithelium, more than one cell type expresses CFTR and the molecular mechanisms controlling its transcription are likely divergent between them. Here, we determine how two extragenic CREs that are prominent in epithelial cells in the lung, regulate expression of the gene. We showed earlier that these CREs, located at -44 and -35 kb upstream of the promoter, have strong cell-type-selective enhancer function. They are also responsive to inflammatory mediators and to oxidative stress, consistent with a key role in CF lung disease. Here, we use CRISPR/Cas9 technology to remove these CREs from the endogenous locus in human bronchial epithelial cells. Loss of either site extinguished CFTR expression and abolished long-range interactions between these sites and the gene promoter, suggesting non-redundant enhancers. The deletions also greatly reduced promoter interactions with the 5' TAD boundary. We show substantial recruitment of RNAPII to the -35 kb element and identify CEBPβ as a key activator of airway expression of CFTR, likely through occupancy at this CRE and the gene promoter.

BRD4 bimodal binding at promoters and drug-induced displacement at Pol II pause sites associates with I-BET sensitivity

Background

Deregulated transcription is a major driver of diseases such as cancer. Bromodomain and extra-terminal (BET) proteins (BRD2, BRD3, BRD4 and BRDT) are chromatin readers essential for maintaining proper gene transcription by specifically binding acetylated lysine residues. Targeted displacement of BET proteins from chromatin, using BET inhibitors (I-BETs), is a promising therapy, especially for acute myeloid leukemia (AML), and evaluation of resistance mechanisms is necessary to optimize the clinical efficacy of these drugs.

Results

To uncover mechanisms of intrinsic I-BET resistance, we quantified chromatin binding and displacement for BRD2, BRD3 and BRD4 after dose response treatment with I-BET151, in sensitive and resistant in vitro models of leukemia, and mapped BET proteins/I-BET interactions genome wide using antibody- and compound-affinity capture methods followed by deep sequencing. The genome-wide map of BET proteins sensitivity to I-BET revealed a bimodal pattern of binding flanking transcription start sites (TSSs), in which drug-mediated displacement from chromatin primarily affects BRD4 downstream of the TSS and prolongs the pausing of RNA Pol II. Correlation of BRD4 binding and drug-mediated displacement at RNA Pol II pause sites with gene expression revealed a differential behavior of sensitive and resistant tumor cells to I-BET and identified a BRD4 signature at promoters of sensitive coding and non-coding genes.

Conclusions

We provide evidence that I-BET-induced shift of Pol II pausing at promoters via displacement of BRD4 is a determinant of intrinsic I-BET sensitivity. This finding may guide pharmacological treatment to enhance the clinical utility of such targeted therapies in AML and potentially other BET proteins-driven diseases.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0286-5) contains supplementary material, which is available to authorized users.

Genome-wide analysis reveals that exon methylation facilitates its selective usage in the human transcriptome

DNA methylation, especially in promoter regions, is a well-characterized epigenetic marker related to gene expression regulation in eukaryotes. However, the role of intragenic DNA methylation in the usage of corresponding exons still remains elusive. In this study, we described the DNA methylome across 10 human tissues. The human genome showed both conserved and varied methylation levels among these tissues. We found that the methylation densities in promoters and first exons were negatively correlated with the corresponding gene expression level. Nevertheless, the methylation densities within introns, internal exons and down 1 kb regions showed weak correlation with gene expression levels. Importantly, we observed a remarkably positive relationship between methylation density and exon expression level of intragenic exons. Notably, skip-in exons were much more methylated than skip-out exons. We also identified 260 exons that showed both differential methylation levels and differential expression levels in lung cancer. Genes harboring these differentially regulated exons were significantly enriched in the cancer hallmark-related biological process. Moreover, a 10-exon signature was identified as a promising prognostic predictor for lung cancer. Our study illuminates the DNA methylome, describes its relationship with gene expression across human tissues and provides new insights into intragenic DNA methylation and exon usage during the transcriptional/alternative splicing process and in cancer.

Metastasis Suppressor NME1 Directly Activates Transcription of the ALDOC Gene in Melanoma Cells

BACKGROUND/AIM

NME/NM23 nucleoside diphosphate kinase 1 (NME1) is a metastasis suppressor gene, exhibiting reduced expression in metastatic cancers and the ability to suppress metastatic activity of cancer cells. We previously identified NME1-regulated genes with prognostic value in human melanoma. This study was conducted in melanoma cell lines aiming to elucidate the mechanism through which NME regulates one of these genes, aldolase C (ALDOC).

MATERIALS AND METHODS

ALDOC mRNA and protein expression was measured using qRT-PCR and immunoblot analyses. Promoter-luciferase constructs and chromatin immunoprecipitation were employed to measure the impact of NME1 on ALDOC transcription.

RESULTS

NME1 enhanced ALDOC transcription, evidenced by increased expression of ALDOC pre-mRNA and activity of an ALDOC promoter-luciferase module. NME1 was detected at the ALDOC promoter, and forced NME1 expression resulted in enhanced occupancy of the promoter by NME1, increased presence of epigenetic activation markers (H3K4me3 and H3K27ac), and recruitment of RNA polymerase II.

CONCLUSION

This is the first study to indicate that NME1 induces transcription through its direct binding to the promoter region of a target gene.

Dynamic turnover of paused Pol II complexes at human promoters

Paused RNA polymerase II (Pol II) that piles up near most human promoters is the target of mechanisms that control entry into productive elongation. Whether paused Pol II is a stable or dynamic target remains unresolved. We report that most 5' paused Pol II throughout the genome is turned over within 2 min. This process is revealed under hypertonic conditions that prevent Pol II recruitment to promoters. This turnover requires cell viability but is not prevented by inhibiting transcription elongation, suggesting that it is mediated at the level of termination. When initiation was prevented by triptolide during recovery from high salt, a novel preinitiated state of Pol II lacking the pausing factor Spt5 accumulated at transcription start sites. We propose that Pol II occupancy near 5' ends is governed by a cycle of ongoing assembly of preinitiated complexes that transition to pause sites followed by eviction from the DNA template. This model suggests that mechanisms regulating the transition to productive elongation at pause sites operate on a dynamic population of Pol II that is turning over at rates far higher than previously suspected. We suggest that a plausible alternative to elongation control via escape from a stable pause is by escape from premature termination.

Rapid Response and Slow Recovery of the H3K4me3 Epigenomic Marker in the Liver after Light-mediated Phase Advances of the Circadian Clock

Mammalian tissues display circadian rhythms in transcription, translation, and histone modifications. Here we asked how an advance of the light-dark cycle alters daily rhythms in the liver epigenome at the H3K4me3 (trimethylation of lysine 4 on histone 3) modification, which is found at active and poised gene promoters. H3K4me3 levels were first measured at 4 time points (zeitgeber time [ZT] 3, 8, 15, and 20) during a normal 12L:12D light-dark cycle. Peak levels were observed during the early dark phase at ZT15 and dropped to low levels around lights-on (ZT0) between ZT20 and ZT3. A 6-h phase advance at ZT18 (new lights-on after only 6 h of darkness) led to a transient extension of peak H3K4me3 levels. Although locomotor activity reentrained within a week after the phase advance, H3K4me3 rhythms failed to do so, with peak levels remaining in the light phase at the 1-week recovery time point. Eight weekly phase advances, with 1-week recovery times between each phase advance, further disrupted the H3K4me3 rhythms. Finally, we used the mPer2^Luc knockin mouse to determine whether the phase advance also disrupted Per2 protein expression. Similar to the results from the histone work, we found both a rapid response to the phase advance and a delayed recovery, the latter in sync with H3K4me3 levels. A model to explain these results is offered.

Identification of active miRNA promoters from nuclear run-on RNA sequencing

The genome-wide identification of microRNA transcription start sites (miRNA TSSs) is essential for understanding how miRNAs are regulated in development and disease. In this study, we developed mirSTP (mirna transcription Start sites Tracking Program), a probabilistic model for identifying active miRNA TSSs from nascent transcriptomes generated by global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq). MirSTP takes advantage of characteristic bidirectional transcription signatures at active TSSs in GRO/PRO-seq data, and provides accurate TSS prediction for human intergenic miRNAs at a high resolution. MirSTP performed better than existing generalized and experiment specific methods, in terms of the enrichment of various promoter-associated marks. MirSTP analysis of 27 human cell lines in 183 GRO-seq and 28 PRO-seq experiments identified TSSs for 480 intergenic miRNAs, indicating a wide usage of alternative TSSs. By integrating predicted miRNA TSSs with matched ENCODE transcription factor (TF) ChIP-seq data, we connected miRNAs into the transcriptional circuitry, which provides a valuable source for understanding the complex interplay between TF and miRNA. With mirSTP, we not only predicted TSSs for 72 miRNAs, but also identified 12 primary miRNAs with significant RNA polymerase pausing alterations after JQ1 treatment; each miRNA was further validated through BRD4 binding to its predicted promoter. MirSTP is available at http://bioinfo.vanderbilt.edu/mirSTP/.

Retrieving Chromatin Patterns from Deep Sequencing Data Using Correlation Functions

Epigenetic modifications and other chromatin features partition the genome on multiple length scales. They define chromatin domains with distinct biological functions that come in sizes ranging from single modified DNA bases to several megabases in the case of heterochromatic histone modifications. Due to chromatin folding, domains that are well separated along the linear nucleosome chain can form long-range interactions in three-dimensional space. It has now become a routine task to map epigenetic marks and chromatin structure by deep sequencing methods. However, assessing and comparing the properties of chromatin domains and their positional relationships across data sets without a priori assumptions remains challenging. Here, we introduce multiscale correlation evaluation (MCORE), which uses the fluctuation spectrum of mapped sequencing reads to quantify and compare chromatin patterns over a broad range of length scales in a model-independent manner. We applied MCORE to map the chromatin landscape in mouse embryonic stem cells and differentiated neural cells. We integrated sequencing data from chromatin immunoprecipitation, RNA expression, DNA methylation, and chromosome conformation capture experiments into network models that reflect the positional relationships among these features on different genomic scales. Furthermore, we used MCORE to compare our experimental data to models for heterochromatin reorganization during differentiation. The application of correlation functions to deep sequencing data complements current evaluation schemes and will support the development of quantitative descriptions of chromatin networks.

PHF13: A new player involved in RNA polymerase II transcriptional regulation and co-transcriptional splicing

We recently identified PHF13 as an H3K4me2/3 chromatin reader and transcriptional co-regulator. We found that PHF13 interacts with RNAPIIS5P and PRC2 stabilizing their association with active and bivalent promoters. Furthermore, mass spectrometry analysis identified ∼50 spliceosomal proteins in PHF13s interactome. Here, we will discuss the potential role of PHF13 in RNAPII pausing and co-transcriptional splicing.

Regulatory cis-elements are located in accessible promoter regions of the CAT2 promoter and affect activating histone modifications in Arabidopsis thaliana

Catalase 2 (CAT2) plays an important role in the detoxification of hydrogen peroxide released either during photorespiration or as a consequence of biotic and abiotic stress as well as in the initiation of senescence. To date, our understanding of the regulation of CAT2 gene expression is rather poor. Chromatin immunoprecipitation experiments revealed that a wide region of the CAT2 promoter is nucleosome depleted, reflecting the ability to rapidly respond to changing environmental and stress conditions and, thus, adjusting the transcript levels of CAT2. The lowest nucleosome density was found in the region of -900 bp relative to the transcription initiation start (TIS) where two regulatory elements are located. The distance of the nucleosome depleted region to the TIS is quite unusual because the majority of nucleosome free regions are generally located in close vicinity to the 5' untranslated region. The analysis of transgenic 5' upstream deletion::gusA Arabidopsis lines showed that this region is important for the regulation of CAT2 promoter activity. To evaluate the function of the two motifs, the contribution of each element to CAT2 promoter activity was analyzed by site directed mutagenesis. The data revealed that the CAT2 promoter is regulated by the ACGT motif (Box2) rather than by the G-Box binding motif (Box1) in the vegetative phase of development. Furthermore, the presence of both Box1 and Box2 positively affected the abundance of activating histone modifications.

Genome-wide uniformity of human 'open' pre-initiation complexes

Transcription of protein-coding and noncoding DNA occurs pervasively throughout the mammalian genome. Their sites of initiation are generally inferred from transcript 5' ends and are thought to be either locally dispersed or focused. How these two modes of initiation relate is unclear. Here, we apply permanganate treatment and chromatin immunoprecipitation (PIP-seq) of initiation factors to identify the precise location of melted DNA separately associated with the preinitiation complex (PIC) and the adjacent paused complex (PC). This approach revealed the two known modes of transcription initiation. However, in contrast to prevailing views, they co-occurred within the same promoter region: initiation originating from a focused PIC, and broad nucleosome-linked initiation. PIP-seq allowed transcriptional orientation of Pol II to be determined, which may be useful near promoters where sufficient sense/anti-sense transcript mapping information is lacking. PIP-seq detected divergently oriented Pol II at both coding and noncoding promoters, as well as at enhancers. Their occupancy levels were not necessarily coupled in the two orientations. DNA sequence and shape analysis of initiation complex sites suggest that both sequence and shape contribute to specificity, but in a context-restricted manner. That is, initiation sites have the locally "best" initiator (INR) sequence and/or shape. These findings reveal a common core to pervasive Pol II initiation throughout the human genome.

Linking Core Promoter Classes to Circadian Transcription

Circadian rhythms in transcription are generated by rhythmic abundances and DNA binding activities of transcription factors. Propagation of rhythms to transcriptional initiation involves the core promoter, its chromatin state, and the basal transcription machinery. Here, I characterize core promoters and chromatin states of genes transcribed in a circadian manner in mouse liver and in Drosophila. It is shown that the core promoter is a critical determinant of circadian mRNA expression in both species. A distinct core promoter class, strong circadian promoters (SCPs), is identified in mouse liver but not Drosophila. SCPs are defined by specific core promoter features, and are shown to drive circadian transcriptional activities with both high averages and high amplitudes. Data analysis and mathematical modeling further provided evidence for rhythmic regulation of both polymerase II recruitment and pause release at SCPs. The analysis provides a comprehensive and systematic view of core promoters and their link to circadian mRNA expression in mouse and Drosophila, and thus reveals a crucial role for the core promoter in regulated, dynamic transcription.

Methylation of RNA polymerase II non-consensus Lysine residues marks early transcription in mammalian cells

Dynamic post-translational modification of RNA polymerase II (RNAPII) coordinates the co-transcriptional recruitment of enzymatic complexes that regulate chromatin states and processing of nascent RNA. Extensive phosphorylation of serine residues at the largest RNAPII subunit occurs at its structurally-disordered C-terminal domain (CTD), which is composed of multiple heptapeptide repeats with consensus sequence Y1-S2-P3-T4-S5-P6-S7. Serine-5 and Serine-7 phosphorylation mark transcription initiation, whereas Serine-2 phosphorylation coincides with productive elongation. In vertebrates, the CTD has eight non-canonical substitutions of Serine-7 into Lysine-7, which can be acetylated (K7ac). Here, we describe mono- and di-methylation of CTD Lysine-7 residues (K7me1 and K7me2). K7me1 and K7me2 are observed during the earliest transcription stages and precede or accompany Serine-5 and Serine-7 phosphorylation. In contrast, K7ac is associated with RNAPII elongation, Serine-2 phosphorylation and mRNA expression. We identify an unexpected balance between RNAPII K7 methylation and acetylation at gene promoters, which fine-tunes gene expression levels.

Uncovering drug-responsive regulatory elements

Nucleotide changes in gene regulatory elements can have a major effect on interindividual differences in drug response. For example, by reviewing all published pharmacogenomic genome-wide association studies, we show here that 96.4% of the associated single nucleotide polymorphisms reside in noncoding regions. We discuss how sequencing technologies are improving our ability to identify drug response-associated regulatory elements genome-wide and to annotate nucleotide variants within them. We highlight specific examples of how nucleotide changes in these elements can affect drug response and illustrate the techniques used to find them and functionally characterize them. Finally, we also discuss challenges in the field of drug-responsive regulatory elements that need to be considered in order to translate these findings into the clinic.

Analysis of paired end Pol II ChIP-seq and short capped RNA-seq in MCF-7 cells

While a role of promoter-proximal RNA Polymerase II (Pol II) pausing in regulation of eukaryotic gene expression is implied, the mechanisms and dynamics of this process are poorly understood. We performed genome-wide analysis of short capped RNAs (scRNAs) and Pol II chromatin immunoprecipitation sequencing (ChIP-seq) in human breast cancer MCF-7 cells to better understand Pol II pausing (Samarakkody, A., Abbas, A., Scheidegger, A., Warns, J., Nnoli, O., Jokinen, B., Zarns, K., Kubat, B., Dhasarathy, A. and Nechaev, S. (2015) RNA polymerase II pausing can be retained or acquired during activation of genes involved in the epithelial to mesenchymal transition. Nucleic Acids Res43, 3938–3949). The data are available at the NCBI Gene Expression Omnibus under accession number GSE67041. For both ChIP and scRNA samples, we used paired end sequencing on the Illumina MiSeq instrument. For ChIP-seq, the use of paired end sequencing allowed us to avoid ambiguities in center-read definition. For scRNA seq, this allowed us to identify both the 5′-end and the 3′-end in the same run that represent, respectively, the transcription start sites and the locations of Pol II pausing. The sharpening of Pol II ChIP-seq metagene profiles when aligned against 5′-ends of scRNAs indicates that these RNAs can be used to define the start sites for the majority of mRNA transcription events.

The core promoter: At the heart of gene expression

The identities of different cells and tissues in multicellular organisms are determined by tightly controlled transcriptional programs that enable accurate gene expression. The mechanisms that regulate gene expression comprise diverse multiplayer molecular circuits of multiple dedicated components. The RNA polymerase II (Pol II) core promoter establishes the center of this spatiotemporally orchestrated molecular machine. Here, we discuss transcription initiation, diversity in core promoter composition, interactions of the basal transcription machinery with the core promoter, enhancer-promoter specificity, core promoter-preferential activation, enhancer RNAs, Pol II pausing, transcription termination, Pol II recycling and translation. We further discuss recent findings indicating that promoters and enhancers share similar features and may not substantially differ from each other, as previously assumed. Taken together, we review a broad spectrum of studies that highlight the importance of the core promoter and its pivotal role in the regulation of metazoan gene expression and suggest future research directions and challenges.

RNA polymerase II pausing can be retained or acquired during activation of genes involved in the epithelial to mesenchymal transition

Promoter-proximal RNA polymerase II (Pol II) pausing is implicated in the regulation of gene transcription. However, the mechanisms of pausing including its dynamics during transcriptional responses remain to be fully understood. We performed global analysis of short capped RNAs and Pol II Chromatin Immunoprecipitation sequencing in MCF-7 breast cancer cells to map Pol II pausing across the genome, and used permanganate footprinting to specifically follow pausing during transcriptional activation of several genes involved in the epithelial to mesenchymal transition (EMT). We find that the gene for EMT master regulator Snail (SNAI1), but not Slug (SNAI2), shows evidence of Pol II pausing before activation. Transcriptional activation of the paused SNAI1 gene is accompanied by a further increase in Pol II pausing signal, whereas activation of non-paused SNAI2 gene results in the acquisition of a typical pausing signature. The increase in pausing signal reflects increased transcription initiation without changes in Pol II pausing. Activation of the heat shock HSP70 gene involves pausing release that speeds up Pol II turnover, but does not change pausing location. We suggest that Pol II pausing is retained during transcriptional activation and can further undergo regulated release in a signal-specific manner.