Search

RNA surveillance pathways detect and degrade defective transcripts to ensure RNA fidelity. We found that disrupted nuclear RNA surveillance is oncogenic. Cyclin-dependent kinase 13 ( CDK13 ) is mutated in melanoma, and patient-mutated CDK13 accelerates zebrafish melanoma. CDK13 mutation causes aberrant RNA stabilization. CDK13 is required for ZC3H14 phosphorylation, which is necessary and sufficient to promote nuclear RNA degradation. Mutant CDK13 fails to activate nuclear RNA surveillance, causing aberrant protein-coding transcripts to be stabilized and translated. Forced aberrant RNA expression accelerates melanoma in zebrafish. We found recurrent mutations in genes encoding nuclear RNA surveillance components in many malignancies, establishing nuclear RNA surveillance as a tumor-suppressive pathway. Activating nuclear RNA surveillance is crucial to avoid accumulation of aberrant RNAs and their ensuing consequences in development and disease.

Supplementary Data from EP300 Selectively Controls the Enhancer Landscape of MYCN-Amplified Neuroblastoma

Gene expression is regulated by promoters and enhancers marked by histone H3 lysine 27 acetylation (H3K27ac), which is established by the paralogous histone acetyltransferases (HAT) EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma (NB) depends on EP300, whereas CBP has a limited role. EP300 controls enhancer acetylation by interacting with TFAP2β, a transcription factor member of the lineage-defining transcriptional core regulatory circuitry (CRC) in NB. To disrupt EP300, we developed a proteolysis-targeting chimera (PROTAC) compound termed “JQAD1” that selectively targets EP300 for degradation. JQAD1 treatment causes loss of H3K27ac at CRC enhancers and rapid NB apoptosis, with limited toxicity to untransformed cells where CBP may compensate. Furthermore, JQAD1 activity is critically determined by cereblon (CRBN) expression across NB cells.Significance:EP300, but not CBP, controls oncogenic CRC-driven transcription in high-risk NB by binding TFAP2β. We developed JQAD1, a CRBN-dependent PROTAC degrader with preferential activity against EP300 and demonstrated its activity in NB. JQAD1 has limited toxicity to untransformed cells and is effective in vivo in a CRBN-dependent manner.This article is highlighted in the In This Issue feature, p. 587

Table S4

Supplementary Figures 1-10 and corresponding legends

Supplementary Data

ARACNe-AP transcriptional network from gene expression profiles

High-risk neuroblastomas show a paucity of recurrent somatic mutations at diagnosis. As a result, the molecular basis for this aggressive phenotype remains elusive. Recent progress in regulatory network analysis helped us elucidate disease-driving mechanisms downstream of genomic alterations, including recurrent chromosomal alterations. Our analysis identified three molecular subtypes of high-risk neuroblastomas, consistent with chromosomal alterations, and identified subtype-specific master regulator proteins that were conserved across independent cohorts. A 10-protein transcriptional module—centered around a TEAD4–MYCN positive feedback loop—emerged as the regulatory driver of the high-risk subtype associated with MYCN amplification. Silencing of either gene collapsed MYCN-amplified (MYCNAmp) neuroblastoma transcriptional hallmarks and abrogated viability in vitro and in vivo. Consistently, TEAD4 emerged as a robust prognostic marker of poor survival, with activity independent of the canonical Hippo pathway transcriptional coactivators YAP and TAZ. These results suggest novel therapeutic strategies for the large subset of MYCN-deregulated neuroblastomas.Significance: Despite progress in understanding of neuroblastoma genetics, little progress has been made toward personalized treatment. Here, we present a framework to determine the downstream effectors of the genetic alterations sustaining neuroblastoma subtypes, which can be easily extended to other tumor types. We show the critical effect of disrupting a 10-protein module centered around a YAP/TAZ-independent TEAD4–MYCN positive feedback loop in MYCNAmp neuroblastomas, nominating TEAD4 as a novel candidate for therapeutic intervention. Cancer Discov; 8(5); 582–99. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 517

Supplementary Experimental Procedures; Supplementary Figures S1-S10; Supplementary Table Legends S1-S9; Supplementary Tables S8-S9; Additional References

Anaplastic large cell lymphoma (ALCL) is an aggressive, CD30+ T-cell lymphoma of children and adults. ALK fusion transcripts or mutations in the JAK-STAT pathway are observed in most ALCL tumors, but the mechanisms underlying tumorigenesis are not fully understood. Here we show that dysregulated STAT3, together with a core transcriptional regulatory circuit consisting of BATF3–IRF4– IKZF1, co-occupies gene enhancers to establish an oncogenic transcription program and maintain the malignant state of ALCL. Critical downstream targets of this network in ALCL cells include the proto-oncogene MYC, which requires active STAT3 to facilitate high levels of MYC transcription. The activity of this auto-regulatory transcription loop is reinforced by MYC binding to the enhancer regions associated with STAT3 and each of the core regulatory transcription factors. These findings provide new insights for understanding how dysregulated signaling pathways hijack cell-type-specific transcriptional machinery to drive tumorigenesis and create therapeutic vulnerabilities in genetically defined tumors.

