Lineage specific transcription factor waves reprogram neuroblastoma from self-renewal to differentiation.

Temporal regulation of super-enhancer (SE) driven transcription factors (TFs) underlies normal developmental programs. Neuroblastoma (NB) arises from an inability of sympathoadrenal progenitors to exit a self-renewal program and terminally differentiate. To identify SEs driving TF regulators, we use all-trans retinoic acid (ATRA) to induce NB growth arrest and differentiation. Time-course H3K27ac ChIP-seq and RNA-seq reveal ATRA coordinated SE waves. SEs that decrease with ATRA link to stem cell development (MYCN, GATA3, SOX11). CRISPR-Cas9 and siRNA verify SOX11 dependency, in vitro and in vivo. Silencing the SOX11 SE using dCAS9-KRAB decreases SOX11 mRNA and inhibits cell growth. Other TFs activate in sequential waves at 2, 4 and 8 days of ATRA treatment that regulate neural development (GATA2 and SOX4). Silencing the gained SOX4 SE using dCAS9-KRAB decreases SOX4 expression and attenuates ATRA-induced differentiation genes. Our study identifies oncogenic lineage drivers of NB self-renewal and TFs critical for implementing a differentiation program.

KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma.

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

Aberrant activation of wound healing programs within the metastatic niche facilitates lung colonization by osteosarcoma cells.

Purpose: Lung metastasis is responsible for nearly all deaths caused by osteosarcoma, the most common pediatric bone tumor. How malignant bone cells coerce the lung microenvironment to support metastatic growth is unclear. This study delineates how osteosarcoma cells educate the lung microenvironment during metastatic progression. Experimental design: Using single-cell transcriptomics (scRNA-seq), we characterized genome- and tissue-wide molecular changes induced within lung tissues by disseminated osteosarcoma cells in both immunocompetent murine models of metastasis and patient samples. We confirmed transcriptomic findings at the protein level and determined spatial relationships with multi-parameter immunofluorescence. We evaluated the ability of nintedanib to impair metastatic colonization and prevent osteosarcoma-induced education of the lung microenvironment in both immunocompetent murine osteosarcoma and immunodeficient human xenograft models. Results: Osteosarcoma cells induced acute alveolar epithelial injury upon lung dissemination. scRNA-seq demonstrated that the surrounding lung stroma adopts a chronic, non-resolving wound-healing phenotype similar to that seen in other models of lung injury. Accordingly, metastasis-associated lung demonstrated marked fibrosis, likely due to the accumulation of pathogenic, pro-fibrotic, partially-differentiated epithelial intermediates. Inhibition of fibrotic pathways with nintedanib prevented metastatic progression in multiple murine and human xenograft models. Conclusions: Our work demonstrates that osteosarcoma cells co-opt fibrosis to promote metastatic outgrowth. When harmonized with data from adult epithelial cancers, our results support a generalized model wherein aberrant mesenchymal-epithelial interactions are critical for promoting lung metastasis. Adult epithelial carcinomas induce fibrotic pathways in normal lung fibroblasts, whereas osteosarcoma, a pediatric mesenchymal tumor, exhibits fibrotic reprogramming in response to the aberrant wound-healing behaviors of an otherwise normal lung epithelium, which are induced by tumor cell interactions. Statement of translational relevance: Therapies that block metastasis have the potential to save the majority of lives lost due to solid tumors. Disseminated tumor cells must educate the foreign, inhospitable microenvironments they encounter within secondary organs to facilitate metastatic colonization. Our study elucidated that disseminated osteosarcoma cells survive within the lung by co-opting and amplifying the lung's endogenous wound healing response program. More broadly, our results support a model wherein mesenchymal-epithelial cooperation is a key driver of lung metastasis. Osteosarcoma, a pediatric mesenchymal tumor, undergoes lung epithelial induced fibrotic activation while also transforming normal lung epithelial cells towards a fibrosis promoting phenotype. Conversely, adult epithelial carcinomas activate fibrotic signaling in normal lung mesenchymal fibroblasts. Our data implicates fibrosis and abnormal wound healing as key drivers of lung metastasis across multiple tumor types that can be targeted therapeutically to disrupt metastasis progression.

MYOD-SKP2 axis boosts tumorigenesis in fusion negative rhabdomyosarcoma by preventing differentiation through p57Kip2 targeting.

