The Ground State and Evolution of Promoter Region Directionality.

Eukaryotic promoter regions are frequently divergently transcribed in vivo, but it is unknown whether the resultant antisense RNAs are a mechanistic by-product of RNA polymerase II (Pol II) transcription or biologically meaningful. Here, we use a functional evolutionary approach that involves nascent transcript mapping in S. cerevisiae strains containing foreign yeast DNA. Promoter regions in foreign environments lose the directionality they have in their native species. Strikingly, fortuitous promoter regions arising in foreign DNA produce equal transcription in both directions, indicating that divergent transcription is a mechanistic feature that does not imply a function for these transcripts. Fortuitous promoter regions arising during evolution promote bidirectional transcription and over time are purged through mutation or retained to enable new functionality. Similarly, human transcription is more bidirectional at newly evolved enhancers and promoter regions. Thus, promoter regions are intrinsically bidirectional and are shaped by evolution to bias transcription toward coding versus non-coding RNAs.

Intrinsically disordered regions (IDRs): A vague and confusing concept for protein function.

The term "intrinsically disordered region" (IDR) in proteins has been used in numerous publications. However, most proteins contain IDRs, the term refers to very different types of structures and functions, and many IDRs become structured upon interaction with other biomolecules. Thus, IDR is an unnecessary, vague, and ultimately confusing concept.

Functional analysis of a random-sequence chromosome reveals a high level and the molecular nature of transcriptional noise in yeast cells.

We measure transcriptional noise in yeast by analyzing chromatin structure and transcription of an 18-kb region of DNA whose sequence was randomly generated. Nucleosomes fully occupy random-sequence DNA, but nucleosome-depleted regions (NDRs) are much less frequent, and there are fewer well-positioned nucleosomes and shorter nucleosome arrays. Steady-state levels of random-sequence RNAs are comparable to yeast mRNAs, although transcription and decay rates are higher. Transcriptional initiation from random-sequence DNA occurs at numerous sites, indicating very low intrinsic specificity of the RNA Pol II machinery. In contrast, poly(A) profiles of random-sequence RNAs are roughly comparable to those of yeast mRNAs, suggesting limited evolutionary restraints on poly(A) site choice. Random-sequence RNAs show higher cell-to-cell variability than yeast mRNAs, suggesting that functional elements limit variability. These observations indicate that transcriptional noise occurs at high levels in yeast, and they provide insight into how chromatin and transcription patterns arise from the evolved yeast genome.

Global analysis of mRNA isoform half-lives reveals stabilizing and destabilizing elements in yeast.

We measured half-lives of 21,248 mRNA 3' isoforms in yeast by rapidly depleting RNA polymerase II from the nucleus and performing direct RNA sequencing throughout the decay process. Interestingly, half-lives of mRNA isoforms from the same gene, including nearly identical isoforms, often vary widely. Based on clusters of isoforms with different half-lives, we identify hundreds of sequences conferring stabilization or destabilization upon mRNAs terminating downstream. One class of stabilizing element is a polyU sequence that can interact with poly(A) tails, inhibit the association of poly(A)-binding protein, and confer increased stability upon introduction into ectopic transcripts. More generally, destabilizing and stabilizing elements are linked to the propensity of the poly(A) tail to engage in double-stranded structures. Isoforms engineered to fold into 3' stem-loop structures not involving the poly(A) tail exhibit even longer half-lives. We suggest that double-stranded structures at 3' ends are a major determinant of mRNA stability.

Non-canonical functions of enhancers: regulation of RNA polymerase III transcription, DNA replication, and V(D)J recombination.

Enhancers are the key regulators of other DNA-based processes by virtue of their unique ability to generate nucleosome-depleted regions in a highly regulated manner. Enhancers regulate cell-type-specific transcription of tRNA genes by RNA polymerase III (Pol III). They are also responsible for the binding of the origin replication complex (ORC) to DNA replication origins, thereby regulating origin utilization, replication timing, and replication-dependent chromosome breaks. Additionally, enhancers regulate V(D)J recombination by increasing access of the recombination-activating gene (RAG) recombinase to target sites and by generating non-coding enhancer RNAs and localized regions of trimethylated histone H3-K4 recognized by the RAG2 PHD domain. Thus, enhancers represent the first step in decoding the genome, and hence they regulate biological processes that, unlike RNA polymerase II (Pol II) transcription, do not have dedicated regulatory proteins.

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate.

