Cis-regulatory control of transcriptional timing and noise in response to estrogen.

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Psi promotes Drosophila wing growth via direct transcriptional activation of cell cycle targets and repression of growth inhibitors.

The first characterised FUSE Binding Protein family member, FUBP1, binds single-stranded DNA to activate MYC transcription. Psi, the sole FUBP protein in Drosophila, binds RNA to regulate P-element and mRNA splicing. Our previous work revealed pro-growth functions for Psi, which depend, in part, on transcriptional activation of Myc. Genome-wide functions for FUBP family proteins in transcriptional control remain obscure. Here, through the first genome-wide binding and expression profiles obtained for a FUBP family protein, we demonstrate that, in addition to being required to activate Myc to promote cell growth, Psi also directly binds and activates stg to couple growth and cell division. Thus, Psi knockdown results in reduced cell division in the wing imaginal disc. In addition to activating these pro-proliferative targets, Psi directly represses transcription of the growth inhibitor tolkin (tok, a metallopeptidase implicated in TGFß signalling). We further demonstrate tok overexpression inhibits proliferation, while tok loss of function increases mitosis alone and suppresses impaired cell division caused by Psi knockdown. Thus, Psi orchestrates growth through concurrent transcriptional activation of the pro-proliferative genes Myc and stg, in combination with repression of the growth inhibitor tok.

Stress-induced transcriptional memory accelerates promoter-proximal pause release and decelerates termination over mitotic divisions.

Heat shock instantly reprograms transcription. Whether gene and enhancer transcription fully recover from stress and whether stress establishes a memory by provoking transcription regulation that persists through mitosis remained unknown. Here, we measured nascent transcription and chromatin accessibility in unconditioned cells and in the daughters of stress-exposed cells. Tracking transcription genome-wide at nucleotide-resolution revealed that cells precisely restored RNA polymerase II (Pol II) distribution at gene bodies and enhancers upon recovery from stress. However, a single heat exposure in embryonic fibroblasts primed a faster gene induction in their daughter cells by increasing promoter-proximal Pol II pausing and by accelerating the pause release. In K562 erythroleukemia cells, repeated stress refined basal and heat-induced transcription over mitotic division and decelerated termination-coupled pre-mRNA processing. The slower termination retained transcripts on the chromatin and reduced recycling of Pol II. These results demonstrate that heat-induced transcriptional memory acts through promoter-proximal pause release and pre-mRNA processing at transcription termination.

EmPC-seq: Accurate RNA-sequencing and Bioinformatics Platform to Map RNA Polymerases and Remove Background Error.

Transcription errors can substantially affect metabolic processes in organisms by altering the epigenome and causing misincorporations in mRNA, which is translated into aberrant mutant proteins. Moreover, within eukaryotic genomes there are specific Transcription Error-Enriched genomic Loci (TEELs) which are transcribed by RNA polymerases with significantly higher error rates and hypothesized to have implications in cancer, aging, and diseases such as Down syndrome and Alzheimer's. Therefore, research into transcription errors is of growing importance within the field of genetics. Nevertheless, methodological barriers limit the progress in accurately identifying transcription errors. Pro-Seq and NET-Seq can purify nascent RNA and map RNA polymerases along the genome but cannot be used to identify transcriptional mutations. Here we present background Error Model-coupled Precision nuclear run-on Circular-sequencing (EmPC-seq), a method combining a nuclear run-on assay and circular sequencing with a background error model to precisely detect nascent transcription errors and effectively discern TEELs within the genome.

Pioneer-like factor GAF cooperates with PBAP (SWI/SNF) and NURF (ISWI) to regulate transcription.

Transcriptionally silent genes must be activated throughout development. This requires nucleosomes be removed from promoters and enhancers to allow transcription factor (TF) binding and recruitment of coactivators and RNA polymerase II (Pol II). Specialized pioneer TFs bind nucleosome-wrapped DNA to perform this chromatin opening by mechanisms that remain incompletely understood. Here, we show that GAGA factor (GAF), a Drosophila pioneer-like factor, functions with both SWI/SNF and ISWI family chromatin remodelers to allow recruitment of Pol II and entry to a promoter-proximal paused state, and also to promote Pol II's transition to productive elongation. We found that GAF interacts with PBAP (SWI/SNF) to open chromatin and allow Pol II to be recruited. Importantly, this activity is not dependent on NURF as previously proposed; however, GAF also synergizes with NURF downstream from this process to ensure efficient Pol II pause release and transition to productive elongation, apparently through its role in precisely positioning the +1 nucleosome. These results demonstrate how a single sequence-specific pioneer TF can synergize with remodelers to activate sets of genes. Furthermore, this behavior of remodelers is consistent with findings in yeast and mice, and likely represents general, conserved mechanisms found throughout eukarya.

