Epigenetic determinants of fusion-driven sarcomas: paradigms and challenges.

We describe exciting recent advances in fusion-driven sarcoma etiology, from an epigenetics perspective. By exploring the current state of the field, we identify and describe the central mechanisms that determine sarcomagenesis. Further, we discuss seminal studies in translational genomics, which enabled epigenetic characterization of fusion-driven sarcomas. Important context for epigenetic mechanisms include, but are not limited to, cell cycle and metabolism, core regulatory circuitry, 3-dimensional chromatin architectural dysregulation, integration with ATP-dependent chromatin remodeling, and translational animal modeling. Paradoxically, while the genetic requirements for oncogenic transformation are highly specific for the fusion partners, the epigenetic mechanisms we as a community have uncovered are categorically very broad. This dichotomy prompts the question of whether the investigation of rare disease epigenomics should prioritize studying individual cell populations, thereby examining whether the mechanisms of chromatin dysregulation are specific to a particular tumor. We review recent advances focusing on rhabdomyosarcoma, synovial sarcoma, alveolar soft part sarcoma, clear cell sarcoma, undifferentiated round cell sarcoma, Ewing sarcoma, myxoid/round liposarcoma, epithelioid hemangioendothelioma and desmoplastic round cell tumor. The growing number of groundbreaking discoveries in the field, motivated us to anticipate further exciting advances in the area of mechanistic epigenomics and direct targeting of fusion transcription factors in the years ahead.

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.

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.

The 3D chromatin landscape of rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric soft tissue cancer with a lack of precision therapy options for patients. We hypothesized that with a general paucity of known mutations in RMS, chromatin structural driving mechanisms are essential for tumor proliferation. Thus, we carried out high-depth in situ Hi-C in representative cell lines and patient-derived xenografts (PDXs) to define chromatin architecture in each major RMS subtype. We report a comprehensive 3D chromatin structural analysis and characterization of fusion-positive (FP-RMS) and fusion-negative RMS (FN-RMS). We have generated spike-in in situ Hi-C chromatin interaction maps for the most common FP-RMS and FN-RMS cell lines and compared our data with PDX models. In our studies, we uncover common and distinct structural elements in large Mb-scale chromatin compartments, tumor-essential genes within variable topologically associating domains and unique patterns of structural variation. Our high-depth chromatin interactivity maps and comprehensive analyses provide context for gene regulatory events and reveal functional chromatin domains in RMS.

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.

Chromatin structure in cancer.

In the past decade, we have seen the emergence of sequence-based methods to understand chromosome organization. With the confluence of in situ approaches to capture information on looping, topological domains, and larger chromatin compartments, understanding chromatin-driven disease is becoming feasible. Excitingly, recent advances in single molecule imaging with capacity to reconstruct "bulk-cell" features of chromosome conformation have revealed cell-to-cell chromatin structural variation. The fundamental question motivating our analysis of the literature is, can altered chromatin structure drive tumorigenesis? As our community learns more about rare disease, including low mutational frequency cancers, understanding "chromatin-driven" pathology will illuminate the regulatory structures of the genome. We describe recent insights into altered genome architecture in human cancer, highlighting multiple pathways toward disruptions of chromatin structure, including structural variation, noncoding mutations, metabolism, and de novo mutations to architectural regulators themselves. Our analysis of the literature reveals that deregulation of genome structure is characteristic in distinct classes of chromatin-driven tumors. As we begin to integrate the findings from single cell imaging studies and chromatin structural sequencing, we will be able to understand the diversity of cells within a common diagnosis, and begin to define structure-function relationships of the misfolded genome.

Succinate dehydrogenase/complex II is critical for metabolic and epigenetic regulation of T cell proliferation and inflammation.

Effective T cell-mediated immune responses require the proper allocation of metabolic resources to sustain growth, proliferation, and cytokine production. Epigenetic control of the genome also governs T cell transcriptome and T cell lineage commitment and maintenance. Cellular metabolic programs interact with epigenetic regulation by providing substrates for covalent modifications of chromatin. By using complementary genetic, epigenetic, and metabolic approaches, we revealed that tricarboxylic acid (TCA) cycle flux fueled biosynthetic processes while controlling the ratio of succinate/α-ketoglutarate (α-KG) to modulate the activities of dioxygenases that are critical for driving T cell inflammation. In contrast to cancer cells, where succinate dehydrogenase (SDH)/complex II inactivation drives cell transformation and growth, SDH/complex II deficiency in T cells caused proliferation and survival defects when the TCA cycle was truncated, blocking carbon flux to support nucleoside biosynthesis. Replenishing the intracellular nucleoside pool partially relieved the dependence of T cells on SDH/complex II for proliferation and survival. SDH deficiency induced a proinflammatory gene signature in T cells and promoted T helper 1 and T helper 17 lineage differentiation. An increasing succinate/α-KG ratio in SDH-deficient T cells promoted inflammation by changing the pattern of the transcriptional and chromatin accessibility signatures and consequentially increasing the expression of the transcription factor, PR domain zinc finger protein 1. Collectively, our studies revealed a role of SDH/complex II in allocating carbon resources for anabolic processes and epigenetic regulation in T cell proliferation and inflammation.

