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.

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.

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.

An Old Dog Has a New Trick: Somatic Exonic Deletions in RUNX1 Are Frequent in AML.

Somatic loss-of-function RUNX1 mutations in acute myeloid leukemia (AML) include missense, nonsense, and frameshift mutations, whereas germline RUNX1 variants in RUNX1-FPDMM also include large exonic deletions. Alternative variant detection approaches revealed that large exonic deletions in RUNX1 are also common in sporadic AML, which has implications for patient stratification and therapeutic decision-making. See related article by Eriksson et al., p. 2826.

Emerging insights into enhancer biology and function.

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

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.

Over-Generalizing About GC (Hypoxia): Pitfalls of Limiting Breadth of Experimental Systems and Analyses in Framing Informatics Conclusions.

Accumulating evidence suggests that many immune responses are influenced by local nutrient concentrations in addition to the programming of intermediary metabolism within immune cells. Humoral immunity and germinal centers (GC) are settings in which these factors are under active investigation. Hypoxia is an example of how a particular nutrient is distributed in lymphoid follicles during an antibody response, and how oxygen sensors may impact the qualities of antibody output after immunization. Using exclusively a bio-informatic analysis of mRNA levels in GC and other B cells, recent work challenged the concept that there is any hypoxia or that it has any influence. To explore this proposition, we performed new analyses of published genomics data, explored potential sources of disparity, and elucidated aspects of the apparently conflicting conclusions. Specifically, replicability and variance among data sets derived from different naïve as well as GC B cells were considered. The results highlight broader issues that merit consideration, especially at a time of heightened focus on scientific reports in the realm of immunity and antibody responses. Based on these analyses, a standard is proposed under which the relationship of new data sets should be compared to prior "fingerprints" of cell types and reported transparently to referees and readers. In light of independent evidence of diversity within and among GC elicited by protein immunization, avoidance of overly broad conclusions about germinal centers in general when experimental systems are subject to substantial constraints imposed by technical features also is warranted.

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.

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.

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.

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.

Deacetylase activity of histone deacetylase 3 is required for productive VDJ recombination and B-cell development.

Histone deacetylase 3 (HDAC3) is the catalytic component of NCoR/SMRT corepressor complexes that mediate the actions of transcription factors implicated in the regulation of B-cell development and function. We crossed Hdac3 conditional knockout mice with Mb1-Cre knockin animals to delete Hdac3 in early progenitor B cells. The spleens of Hdac3F/-Mb1-Cre+/- mice were virtually devoid of mature B cells, and B220+CD43+ B-cell progenitors accumulated within the bone marrow. Quantitative deep sequencing of the Ig heavy chain locus from B220+CD43+ populations identified a defect in VHDJH recombination with a severe reduction in productive rearrangements, which directly corresponded to the loss of pre-B cells from Hdac3Δ/- bone marrow. For Hdac3Δ/- B cells that did show productive VDJ rearrangement, there was significant skewing toward the incorporation of proximal VH gene segments and a corresponding reduction in distal VH gene segment use. Although transcriptional effects within these loci were modest, Hdac3Δ/- progenitor cells displayed global changes in chromatin structure that likely hindered effective distal V-DJ recombination. Reintroduction of wild-type Hdac3 restored normal B-cell development, whereas an Hdac3 point mutant lacking deacetylase activity failed to complement this defect. Thus, the deacetylase activity of Hdac3 is required for the generation of mature B cells.

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/.

CREBBP Inactivation Promotes the Development of HDAC3-Dependent Lymphomas.

Somatic mutations in CREBBP occur frequently in B-cell lymphoma. Here, we show that loss of CREBBP facilitates the development of germinal center (GC)-derived lymphomas in mice. In both human and murine lymphomas, CREBBP loss-of-function resulted in focal depletion of enhancer H3K27 acetylation and aberrant transcriptional silencing of genes that regulate B-cell signaling and immune responses, including class II MHC. Mechanistically, CREBBP-regulated enhancers are counter-regulated by the BCL6 transcriptional repressor in a complex with SMRT and HDAC3, which we found to bind extensively to MHC class II loci. HDAC3 loss-of-function rescued repression of these enhancers and corresponding genes, including MHC class II, and more profoundly suppressed CREBBP-mutant lymphomas in vitro and in vivo Hence, CREBBP loss-of-function contributes to lymphomagenesis by enabling unopposed suppression of enhancers by BCL6/SMRT/HDAC3 complexes, suggesting HDAC3-targeted therapy as a precision approach for CREBBP-mutant lymphomas. SIGNIFICANCE: Our findings establish the tumor suppressor function of CREBBP in GC lymphomas in which CREBBP mutations disable acetylation and result in unopposed deacetylation by BCL6/SMRT/HDAC3 complexes at enhancers of B-cell signaling and immune response genes. Hence, inhibition of HDAC3 can restore the enhancer histone acetylation and may serve as a targeted therapy for CREBBP-mutant lymphomas. Cancer Discov; 7(1); 38-53. ©2016 AACR.See related commentary by Höpken, p. 14This article is highlighted in the In This Issue feature, p. 1.