PAX3-FOXO1 coordinates enhancer architecture, eRNA transcription, and RNA polymerase pause release at select gene targets

PAX3-FOXO1 drives targetable cell state-dependent metabolic vulnerabilities in rhabdomyosarcoma

PAX3-FOXO1, an oncogenic transcription factor, drives a particularly aggressive subtype of rhabdomyosarcoma (RMS) by enforcing gene expression programs that support malignant cell states. Here we show that PAX3-FOXO1 + RMS cells exhibit altered pyrimidine metabolism and increased dependence on enzymes involved in de novo pyrimidine synthesis, including dihydrofolate reductase (DHFR). Consequently, PAX3-FOXO1 + cells display increased sensitivity to inhibition of DHFR by the chemotherapeutic drug methotrexate, and this dependence is rescued by provision of pyrimidine nucleotides. Methotrexate treatment mimics the metabolic and transcriptional impact of PAX3-FOXO1 silencing, reducing expression of genes related to PAX3-FOXO1-driven malignant cell states. Accordingly, methotrexate treatment slows growth of multiple PAX3-FOXO1 + tumor xenograft models, but not fusion-negative counterparts. Taken together, these data demonstrate that PAX3-FOXO1 induces cell states characterized by altered pyrimidine dependence and nominate methotrexate as an addition to the current therapeutic arsenal for treatment of these malignant pediatric tumors.

Developmental Heterogeneity of Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric embryonal solid tumor and the most common pediatric soft tissue sarcoma. The histology and transcriptome of RMS resemble skeletal muscle progenitor cells that have failed to terminally differentiate. Thus, RMS is typically thought to arise from corrupted skeletal muscle progenitor cells during development. However, RMS can occur in body regions devoid of skeletal muscle, suggesting the potential for nonmyogenic cells of origin. Here, we discuss the interplay between RMS driver mutations and cell(s) of origin with an emphasis on driving location specificity. Additionally, we discuss the mechanisms governing RMS transformation events and tumor heterogeneity through the lens of transcriptional networks and epigenetic control. Finally, we reimagine Waddington's developmental landscape to include a plane of transformation connecting distinct lineage landscapes to more accurately reflect the phenomena observed in pediatric cancers.

Dependence of PAX3-FOXO1 chromatin occupancy on ETS1 at important disease-promoting genes exposes new targetable vulnerability in Fusion-Positive Rhabdomyosarcoma

Rhabdomyosarcoma (RMS), a malignancy of impaired myogenic differentiation, is the most common soft tissue pediatric cancer. PAX3-FOXO1 oncofusions drive the majority of the clinically more aggressive fusion-positive rhabdomyosarcoma (FP-RMS). Recent studies have established an epigenetic basis for PAX3-FOXO1-driven oncogenic processes. However, details of PAX3-FOXO1 epigenetic mechanisms, including interactions with, and dependence on, other chromatin and transcription factors, are incompletely understood. We previously identified a novel disease-promoting epigenetic axis in RMS, involving the histone demethylase KDM3A and the ETS1 transcription factor, and demonstrated that this epigenetic axis interfaces with PAX3-FOXO1 both phenotypically and transcriptomically, including co-regulation of biological processes and genes important to FP-RMS progression. In this study, we demonstrate that KDM3A and ETS1 colocalize with PAX3-FOXO1 to enhancers of important disease-promoting genes in FP-RMS, including FGF8, IL4R, and MEST, as well as PODXL, which we define herein as a new FP-RMS-promoting gene. We show that ETS1, which is induced by both PAX3-FOXO1 and KDM3A, exists in complex with PAX3-FOXO1, and augments PAX3-FOXO1 chromatin occupancy. We further show that the PAX3-FOXO1/ETS1 complex can be disrupted by the clinically relevant small molecule inhibitor YK-4-279. YK-4-279 displaces PAX3-FOXO1 from chromatin and interferes with PAX3-FOXO1-dependent gene regulation, resulting in potent inhibition of growth and invasive properties in FP-RMS, along with downregulation of FGF8, IL4R, MEST and PODXL expression. We additionally show that, in some FP-RMS, KDM3A also increases PAX3-FOXO1 levels. Together, our studies illuminate mechanisms of action of the KDM3A/ETS1 regulatory module, and reveal novel targetable mechanisms of PAX3-FOXO1 chromatin complex regulation, in FP-RMS.

Cross-species regulatory network analysis identifies FOXO1 as a driver of ovarian follicular recruitment

