Specific oncogene activation of the cell of origin in mucosal melanoma.

Mucosal melanoma (MM) is a deadly cancer derived from mucosal melanocytes. To test the consequences of MM genetics, we developed a zebrafish model in which all melanocytes experienced CCND1 expression and loss of PTEN and TP53. Surprisingly, melanoma only developed from melanocytes lining internal organs, analogous to the location of patient MM. We found that zebrafish MMs had a unique chromatin landscape from cutaneous melanoma. Internal melanocytes could be labeled using a MM-specific transcriptional enhancer. Normal zebrafish internal melanocytes shared a gene expression signature with MMs. Patient and zebrafish MMs have increased migratory neural crest gene and decreased antigen presentation gene expression, consistent with the increased metastatic behavior and decreased immunotherapy sensitivity of MM. Our work suggests the cell state of the originating melanocyte influences the behavior of derived melanomas. Our animal model phenotypically and transcriptionally mimics patient tumors, allowing this model to be used for MM therapeutic discovery.

Bone morphogenetic protein (BMP) signaling determines neuroblastoma cell fate and sensitivity to retinoic acid.

Retinoic acid (RA) is a standard-of-care neuroblastoma drug thought to be effective by inducing differentiation. Curiously, RA has little effect on primary human tumors during upfront treatment but can eliminate neuroblastoma cells from the bone marrow during post-chemo consolidation therapy-a discrepancy that has never been explained. To investigate this, we treated a large cohort of neuroblastoma cell lines with RA and observed that the most RA-sensitive cells predominantly undergo apoptosis or senescence, rather than differentiation. We conducted genome-wide CRISPR knockout screens under RA treatment, which identified BMP signaling as controlling the apoptosis/senescence vs differentiation cell fate decision and determining RA's overall potency. We then discovered that BMP signaling activity is markedly higher in neuroblastoma patient samples at bone marrow metastatic sites, providing a plausible explanation for RA's ability to clear neuroblastoma cells specifically from the bone marrow, seemingly mimicking interactions between BMP and RA during normal development.

Group 3 medulloblastoma transcriptional networks collapse under domain specific EP300/CBP inhibition.

Chemical discovery efforts commonly target individual protein domains. Many proteins, including the EP300/CBP histone acetyltransferases (HATs), contain several targetable domains. EP300/CBP are critical gene-regulatory targets in cancer, with existing high potency inhibitors of either the catalytic HAT domain or protein-binding bromodomain (BRD). A domain-specific inhibitory approach to multidomain-containing proteins may identify exceptional-responding tumor types, thereby expanding a therapeutic index. Here, we discover that targeting EP300/CBP using the domain-specific inhibitors, A485 (HAT) or CCS1477 (BRD) have different effects in select tumor types. Group 3 medulloblastoma (G3MB) cells are especially sensitive to BRD, compared with HAT inhibition. Structurally, these effects are mediated by the difluorophenyl group in the catalytic core of CCS1477. Mechanistically, bromodomain inhibition causes rapid disruption of genetic dependency networks that are required for G3MB growth. These studies provide a domain-specific structural foundation for drug discovery efforts targeting EP300/CBP and identify a selective role for the EP300/CBP bromodomain in maintaining genetic dependency networks in G3MB.

STAT3 couples activated tyrosine kinase signaling to the oncogenic core transcriptional regulatory circuitry of anaplastic large cell lymphoma.

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 in ALCL cooccupies enhancers with master transcription factors BATF3, IRF4, and IKZF1 to form a core regulatory circuit that establishes and maintains the malignant cell state in ALCL. Critical downstream targets of this network in ALCL cells include the protooncogene MYC, which requires active STAT3 to facilitate high levels of MYC transcription. The core autoregulatory transcriptional circuitry activity is reinforced by MYC binding to the enhancer regions associated with STAT3 and each of the core regulatory transcription factors. Thus, activation of STAT3 provides the crucial link between aberrant tyrosine kinase signaling and the core transcriptional machinery that drives tumorigenesis and creates therapeutic vulnerabilities in ALCL.

PAX3-FOXO1 dictates myogenic reprogramming and rhabdomyosarcoma identity in endothelial progenitors.

