Mutant FOXO1 controls an oncogenic network via enhancer accessibility.

Transcriptional dysregulation is a hallmark of diffuse large B cell lymphoma (DLBCL), as transcriptional regulators are frequently mutated. However, our mechanistic understanding of how normal transcriptional programs are co-opted in DLBCL has been hindered by a lack of methodologies that provide the temporal resolution required to separate direct and indirect effects on transcriptional control. We applied a chemical-genetic approach to engineer the inducible degradation of the transcription factor FOXO1, which is recurrently mutated (mFOXO1) in DLBCL. The combination of rapid degradation of mFOXO1, nascent transcript detection, and assessment of chromatin accessibility allowed us to identify the direct targets of mFOXO1. mFOXO1 was required to maintain accessibility at specific enhancers associated with multiple oncogenes, and mFOXO1 degradation impaired RNA polymerase pause-release at some targets. Wild-type FOXO1 appeared to weakly regulate many of the same targets as mFOXO1 and was able to complement the degradation of mFOXO1 in the context of AKT inhibition.

Reducing the stress of revisions: A chat with Scott W. Hiebert.

While reviewers' recommendations can strengthen a manuscript, responding to their concerns can be a tricky and frustrating process for authors. For our special issue on stress, Molecular Cell speaks with Scott Hiebert about his advice for responding to reviewers, crafting an effective revision plan, and involving trainees in the process.

Histone deacetylases maintain expression of the pluripotent gene network via recruitment of RNA polymerase II to coding and noncoding loci.

Histone acetylation is a dynamic modification regulated by the opposing actions of histone acetyltransferases (HATs) and histone deacetylases (HDACs). Deacetylation of histone tails results in chromatin tightening, and therefore, HDACs are generally regarded as transcriptional repressors. Counterintuitively, simultaneous deletion of Hdac1 and Hdac2 in embryonic stem cells (ESCs) reduces expression of the pluripotency-associated transcription factors Pou5f1, Sox2, and Nanog (PSN). By shaping global histone acetylation patterns, HDACs indirectly regulate the activity of acetyl-lysine readers, such as the transcriptional activator BRD4. Here, we use inhibitors of HDACs and BRD4 (LBH589 and JQ1, respectively) in combination with precision nuclear run-on and sequencing (PRO-seq) to examine their roles in defining the ESC transcriptome. Both LBH589 and JQ1 cause a marked reduction in the pluripotent gene network. However, although JQ1 treatment induces widespread transcriptional pausing, HDAC inhibition causes a reduction in both paused and elongating polymerase, suggesting an overall reduction in polymerase recruitment. Using enhancer RNA (eRNA) expression to measure enhancer activity, we find that LBH589-sensitive eRNAs are preferentially associated with superenhancers and PSN binding sites. These findings suggest that HDAC activity is required to maintain pluripotency by regulating the PSN enhancer network via the recruitment of RNA polymerase II.

Chemical-genetics refines transcription factor regulatory circuits.

Transcriptional dysregulation is a key step in oncogenesis, but our understanding of transcriptional control has relied on genetic approaches that are slow and allow for compensation. Chemical-genetic approaches have shortened the time frame for the analysis of transcription factors from days or weeks to minutes. These studies show that while DNA-binding proteins bind to thousands of sites, they are directly required to regulate only a small cadre of genes. Moreover, these transcriptional control networks are far more distinct, with much less overlap and interconnectivity than predicted from DNA binding. The identified direct targets can then be used to dissect the mechanism of action of these factors, which could identify ways to therapeutically manipulate these oncogenic transcriptional control networks.

Myo-differentiation reporter screen reveals NF-Y as an activator of PAX3-FOXO1 in rhabdomyosarcoma.