The transcriptional regulation of gene expression in the latter stages of follicular development in laying hen ovarian follicles is not well understood. Although differentially expressed genes (DEGs) have been identified in pre-recruitment and pre-ovulatory stages, the master regulators driving these DEGs remain unknown. This study addresses this knowledge gap by utilizing Master Regulator Analysis (MRA) combined with the Algorithm for the Reconstruction of Accurate Cellular Networks (ARACNe) for the first time in laying hen research to identify master regulators that are controlling DEGs in pre-recruitment and pre-ovulatory phases. The constructed ARACNe network included 10,466 nodes and 292,391 edges. The ARACNe network was then used in conjunction with the Virtual Inference of Protein-activity by Enriched Regulon (VIPER) for the MRA to identify top up- and down-regulated master regulators. VIPER analysis revealed FOXO1 as a master regulator, influencing 275 DEGs and impacting pathways related to apoptosis, proliferation, and hormonal regulation. Additionally, CLOCK, known as a crucial regulator of circadian rhythm, emerged as an upregulated master regulator in the pre-ovulatory stage. These findings provide new insights into the transcriptional landscape of laying hen ovarian follicles, offering a foundation for further exploration of follicle development and enhancing reproductive efficiency in avian species.

The Capicua C1 Domain Is Required for Full Activity of the CIC::DUX4 Fusion Oncoprotein

SIGNIFICANCE

We show in mammalian settings that the capicua C1 functional domain is a supercharger for CIC::DUX4, a poorly studied fusion oncoprotein which drives a rare sarcoma with dismal outcomes.

Pharmacological restriction of genomic binding sites redirects PU.1 pioneer transcription factor activity

Transcription factor (TF) DNA-binding dynamics govern cell fate and identity. However, our ability to pharmacologically control TF localization is limited. Here we leverage chemically driven binding site restriction leading to robust and DNA-sequence-specific redistribution of PU.1, a pioneer TF pertinent to many hematopoietic malignancies. Through an innovative technique, 'CLICK-on-CUT&Tag', we characterize the hierarchy of de novo PU.1 motifs, predicting occupancy in the PU.1 cistrome under binding site restriction. Temporal and single-molecule studies of binding site restriction uncover the pioneering dynamics of native PU.1 and identify the paradoxical activation of an alternate target gene set driven by PU.1 localization to second-tier binding sites. These transcriptional changes were corroborated by genetic blockade and site-specific reporter assays. Binding site restriction and subsequent PU.1 network rewiring causes primary human leukemia cells to differentiate. In summary, pharmacologically induced TF redistribution can be harnessed to govern TF localization, actuate alternate gene networks and direct cell fate.

The Capicua C1 Domain is Required for Full Activity of the CIC::DUX4 Fusion Oncoprotein

Rearrangements between genes can yield neomorphic fusions that drive oncogenesis. Fusion oncogenes are made up of fractional segments of the partner genes that comprise them, with each partner potentially contributing some of its own function to the nascent fusion oncoprotein. Clinically, fusion oncoproteins driving one diagnostic entity are typically clustered into a single molecular subset and are often treated a similar fashion. However, knowledge of where specific fusion breakpoints occur in partner genes, and the resulting retention of functional domains in the fusion, is an important determinant of fusion oncoprotein activity and may differ between patients. This study investigates this phenomena through the example of CIC::DUX4, a fusion between the transcriptional repressor capicua (CIC) and the double homeobox 4 gene (DUX4), which drives an aggressive subset of undifferentiated round cell sarcoma. Using a harmonized dataset of over 100 patient fusion breakpoints from the literature, we show that most bona fide CIC::DUX4 fusions retain the C1 domain, which is known to contribute to DNA binding by wild type CIC. Mechanistically, deletion or mutation of the C1 domain reduces, but does not eliminate, activation of CIC target genes by CIC::DUX4. We also find that expression of C1-deleted CIC::DUX4 is capable of exerting intermediate transformation-related phenotypes compared with those imparted by full-length CIC::DUX4, but was not sufficient for tumorigenesis in a subcutaneous mouse model. In summary, our results suggest a supercharging role for the C1 domain in the activity of CIC::DUX4.

Mutant FOXO1 controls an oncogenic network via enhancer accessibility

Transcriptional dysregulation is a hallmark of diffuse large B cell lymphoma (DLBCL), as transcriptional regulators are frequently mutated. However, our mechanistic understanding of how normal transcriptional programs are co-opted in DLBCL has been hindered by a lack of methodologies that provide the temporal resolution required to separate direct and indirect effects on transcriptional control. We applied a chemical-genetic approach to engineer the inducible degradation of the transcription factor FOXO1, which is recurrently mutated (mFOXO1) in DLBCL. The combination of rapid degradation of mFOXO1, nascent transcript detection, and assessment of chromatin accessibility allowed us to identify the direct targets of mFOXO1. mFOXO1 was required to maintain accessibility at specific enhancers associated with multiple oncogenes, and mFOXO1 degradation impaired RNA polymerase pause-release at some targets. Wild-type FOXO1 appeared to weakly regulate many of the same targets as mFOXO1 and was able to complement the degradation of mFOXO1 in the context of AKT inhibition.