RNA aptamer capture of macromolecular complexes for mass spectrometry analysis.

Specific genomic functions are dictated by macromolecular complexes (MCs) containing multiple proteins. Affinity purification of these complexes, often using antibodies, followed by mass spectrometry (MS) has revolutionized our ability to identify the composition of MCs. However, conventional immunoprecipitations suffer from contaminating antibody/serum-derived peptides that limit the sensitivity of detection for low-abundant interacting partners using MS. Here, we present AptA-MS (aptamer affinity-mass spectrometry), a robust strategy primarily using a specific, high-affinity RNA aptamer against Green Fluorescent Protein (GFP) to identify interactors of a GFP-tagged protein of interest by high-resolution MS. Utilizing this approach, we have identified the known molecular chaperones that interact with human Heat Shock Factor 1 (HSF1), and observed an increased association with several proteins upon heat shock, including translation elongation factors and histones. HSF1 is known to be regulated by multiple post-translational modifications (PTMs), and we observe both known and new sites of modifications on HSF1. We show that AptA-MS provides a dramatic target enrichment and detection sensitivity in evolutionarily diverse organisms and allows identification of PTMs without the need for modification-specific enrichments. In combination with the expanding libraries of GFP-tagged cell lines, this strategy offers a general, inexpensive, and high-resolution alternative to conventional approaches for studying MCs.

An AR-ERG transcriptional signature defined by long-range chromatin interactomes in prostate cancer cells.

The aberrant activities of transcription factors such as the androgen receptor (AR) underpin prostate cancer development. While the AR cis-regulation has been extensively studied in prostate cancer, information pertaining to the spatial architecture of the AR transcriptional circuitry remains limited. In this paper, we propose a novel framework to profile long-range chromatin interactions associated with AR and its collaborative transcription factor, erythroblast transformation-specific related gene (ERG), using chromatin interaction analysis by paired-end tag (ChIA-PET). We identified ERG-associated long-range chromatin interactions as a cooperative component in the AR-associated chromatin interactome, acting in concert to achieve coordinated regulation of a subset of AR target genes. Through multifaceted functional data analysis, we found that AR-ERG interaction hub regions are characterized by distinct functional signatures, including bidirectional transcription and cotranscription factor binding. In addition, cancer-associated long noncoding RNAs were found to be connected near protein-coding genes through AR-ERG looping. Finally, we found strong enrichment of prostate cancer genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) at AR-ERG co-binding sites participating in chromatin interactions and gene regulation, suggesting GWAS target genes identified from chromatin looping data provide more biologically relevant findings than using the nearest gene approach. Taken together, our results revealed an AR-ERG-centric higher-order chromatin structure that drives coordinated gene expression in prostate cancer progression and the identification of potential target genes for therapeutic intervention.

Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme.

The human genome encodes a variety of poorly understood RNA species that remain challenging to identify using existing genomic tools. We developed chromatin run-on and sequencing (ChRO-seq) to map the location of RNA polymerase for almost any input sample, including samples with degraded RNA that are intractable to RNA sequencing. We used ChRO-seq to map nascent transcription in primary human glioblastoma (GBM) brain tumors. Enhancers identified in primary GBMs resemble open chromatin in the normal human brain. Rare enhancers that are activated in malignant tissue drive regulatory programs similar to the developing nervous system. We identified enhancers that regulate groups of genes that are characteristic of each known GBM subtype and transcription factors that drive them. Finally we discovered a core group of transcription factors that control the expression of genes associated with clinical outcomes. This study characterizes the transcriptional landscape of GBM and introduces ChRO-seq as a method to map regulatory programs that contribute to complex diseases.

Nascent RNA sequencing reveals a dynamic global transcriptional response at genes and enhancers to the natural medicinal compound celastrol.