The bromodomain and extraterminal motif (BET) proteins are critical drug targets for diseases. The precise functions and relationship of BRD2 with other BET proteins remain elusive mechanistically. Here, we used acute protein degradation and quantitative genomic and proteomic approaches to investigate the primary functions of BRD2 in transcription. We report that BRD2 is required for TAF3-mediated Pol II initiation at promoters with low levels of H3K4me3 and for R-loop suppression during Pol II elongation. Single and double depletion revealed that BRD2 and BRD3 function additively, independently, or perhaps antagonistically in Pol II transcription at different promoters. Furthermore, we found that BRD2 regulates the expression of different genes during embryonic body differentiation processes by promoter priming in embryonic stem cells. Therefore, our results suggest complex interconnections between BRD2 and BRD3 at promoters to fine-tune Pol II initiation and elongation for control of cell state.

Metabolic pathway regulates cell fate Lineage-specific regulators direct cell fate decisions, but the precise mechanisms are not well known. Using an in vivo chemical suppressor screen of a bloodless zebrafish mutant, Rossmann et al. show that the lineage-specific chromatin factor tif1γ directly regulates mitochondrial genes to drive red blood cell differentiation. Loss of tif1γ reduces coenzyme Q synthesis and function, impeding mitochondrial respiration and leading to epigenetic alterations and repression of erythropoiesis. The loss of blood in the mutant fish can be rescued by the addition of coenzyme Q. This work establishes a mechanism by which a chromatin factor tunes a metabolic pathway in a tissue-specific manner. Science , this issue p. 716

SummaryThe capacity to generate functional hepatocytes from renewable human pluripotent stem cells (hPSCs) could address limited supplies of primary human hepatocytes. However, hepatocytes differentiated from hPSCs in vitro are functionally immature. To understand mechanisms regulating maturation of in vitro derived hepatocytes, we developed a 3D spheroid differentiation system and compared gene regulatory elements in uncultured human primary hepatocytes with those in hepatocytes that were differentiated in 2D or 3D conditions from human PSCs by RNA-seq, ATAC-seq, and H3K27Ac ChIP-seq. Three-dimensional differentiation improved enhancer activity and expression of transcription factor ONECUT1, but was insufficient to upregulate human-specific mature hepatocytes marker gene CYP3A4 or super-enhancer regulated transcription factor gene NFIC. Regulome comparisons showed reduced enrichment of thyroid receptor THRB motifs in accessible chromatin and in active enhancers without reduced transcription of THRB, suggesting the regulation at the level of THRB ligands in PSC-differentiated hepatocytes. Addition of thyroid hormone T3 to the PSC-differentiated hepatocytes increased CYP3A4 expression. T3 increased binding of THRB to the CYP3A4 proximal enhancer and restored the super-enhancer status and gene expression of NFIC and reduced expression of AFP. The resultant hPSC-hepatocytes showed gene expression, epigenetic status and super-enhancer landscape closer to primary hepatocytes and activated regulatory regions including non-coding SNPs associated with liver-related diseases. Transplanting the 3D PSC-hepatocytes into immunocompromised mice resulted in engraftment of human hepatocytes in the mouse liver parenchyma without disrupting normal liver histology at 6 months after transplantation. This work provides insights into the functions of nuclear receptor THRB and highlights the importance of the environmental factors-nuclear receptors axis in regulating maturation of human PSC-differentiated cell types.

Genome-wide association studies identify genomic variants associated with human traits and diseases. Most trait-associated variants are located within cell-type-specific enhancers, but the molecular mechanisms governing phenotypic variation are less well understood. Here, we show that many enhancer variants associated with red blood cell (RBC) traits map to enhancers that are co-bound by lineage-specific master transcription factors (MTFs) and signaling transcription factors (STFs) responsive to extracellular signals. The majority of enhancer variants reside on STF and not MTF motifs, perturbing DNA binding by various STFs (BMP/TGF-β-directed SMADs or WNT-induced TCFs) and affecting target gene expression. Analyses of engineered human blood cells and expression quantitative trait loci verify that disrupted STF binding leads to altered gene expression. Our results propose that the majority of the RBC-trait-associated variants that reside on transcription-factor-binding sequences fall in STF target sequences, suggesting that the phenotypic variation of RBC traits could stem from altered responsiveness to extracellular stimuli.