eRNAbase: a comprehensive database for decoding the regulatory eRNAs in human and mouse

Enhancer RNAs (eRNAs) transcribed from distal active enhancers serve as key regulators in gene transcriptional regulation. The accumulation of eRNAs from multiple sequencing assays has led to an urgent need to comprehensively collect and process these data to illustrate the regulatory landscape of eRNAs. To address this need, we developed the eRNAbase (http://bio.liclab.net/eRNAbase/index.php) to store the massive available resources of human and mouse eRNAs and provide comprehensive annotation and analyses for eRNAs. The current version of eRNAbase cataloged 10 399 928 eRNAs from 1012 samples, including 858 human samples and 154 mouse samples. These eRNAs were first identified and uniformly processed from 14 eRNA-related experiment types manually collected from GEO/SRA and ENCODE. Importantly, the eRNAbase provides detailed and abundant (epi)genetic annotations in eRNA regions, such as super enhancers, enhancers, common single nucleotide polymorphisms, expression quantitative trait loci, transcription factor binding sites, CRISPR/Cas9 target sites, DNase I hypersensitivity sites, chromatin accessibility regions, methylation sites, chromatin interactions regions, topologically associating domains and RNA spatial interactions. Furthermore, the eRNAbase provides users with three novel analyses including eRNA-mediated pathway regulatory analysis, eRNA-based variation interpretation analysis and eRNA-mediated TF-target gene analysis. Hence, eRNAbase is a powerful platform to query, browse and visualize regulatory cues associated with eRNAs.

Chemical-genetics refines transcription factor regulatory circuits

Transcriptional dysregulation is a key step in oncogenesis, but our understanding of transcriptional control has relied on genetic approaches that are slow and allow for compensation. Chemical-genetic approaches have shortened the time frame for the analysis of transcription factors from days or weeks to minutes. These studies show that while DNA-binding proteins bind to thousands of sites, they are directly required to regulate only a small cadre of genes. Moreover, these transcriptional control networks are far more distinct, with much less overlap and interconnectivity than predicted from DNA binding. The identified direct targets can then be used to dissect the mechanism of action of these factors, which could identify ways to therapeutically manipulate these oncogenic transcriptional control networks.

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.

PROTAC-Based Protein Degradation as a Promising Strategy for Targeted Therapy in Sarcomas

Sarcomas are heterogeneous bone and soft tissue cancers representing the second most common tumor type in children and adolescents. Histology and genetic profiling discovered more than 100 subtypes, which are characterized by peculiar molecular vulnerabilities. However, limited therapeutic options exist beyond standard therapy and clinical benefits from targeted therapies were observed only in a minority of patients with sarcomas. The rarity of these tumors, paucity of actionable mutations, and limitations in the chemical composition of current targeted therapies hindered the use of these approaches in sarcomas. Targeted protein degradation (TPD) is an innovative pharmacological modality to directly alter protein abundance with promising clinical potential in cancer, even for undruggable proteins. TPD is based on the use of small molecules called degraders or proteolysis-targeting chimeras (PROTACs), which trigger ubiquitin-dependent degradation of protein of interest. In this review, we will discuss major features of PROTAC and PROTAC-derived genetic systems for target validation and cancer treatment and focus on the potential of these approaches to overcome major issues connected to targeted therapies in sarcomas, including drug resistance, target specificity, and undruggable targets. A deeper understanding of these strategies might provide new fuel to drive molecular and personalized medicine to sarcomas.

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.

Differentiation-dependent chromosomal organization changes in normal myogenic cells are absent in rhabdomyosarcoma cells

Myogenesis, the progression of proliferating skeletal myoblasts to terminally differentiated myotubes, regulates thousands of target genes. Uninterrupted linear arrays of such genes are differentially associated with specific chromosomes, suggesting chromosome specific regulatory roles in myogenesis. Rhabdomyosarcoma (RMS), a tumor of skeletal muscle, shares common features with normal muscle cells. We hypothesized that RMS and myogenic cells possess differences in chromosomal organization related to myogenic gene arrangement. We compared the organizational characteristics of chromosomes 2 and 18, chosen for their difference in myogenic gene arrangement, in cultured RMS cell lines and normal myoblasts and myotubes. We found chromosome-specific differences in organization during normal myogenesis, with increased area occupied and a shift in peripheral localization specifically for chromosome 2. Most strikingly, we found a differentiation-dependent difference in positioning of chromosome 2 relative to the nuclear axis, with preferential positioning along the major nuclear axis present only in myotubes. RMS cells demonstrated no preference for such axial positioning, but induced differentiation through transfection of the pro-myogenic miRNA miR-206 resulted in an increase of major axial positioning of chromosome 2. Our findings identify both a differentiation-dependent, chromosome-specific change in organization in normal myogenesis, and highlight the role of chromosomal spatial organization in myogenic differentiation.

Emerging insights into enhancer biology and function

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

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.

Targeting epigenetic aberrations of sarcoma in CRISPR era

Sarcomas are rare malignancies that exhibit diverse biological, genetic, morphological, and clinical characteristics. Genetic alterations, such as gene fusions, mutations in transcriptional machinery components, histones, and DNA methylation regulatory molecules, play an essential role in sarcomagenesis. These mutations induce and/or cooperate with specific epigenetic aberrations required for the growth and maintenance of sarcomas. Appropriate mouse models have been developed to clarify the significance of genetic and epigenetic interactions in sarcomas. Studies using the mouse models for human sarcomas have demonstrated major advances in our understanding the developmental processes as well as tumor microenvironment of sarcomas. Recent technological progresses in epigenome editing will not only improve the studies using animal models but also provide a direct clue for epigenetic therapies. In this manuscript, we review important epigenetic aberrations in sarcomas and their representative mouse models, current methods of epigenetic editing using CRISPR/dCas9 systems, and potential applications in sarcoma studies and therapeutics.

Differentiation-dependent chromosomal organization changes in normal myogenic cells are absent in rhabdomyosarcoma cells

Myogenesis, the progression of proliferating skeletal myoblasts to terminally differentiated myotubes, regulates thousands of target genes. Uninterrupted linear arrays of such genes are differentially associated with specific chromosomes, suggesting chromosome specific regulatory roles in myogenesis. Rhabdomyosarcoma (RMS), a tumor of skeletal muscle, shares common features with normal muscle cells. We hypothesized that RMS and myogenic cells possess differences in chromosomal organization related to myogenic gene arrangement. We compared the organizational characteristics of chromosomes 2 and 18, chosen for their difference in myogenic gene arrangement, in cultured RMS cell lines and normal myoblasts and myotubes. We found chromosome-specific differences in organization during normal myogenesis, with increased area occupied and a shift in peripheral localization specifically for chromosome 2. Most strikingly, we found a differentiation-dependent difference in positioning of chromosome 2 relative to the nuclear axis, with preferential positioning along the major nuclear axis present only in myotubes. RMS cells demonstrated no preference for such axial positioning, but induced differentiation through transfection of the pro-myogenic miRNA miR-206 resulted in an increase of major axial positioning of chromosome 2. Our findings identify both a differentiation-dependent, chromosome-specific change in organization in normal myogenesis, and highlight the role of chromosomal spatial organization in myogenic differentiation.

Super-enhancer landscape rewiring in cancer: The epigenetic control at distal sites

Super-enhancers evolve as elements at the top of the hierarchical control of gene expression. They are important end-gatherers of signaling pathways that control stemness, differentiation or adaptive responses. Many epigenetic regulations focus on these regions, and not surprisingly, during the process of tumorigenesis, various alterations can account for their dysfunction. Super-enhancers are emerging as key drivers of the aberrant gene expression landscape that sustain the aggressiveness of cancer cells. In this review, we will describe and discuss about the structure of super-enhancers, their epigenetic regulation, and the major changes affecting their functionality in cancer.

Rapid-kinetics degron benchmarking reveals off-target activities and mixed agonism-antagonism of MYB inhibitors

Attenuating aberrant transcriptional circuits holds great promise for the treatment of numerous diseases, including cancer. However, development of transcriptional inhibitors is hampered by the lack of a generally accepted functional cellular readout to characterize their target specificity and on-target activity. We benchmarked the direct gene-regulatory signatures of six agents reported as inhibitors of the oncogenic transcription factor MYB against targeted MYB degradation in a nascent transcriptomics assay. The inhibitors demonstrated partial specificity for MYB target genes but displayed significant off-target activity. Unexpectedly, the inhibitors displayed bimodal on-target effects, acting as mixed agonists-antagonists. Our data uncover unforeseen agonist effects of small molecules originally developed as TF inhibitors and argue that rapid-kinetics benchmarking against degron models should be used for functional characterization of transcriptional modulators.