Most studies of responses to transcriptional stimuli measure changes in cellular mRNA concentrations. By sequencing nascent RNA instead, it is possible to detect changes in transcription in minutes rather than hours and thereby distinguish primary from secondary responses to regulatory signals. Here, we describe the use of PRO-seq to characterize the immediate transcriptional response in human cells to celastrol, a compound derived from traditional Chinese medicine that has potent anti-inflammatory, tumor-inhibitory, and obesity-controlling effects. Celastrol is known to elicit a cellular stress response resembling the response to heat shock, but the transcriptional basis of this response remains unclear. Our analysis of PRO-seq data for K562 cells reveals dramatic transcriptional effects soon after celastrol treatment at a broad collection of both coding and noncoding transcription units. This transcriptional response occurred in two major waves, one within 10 min, and a second 40-60 min after treatment. Transcriptional activity was generally repressed by celastrol, but one distinct group of genes, enriched for roles in the heat shock response, displayed strong activation. Using a regression approach, we identified key transcription factors that appear to drive these transcriptional responses, including members of the E2F and RFX families. We also found sequence-based evidence that particular transcription factors drive the activation of enhancers. We observed increased polymerase pausing at both genes and enhancers, suggesting that pause release may be widely inhibited during the celastrol response. Our study demonstrates that a careful analysis of PRO-seq time-course data can disentangle key aspects of a complex transcriptional response, and it provides new insights into the activity of a powerful pharmacological agent.

High-Resolution Mapping of RNA Polymerases Identifies Mechanisms of Sensitivity and Resistance to BET Inhibitors in t(8;21) AML.

Bromodomain and extra-terminal domain (BET) family inhibitors offer an approach to treating hematological malignancies. We used precision nuclear run-on transcription sequencing (PRO-seq) to create high-resolution maps of active RNA polymerases across the genome in t(8;21) acute myeloid leukemia (AML), as these polymerases are exceptionally sensitive to BET inhibitors. PRO-seq identified over 1,400 genes showing impaired release of promoter-proximal paused RNA polymerases, including the stem cell factor receptor tyrosine kinase KIT that is mutated in t(8;21) AML. PRO-seq also identified an enhancer 3' to KIT. Chromosome conformation capture confirmed contacts between this enhancer and the KIT promoter, while CRISPRi-mediated repression of this enhancer impaired cell growth. PRO-seq also identified microRNAs, including MIR29C and MIR29B2, that target the anti-apoptotic factor MCL1 and were repressed by BET inhibitors. MCL1 protein was upregulated, and inhibition of BET proteins sensitized t(8:21)-containing cells to MCL1 inhibition, suggesting a potential mechanism of resistance to BET-inhibitor-induced cell death.

Getting up to speed with transcription elongation by RNA polymerase II.

Recent advances in sequencing techniques that measure nascent transcripts and that reveal the positioning of RNA polymerase II (Pol II) have shown that the pausing of Pol II in promoter-proximal regions and its release to initiate a phase of productive elongation are key steps in transcription regulation. Moreover, after the release of Pol II from the promoter-proximal region, elongation rates are highly dynamic throughout the transcription of a gene, and vary on a gene-by-gene basis. Interestingly, Pol II elongation rates affect co-transcriptional processes such as splicing, termination and genome stability. Increasing numbers of factors and regulatory mechanisms have been associated with the steps of transcription elongation by Pol II, revealing that elongation is a highly complex process. Elongation is thus now recognized as a key phase in the regulation of transcription by Pol II.

Targeted H3R26 deimination specifically facilitates estrogen receptor binding by modifying nucleosome structure.

Transcription factor binding to DNA in vivo causes the recruitment of chromatin modifiers that can cause changes in chromatin structure, including the modification of histone tails. We previously showed that estrogen receptor (ER) target gene activation is facilitated by peptidylarginine deiminase 2 (PAD2)-catalyzed histone H3R26 deimination (H3R26Cit). Here we report that the genomic distributions of ER and H3R26Cit in breast cancer cells are strikingly coincident, linearly correlated, and observed as early as 2 minutes following estradiol treatment. The H3R26Cit profile is unlike that of previously described histone modifications and is characterized by sharp, narrow peaks. Paired-end MNase ChIP-seq indicates that the charge-neutral H3R26Cit modification facilitates ER binding to DNA by altering the fine structure of the nucleosome. Clinically, we find that PAD2 and H3R26Cit levels correlate with ER expression in breast tumors and that high PAD2 expression is associated with increased survival in ER+ breast cancer patients. These findings provide insight into how transcription factors gain access to nucleosomal DNA and implicate PAD2 as a novel therapeutic target for ER+ breast cancer.