Rhabdomyosarcomas (RMS) are pediatric mesenchymal-derived malignancies encompassing PAX3/7-FOXO1 Fusion Positive (FP)-RMS, and Fusion Negative (FN)-RMS with frequent RAS pathway mutations. RMS express the master myogenic transcription factor MYOD that, whilst essential for survival, cannot support differentiation. Here we discover SKP2, an oncogenic E3-ubiquitin ligase, as a critical pro-tumorigenic driver in FN-RMS. We show that SKP2 is overexpressed in RMS through the binding of MYOD to an intronic enhancer. SKP2 in FN-RMS promotes cell cycle progression and prevents differentiation by directly targeting p27Kip1 and p57Kip2, respectively. SKP2 depletion unlocks a partly MYOD-dependent myogenic transcriptional program and strongly affects stemness and tumorigenic features and prevents in vivo tumor growth. These effects are mirrored by the investigational NEDDylation inhibitor MLN4924. Results demonstrate a crucial crosstalk between transcriptional and post-translational mechanisms through the MYOD-SKP2 axis that contributes to tumorigenesis in FN-RMS. Finally, NEDDylation inhibition is identified as a potential therapeutic vulnerability in FN-RMS.

PAX3-FOXO1 uses its activation domain to recruit CBP/P300 and shape RNA Pol2 cluster distribution.

Activation of oncogenic gene expression from long-range enhancers is initiated by the assembly of DNA-binding transcription factors (TF), leading to recruitment of co-activators such as CBP/p300 to modify the local genomic context and facilitate RNA-Polymerase 2 (Pol2) binding. Yet, most TF-to-coactivator recruitment relationships remain unmapped. Here, studying the oncogenic fusion TF PAX3-FOXO1 (P3F) from alveolar rhabdomyosarcoma (aRMS), we show that a single cysteine in the activation domain (AD) of P3F is important for a small alpha helical coil that recruits CBP/p300 to chromatin. P3F driven transcription requires both this single cysteine and CBP/p300. Mutants of the cysteine reduce aRMS cell proliferation and induce cellular differentiation. Furthermore, we discover a profound dependence on CBP/p300 for clustering of Pol2 loops that connect P3F to its target genes. In the absence of CBP/p300, Pol2 long range enhancer loops collapse, Pol2 accumulates in CpG islands and fails to exit the gene body. These results reveal a potential novel axis for therapeutic interference with P3F in aRMS and clarify the molecular relationship of P3F and CBP/p300 in sustaining active Pol2 clusters essential for oncogenic transcription.

Temporal chromatin accessibility changes define transcriptional states essential for osteosarcoma metastasis.

The metastasis-invasion cascade describes the series of steps required for a cancer cell to successfully spread from its primary tumor and ultimately grow within a secondary organ. Despite metastasis being a dynamic, multistep process, most omics studies to date have focused on comparing primary tumors to the metastatic deposits that define end-stage disease. This static approach means we lack information about the genomic and epigenomic changes that occur during the majority of tumor progression. One particularly understudied phase of tumor progression is metastatic colonization, during which cells must adapt to the new microenvironment of the secondary organ. Through temporal profiling of chromatin accessibility and gene expression in vivo, we identify dynamic changes in the epigenome that occur as osteosarcoma tumors form and grow within the lung microenvironment. Furthermore, we show through paired in vivo and in vitro CRISPR drop-out screens and pharmacological validation that the upstream transcription factors represent a class of metastasis-specific dependency genes. While current models depict lung colonization as a discrete step within the metastatic cascade, our study shows it is a defined trajectory through multiple epigenetic states, revealing new therapeutic opportunities undetectable with standard approaches.

A SNAI2/CTCF Interaction is Required for NOTCH1 Expression in Rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of the muscle with characteristics of cells blocked in differentiation. NOTCH1 is an oncogene that promotes self-renewal and blocks differentiation in the fusion negative-RMS sub-type. However, how NOTCH1 expression is transcriptionally maintained in tumors is unknown. Analyses of SNAI2 and CTCF chromatin binding and HiC analyses revealed a conserved SNAI2/CTCF overlapping peak downstream of the NOTCH1 locus marking a sub-topologically associating domain (TAD) boundary. Deletion of the SNAI2-CTCF peak showed that it is essential for NOTCH1 expression and viability of FN-RMS cells. Reintroducing constitutively activated NOTCH1-ΔE in cells with the SNAI2-CTCF peak deleted restored cell-viability. Ablation of SNAI2 using CRISPR/Cas9 reagents resulted in the loss of majority of RD and SMS-CTR FN-RMS cells. However, the few surviving clones that repopulate cultures have recovered NOTCH1. Cells that re-establish NOTCH1 expression after SNAI2 ablation are unable to differentiate robustly as SNAI2 shRNA knockdown cells; yet, SNAI2-ablated cells continued to be exquisitely sensitive to ionizing radiation. Thus, we have uncovered a novel mechanism by which SNAI2 and CTCF maintenance of a sub-TAD boundary promotes rather than represses NOTCH1 expression. Further, we demonstrate that SNAI2 suppression of apoptosis post-radiation is independent of SNAI2/NOTCH1 effects on self-renewal and differentiation.