BAF complexes drive proliferation and block myogenic differentiation in fusion-positive rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of skeletal muscle lineage. The aggressive alveolar subtype is characterized by t(2;13) or t(1;13) translocations encoding for PAX3- or PAX7-FOXO1 chimeric transcription factors, respectively, and are referred to as fusion positive RMS (FP-RMS). The fusion gene alters the myogenic program and maintains the proliferative state while blocking terminal differentiation. Here, we investigated the contributions of chromatin regulatory complexes to FP-RMS tumor maintenance. We define the mSWI/SNF functional repertoire in FP-RMS. We find that SMARCA4 (encoding BRG1) is overexpressed in this malignancy compared to skeletal muscle and is essential for cell proliferation. Proteomic studies suggest proximity between PAX3-FOXO1 and BAF complexes, which is further supported by genome-wide binding profiles revealing enhancer colocalization of BAF with core regulatory transcription factors. Further, mSWI/SNF complexes localize to sites of de novo histone acetylation. Phenotypically, interference with mSWI/SNF complex function induces transcriptional activation of the skeletal muscle differentiation program associated with MYCN enhancer invasion at myogenic target genes, which is recapitulated by BRG1 targeting compounds. We conclude that inhibition of BRG1 overcomes the differentiation blockade of FP-RMS cells and may provide a therapeutic strategy for this lethal childhood tumor.

Epigenetic regulation of nuclear lamina-associated heterochromatin by HAT1 and the acetylation of newly synthesized histones.

A central component of the epigenome is the pattern of histone post-translational modifications that play a critical role in the formation of specific chromatin states. Following DNA replication, nascent chromatin is a 1:1 mixture of parental and newly synthesized histones and the transfer of modification patterns from parental histones to new histones is a fundamental step in epigenetic inheritance. Here we report that loss of HAT1, which acetylates lysines 5 and 12 of newly synthesized histone H4 during replication-coupled chromatin assembly, results in the loss of accessibility of large domains of heterochromatin, termed HAT1-dependent Accessibility Domains (HADs). HADs are mega base-scale domains that comprise ∼10% of the mouse genome. HAT1 globally represses H3 K9 me3 levels and HADs correspond to the regions of the genome that display HAT1-dependent increases in H3 K9me3 peak density. HADs display a high degree of overlap with a subset of Lamin-Associated Domains (LADs). HAT1 is required to maintain nuclear structure and integrity. These results indicate that HAT1 and the acetylation of newly synthesized histones may be critical regulators of the epigenetic inheritance of heterochromatin and suggest a new mechanism for the epigenetic regulation of nuclear lamina-heterochromatin interactions.

Pioneer factors in development and cancer.

Transcription factors (TFs) are essential mediators of epigenetic regulation and modifiers of penetrance. Studies from the past decades have revealed a sub-class of TF that is capable of remodeling closed chromatin states through targeting nucleosomal motifs. This pioneer factor (PF) class of chromatin remodeler is ATP independent in its roles in epigenetic initiation, with nucleosome-motif recognition and association with repressive chromatin regions. Increasing evidence suggests that the fundamental properties of PFs can be coopted in human cancers. We explore the role of PFs in the larger context of tissue-specific epigenetic regulation. Moreover, we highlight an emerging class of chimeric PF derived from translocation partners in human disease and PFs associated with rare tumors. In the age of site-directed genome editing and targeted protein degradation, increasing our understanding of PFs will provide access to next-generation therapy for human disease driven from altered transcriptional circuitry.

SNAI2-Mediated Repression of BIM Protects Rhabdomyosarcoma from Ionizing Radiation.