Recurrent chromosomal rearrangements found in rhabdomyosarcoma (RMS) produce the PAX3-FOXO1 fusion protein, which is an oncogenic driver and a dependency in this disease. One important function of PAX3-FOXO1 is to arrest myogenic differentiation, which is linked to the ability of RMS cells to gain an unlimited proliferation potential. Here, we developed a phenotypic screening strategy for identifying factors that collaborate with PAX3-FOXO1 to block myo-differentiation in RMS. Unlike most genes evaluated in our screen, we found that loss of any of the three subunits of the Nuclear Factor Y (NF-Y) complex leads to a myo-differentiation phenotype that resembles the effect of inactivating PAX3-FOXO1. While the transcriptomes of NF-Y- and PAX3-FOXO1-deficient RMS cells bear remarkable similarity to one another, we found that these two transcription factors occupy nonoverlapping sites along the genome: NF-Y preferentially occupies promoters, whereas PAX3-FOXO1 primarily binds to distal enhancers. By integrating multiple functional approaches, we map the PAX3 promoter as the point of intersection between these two regulators. We show that NF-Y occupies CCAAT motifs present upstream of PAX3 to function as a transcriptional activator of PAX3-FOXO1 expression in RMS. These findings reveal a critical upstream role of NF-Y in the oncogenic PAX3-FOXO1 pathway, highlighting how a broadly essential transcription factor can perform tumor-specific roles in governing cellular state.

Human SMARCA5 is continuously required to maintain nucleosome spacing.

Genetic models suggested that SMARCA5 was required for DNA-templated events including transcription, DNA replication, and DNA repair. We engineered a degron tag into the endogenous alleles of SMARCA5, a catalytic component of the imitation switch complexes in three different human cell lines to define the effects of rapid degradation of this key regulator. Degradation of SMARCA5 was associated with a rapid increase in global nucleosome repeat length, which may allow greater chromatin compaction. However, there were few changes in nascent transcription within the first 6 h of degradation. Nevertheless, we demonstrated a requirement for SMARCA5 to control nucleosome repeat length at G1/S and during the S phase. SMARCA5 co-localized with CTCF and H2A.Z, and we found a rapid loss of CTCF DNA binding and disruption of nucleosomal phasing around CTCF binding sites. This spatiotemporal analysis indicates that SMARCA5 is continuously required for maintaining nucleosomal spacing.

PAX3-FOXO1 coordinates enhancer architecture, eRNA transcription, and RNA polymerase pause release at select gene targets.

Transcriptional control is a highly dynamic process that changes rapidly in response to various cellular and extracellular cues, making it difficult to define the mechanism of transcription factor function using slow genetic methods. We used a chemical-genetic approach to rapidly degrade a canonical transcriptional activator, PAX3-FOXO1, to define the mechanism by which it regulates gene expression programs. By coupling rapid protein degradation with the analysis of nascent transcription over short time courses and integrating CUT&RUN, ATAC-seq, and eRNA analysis with deep proteomic analysis, we defined PAX3-FOXO1 function at a small network of direct transcriptional targets. PAX3-FOXO1 degradation impaired RNA polymerase pause release and transcription elongation at most regulated gene targets. Moreover, the activity of PAX3-FOXO1 at enhancers controlling this core network was surprisingly selective, affecting single elements in super-enhancers. This combinatorial analysis indicated that PAX3-FOXO1 was continuously required to maintain chromatin accessibility and enhancer architecture at regulated enhancers.

MTG16 regulates colonic epithelial differentiation, colitis, and tumorigenesis by repressing E protein transcription factors.