The RNA polymerase II (Pol II) elongation rate influences poly(A) site selection, with slow and fast Pol II derivatives causing upstream and downstream shifts, respectively, in poly(A) site utilization. In yeast, depletion of either of the histone chaperones FACT or Spt6 causes an upstream shift of poly(A) site use that strongly resembles the poly(A) profiles of slow Pol II mutant strains. Like slow Pol II mutant strains, FACT- and Spt6-depleted cells exhibit Pol II processivity defects, indicating that both Spt6 and FACT stimulate the Pol II elongation rate. Poly(A) profiles of some genes show atypical downstream shifts; this subset of genes overlaps well for FACT- or Spt6-depleted strains but is different from the atypical genes in Pol II speed mutant strains. In contrast, depletion of histone H3 or H4 causes a downstream shift of poly(A) sites for most genes, indicating that nucleosomes inhibit the Pol II elongation rate in vivo. Thus, chromatin-based control of the Pol II elongation rate is a potential mechanism, distinct from direct effects on the cleavage/polyadenylation machinery, to regulate alternative polyadenylation in response to genetic or environmental changes.

An HNF4α-miRNA inflammatory feedback circuit regulates hepatocellular oncogenesis.

Hepatocyte nuclear factor 4α (HNF4α) is essential for liver development and hepatocyte function. Here, we show that transient inhibition of HNF4α initiates hepatocellular transformation through a microRNA-inflammatory feedback loop circuit consisting of miR-124, IL6R, STAT3, miR-24, and miR-629. Moreover, we show that, once this circuit is activated, it maintains suppression of HNF4α and sustains oncogenesis. Systemic administration of miR-124, which modulates inflammatory signaling, prevents and suppresses hepatocellular carcinogenesis by inducing tumor-specific apoptosis without toxic side effects. As we also show that this HNF4α circuit is perturbed in human hepatocellular carcinomas, our data raise the possibility that manipulation of this microRNA feedback-inflammatory loop has therapeutic potential for treating liver cancer.

Extensive Structural Differences of Closely Related 3' mRNA Isoforms: Links to Pab1 Binding and mRNA Stability.

Alternative polyadenylation generates numerous 3' mRNA isoforms that can vary in biological properties, such as stability and localization. We developed methods to obtain transcriptome-scale structural information and protein binding on individual 3' mRNA isoforms in vivo. Strikingly, near-identical mRNA isoforms can possess dramatically different structures throughout the 3' UTR. Analyses of identical mRNAs in different species or refolded in vitro indicate that structural differences in vivo are often due to trans-acting factors. The level of Pab1 binding to poly(A)-containing isoforms is surprisingly variable, and differences in Pab1 binding correlate with the extent of structural variation for closely spaced isoforms. A pattern encompassing single-strandedness near the 3' terminus, double-strandedness of the poly(A) tail, and low Pab1 binding is associated with mRNA stability. Thus, individual 3' mRNA isoforms can be remarkably different physical entities in vivo. Sequences responsible for isoform-specific structures, differential Pab1 binding, and mRNA stability are evolutionarily conserved, indicating biological function.

Condition-specific 3' mRNA isoform half-lives and stability elements in yeast.