Inhibiting heat shock factor 1 in human cancer cells with a potent RNA aptamer.

Heat shock factor 1 (HSF1) is a master regulator that coordinates chaperone protein expression to enhance cellular survival in the face of heat stress. In cancer cells, HSF1 drives a transcriptional program distinct from heat shock to promote metastasis and cell survival. Its strong association with the malignant phenotype implies that HSF1 antagonists may have general and effective utilities in cancer therapy. For this purpose, we had identified an avid RNA aptamer for HSF1 that is portable among different model organisms. Extending our previous work in yeast and Drosophila, here we report the activity of this aptamer in human cancer cell lines. When delivered into cells using a synthetic gene and strong promoter, this aptamer was able to prevent HSF1 from binding to its DNA regulation elements. At the cellular level, expression of this aptamer induced apoptosis and abolished the colony-forming capability of cancer cells. At the molecular level, it reduced chaperones and attenuated the activation of the MAPK signaling pathway. Collectively, these data demonstrate the advantage of aptamers in drug target validation and support the hypothesis that HSF1 DNA binding activity is a potential target for controlling oncogenic transformation and neoplastic growth.

The metazoan-specific mediator subunit 26 (Med26) is essential for viability and is found at both active genes and pericentric heterochromatin in Drosophila melanogaster.

Human MED26 was originally purified in the cofactor required for the Sp1 activation complex (CRSP) as a 70-kDa component named CRSP70. This polypeptide was specific to metazoans and the "small" form of the Mediator complex. We report here that a Drosophila melanogaster homologue of MED26 similarly interacts with other components of the core Drosophila Mediator complex but not with the kinase module and is recruited to genes upon activation. Using a null allele of Med26, we show that Med26 is required for organismal viability but not for cell proliferation or survival. Clones lacking Med26 in the wing disc lead to loss of the adult wing margin and reduced expression of genes involved in wing margin formation. Surprisingly, when polytene chromosomes from the salivary gland were examined using antibodies to Med26, it was apparent that a fraction of the protein was associated with the chromocenter, which contains pericentric heterochromatin. This staining colocalizes with heterochromatin protein 1 (HP1). Immunoprecipitation experiments show that Med26 interacts with HP1. The interaction is mediated through the chromoshadow domain of HP1 and through the conserved motif in the carboxy terminus of the Med26 protein. This work is the first characterization of the metazoan-specific Mediator subunit in an animal model.

Discovering aptamers by cell-SELEX against human soluble growth factors ectopically expressed on yeast cell surface.

SELEX, the process of selecting aptamers, is often hampered by the difficulty of preparing target molecules in their native forms and by a lack of a simple yet quantitative assay for monitoring enrichment and affinity of reactive aptamers. In this study, we sought to discover DNA aptamers against human serum markers for potential therapeutic and diagnostic applications. To circumvent soluble expression and immobilization for performing SELEX, we ectopically expressed soluble growth factors on the surface of yeast cells to enable cell-SELEX and devised a flow cytometry-based method to quantitatively monitor progressive enrichment of specific aptamers. High-throughput sequencing of selected pools revealed that the emergence of highly enriched sequences concurred with the increase in the percentage of reactive aptamers shown by flow cytometry. Particularly, selected DNA aptamers against VEGF were specific and of high affinity (K(D)  = ∼ 1 nM) and demonstrated a potent inhibition of capillary tube formation of endothelial cells, comparable to the effect of a clinically approved anti-VEGF antibody drug, bevacizumab. Considering the fact that many mammalian secretory proteins have been functionally expressed in yeast, the strategy of implementing cell-SELEX and quantitative binding assay can be extended to discover aptamers against a broad array of soluble antigens.

Condensin controls recruitment of RNA polymerase II to achieve nematode X-chromosome dosage compensation.