Aberrant epithelial differentiation and regeneration contribute to colon pathologies, including inflammatory bowel disease (IBD) and colitis-associated cancer (CAC). Myeloid translocation gene 16 (MTG16, also known as CBFA2T3) is a transcriptional corepressor expressed in the colonic epithelium. MTG16 deficiency in mice exacerbates colitis and increases tumor burden in CAC, though the underlying mechanisms remain unclear. Here, we identified MTG16 as a central mediator of epithelial differentiation, promoting goblet and restraining enteroendocrine cell development in homeostasis and enabling regeneration following dextran sulfate sodium-induced (DSS-induced) colitis. Transcriptomic analyses implicated increased Ephrussi box-binding transcription factor (E protein) activity in MTG16-deficient colon crypts. Using a mouse model with a point mutation that attenuates MTG16:E protein interactions (Mtg16P209T), we showed that MTG16 exerts control over colonic epithelial differentiation and regeneration by repressing E protein-mediated transcription. Mimicking murine colitis, MTG16 expression was increased in biopsies from patients with active IBD compared with unaffected controls. Finally, uncoupling MTG16:E protein interactions partially phenocopied the enhanced tumorigenicity of Mtg16-/- colon in the azoxymethane/DSS-induced model of CAC, indicating that MTG16 protects from tumorigenesis through additional mechanisms. Collectively, our results demonstrate that MTG16, via its repression of E protein targets, is a key regulator of cell fate decisions during colon homeostasis, colitis, and cancer.

A protocol for rapid degradation of endogenous transcription factors in mammalian cells and identification of direct regulatory targets.

Transcriptional changes happen within minutes; however, RNAi or genetic deletion requires days to weeks before transcription networks can be analyzed. This limitation has made it challenging to distinguish direct from indirect targets of sequence-specific transcription factors. This inability to define direct transcriptional targets hinders detailed studies of transcriptional mechanisms. This protocol combines rapid degradation of endogenous transcription factors with nascent transcript analysis to define the earliest, and likely direct, regulatory targets of transcription factors. For complete details on the use and execution of this protocol, please refer to Stengel et al., 2021).

Definition of a small core transcriptional circuit regulated by AML1-ETO.

Transcription factors regulate gene networks controlling normal hematopoiesis and are frequently deregulated in acute myeloid leukemia (AML). Critical to our understanding of the mechanism of cellular transformation by oncogenic transcription factors is the ability to define their direct gene targets. However, gene network cascades can change within minutes to hours, making it difficult to distinguish direct from secondary or compensatory transcriptional changes by traditional methodologies. To overcome this limitation, we devised cell models in which the AML1-ETO protein could be quickly degraded upon addition of a small molecule. The rapid kinetics of AML1-ETO removal, when combined with analysis of transcriptional output by nascent transcript analysis and genome-wide AML1-ETO binding by CUT&RUN, enabled the identification of direct gene targets that constitute a core AML1-ETO regulatory network. Moreover, derepression of this gene network was associated with RUNX1 DNA binding and triggered a transcription cascade ultimately resulting in myeloid differentiation.

BET Inhibition Enhances the Antileukemic Activity of Low-dose Venetoclax in Acute Myeloid Leukemia.

PURPOSE: The BCL2 inhibitor, venetoclax, has transformed clinical care in acute myeloid leukemia (AML). However, subsets of patients do not respond or eventually acquire resistance. Venetoclax-based regimens can lead to considerable marrow suppression in some patients. Bromodomain and extraterminal inhibitors (BETi) are potential treatments for AML, as regulators of critical AML oncogenes. We tested the efficacy of novel BET inhibitor INCB054329, and its synergy with venetoclax to reduce AML without induction of hematopoietic toxicity. EXPERIMENTAL DESIGN: INCB054329 efficacy was assessed by changes in cell cycle and apoptosis in treated AML cell lines. In vivo efficacy was assessed by tumor reduction in MV-4-11 cell line-derived xenografts. Precision run-on and sequencing (PRO-seq) evaluated effects of INCB054329. Synergy between low-dose BETi and venetoclax was assessed in cell lines and patient samples in vitro and in vivo while efficacy and toxicity was assessed in patient-derived xenograft (PDX) models. RESULTS: INCB054329 induced dose-dependent apoptosis and quiescence in AML cell lines. PRO-seq analysis evaluated the effects of INCB054329 on transcription and confirmed reduced transcriptional elongation of key oncogenes, MYC and BCL2, and genes involved in the cell cycle and metabolism. Combinations of BETi and venetoclax led to reduced cell viability in cell lines and patient samples. Low-dose combinations of INCB054329 and venetoclax in cell line and PDX models reduced AML burden, regardless of the sensitivity to monotherapy without development of toxicity. CONCLUSIONS: Our findings suggest low dose combinations of venetoclax and BETi may be more efficacious for patients with AML than either monotherapy, potentially providing a longer, more tolerable dosing regimen.