Fusion-positive rhabdomyosarcoma (FP-RMS) driven by the expression of the PAX3-FOXO1 (P3F) fusion oncoprotein is an aggressive subtype of pediatric rhabdomyosarcoma. FP-RMS histologically resembles developing muscle yet occurs throughout the body in areas devoid of skeletal muscle highlighting that FP-RMS is not derived from an exclusively myogenic cell of origin. Here we demonstrate that P3F reprograms mouse and human endothelial progenitors to FP-RMS. We show that P3F expression in aP2-Cre expressing cells reprograms endothelial progenitors to functional myogenic stem cells capable of regenerating injured muscle fibers. Further, we describe a FP-RMS mouse model driven by P3F expression and Cdkn2a loss in endothelial cells. Additionally, we show that P3F expression in TP53-null human iPSCs blocks endothelial-directed differentiation and guides cells to become myogenic cells that form FP-RMS tumors in immunocompromised mice. Together these findings demonstrate that FP-RMS can originate from aberrant development of non-myogenic cells driven by P3F.

Regulatory architecture of housekeeping genes is driven by promoter assemblies.

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

Genetic predisposition to neuroblastoma results from a regulatory polymorphism that promotes the adrenergic cell state.

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like mesenchymal cell state and a more differentiated sympathetic adrenergic cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional cofactor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism, G/T, that affects a GATA motif in the first intron of LMO1. Here, we showed that WT zebrafish with the GATA genotype developed adrenergic neuroblastoma, while knock-in of the protective TATA allele at this locus reduced the penetrance of MYCN-driven tumors, which were restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrated that TATA/TATA tumors also exhibited a mesenchymal cell state and were low risk at diagnosis. Thus, conversion of the regulatory GATA to a TATA allele in the first intron of LMO1 reduced the neuroblastoma-initiation rate by preventing formation of the adrenergic cell state. This mechanism was conserved over 400 million years of evolution, separating zebrafish and humans.

Oncogenic CDK13 mutations impede nuclear RNA surveillance.

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.

Genetic Predisposition to Neuroblastoma Results from a Regulatory Polymorphism that Promotes the Adrenergic Cell State.

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like "mesenchymal" cell state and a more differentiated sympathetic "adrenergic" cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional co-factor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism G/T that affects a G ATA motif in the first intron of LMO1. Here we show that wild-type zebrafish with the G ATA genotype develop adrenergic neuroblastoma, while knock-in of the protective T ATA allele at this locus reduces the penetrance of MYCN-driven tumors, which are restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrates that T ATA/ T ATA tumors also exhibit a mesenchymal cell state and are low risk at diagnosis. Thus, conversion of the regulatory G ATA to a T ATA allele in the first intron of LMO1 reduces the neuroblastoma initiation rate by preventing formation of the adrenergic cell state, a mechanism that is conserved over 400 million years of evolution separating zebrafish and humans.

BRD2 interconnects with BRD3 to facilitate Pol II transcription initiation and elongation to prime promoters for cell differentiation.

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.

Intrinsic transcriptional heterogeneity in neuroblastoma guides mechanistic and therapeutic insights.

Cell state is controlled by master transcription factors (mTFs) that determine the cellular gene expression program. Cancer cells acquire dysregulated gene expression programs by mutational and non-mutational processes. Intratumoral heterogeneity can result from cells displaying distinct mTF-regulated cell states, which co-exist within the tumor. One archetypal tumor associated with transcriptionally regulated heterogeneity is high-risk neuroblastoma (NB). Patients with NB have poor overall survival despite intensive therapies, and relapsed patients are commonly refractory to treatment. The cellular populations that comprise NB are marked by different cohorts of mTFs and differential sensitivity to conventional therapies. Recent studies have highlighted mechanisms by which NB cells dynamically shift the cell state with treatment, revealing new opportunities to control the cellular response to treatment by manipulating cell-state-defining transcriptional programs. Here, we review recent advances in understanding transcriptionally defined cancer heterogeneity. We offer challenges to the field to encourage translation of basic science into clinical benefit.

The nuclear receptor THRB facilitates differentiation of human PSCs into more mature hepatocytes.