Preclinical development of a chimeric antigen receptor T cell therapy targeting FGFR4 in rhabdomyosarcoma.

Pediatric patients with relapsed or refractory rhabdomyosarcoma (RMS) have dismal cure rates, and effective therapy is urgently needed. The oncogenic receptor tyrosine kinase fibroblast growth factor receptor 4 (FGFR4) is highly expressed in RMS and lowly expressed in healthy tissues. Here, we describe a second-generation FGFR4-targeting chimeric antigen receptor (CAR), based on an anti-human FGFR4-specific murine monoclonal antibody 3A11, as an adoptive T cell treatment for RMS. The 3A11 CAR T cells induced robust cytokine production and cytotoxicity against RMS cell lines in vitro. In contrast, a panel of healthy human primary cells failed to activate 3A11 CAR T cells, confirming the selectivity of 3A11 CAR T cells against tumors with high FGFR4 expression. Finally, we demonstrate that 3A11 CAR T cells are persistent in vivo and can effectively eliminate RMS tumors in two metastatic and two orthotopic models. Therefore, our study credentials CAR T cell therapy targeting FGFR4 to treat patients with RMS.

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.

BET Bromodomain Degradation Disrupts Function but Not 3D Formation of RNA Pol2 Clusters.

Fusion-positive rhabdomyosarcoma (FP-RMS) is driven by a translocation that creates the chimeric transcription factor PAX3-FOXO1 (P3F), which assembles de novo super enhancers to drive high levels of transcription of other core regulatory transcription factors (CRTFs). P3F recruits co-regulatory factors to super enhancers such as BRD4, which recognizes acetylated lysines via BET bromodomains. In this study, we demonstrate that inhibition or degradation of BRD4 leads to global decreases in transcription, and selective downregulation of CRTFs. We also show that the BRD4 degrader ARV-771 halts transcription while preserving RNA Polymerase II (Pol2) loops between super enhancers and their target genes, and causes the removal of Pol2 only past the transcriptional end site of CRTF genes, suggesting a novel effect of BRD4 on Pol2 looping. We finally test the most potent molecule, inhibitor BMS-986158, in an orthotopic PDX mouse model of FP-RMS with additional high-risk mutations, and find that it is well tolerated in vivo and leads to an average decrease in tumor size. This effort represents a partnership with an FP-RMS patient and family advocates to make preclinical data rapidly accessible to the family, and to generate data to inform future patients who develop this disease.

High-throughput approaches to uncover synergistic drug combinations in leukemia.

We report a comprehensive drug synergy study in acute myeloid leukemia (AML). In this work, we investigate a panel of cell lines spanning both MLL-rearranged and non-rearranged subtypes. The work comprises a resource for the community, with many synergistic drug combinations that could not have been predicted a priori, and open source code for automation and analyses. We base our definitions of drug synergy on the Chou-Talalay method, which is useful for visualizations of synergy experiments in isobolograms, and median-effects plots, among other representations. Our key findings include drug synergies affecting the chromatin state, specifically in the context of regulation of the modification state of histone H3 lysine-27. We report open source high throughput methodology such that multidimensional drug screening can be accomplished with equipment that is accessible to most laboratories. This study will enable preclinical investigation of new drug combinations in a lethal blood cancer, with data analysis and automation workflows freely available to the community.

A stem cell epigenome is associated with primary nonresponse to CD19 CAR T cells in pediatric acute lymphoblastic leukemia.