Ionizing radiation (IR) and chemotherapy are mainstays of treatment for patients with rhabdomyosarcoma, yet the molecular mechanisms that underlie the success or failure of radiotherapy remain unclear. The transcriptional repressor SNAI2 was previously identified as a key regulator of IR sensitivity in normal and malignant stem cells through its repression of the proapoptotic BH3-only gene PUMA/BBC3. Here, we demonstrate a clear correlation between SNAI2 expression levels and radiosensitivity across multiple rhabdomyosarcoma cell lines. Modulating SNAI2 levels in rhabdomyosarcoma cells through its overexpression or knockdown altered radiosensitivity in vitro and in vivo. SNAI2 expression reliably promoted overall cell growth and inhibited mitochondrial apoptosis following exposure to IR, with either variable or minimal effects on differentiation and senescence, respectively. Importantly, SNAI2 knockdown increased expression of the proapoptotic BH3-only gene BIM, and chromatin immunoprecipitation sequencing experiments established that SNAI2 is a direct repressor of BIM/BCL2L11. Because the p53 pathway is nonfunctional in the rhabdomyosarcoma cells used in this study, we have identified a new, p53-independent SNAI2/BIM signaling axis that could potentially predict clinical responses to IR treatment and be exploited to improve rhabdomyosarcoma therapy. SIGNIFICANCE: SNAI2 is identified as a major regulator of radiation-induced apoptosis in rhabdomyosarcoma through previously unknown mechanisms independent of p53.

Evidence of pioneer factor activity of an oncogenic fusion transcription factor.

Recent characterizations of pioneer transcription factors provide insights into their structures and patterns of chromatin recognition associated with their roles in cell fate commitment and transformation. Intersecting with these basic science concepts, identification of pioneer factors (PFs) fused together as driver translocations in childhood cancers raises questions of whether these fusions retain the fundamental ability to invade repressed chromatin, consistent with their monomeric PF constituents. This study defines the cellular and chromatin localization of PAX3-FOXO1, an oncogenic driver of childhood rhabdomyosarcoma (RMS), derived from a fusion of PFs. To quantitatively define its chromatin-targeting functions and capacity to drive epigenetic reprogramming, we developed a ChIP-seq workflow with per-cell normalization (pc-ChIP-seq). Our quantitative localization studies address structural variation in RMS genomes and reveal insights into inactive chromatin localization of PAX3-FOXO1. Taken together, our studies are consistent with pioneer function for a driver oncoprotein in RMS, with repressed chromatin binding and nucleosome-motif targeting.

Inducible Protein Degradation to Understand Genome Architecture.

We review exciting recent advances in protein degradation, with a focus on chromatin structure. In our analysis of the literature, we highlight studies of kinetic control of protein stability for cohesin, condensin, ATP-dependent chromatin remodeling, and pioneer transcription factors. With new connections emerging between chromatin remodeling and genome structure, we anticipate exciting developments at the intersection of these topics to be revealed in the coming years. Moreover, we pay special attention to the 20-year anniversary of PROTACs, with an overview of E3 ligase/target pairings and central questions that might lead to the next generation of PROTACs with an expanded scope and generality. While steady-state experimental measurements with constitutive genome editing are impactful, we highlight complementary approaches for rapid kinetic protein degradation to uncover early targeting functions and to understand the central determinants of genome structure-function relationships.

mSWI/SNF promotes Polycomb repression both directly and through genome-wide redistribution.

The mammalian SWI/SNF complex, or BAF complex, has a conserved and direct role in antagonizing Polycomb-mediated repression. Yet, BAF also promotes repression by Polycomb in stem cells and cancer. How BAF both antagonizes and promotes Polycomb-mediated repression remains unknown. Here, we utilize targeted protein degradation to dissect the BAF-Polycomb axis in mouse embryonic stem cells on short timescales. We report that rapid BAF depletion redistributes Polycomb repressive complexes PRC1 and PRC2 from highly occupied domains, like Hox clusters, to weakly occupied sites normally opposed by BAF. Polycomb redistribution from highly repressed domains results in their decompaction, gain of active epigenomic features and transcriptional derepression. Surprisingly, through dose-dependent degradation of PRC1 and PRC2, we identify a conventional role for BAF in Polycomb-mediated repression, in addition to global Polycomb redistribution. These findings provide new mechanistic insight into the highly dynamic state of the Polycomb-Trithorax axis.

The FLI portion of EWS/FLI contributes a transcriptional regulatory function that is distinct and separable from its DNA-binding function in Ewing sarcoma.