Alternative polyadenylation generates numerous 3' mRNA isoforms that can differ in their stability, structure, and function. These isoforms can be used to map mRNA stabilizing and destabilizing elements within 3' untranslated regions (3'UTRs). Here, we examine how environmental conditions affect 3' mRNA isoform turnover and structure in yeast cells on a transcriptome scale. Isoform stability broadly increases when cells grow more slowly, with relative half-lives of most isoforms being well correlated across multiple conditions. Surprisingly, dimethyl sulfate probing reveals that individual 3' isoforms have similar structures across different conditions, in contrast to the extensive structural differences that can exist between closely related isoforms in an individual condition. Unexpectedly, most mRNA stabilizing and destabilizing elements function only in a single growth condition. The genes associated with some classes of condition-specific stability elements are enriched for different functional categories, suggesting that regulated mRNA stability might contribute to adaptation to different growth environments. Condition-specific stability elements do not result in corresponding condition-specific changes in steady-state mRNA isoform levels. This observation is consistent with a compensatory mechanism between polyadenylation and stability, and it suggests that condition-specific mRNA stability elements might largely reflect condition-specific regulation of mRNA 3' end formation.

A compensatory link between cleavage/polyadenylation and mRNA turnover regulates steady-state mRNA levels in yeast.

Cells have compensatory mechanisms to coordinate the rates of major biological processes, thereby permitting growth in a wide variety of conditions. Here, we uncover a compensatory link between cleavage/polyadenylation in the nucleus and messenger RNA (mRNA) turnover in the cytoplasm. On a global basis, same-gene 3' mRNA isoforms with twofold or greater differences in half-lives have steady-state mRNA levels that differ by significantly less than a factor of 2. In addition, increased efficiency of cleavage/polyadenylation at a specific site is associated with reduced stability of the corresponding 3' mRNA isoform. This inverse relationship between cleavage/polyadenylation and mRNA isoform half-life reduces the variability in the steady-state levels of mRNA isoforms, and it occurs in all four growth conditions tested. These observations suggest that during cleavage/polyadenylation in the nucleus, mRNA isoforms are marked in a manner that persists upon translocation to the cytoplasm and affects the activity of mRNA degradation machinery, thus influencing mRNA stability.

Elongation rate of RNA polymerase II affects pausing patterns across 3' UTRs.

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located < 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

How is polyadenylation restricted to 3'-untranslated regions?

Polyadenylation occurs at numerous sites within 3'-untranslated regions (3'-UTRs) but rarely within coding regions. How does Pol II travel through long coding regions without generating poly(A) sites, yet then permits promiscuous polyadenylation once it reaches the 3'-UTR? The cleavage/polyadenylation (CpA) machinery preferentially associates with 3'-UTRs, but it is unknown how its recruitment is restricted to 3'-UTRs during Pol II elongation. Unlike coding regions, 3'-UTRs have long AT-rich stretches of DNA that may be important for restricting polyadenylation to 3'-UTRs. Recognition of the 3'-UTR could occur at the DNA (AT-rich), RNA (AU-rich), or RNA:DNA hybrid (rU:dA- and/or rA:dT-rich) level. Based on the nucleic acid critical for 3'-UTR recognition, there are three classes of models, not mutually exclusive, for how the CpA machinery is selectively recruited to 3'-UTRs, thereby restricting where polyadenylation occurs: (1) RNA-based models suggest that the CpA complex directly (or indirectly through one or more intermediary proteins) binds long AU-rich stretches that are exposed after Pol II passes through these regions. (2) DNA-based models suggest that the AT-rich sequence affects nucleosome depletion or the elongating Pol II machinery, resulting in dissociation of some elongation factors and subsequent recruitment of the CpA machinery. (3) RNA:DNA hybrid models suggest that preferential destabilization of the Pol II elongation complex at rU:dA- and/or rA:dT-rich duplexes bridging the nucleotide addition and RNA exit sites permits preferential association of the CpA machinery with 3'-UTRs. Experiments to provide evidence for one or more of these models are suggested.

An epigenetic switch involving NF-kappaB, Lin28, Let-7 MicroRNA, and IL6 links inflammation to cell transformation.

Inflammation is linked clinically and epidemiologically to cancer, and NF-kappaB appears to play a causative role, but the mechanisms are poorly understood. We show that transient activation of Src oncoprotein can mediate an epigenetic switch from immortalized breast cells to a stably transformed line that forms self-renewing mammospheres that contain cancer stem cells. Src activation triggers an inflammatory response mediated by NF-kappaB that directly activates Lin28 transcription and rapidly reduces let-7 microRNA levels. Let-7 directly inhibits IL6 expression, resulting in higher levels of IL6 than achieved by NF-kappaB activation. IL6-mediated activation of the STAT3 transcription factor is necessary for transformation, and IL6 activates NF-kappaB, thereby completing a positive feedback loop. This regulatory circuit operates in other cancer cells lines, and its transcriptional signature is found in human cancer tissues. Thus, inflammation activates a positive feedback loop that maintains the epigenetic transformed state for many generations in the absence of the inducing signal.