CD19 chimeric antigen receptor T-cell therapy (CD19-CAR) has changed the treatment landscape and outcomes for patients with pre-B-cell acute lymphoblastic leukemia (B-ALL). Unfortunately, primary nonresponse (PNR), sustained CD19+ disease, and concurrent expansion of CD19-CAR occur in 20% of the patients and is associated with adverse outcomes. Although some failures may be attributable to CD19 loss, mechanisms of CD19-independent, leukemia-intrinsic resistance to CD19-CAR remain poorly understood. We hypothesize that PNR leukemias are distinct compared with primary sensitive (PS) leukemias and that these differences are present before treatment. We used a multiomic approach to investigate this in 14 patients (7 with PNR and 7 with PS) enrolled in the PLAT-02 trial at Seattle Children's Hospital. Long-read PacBio sequencing helped identify 1 PNR in which 47% of CD19 transcripts had exon 2 skipping, but other samples lacked CD19 transcript abnormalities. Epigenetic profiling discovered DNA hypermethylation at genes targeted by polycomb repressive complex 2 (PRC2) in embryonic stem cells. Similarly, assays of transposase-accessible chromatin-sequencing revealed reduced accessibility at these PRC2 target genes, with a gain in accessibility of regions characteristic of hematopoietic stem cells and multilineage progenitors in PNR. Single-cell RNA sequencing and cytometry by time of flight analyses identified leukemic subpopulations expressing multilineage markers and decreased antigen presentation in PNR. We thus describe the association of a stem cell epigenome with primary resistance to CD19-CAR therapy. Future trials incorporating these biomarkers, with the addition of multispecific CAR T cells targeting against leukemic stem cell or myeloid antigens, and/or combined epigenetic therapy to disrupt this distinct stem cell epigenome may improve outcomes of patients with B-ALL.

Dynamic Opposition of Histone Modifications.

Gene transcription does not only require writers of active histone modifications; on-site opposition by erasers is essential for many genes. Here, we propose the concept of dynamic opposition of histone modifications to explain this conundrum. We highlight the requirement of HDACs for acetylation balance at superenhancers, and the requirement of KDM5A for H4K3me3 recycling at highly active gene promoters. We propose that histone post-translational modifications regulate charge balance for biomolecular condensate formation and nucleosome turnover and form a short-term memory that informs lock-and-step checkpoints for chromatin engagement by RNA polymerase II.

A local tumor microenvironment acquired super-enhancer induces an oncogenic driver in colorectal carcinoma.

Tumors exhibit enhancer reprogramming compared to normal tissue. The etiology is largely attributed to cell-intrinsic genomic alterations. Here, using freshly resected primary CRC tumors and patient-matched adjacent normal colon, we find divergent epigenetic landscapes between CRC tumors and cell lines. Intriguingly, this phenomenon extends to highly recurrent aberrant super-enhancers gained in CRC over normal. We find one such super-enhancer activated in epithelial cancer cells due to surrounding inflammation in the tumor microenvironment. We restore this super-enhancer and its expressed gene, PDZK1IP1, following treatment with cytokines or xenotransplantation into nude mice, thus demonstrating cell-extrinsic etiology. We demonstrate mechanistically that PDZK1IP1 enhances the reductive capacity CRC cancer cells via the pentose phosphate pathway. We show this activation enables efficient growth under oxidative conditions, challenging the previous notion that PDZK1IP1 acts as a tumor suppressor in CRC. Collectively, these observations highlight the significance of epigenomic profiling on primary specimens.

Dentithecamides A-H, Diacylated Zoanthoxanthin Derivatives with PAX3-FOXO1 Inhibitory Activity from the Hydroid Dentitheca habereri.

Chemical investigation of the marine hydroid Dentitheca habereri led to the identification of eight new diacylated zoanthoxanthin alkaloids, named dentithecamides A-H (1-8), along with three previously reported analogues, zoamides B-D (9-11). The structures of compounds 1-11 were elucidated by spectroscopic and spectrometric analyses, including IR, HRESIMS, and NMR experiments, and by comparison with literature data. Compounds 1-11 are the first zoanthoxanthin alkaloids to be reported from a hydroid. Dentithecamides A (1) and B (2) along with zoamides B-D (9-11), which all share a conformationally mobile cycloheptadiene core, inhibited PAX3-FOXO1 regulated transcriptional activity and thus provided a structural framework for the potential development of more potent PAX3-FOXO1 inhibitors.

Convergence of case-specific epigenetic alterations identify a confluence of genetic vulnerabilities tied to opioid overdose.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Chromatin Mechanisms Driving Cancer.

The change in cell state from normal to malignant is driven fundamentally by oncogenic mutations in cooperation with epigenetic alterations of chromatin. These alterations in chromatin can be a consequence of environmental stressors or germline and/or somatic mutations that directly alter the structure of chromatin machinery proteins, their levels, or their regulatory function. These changes can result in an inability of the cell to differentiate along a predefined lineage path, or drive a hyperactive, highly proliferative state with addiction to high levels of transcriptional output. We discuss how these genetic alterations hijack the chromatin machinery for the oncogenic process to reveal unique vulnerabilities and novel targets for cancer therapy.