Ewing sarcoma is an aggressive bone cancer of children and young adults defined by the presence of a chromosomal translocation: t(11;22)(q24;q12). The encoded protein, EWS/FLI, fuses the amino-terminal domain of EWS to the carboxyl-terminus of FLI. The EWS portion is an intrinsically disordered transcriptional regulatory domain, while the FLI portion contains an ETS DNA-binding domain and two flanking regions of unknown function. Early studies using non-Ewing sarcoma models provided conflicting information on the roles of each domain of FLI in EWS/FLI oncogenic function. We therefore sought to define the specific contributions of each FLI domain to EWS/FLI activity in a well-validated Ewing sarcoma model and, in doing so, to better understand Ewing sarcoma development mediated by the fusion protein. We analyzed a series of engineered EWS/FLI mutants with alterations in the FLI portion using a variety of assays. Fluorescence anisotropy, CUT&RUN, and ATAC-sequencing experiments revealed that the isolated ETS domain is sufficient to maintain the normal DNA-binding and chromatin accessibility function of EWS/FLI. In contrast, RNA-sequencing and soft agar colony formation assays revealed that the ETS domain alone was insufficient for transcriptional regulatory and oncogenic transformation functions of the fusion protein. We found that an additional alpha-helix immediately downstream of the ETS domain is required for full transcriptional regulation and EWS/FLI-mediated oncogenesis. These data demonstrate a previously unknown role for FLI in transcriptional regulation that is distinct from its DNA-binding activity. This activity is critical for the cancer-causing function of EWS/FLI and may lead to novel therapeutic approaches.

Interaction between SNAI2 and MYOD enhances oncogenesis and suppresses differentiation in Fusion Negative Rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is an aggressive pediatric malignancy of the muscle, that includes Fusion Positive (FP)-RMS harboring PAX3/7-FOXO1 and Fusion Negative (FN)-RMS commonly with RAS pathway mutations. RMS express myogenic master transcription factors MYOD and MYOG yet are unable to terminally differentiate. Here, we report that SNAI2 is highly expressed in FN-RMS, is oncogenic, blocks myogenic differentiation, and promotes growth. MYOD activates SNAI2 transcription via super enhancers with striped 3D contact architecture. Genome wide chromatin binding analysis demonstrates that SNAI2 preferentially binds enhancer elements and competes with MYOD at a subset of myogenic enhancers required for terminal differentiation. SNAI2 also suppresses expression of a muscle differentiation program modulated by MYOG, MEF2, and CDKN1A. Further, RAS/MEK-signaling modulates SNAI2 levels and binding to chromatin, suggesting that the differentiation blockade by oncogenic RAS is mediated in part by SNAI2. Thus, an interplay between SNAI2, MYOD, and RAS prevents myogenic differentiation and promotes tumorigenesis.

Miswired Enhancer Logic Drives a Cancer of the Muscle Lineage.

Core regulatory transcription factors (CR TFs) establish enhancers with logical ordering during embryogenesis and development. Here we report that in fusion-positive rhabdomyosarcoma, a cancer of the muscle lineage, the chief oncogene PAX3-FOXO1 is driven by a translocated FOXO1 super enhancer (SE) restricted to a late stage of myogenesis. Using chromatin conformation capture techniques, we demonstrate that the extensive FOXO1 cis-regulatory domain interacts with PAX3. Furthermore, RNA sequencing and chromatin immunoprecipitation sequencing data in tumors bearing rare PAX translocations implicate enhancer miswiring across all fusion-positive tumors. HiChIP of H3K27ac showed connectivity between the FOXO1 SE, additional intra-domain enhancers, and the PAX3 promoter. We show that PAX3-FOXO1 transcription is diminished when this network of enhancers is ablated by CRISPR. Our data reveal a hijacked enhancer network that disrupts the stepwise CR TF logic of normal skeletal muscle development (PAX3 to MYOD to MYOG), replacing it with an "infinite loop" enhancer logic that locks rhabdomyosarcoma in an undifferentiated stage.

Measurement of differential chromatin interactions with absolute quantification of architecture (AQuA-HiChIP).

Methods developed to capture protein-anchored chromatin interactions (chromatin interaction analysis by paired-end tag sequencing and HiChIP) have yielded tremendous insights into the 3D folding principles of the genome, but are normalized by sequencing depth and therefore unable to accurately measure global changes in chromatin interactions and contact domain organization. We herein describe the protocol for absolute quantification of chromatin architecture (AQuA)-HiChIP, an advance that allows the absolute differences in protein-anchored chromatin interactions between samples to be determined. With our method, defined ratios of mouse and human fixed nuclei are mixed and subjected to endonuclease digestion. Chromatin contacts are captured by biotin-dATP incorporation and proximity ligation, followed by gentle shearing, ChIP, biotin capture and paired-end sequencing. 3D contacts are counted from paired-end tags (PETs) from the human genome and are normalized to the total PETs from the mouse genome. As orthogonal normalization allows observation of global changes, the approach will enable more quantitative insights into the topological determinants of transcriptional control and tissue-specific epigenetic memory. With our approach, we have discovered that rapid histone deacetylase inhibition disrupts super enhancer function by creating many new aberrant contacts. The code for data analysis is available at https://github.com/GryderArt/AQuA-HiChIP. This protocol reports both experimental and bioinformatic details to perform AQuA-HiChIP, going from cell culture to ranking chromatin interactions within 6 d.