Mediator Undergoes a Compositional Change during Transcriptional Activation.

Mediator is a transcriptional co-activator recruited to enhancers by DNA-binding activators, and it also interacts with RNA polymerase (Pol) II as part of the preinitiation complex (PIC). We demonstrate that a single Mediator complex associates with the enhancer and core promoter in vivo, indicating that it can physically bridge these transcriptional elements. However, the Mediator kinase module associates strongly with the enhancer, but not with the core promoter, and it dissociates from the enhancer upon depletion of the TFIIH kinase. Severing the kinase module from Mediator by removing the connecting subunit Med13 does not affect Mediator association at the core promoter but increases occupancy at enhancers. Thus, Mediator undergoes a compositional change in which the kinase module, recruited via Mediator to the enhancer, dissociates from Mediator to permit association with Pol II and the PIC. As such, Mediator acts as a dynamic bridge between the enhancer and core promoter.

Different SP1 binding dynamics at individual genomic loci in human cells.

Using a tamoxifen-inducible time-course ChIP-sequencing (ChIP-seq) approach, we show that the ubiquitous transcription factor SP1 has different binding dynamics at its target sites in the human genome. SP1 very rapidly reaches maximal binding levels at some sites, but binding kinetics at other sites is biphasic, with rapid half-maximal binding followed by a considerably slower increase to maximal binding. While ∼70% of SP1 binding sites are located at promoter regions, loci with slow SP1 binding kinetics are enriched in enhancer and Polycomb-repressed regions. Unexpectedly, SP1 sites with fast binding kinetics tend to have higher quality and more copies of the SP1 sequence motif. Different cobinding factors associate near SP1 binding sites depending on their binding kinetics and on their location at promoters or enhancers. For example, NFY and FOS are preferentially associated near promoter-bound SP1 sites with fast binding kinetics, whereas DNA motifs of ETS and homeodomain proteins are preferentially observed at sites with slow binding kinetics. At promoters but not enhancers, proteins involved in sumoylation and PML bodies associate more strongly with slow SP1 binding sites than with the fast binding sites. The speed of SP1 binding is not associated with nucleosome occupancy, and it is not necessarily coupled to higher transcriptional activity. These results with SP1 are in contrast to those of human TBP, indicating that there is no common mechanism affecting transcription factor binding kinetics. The biphasic kinetics at some SP1 target sites suggest the existence of distinct chromatin states at these loci in different cells within the overall population.

Pheno-RNA, a method to associate genes with a specific phenotype, identifies genes linked to cellular transformation.

Cellular transformation is associated with dramatic changes in gene expression, but it is difficult to determine which regulated genes are oncogenically relevant. Here we describe Pheno-RNA, a general approach to identifying candidate genes associated with a specific phenotype. Specifically, we generate a "phenotypic series" by treating a nontransformed breast cell line with a wide variety of molecules that induce cellular transformation to various extents. By performing transcriptional profiling across this phenotypic series, the expression profile of every gene can be correlated with the strength of the transformed phenotype. We identify ∼200 genes whose expression profiles are very highly correlated with the transformation phenotype, strongly suggesting their importance in transformation. Within biological categories linked to cancer, some genes show high correlations with the transformed phenotype, but others do not. Many genes whose expression profiles are highly correlated with transformation have never been associated with cancer, suggesting the involvement of heretofore unknown genes in cancer.

Promoter-specific dynamics of TATA-binding protein association with the human genome.

Transcription factor binding to target sites in vivo is a dynamic process that involves cycles of association and dissociation, with individual proteins differing in their binding dynamics. The dynamics at individual sites on a genomic scale have been investigated in yeast cells, but comparable experiments have not been done in multicellular eukaryotes. Here, we describe a tamoxifen-inducible, time-course ChIP-seq approach to measure transcription factor binding dynamics at target sites throughout the human genome. As observed in yeast cells, the TATA-binding protein (TBP) typically displays rapid turnover at RNA polymerase (Pol) II-transcribed promoters, slow turnover at Pol III promoters, and very slow turnover at the Pol I promoter. Turnover rates vary widely among Pol II promoters in a manner that does not correlate with the level of TBP occupancy. Human Pol II promoters with slow TBP dissociation preferentially contain a TATA consensus motif, support high transcriptional activity of downstream genes, and are linked with specific activators and chromatin remodelers. These properties of human promoters with slow TBP turnover differ from those of yeast promoters with slow turnover. These observations suggest that TBP binding dynamics differentially affect promoter function and gene expression, possibly at the level of transcriptional reinitiation/bursting.