Significance LIN28B is well known as a RNA-binding protein and a suppressor of microRNA biogenesis by selectively blocking the processing of let-7 precursors. However, little is known about let-7–independent roles of LIN28B. Here, we show that LIN28B is recruited to active promoters by binding to the zinc-finger transcription factor ZNF143. LIN28B acts as a cofactor to upregulate expression of a subset of downstream target genes that are essential for neuroblastoma cell survival and migration. Our paper reveals an unexpected role of LIN28B in transcriptional regulation that is independent of let-7 during neuroblastoma pathogenesis., LIN28B is highly expressed in neuroblastoma and promotes tumorigenesis, at least, in part, through inhibition of let-7 microRNA biogenesis. Here, we report that overexpression of either wild-type (WT) LIN28B or a LIN28B mutant that is unable to inhibit let-7 processing increases the penetrance of MYCN-induced neuroblastoma, potentiates the invasion and migration of transformed sympathetic neuroblasts, and drives distant metastases in vivo. Genome-wide chromatin immunoprecipitation coupled with massively parallel DNA sequencing (ChIP-seq) and coimmunoprecipitation experiments show that LIN28B binds active gene promoters in neuroblastoma cells through protein–protein interaction with the sequence-specific zinc-finger transcription factor ZNF143 and activates the expression of downstream targets, including transcription factors forming the adrenergic core regulatory circuitry that controls the malignant cell state in neuroblastoma as well as GSK3B and L1CAM that are involved in neuronal cell adhesion and migration. These findings reveal an unexpected let-7–independent function of LIN28B in transcriptional regulation during neuroblastoma pathogenesis.

The availability of nucleotides has a direct impact on transcription. The inhibition of dihydroorotate dehydrogenase (DHODH) with leflunomide impacts nucleotide pools by reducing pyrimidine levels. Leflunomide abrogates the effective transcription elongation of genes required for neural crest development and melanoma growth in vivo1. To define the mechanism of action, we undertook an in vivo chemical suppressor screen for restoration of neural crest after leflunomide treatment. Surprisingly, we found that alterations in progesterone and progesterone receptor (Pgr) signalling strongly suppressed leflunomide-mediated neural crest effects in zebrafish. In addition, progesterone bypasses the transcriptional elongation block resulting from Paf complex deficiency, rescuing neural crest defects in ctr9 morphant and paf1(alnz24) mutant embryos. Using proteomics, we found that Pgr binds the RNA helicase protein Ddx21. ddx21-deficient zebrafish show resistance to leflunomide-induced stress. At a molecular level, nucleotide depletion reduced the chromatin occupancy of DDX21 in human A375 melanoma cells. Nucleotide supplementation reversed the gene expression signature and DDX21 occupancy changes prompted by leflunomide. Together, our results show that DDX21 acts as a sensor and mediator of transcription during nucleotide stress. Santoriello, Sporrij et al. show that the progesterone receptor associates with RNA helicase DDX21 during nucleotide depletion, promotes its binding on chromatin and rescues efficient transcription in melanoma cells.

A heritable polymorphism within regulatory sequences of the LMO1 gene is associated with its elevated expression and increased susceptibility to develop neuroblastoma, but the oncogenic pathways downstream of the LMO1 transcriptional co-regulatory protein are unknown. Our ChIP-seq and RNA-seq analyses reveal that a key gene directly regulated by LMO1 and MYCN is ASCL1, which encodes a basic helix-loop-helix transcription factor. Regulatory elements controlling ASCL1 expression are bound by LMO1, MYCN and the transcription factors GATA3, HAND2, PHOX2B, TBX2 and ISL1—all members of the adrenergic (ADRN) neuroblastoma core regulatory circuitry (CRC). ASCL1 is required for neuroblastoma cell growth and arrest of differentiation. ASCL1 and LMO1 directly regulate the expression of CRC genes, indicating that ASCL1 is a member and LMO1 is a coregulator of the ADRN neuroblastoma CRC., Polymorphisms in LMO1 are associated with increased susceptibility to develop neuroblastoma. Here, the authors show that LMO1 directly induces the transcription factor ASCL1, which regulates the differentiation of neurons, demonstrating that ASCL1 is part of the adrenergic neuroblastoma core regulatory circuit.

Gene expression is regulated by promoters and enhancers marked by histone H3 lysine 27 acetylation (H3K27ac), which is established by the paralogous, multidomain-containing histone acetyltransferases (HAT) EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma (NB) depends on EP300, whereas CBP has a limited role. To disrupt EP300, we developed a proteolysis-targeting chimera (PROTAC) compound termed “JQAD1” that selectively targets EP300 for degradation. JQAD1 treatment causes loss of H3K27ac at the enhancers that regulate NB master transcription factors, and drives rapid NB apoptosis, with limited toxicity to untransformed cells where CBP may compensate. In parallel, we demonstrate that EP300 and CBP have subdomain-specific functions, with enriched activity for individual domains in distinct tumor lineages. These data provide a foundation for interrogation of these histone acetyltransferases in distinct tumor types and new strategies for therapeutic disruption of subdomain and scaffolding activities of these multidomain proteins in cancer. Citation Format: Adam D. Durbin, Noha Shendy, Audrey Mercier, Yang Zhang, Melissa J. Bikowitz, Logan H. Sigua, Sarah Robinson, Tingjian Wang, Barbara Jonchere, A. Thomas Look, Mark W. Zimmerman, Martine Roussel, Brian J. Abraham, Ernst Schonbrunn, Kimberly Stegmaier, Jun Qi. Targeting EP300 and CBP for therapeutic benefit in pediatric solid tumors. [abstract]. In: Proceedings of the AACR Special Conference: Cancer Epigenomics; 2022 Oct 6-8; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_2):Abstract nr PR005.