To understand the mechanisms regulating the in vitro maturation of hPSC-derived hepatocytes, we developed a 3D differentiation system and compared gene regulatory elements in human primary hepatocytes with those in hPSC-hepatocytes that were differentiated in 2D or 3D conditions by RNA-seq, ATAC-seq, and H3K27Ac ChIP-seq. Regulome comparisons showed a reduced enrichment of thyroid receptor THRB motifs in accessible chromatin and active enhancers without a reduced transcription of THRB. The addition of thyroid hormone T3 increased the binding of THRB to the CYP3A4 proximal enhancer, restored the super-enhancer status and gene expression of NFIC, and reduced the 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 hPSC-hepatocytes resulted in the engraftment of human hepatocytes into the mouse liver without disrupting normal liver histology. This work implicates the environmental factor-nuclear receptor axis in regulating the maturation of hPSC-hepatocytes.

SEAseq: a portable and cloud-based chromatin occupancy analysis suite.

BACKGROUND: 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 (CUT&RUN) 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. RESULTS: We present a comprehensive, infrastructure-independent, computational pipeline called SEAseq, which leverages field-standard, open-source tools for processing and analyzing ChIP-Seq/CUT&RUN 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. CONCLUSIONS: 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 .

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.

Predicting master transcription factors from pan-cancer expression data.

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.

Retinoic acid rewires the adrenergic core regulatory circuitry of childhood neuroblastoma.

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.

MEIS2 Is an Adrenergic Core Regulatory Transcription Factor Involved in Early Initiation of TH-MYCN-Driven Neuroblastoma Formation.

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.

Super-enhancer-based identification of a BATF3/IL-2R-module reveals vulnerabilities in anaplastic large cell lymphoma.

Anaplastic large cell lymphoma (ALCL), an aggressive CD30-positive T-cell lymphoma, comprises systemic anaplastic lymphoma kinase (ALK)-positive, and ALK-negative, primary cutaneous and breast implant-associated ALCL. Prognosis of some ALCL subgroups is still unsatisfactory, and already in second line effective treatment options are lacking. To identify genes defining ALCL cell state and dependencies, we here characterize super-enhancer regions by genome-wide H3K27ac ChIP-seq. In addition to known ALCL key regulators, the AP-1-member BATF3 and IL-2 receptor (IL2R)-components are among the top hits. Specific and high-level IL2R expression in ALCL correlates with BATF3 expression. Confirming a regulatory link, IL-2R-expression decreases following BATF3 knockout, and BATF3 is recruited to IL2R regulatory regions. Functionally, IL-2, IL-15 and Neo-2/15, a hyper-stable IL-2/IL-15 mimic, accelerate ALCL growth and activate STAT1, STAT5 and ERK1/2. In line, strong IL-2Rα-expression in ALCL patients is linked to more aggressive clinical presentation. Finally, an IL-2Rα-targeting antibody-drug conjugate efficiently kills ALCL cells in vitro and in vivo. Our results highlight the importance of the BATF3/IL-2R-module for ALCL biology and identify IL-2Rα-targeting as a promising treatment strategy for ALCL.

Cell-specific transcriptional control of mitochondrial metabolism by TIF1γ drives erythropoiesis.

Transcription and metabolism both influence cell function, but dedicated transcriptional control of metabolic pathways that regulate cell fate has rarely been defined. We discovered, using a chemical suppressor screen, that inhibition of the pyrimidine biosynthesis enzyme dihydroorotate dehydrogenase (DHODH) rescues erythroid differentiation in bloodless zebrafish moonshine (mon) mutant embryos defective for transcriptional intermediary factor 1 gamma (tif1γ). This rescue depends on the functional link of DHODH to mitochondrial respiration. The transcription elongation factor TIF1γ directly controls coenzyme Q (CoQ) synthesis gene expression. Upon tif1γ loss, CoQ levels are reduced, and a high succinate/α-ketoglutarate ratio leads to increased histone methylation. A CoQ analog rescues mon's bloodless phenotype. These results demonstrate that mitochondrial metabolism is a key output of a lineage transcription factor that drives cell fate decisions in the early blood lineage.