The X-chromosome gene regulatory process called dosage compensation ensures that males (1X) and females (2X) express equal levels of X-chromosome transcripts. The mechanism in Caenorhabditis elegans has been elusive due to improperly annotated transcription start sites (TSSs). Here we define TSSs and the distribution of transcriptionally engaged RNA polymerase II (Pol II) genome-wide in wild-type and dosage-compensation-defective animals to dissect this regulatory mechanism. Our TSS-mapping strategy integrates GRO-seq, which tracks nascent transcription, with a new derivative of this method, called GRO-cap, which recovers nascent RNAs with 5' caps prior to their removal by co-transcriptional processing. Our analyses reveal that promoter-proximal pausing is rare, unlike in other metazoans, and promoters are unexpectedly far upstream from the 5' ends of mature mRNAs. We find that C. elegans equalizes X-chromosome expression between the sexes, to a level equivalent to autosomes, by reducing Pol II recruitment to promoters of hermaphrodite X-linked genes using a chromosome-restructuring condensin complex. DOI:http://dx.doi.org/10.7554/eLife.00808.001.

Density-dependent cooperative non-specific binding in solid-phase SELEX affinity selection.

The non-specific binding of undesired ligands to a target is the primary factor limiting the enrichment of tight-binding ligands in affinity selection. Solution-phase non-specific affinity is determined by the free-energy of ligand binding to a single target. However, the solid-phase affinity might be higher if a ligand bound concurrently to multiple adjacent immobilized targets in a cooperative manner. Cooperativity could emerge in this case as a simple consequence of the relationship between the free energy of binding, localization entropy and the spatial distribution of the immobilized targets. We tested this hypothesis using a SELEX experimental design and found that non-specific RNA aptamer ligands can concurrently bind up to four bead-immobilized peptide targets, and that this can increase their effective binding affinity by two orders-of-magnitude. Binding curves were quantitatively explained by a new statistical mechanical model of density-dependent cooperative binding, which relates cooperative binding to both the target concentration and the target surface density on the immobilizing substrate. Target immobilization plays a key role in SELEX and other ligand enrichment methods, particularly in new multiplexed microfluidic purification devices, and these results have strong implications for optimizing their performance.

Signaling pathways differentially affect RNA polymerase II initiation, pausing, and elongation rate in cells.

RNA polymerase II (Pol II) transcribes hundreds of kilobases of DNA, limiting the production of mRNAs and lncRNAs. We used global run-on sequencing (GRO-seq) to measure the rates of transcription by Pol II following gene activation. Elongation rates vary as much as 4-fold at different genomic loci and in response to two distinct cellular signaling pathways (i.e., 17ß-estradiol [E2] and TNF-α). The rates are slowest near the promoter and increase during the first ~15 kb transcribed. Gene body elongation rates correlate with Pol II density, resulting in systematically higher rates of transcript production at genes with higher Pol II density. Pol II dynamics following short inductions indicate that E2 stimulates gene expression by increasing Pol II initiation, whereas TNF-α reduces Pol II residence time at pause sites. Collectively, our results identify previously uncharacterized variation in the rate of transcription and highlight elongation as an important, variable, and regulated rate-limiting step during transcription.

Peptidylarginine deiminase 2-catalyzed histone H3 arginine 26 citrullination facilitates estrogen receptor α target gene activation.

Cofactors for estrogen receptor α (ERα) can modulate gene activity by posttranslationally modifying histone tails at target promoters. Here, we found that stimulation of ERα-positive cells with 17ß-estradiol (E2) promotes global citrullination of histone H3 arginine 26 (H3R26) on chromatin. Additionally, we found that the H3 citrulline 26 (H3Cit26) modification colocalizes with ERα at decondensed chromatin loci surrounding the estrogen-response elements of target promoters. Surprisingly, we also found that citrullination of H3R26 is catalyzed by peptidylarginine deiminase (PAD) 2 and not by PAD4 (which citrullinates H4R3). Further, we showed that PAD2 interacts with ERα after E2 stimulation and that inhibition of either PAD2 or ERα strongly suppresses E2-induced H3R26 citrullination and ERα recruitment at target gene promoters. Collectively, our data suggest that E2 stimulation induces the recruitment of PAD2 to target promoters by ERα, whereby PAD2 then citrullinates H3R26, which leads to local chromatin decondensation and transcriptional activation.