Genome-wide protein-DNA binding is popularly assessed using specific antibody pulldown in Chromatin Immunoprecipitation Sequencing (ChIP-Seq) or Cleavage Under Targets and Release Using Nuclease (CUTRUN) sequencing experiments. These technologies generate high-throughput sequencing data that necessitate the use of multiple sophisticated, computationally intensive genomic tools to make discoveries, but these genomic tools often have a high barrier to use because of computational resource constraints.We present a comprehensive, infrastructure-independent, computational pipeline called SEAseq, which leverages field-standard, open-source tools for processing and analyzing ChIP-Seq/CUTRUN data. SEAseq performs extensive analyses from the raw output of the experiment, including alignment, peak calling, motif analysis, promoters and metagene coverage profiling, peak annotation distribution, clustered/stitched peaks (e.g. super-enhancer) identification, and multiple relevant quality assessment metrics, as well as automatic interfacing with data in GEO/SRA. SEAseq enables rapid and cost-effective resource for analysis of both new and publicly available datasets as demonstrated in our comparative case studies.The easy-to-use and versatile design of SEAseq makes it a reliable and efficient resource for ensuring high quality analysis. Its cloud implementation enables a broad suite of analyses in environments with constrained computational resources. SEAseq is platform-independent and is aimed to be usable by everyone with or without programming skills. It is available on the cloud at https://platform.stjude.cloud/workflows/seaseq and can be locally installed from the repository at https://github.com/stjude/seaseq .