TFIIH phosphorylation of the Pol II CTD stimulates mediator dissociation from the preinitiation complex and promoter escape.

The transition between transcriptional initiation and elongation by RNA polymerase (Pol) II is associated with phosphorylation of its C-terminal tail (CTD). Depletion of Kin28, the TFIIH subunit that phosphorylates the CTD, does not affect elongation but causes Pol II occupancy profiles to shift upstream in a FACT-independent manner indicative of a defect in promoter escape. Stronger defects in promoter escape are linked to stronger effects on preinitiation complex formation and transcription, suggesting that impairment in promoter escape results in premature dissociation of general factors and Pol II near the promoter. Kin28 has a stronger effect on genes whose transcription is dependent on SAGA as opposed to TFIID. Strikingly, Kin28 depletion causes a dramatic increase in Mediator at the core promoter. These observations suggest that TFIIH phosphorylation of the CTD causes Mediator dissociation, thereby permitting rapid promoter escape of Pol II from the preinitiation complex.

Inflammatory regulatory network mediated by the joint action of NF-kB, STAT3, and AP-1 factors is involved in many human cancers.

Using an inducible, inflammatory model of breast cellular transformation, we describe the transcriptional regulatory network mediated by STAT3, NF-κB, and AP-1 factors on a genomic scale. These proinflammatory regulators form transcriptional complexes that directly regulate the expression of hundreds of genes in oncogenic pathways via a positive feedback loop. This transcriptional feedback loop and associated network functions to various extents in many types of cancer cells and patient tumors, and it is the basis for a cancer inflammation index that defines cancer types by functional criteria. We identify a network of noninflammatory genes whose expression is well correlated with the cancer inflammatory index. Conversely, the cancer inflammation index is negatively correlated with the expression of genes involved in DNA metabolism, and transformation is associated with genome instability. We identify drugs whose efficacy in cell lines is correlated with the cancer inflammation index, suggesting the possibility of using this index for personalized cancer therapy. Inflammatory tumors are preferentially associated with infiltrating immune cells that might be recruited to the site of the tumor via inflammatory molecules produced by the cancer cells.

Targeted profiling of RNA translation reveals mTOR-4EBP1/2-independent translation regulation of mRNAs encoding ribosomal proteins.

The PI3K-Akt-mTOR signaling pathway is a master regulator of RNA translation. Pharmacological inhibition of this pathway preferentially and coordinately suppresses, in a 4EBP1/2-dependent manner, translation of mRNAs encoding ribosomal proteins. However, it is unclear whether mechanistic target of rapamycin (mTOR)-4EBP1/2 is the exclusive translation regulator of this group of genes, and furthermore, systematic searches for novel translation modulators have been immensely challenging because of difficulties in scaling existing RNA translation profiling assays. Here, we developed a rapid and highly scalable approach for gene-specific quantitation of RNA translation, termed Targeted Profiling of RNA Translation (TPRT). We applied this technique in a chemical screen for translation modulators, and identified numerous preclinical and clinical therapeutic compounds, with diverse nominal targets, that preferentially suppress translation of ribosomal proteins. Surprisingly, some of these compounds act in a manner that bypasses canonical regulation by mTOR-4EBP1/2. Instead, these compounds exert their translation effects in a manner that is dependent on GCN2-eIF2α, a central signaling axis within the integrated stress response. Furthermore, we were also able to identify metabolic perturbations that also suppress ribosomal protein translation in an mTOR-independent manner. Together, we describe a translation assay that is directly applicable to large-scale RNA translation studies, and that enabled us to identify a noncanonical, mTOR-independent mode for translation regulation of ribosomal proteins.