Nascent transcript and single-cell RNA-seq analysis defines the mechanism of action of the LSD1 inhibitor INCB059872 in myeloid leukemia.

Drugs targeting chromatin-modifying enzymes have entered clinical trials for myeloid malignancies, including INCB059872, a selective irreversible inhibitor of Lysine-Specific Demethylase 1 (LSD1). While initial studies of LSD1 inhibitors suggested these compounds may be used to induce differentiation of acute myeloid leukemia (AML), the mechanisms underlying this effect and dose-limiting toxicities are not well understood. Here, we used precision nuclear run-on sequencing (PRO-seq) and ChIP-seq in AML cell lines to probe for the earliest regulatory events associated with INCB059872 treatment. The changes in nascent transcription could be traced back to a loss of CoREST activity and activation of GFI1-regulated genes. INCB059872 is in phase I clinical trials, and we evaluated a pre-treatment bone marrow sample of a patient who showed a clinical response to INCB059872 while being treated with azacitidine. We used single-cell RNA-sequencing (scRNA-seq) to show that INCB059872 caused a shift in gene expression that was again associated with GFI1/GFI1B regulation. Finally, we treated mice with INCB059872 and performed scRNA-seq of lineage-negative bone marrow cells, which showed that INCB059872 triggered accumulation of megakaryocyte early progenitor cells with gene expression hallmarks of stem cells. Accumulation of these stem/progenitor cells may contribute to the thrombocytopenia observed in patients treated with LSD1 inhibitors.

Targeting MYCN-expressing triple-negative breast cancer with BET and MEK inhibitors.

Triple-negative breast cancer (TNBC) is an aggressive form of breast cancer that does not respond to endocrine therapy or human epidermal growth factor receptor 2 (HER2)-targeted therapies. Individuals with TNBC experience higher rates of relapse and shorter overall survival compared to patients with receptor-positive breast cancer subtypes. Preclinical discoveries are needed to identify, develop, and advance new drug targets to improve outcomes for patients with TNBC. Here, we report that MYCN, an oncogene typically overexpressed in tumors of the nervous system or with neuroendocrine features, is heterogeneously expressed within a substantial fraction of primary and recurrent TNBC and is expressed in an even higher fraction of TNBCs that do not display a pathological complete response after neoadjuvant chemotherapy. We performed high-throughput chemical screens on TNBC cell lines with varying amounts of MYCN expression and determined that cells with higher expression of MYCN were more sensitive to bromodomain and extraterminal motif (BET) inhibitors. Combined BET and MEK inhibition resulted in a synergistic decrease in viability, both in vitro and in vivo, using cell lines and patient-derived xenograft (PDX) models. Our preclinical data provide a rationale to advance a combination of BET and MEK inhibitors to clinical investigation for patients with advanced MYCN-expressing TNBC.

Histone deacetylase 3 controls a transcriptional network required for B cell maturation.