Description, The CaCTS algorithm nominates cancer cell master transcription factors and guides a model of ovarian cancer regulatory circuitry., Critical developmental “master transcription factors” (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Simple Summary Neuroblastoma is a pediatric tumor originating from the sympathetic nervous system responsible for 10–15% of all childhood cancer deaths. Half of all neuroblastoma patients present with high-risk disease, of which nearly 50% relapse and die of their disease. In addition, standard therapies cause serious lifelong side effects and increased risk for secondary tumors. Further research is crucial to better understand the molecular basis of neuroblastomas and to identify novel druggable targets. Neuroblastoma tumorigenesis has to this end been modeled in both mice and zebrafish. Here, we present a detailed dissection of the gene expression patterns that underlie tumor formation in the murine TH-MYCN-driven neuroblastoma model. We identified key factors that are putatively important for neuroblastoma tumor initiation versus tumor progression, pinpointed crucial regulators of the observed expression patterns during neuroblastoma development and scrutinized which factors could be innovative and vulnerable nodes for therapeutic intervention. Abstract Roughly half of all high-risk neuroblastoma patients present with MYCN amplification. The molecular consequences of MYCN overexpression in this aggressive pediatric tumor have been studied for decades, but thus far, our understanding of the early initiating steps of MYCN-driven tumor formation is still enigmatic. We performed a detailed transcriptome landscaping during murine TH-MYCN-driven neuroblastoma tumor formation at different time points. The neuroblastoma dependency factor MEIS2, together with ASCL1, was identified as a candidate tumor-initiating factor and shown to be a novel core regulatory circuit member in adrenergic neuroblastomas. Of further interest, we found a KEOPS complex member (gm6890), implicated in homologous double-strand break repair and telomere maintenance, to be strongly upregulated during tumor formation, as well as the checkpoint adaptor Claspin (CLSPN) and three chromosome 17q loci CBX2, GJC1 and LIMD2. Finally, cross-species master regulator analysis identified FOXM1, together with additional hubs controlling transcriptome profiles of MYCN-driven neuroblastoma. In conclusion, time-resolved transcriptome analysis of early hyperplastic lesions and full-blown MYCN-driven neuroblastomas yielded novel components implicated in both tumor initiation and maintenance, providing putative novel drug targets for MYCN-driven neuroblastoma.

The synthesis of pre-mRNA by RNA polymerase II (Pol II) involves the formation of a transcription initiation complex, and a transition to an elongation complex1–4. The large subunit of Pol II contains an intrinsically disordered C-terminal domain that is phosphorylated by cyclin-dependent kinases during the transition from initiation to elongation, thus influencing the interaction of the C-terminal domain with different components of the initiation or the RNA-splicing apparatus5,6. Recent observations suggest that this model provides only a partial picture of the effects of phosphorylation of the C-terminal domain7–12. Both the transcription-initiation machinery and the splicing machinery can form phase-separated condensates that contain large numbers of component molecules: hundreds of molecules of Pol II and mediator are concentrated in condensates at super-enhancers7,8, and large numbers of splicing factors are concentrated in nuclear speckles, some of which occur at highly active transcription sites9–12. Here we investigate whether the phosphorylation of the Pol II C-terminal domain regulates the incorporation of Pol II into phase-separated condensates that are associated with transcription initiation and splicing. We find that the hypophosphorylated C-terminal domain of Pol II is incorporated into mediator condensates and that phosphorylation by regulatory cyclin-dependent kinases reduces this incorporation. We also find that the hyperphosphorylated C-terminal domain is preferentially incorporated into condensates that are formed by splicing factors. These results suggest that phosphorylation of the Pol II C-terminal domain drives an exchange from condensates that are involved in transcription initiation to those that are involved in RNA processing, and implicates phosphorylation as a mechanism that regulates condensate preference. RNA polymerase II with a hypophosphorylated C-terminal domain preferentially incorporates into mediator condensates, and with a hyperphosphorylated C-terminal domain into splicing-factor condensates, revealing phosphorylation as a regulatory mechanism in condensate preference.

Gene regulation is critical for cell identity, and its dysregulation is a defining characteristic of common diseases including cancers. Although promoter activity is a strong predictor of gene expression, activities at other regulatory elements, including enhancers and super-enhancers (SE), are major contributors to gene regulation. For example, distinct transcription factor-regulating SEs define neuroblastoma (NB) subtypes with distinct clinical outcomes. However, technical limitations largely prevent ChIP-seq based enhancer/SE activity profiling in primary patient samples. We previously developed MethylationToActivity (M2A) and demonstrated that high-order DNA methylation (DNAm) features are strong predictors of promoter activity measured by ChIP-seq. This tool is important because genomewide DNAm assays are widely used in clinic to classify tumors and stratify patients on clinical trials. However, it is limited by its reliance on known gene position annotations. Here we present MethylationToRegulation (M2R), a method that utilizes a convolutional neural network (CNN)-based deep learning framework to infer intergenic enhancer/SE activities from DNA methylomes. We obtained paired WGBS and H3K27ac ChIP-seq data for 16 pediatric NB samples profiled in the Pediatric Cancer Genome Project. M2R was trained on 6 samples and tested on the remaining 10 samples. It achieved an average relative prediction accuracy of 84% on test samples when compared to the H3K27ac ChIP-seq replicate consistency from the ENCODE project. We adapted the Rank Ordering of Super-Enhancers algorithm to interpret H3K27ac signals inferred from DNA methylomes to identify SEs. M2R faithfully captured subtype-defining SEs with high specificity, including those associated with master NB transcription factors. Our results demonstrate that M2R accurately quantifies enhancer/SE activities and infers critical epigenetic marks from DNA methylomes. Application of M2R will enable the profiling of epigenetic dysregulation in patient tumor samples, seeking to improve clinical outcomes. Citation Format: Daniel K. Putnam, Brian J. Abraham, Xiang Chen. MethylationToRegulation: A deep-learning approach to infer chromatin properties from DNA methylomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1936.

Gene expression is regulated by promoters and enhancers marked by histone H3 lysine 27 acetylation (H3K27ac), which is established by the paralogous histone acetyltransferases (HAT) EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma (NB) depends on EP300, whereas CBP has a limited role. EP300 controls enhancer acetylation by interacting with TFAP2β, a transcription factor member of the lineage-defining transcriptional core regulatory circuitry (CRC) in NB. To disrupt EP300, we developed a proteolysis-targeting chimera (PROTAC) compound termed “JQAD1” that selectively targets EP300 for degradation. JQAD1 treatment causes loss of H3K27ac at CRC enhancers and rapid NB apoptosis, with limited toxicity to untransformed cells where CBP may compensate. Furthermore, JQAD1 activity is critically determined by cereblon (CRBN) expression across NB cells. Significance: EP300, but not CBP, controls oncogenic CRC-driven transcription in high-risk NB by binding TFAP2β. We developed JQAD1, a CRBN-dependent PROTAC degrader with preferential activity against EP300 and demonstrated its activity in NB. JQAD1 has limited toxicity to untransformed cells and is effective in vivo in a CRBN-dependent manner. This article is highlighted in the In This Issue feature, p. 587

Description, Retinoic acid reprograms neuroblastoma cell identity by remodeling the core transcriptional regulatory circuitry., Neuroblastoma cell identity depends on a core regulatory circuit (CRC) of transcription factors that collaborate with MYCN to drive the oncogenic gene expression program. For neuroblastomas dependent on the adrenergic CRC, treatment with retinoids can inhibit cell growth and induce differentiation. Here, we show that when MYCN-amplified neuroblastoma cells are treated with retinoic acid, histone H3K27 acetylation and methylation become redistributed to decommission super-enhancers driving the expression of PHOX2B and GATA3, together with the activation of new super-enhancers that drive high levels of MEIS1 and SOX4 expression. These findings indicate that treatment with retinoids can reprogram the enhancer landscape, resulting in down-regulation of MYCN expression, while establishing a new retino-sympathetic CRC that causes proliferative arrest and sympathetic differentiation. Thus, we provide mechanisms that account for the beneficial effects of retinoids in high-risk neuroblastoma and explain the rapid down-regulation of expression of MYCN despite massive levels of amplification of this gene.

ARID1A is a tumor suppressor in neuroblastoma, whose loss promotes the adrenergic-to-mesenchymal transition., Mutations in genes encoding SWI/SNF chromatin remodeling complexes are found in approximately 20% of all human cancers, with ARID1A being the most frequently mutated subunit. Here, we show that disruption of ARID1A homologs in a zebrafish model accelerates the onset and increases the penetrance of MYCN-driven neuroblastoma by increasing cell proliferation in the sympathoadrenal lineage. Depletion of ARID1A in human NGP neuroblastoma cells promoted the adrenergic-to-mesenchymal transition with changes in enhancer-mediated gene expression due to alterations in the genomic occupancies of distinct SWI/SNF assemblies, BAF and PBAF. Our findings indicate that ARID1A is a haploinsufficient tumor suppressor in MYCN-driven neuroblastoma, whose depletion enhances tumor development and promotes the emergence of the more drug-resistant mesenchymal cell state.

© 2020 Elsevier Inc. The cyclin-dependent kinase 1 (Cdk1) drives cell division. To uncover additional functions of Cdk1, we generated knockin mice expressing an analog-sensitive version of Cdk1 in place of wild-type Cdk1. In our study, we focused on embryonic stem cells (ESCs), because this cell type displays particularly high Cdk1 activity. We found that in ESCs, a large fraction of Cdk1 substrates is localized on chromatin. Cdk1 phosphorylates many proteins involved in epigenetic regulation, including writers and erasers of all major histone marks. Consistent with these findings, inhibition of Cdk1 altered histone-modification status of ESCs. High levels of Cdk1 in ESCs phosphorylate and partially inactivate Dot1l, the H3K79 methyltransferase responsible for placing activating marks on gene bodies. Decrease of Cdk1 activity during ESC differentiation de-represses Dot1l, thereby allowing coordinated expression of differentiation genes. These analyses indicate that Cdk1 functions to maintain the epigenetic identity of ESCs.

Cutaneous melanoma is an aggressive tumor that accounts for most skin cancer deaths. Among the physiological barriers against therapeutic success is a strong survival program driven by genes such as MITF that specify melanocyte identity, a phenomenon known in melanoma biology as lineage dependency. MITF overexpression is occasionally explained by gene amplification, but here we show that super-enhancers are also important determinants of MITF overexpression in some melanoma cell lines and tumors. Although compounds that directly inhibit MITF are unavailable, a covalent CDK7 inhibitor, THZ1, has recently been shown to potently suppress the growth of various cancers through the depletion of master transcription-regulating oncogenes and the disruption of their attendant super-enhancers. We also show that melanoma cells are highly sensitive to CDK7 inhibition both in vitro and in vivo and that THZ1 can dismantle the super-enhancer apparatus at MITF and SOX10 in some cell lines, thereby extinguishing their intracellular levels. Our results show a dimension to MITF regulation in melanoma cells and point to CDK7 inhibition as a potential strategy to deprive oncogenic transcription and suppress tumor growth in melanoma. Keywords: ChIP-seq; chromatin immunoprecipitation sequencing; CTD; carboxy-terminal domain; GI50; GSEA; gene set enrichment analysis; MITF-hi; high MITF expression level; MITF-lo; low MITF expression level; Pol IIpolymerase II; super-enhancer; Serserine; siRNA; small interfering RNA, National Institutes of Health (U.S.) (Grant K24-CA149202), National Institutes of Health (U.S.) (Grant P01-CA163222), National Institutes of Health (U.S.) (Grant T32-CA071345)

High-risk neuroblastomas show a paucity of recurrent somatic mutations at diagnosis. As a result, the molecular basis for this aggressive phenotype remains elusive. Recent progress in regulatory network analysis helped us elucidate disease-driving mechanisms downstream of genomic alterations, including recurrent chromosomal alterations. Our analysis identified three molecular subtypes of high-risk neuroblastomas, consistent with chromosomal alterations, and identified subtype-specific master regulator proteins that were conserved across independent cohorts. A 10-protein transcriptional module—centered around a TEAD4–MYCN positive feedback loop—emerged as the regulatory driver of the high-risk subtype associated with MYCN amplification. Silencing of either gene collapsed MYCN-amplified (MYCNAmp) neuroblastoma transcriptional hallmarks and abrogated viability in vitro and in vivo. Consistently, TEAD4 emerged as a robust prognostic marker of poor survival, with activity independent of the canonical Hippo pathway transcriptional coactivators YAP and TAZ. These results suggest novel therapeutic strategies for the large subset of MYCN-deregulated neuroblastomas. Significance: Despite progress in understanding of neuroblastoma genetics, little progress has been made toward personalized treatment. Here, we present a framework to determine the downstream effectors of the genetic alterations sustaining neuroblastoma subtypes, which can be easily extended to other tumor types. We show the critical effect of disrupting a 10-protein module centered around a YAP/TAZ-independent TEAD4–MYCN positive feedback loop in MYCNAmp neuroblastomas, nominating TEAD4 as a novel candidate for therapeutic intervention. Cancer Discov; 8(5); 582–99. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 517

Gene expression is regulated by promoters and enhancers marked by histone H3-lysine-27 acetylation (H3K27ac) and established by the paralogous histone acetyltransferases (HATs), EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma depend on EP300, while CBP has a limited role. EP300 controls enhancer acetylation by interacting with TFAP2β, a transcription factor member of the lineage-defining core-regulatory transcriptional circuitry (CRC). To disrupt EP300, we developed a proteolysis-targeted-chimaera (PROTAC) compound “JQAD1” that selectively targets EP300 for degradation. This results in loss of H3K27ac at CRC enhancers and neuroblastoma apoptosis, with limited toxicity to untransformed cells where CBP compensates. EP300 physically interacts with the MYCN oncoprotein, promoting its chromatin localization in neuroblastoma. EP300 degradation therefore capitalizes on EP300 dependency in neuroblastoma by controlling MYCN localization and preventing MYCN from sustaining the malignant cell state. Citation Format: Adam D. Durbin, Tingjian Wang, Virangika K. Wimalasena, Mark W. Zimmerman, Deyao Li, Neekesh V. Dharia, Luca Mariani, Paul M. Park, Logan H. Sigua, Ken Morita, Amy S. Conway, Amanda L. Robichaud, Antonio Perez-Atayde, Melissa Bikowitz, Taylor Quinn, Olaf Wiest, Ernst Schonbrunn, Martha Bulyk, Brian J. Abraham, Kimberly Stegmaier, A. Thomas Look, Jun Qi. EP300 selectively controls the enhancer landscape of high risk neuroblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2318.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignancy of thymocytes. Using a transgenic screen in zebrafish, thymocyte selection–associated high mobility group box protein (TOX) was uncovered as a collaborating oncogenic driver that accelerated T-ALL onset by expanding the initiating pool of transformed clones and elevating genomic instability. TOX is highly expressed in a majority of human T-ALL and is required for proliferation and continued xenograft growth in mice. Using a wide array of functional analyses, we uncovered that TOX binds directly to KU70/80 and suppresses recruitment of this complex to DNA breaks to inhibit nonhomologous end joining (NHEJ) repair. Impaired NHEJ is well known to cause genomic instability, including development of T-cell malignancies in KU70- and KU80-deficient mice. Collectively, our work has uncovered important roles for TOX in regulating NHEJ by elevating genomic instability during leukemia initiation and sustaining leukemic cell proliferation following transformation. Significance: TOX is an HMG box–containing protein that has important roles in T-ALL initiation and maintenance. TOX inhibits the recruitment of KU70/KU80 to DNA breaks, thereby inhibiting NHEJ repair. Thus, TOX is likely a dominant oncogenic driver in a large fraction of human T-ALL and enhances genomic instability. Cancer Discov; 7(11); 1336–53. ©2017 AACR. This article is highlighted in the In This Issue feature, p. 1201

CD47 is a cell surface molecule that inhibits phagocytosis of cells that express it by binding to its receptor, SIRPα, on macrophages and other immune cells. CD47 is expressed at different levels by neoplastic and normal cells. Here, to reveal mechanisms by which different neoplastic cells generate this dominant ‘don’t eat me’ signal, we analyse the CD47 regulatory genomic landscape. We identify two distinct super-enhancers (SEs) associated with CD47 in certain cancer cell types. We show that a set of active constituent enhancers, located within the two CD47 SEs, regulate CD47 expression in different cancer cell types and that disruption of CD47 SEs reduces CD47 gene expression. Finally we report that the TNF-NFKB1 signalling pathway directly regulates CD47 by interacting with a constituent enhancer located within a CD47-associated SE specific to breast cancer. These results suggest that cancers can evolve SE to drive CD47 overexpression to escape immune surveillance.

Oncogenic driver mutations are those that provide a proliferative or survival advantage to neoplastic cells, resulting in clonal selection. Although most cancer-causing mutations have been detected in the protein-coding regions of the cancer genome; driver mutations have recently also been discovered within noncoding genomic sequences. Thus, a current challenge is to gain precise understanding of how these unique genomic elements function in cancer pathogenesis, while clarifying mechanisms of gene regulation and identifying new targets for therapeutic intervention. Here we report a C-to-T single nucleotide transition that occurs as a somatic mutation in noncoding sequences 4 kb upstream of the transcriptional start site of the LMO1 oncogene in primary samples from patients with T-cell acute lymphoblastic leukaemia. This single nucleotide alteration conforms to an APOBEC-like cytidine deaminase mutational signature, and generates a new binding site for the MYB transcription factor, leading to the formation of an aberrant transcriptional enhancer complex that drives high levels of expression of the LMO1 oncogene. Since APOBEC-signature mutations are common in a broad spectrum of human cancers, we suggest that noncoding nucleotide transitions such as the one described here may activate potent oncogenic enhancers not only in T-lymphoid cells but in other cell lineages as well., National Institutes of Health (U.S.) (Grant 5P01CA109901)

The function of critical developmental regulators can be subverted by cancer cells to control expression of oncogenic transcriptional programs. These "master transcription factors" (MTFs) are often essential for cancer cell survival and represent vulnerabilities that can be exploited therapeutically. The current approaches to identify candidate MTFs examine super-enhancer associated transcription factor-encoding genes with high connectivity in network models. This relies on chromatin immunoprecipitation-sequencing (ChIP-seq) data, which is technically challenging to obtain from primary tumors, and is currently unavailable for many cancer types and clinically relevant subtypes. In contrast, gene expression data are more widely available, especially for rare tumors and subtypes where MTFs have yet to be discovered. We have developed a predictive algorithm called CaCTS (Cancer Core Transcription factor Specificity) to identify candidate MTFs using pan-cancer RNA-sequencing data from The Cancer Genome Atlas. The algorithm identified 273 candidate MTFs across 34 tumor types and recovered known tumor MTFs. We also made novel predictions, including for cancer types and subtypes for which MTFs have not yet been characterized. Clustering based on MTF predictions reproduced anatomic groupings of tumors that share 1-2 lineage-specific candidates, but also dictated functional groupings, such as a squamous group that comprised five tumor subtypes sharing 3 common MTFs. PAX8, SOX17, and MECOM were candidate factors in high-grade serous ovarian cancer (HGSOC), an aggressive tumor type where the core regulatory circuit is currently uncharacterized. PAX8, SOX17, and MECOM are required for cell viability and lie proximal to super-enhancers in HGSOC cells. ChIP-seq revealed that these factors co-occupy HGSOC regulatory elements globally and co-bind at critical gene loci including MUC16 (CA-125). Addiction to these factors was confirmed in studies using THZ1 to inhibit transcription in HGSOC cells, suggesting early down-regulation of these genes may be responsible for cytotoxic effects of THZ1 on HGSOC models. Identification of MTFs across 34 tumor types and 140 subtypes, especially for those with limited understanding of transcriptional drivers paves the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Transcriptional Cyclin Dependent Kinases modulate RNA Polymerase II function to impact gene expression. Here, we show that CDK13 is mutated in 4% of patient melanomas and mutation or downregulation is associated with poor overall survival. Mutant CDK13 lacks kinase activity and overexpression in zebrafish leads to accelerated melanoma. CDK13 mutant fish and human melanomas accumulate prematurely terminated RNAs that are translated into truncated proteins. CDK13 binds to and regulates the phosphorylation of ZC3H14, a member of the PolyA eXosome Targeting (PAXT) RNA degradation complex. ZC3H14 phosphorylation recruits the PAXT complex to degrade prematurely terminated polyadenylated transcripts in the nucleus. In the presence of mutant CDK13, ZC3H14 phosphorylation is compromised and consequently fails to recruit the PAXT complex, leading to truncated transcript stabilization. This work establishes a role for CDK13 and the PAXT nuclear RNA degradation complex in cancer and has prognostic significance for melanoma patients with mutated or downregulated CDK13.

The gene expression programs that define the identity of each cell are controlled by master transcription factors (TFs) that bind cell-type-specific enhancers, as well as signaling factors, which bring extracellular stimuli to these enhancers. Recent studies have revealed that master TFs form phase-separated condensates with the Mediator coactivator at super-enhancers. Here, we present evidence that signaling factors for the WNT, TGF-β, and JAK/STAT pathways use their intrinsically disordered regions (IDRs) to enter and concentrate in Mediator condensates at super-enhancers. We show that the WNT coactivator β-catenin interacts both with components of condensates and DNA-binding factors to selectively occupy super-enhancer-associated genes. We propose that the cell-type specificity of the response to signaling is mediated in part by the IDRs of the signaling factors, which cause these factors to partition into condensates established by the master TFs and Mediator at genes with prominent roles in cell identity., National Science Foundation (U.S.). (Grant PHY1743900), National Institutes of Health (U.S.) (Grants GM123511, GM117370 and T32CA009172), German Research Foundation Research Fellowship (DE 3069/1-1)