Histone deacetylase 3 (Hdac3) is a target of the FDA approved HDAC inhibitors, which are used for the treatment of lymphoid malignancies. Here, we used Cd19-Cre to conditionally delete Hdac3 to define its role in germinal center B cells, which represent the cell of origin for many B cell malignancies. Cd19-Cre-Hdac3-/- mice showed impaired germinal center formation along with a defect in plasmablast production. Analysis of Hdac3-/- germinal centers revealed a reduction in dark zone centroblasts and accumulation of light zone centrocytes. RNA-seq revealed a significant correlation between genes up-regulated upon Hdac3 loss and those up-regulated in Foxo1-deleted germinal center B cells, even though Foxo1 typically activates transcription. Therefore, to determine whether gene expression changes observed in Hdac3-/- germinal centers were a result of direct effects of Hdac3 deacetylase activity, we used an HDAC3 selective inhibitor and examined nascent transcription in germinal center-derived cell lines. Transcriptional changes upon HDAC3 inhibition were enriched for light zone gene signatures as observed in germinal centers. Further comparison of PRO-seq data with ChIP-seq/exo data for BCL6, SMRT, FOXO1 and H3K27ac identified direct targets of HDAC3 function including CD86, CD83 and CXCR5 that are likely responsible for driving the light zone phenotype observed in vivo.

Correction: Kaiso is required for MTG16-dependent effects on colitis-associated carcinoma.

In the original version of this article the authors noted that the GEO accession number for the relevant dataset was listed incorrectly as GSE12454.

Kaiso is required for MTG16-dependent effects on colitis-associated carcinoma.

The myeloid translocation gene family member MTG16 is a transcriptional corepressor that relies on the DNA-binding ability of other proteins to determine specificity. One such protein is the ZBTB family member Kaiso, and the MTG16:Kaiso interaction is necessary for repression of Kaiso target genes, such as matrix metalloproteinase-7. Using the azoxymethane and dextran sodium sulfate (AOM/DSS) murine model of colitis-associated carcinoma, we previously determined that MTG16 loss accelerates tumorigenesis and inflammation. However, it was unknown whether this effect was modified by Kaiso-dependent transcriptional repression. To test for a genetic interaction between MTG16 and Kaiso in inflammatory carcinogenesis, we subjected single and double knockout (DKO) mice to the AOM/DSS protocol. Mtg16-/- mice demonstrated increased colitis and tumor burden; in contrast, disease severity in Kaiso-/- mice was equivalent to wild-type controls. Surprisingly, Kaiso deficiency in the context of MTG16 loss reversed injury and pro-tumorigenic responses in the intestinal epithelium following AOM/DSS treatment, and tumor numbers were returned to near to wild-type levels. Transcriptomic analysis of non-tumor colon tissue demonstrated that changes induced by MTG16 loss were widely mitigated by concurrent Kaiso loss, and DKO mice demonstrated downregulation of metabolism and cytokine-associated gene sets with concurrent activation of DNA damage checkpoint pathways as compared with Mtg16-/-. Further, Kaiso knockdown in intestinal enteroids reduced stem- and WNT-associated phenotypes, thus abrogating the induction of these pathways observed in Mtg16-/- samples. Together, these data suggest that Kaiso modifies MTG16-driven inflammation and tumorigenesis and suggests that Kaiso deregulation contributes to MTG16-dependent colitis and CAC phenotypes.

Displacement of WDR5 from Chromatin by a WIN Site Inhibitor with Picomolar Affinity.

The chromatin-associated protein WDR5 is a promising target for pharmacological inhibition in cancer. Drug discovery efforts center on the blockade of the "WIN site" of WDR5, a well-defined pocket that is amenable to small molecule inhibition. Various cancer contexts have been proposed to be targets for WIN site inhibitors, but a lack of understanding of WDR5 target genes and of the primary effects of WIN site inhibitors hampers their utility. Here, by the discovery of potent WIN site inhibitors, we demonstrate that the WIN site links WDR5 to chromatin at a small cohort of loci, including a specific subset of ribosome protein genes. WIN site inhibitors rapidly displace WDR5 from chromatin and decrease the expression of associated genes, causing translational inhibition, nucleolar stress, and p53 induction. Our studies define a mode by which WDR5 engages chromatin and forecast that WIN site blockade could have utility against multiple cancer types.

The CDK7 inhibitor THZ1 alters RNA polymerase dynamics at the 5' and 3' ends of genes.

The t(8;21) is one of the most frequent chromosomal translocations associated with acute myeloid leukemia (AML). We found that t(8;21) AML were extremely sensitive to THZ1, which triggered apoptosis after only 4 h. We used precision nuclear run-on transcription sequencing (PROseq) to define the global effects of THZ1 and other CDK inhibitors on RNA polymerase II dynamics. Inhibition of CDK7 using THZ1 caused wide-spread loss of promoter-proximal paused RNA polymerase. This loss of 5' pausing was associated with accumulation of polymerases in the body of a large number of genes. However, there were modest effects on genes regulated by 'super-enhancers'. At the 3' ends of genes, treatment with THZ1 suppressed RNA polymerase 'read through' at the end of the last exon, which resembled a phenotype associated with a mutant RNA polymerase with slower elongation rates. Consistent with this hypothesis, polyA site-sequencing (PolyA-seq) did not detect differences in poly A sites after THZ1 treatment. PROseq analysis after short treatments with THZ1 suggested that these 3' effects were due to altered CDK7 activity at the 5' end of long genes, and were likely to be due to slower rates of elongation.

BET inhibitors reduce cell size and induce reversible cell cycle arrest in AML.

Inhibitors of the bromodomain and extraterminal domain family (BETi) offer a new approach to treat hematological malignancies, with leukemias containing mixed lineage leukemia rearrangements being especially sensitive due to a reliance on the regulation of transcription elongation. We explored the mechanism of action of BETi in cells expressing the t(8;21), and show that these compounds reduced the size of acute myeloid leukemia cells, triggered a rapid but reversible G0 /G1 arrest, and with time, cause cell death. Meta-analysis of PRO-seq data identified ribosomal genes, which are regulated by MYC, were downregulated within 3 hours of addition of the BETi. This reduction of MYC regulated metabolic genes coincided with the loss of mitochondrial respiration and large reductions in the glycolytic rate. In addition, gene expression analysis showed that transcription of BCL2 was rapidly affected by BETi but this did not cause dramatic increases in cell death. Cell cycle arrest, lowered metabolic activity, and reduced BCL2 levels suggested that a second compound was needed to push these cells over the apoptotic threshold. Indeed, low doses of the BCL2 inhibitor, venetoclax, in combination with the BETi was a potent combination in t(8;21) containing cells. Thus, BET inhibitors that affect MYC and BCL2 expression should be considered for combination therapy with venetoclax.

Nascent RNA sequencing analysis provides insights into enhancer-mediated gene regulation.

BACKGROUND: Enhancers are distal cis-regulatory elements that control gene expression. Despite an increasing appreciation of the importance of enhancers in cellular function and disease, our knowledge of enhancer-regulated transcription is very limited. Nascent RNA sequencing technologies, such as global nuclear run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), not only provide a direct and reliable measurement of enhancer activity, but also allow for quantifying transcription of enhancers and target genes simultaneously, making these technologies extremely useful for exploring enhancer-mediated regulation. RESULTS: Nascent RNA sequencing analysis (NRSA) provides a comprehensive view of enhancer-mediated gene regulation. NRSA not only outperforms existing methods for enhancer identification, but also enables annotation and quantification of active enhancers, and prediction of their target genes. Furthermore, NRSA identifies functionally important enhancers by integrating 1) nascent transcriptional changes in enhancers and their target genes and 2) binding profiles from regulator(s) of interest. Applied to wildtype and histone deacetylase 3 (Hdac3) knockout mouse livers, NRSA showed that HDAC3 regulates RNA polymerase recruitment through both proximal (promoter) and distal (enhancer) regulatory elements. Integrating ChIP-seq with PRO-seq data, NRSA prioritized enhancers based on their potential contribution to mediating HDAC3 regulation. CONCLUSIONS: NRSA will greatly facilitate the usage of nascent RNA sequencing techniques and accelerate the study of enhancer-mediated regulation.