Regulatory landscape fusion in rhabdomyosarcoma through interactions between the PAX3 promoter and FOXO1 regulatory elements

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.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

CRISPR-Cas9-Mediated Bioluminescent Tagging of Endogenous Proteins by Fluorescent Protein-Assisted Cell Sorting

Oncogenic fusion genes are attractive therapeutic targets because of their tumor-specific expression and central "driver" roles in various human cancers. However, oncogenic fusions involving transcription factors such as PAX3-FOXO1 in alveolar fusion gene-positive rhabdomyosarcoma (FP-RMS) have been difficult to inhibit due to the apparent lack of tractable drug-like binding sites comparable to that recognized by Gleevec (imatinib mesylate) on the BCR-ABL1 tyrosine kinase fusion protein. Toward the identification of novel small molecules that selectively target PAX3-FOXO1, we used CRISPR-Cas9-mediated knock-in to append the pro-luminescent HiBiT tag onto the carboxy terminus of the endogenous PAX3-FOXO1 fusion protein in two human FP-RMS cell lines (RH4 and SCMC). HiBiT is an 11-amino acid peptide derived from the NanoLuc luciferase that produces a luminescence signal which is ~100-fold brighter than firefly or Renilla luciferases through high-affinity binding to a complementary NanoLuc peptide fragment called LgBiT. To facilitate single-cell clonal isolation of knock-ins, the homology-directed repair template encoding HiBiT was followed by a P2A self-cleaving peptide for coexpression of an mCherry fluorescent protein as a fluorescence-activated cell sorter (FACS)-selectable marker. HiBiT tagging thus allows highly sensitive luminescence detection of endogenous PAX3-FOXO1 levels permitting quantitative high-throughput screening of large compound libraries for the discovery of PAX3-FOXO1 inhibitors and degraders.

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.

Disruption of regulatory domains and novel transcripts as disease-causing mechanisms

Deletions, duplications, insertions, inversions, and translocations, collectively called structural variations (SVs), affect more base pairs of the genome than any other sequence variant. The recent technological advancements in genome sequencing have enabled the discovery of tens of thousands of SVs per human genome. These SVs primarily affect non-coding DNA sequences, but the difficulties in interpreting their impact limit our understanding of human disease etiology. The functional annotation of non-coding DNA sequences and methodologies to characterize their three-dimensional (3D) organization in the nucleus have greatly expanded our understanding of the basic mechanisms underlying gene regulation, thereby improving the interpretation of SVs for their pathogenic impact. Here, we discuss the various mechanisms by which SVs can result in altered gene regulation and how these mechanisms can result in rare genetic disorders. Beyond changing gene expression, SVs can produce novel gene-intergenic fusion transcripts at the SV breakpoints.

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.

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 regulatory mechanisms of noncoding RNAs in topologically associating domains

Topologically associating domains (TADs) are integral to spatial genome organization, instructing gene expression, and cell fate. Recently, several advances have uncovered roles for noncoding RNAs (ncRNAs) in the regulation of the form and function of mammalian TADs. Phase separation has also emerged as a potential arbiter of ncRNAs in the regulation of TADs. In this review we discuss the implications of these novel findings in relation to how ncRNAs might structurally and functionally regulate TADs from two perspectives: moderating loop extrusion through interactions with architectural proteins, and facilitating TAD phase separation. Additionally, we propose future studies and directions to investigate these phenomena.

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

The PAX-FOXO1s trigger fast trans-differentiation of chick embryonic neural cells into alveolar rhabdomyosarcoma with tissue invasive properties limited by S phase entry inhibition

The chromosome translocations generating PAX3-FOXO1 and PAX7-FOXO1 chimeric proteins are the primary hallmarks of the paediatric fusion-positive alveolar subtype of Rhabdomyosarcoma (FP-RMS). Despite the ability of these transcription factors to remodel chromatin landscapes and promote the expression of tumour driver genes, they only inefficiently promote malignant transformation in vivo. The reason for this is unclear. To address this, we developed an in ovo model to follow the response of spinal cord progenitors to PAX-FOXO1s. Our data demonstrate that PAX-FOXO1s, but not wild-type PAX3 or PAX7, trigger the trans-differentiation of neural cells into FP-RMS-like cells with myogenic characteristics. In parallel, PAX-FOXO1s remodel the neural pseudo-stratified epithelium into a cohesive mesenchyme capable of tissue invasion. Surprisingly, expression of PAX-FOXO1s, similar to wild-type PAX3/7, reduce the levels of CDK-CYCLIN activity and increase the fraction of cells in G1. Introduction of CYCLIN D1 or MYCN overcomes this PAX-FOXO1-mediated cell cycle inhibition and promotes tumour growth. Together, our findings reveal a mechanism that can explain the apparent limited oncogenicity of PAX-FOXO1 fusion transcription factors. They are also consistent with certain clinical reports indicative of a neural origin of FP-RMS.

The fusion-positive subtype of rhabdomyosarcoma (FP-RMS) is a rare malignant paediatric cancer, whose induction and evolution still remain to be deciphered. Out of the gross genetic aberrations found in these cancers, t(2:13) and t(1,13) chromosome translocations are the first to appear and lead to the expression of fusion proteins made of the DNA binding domains of either PAX3 or PAX7 and the transactivation domain of FOXO1. Both PAX3-FOXO1 and PAX7-FOXO1 have a strong impact on gene transcription, yet they only inefficiently promote the transformation of healthy cells into tumorigenic cells. To address this issue, we have used chick embryos to monitor in vivo the early response of cells to PAX-FOXO1 chimeric proteins. We showed that both proteins, but not the normal PAX3 and PAX7, transform neural cells into cells with FP-RMS molecular features. The PAX-FOXO1s also force polarized epithelial neural cells to adopt a mesenchymal phenotype with tissue invasive properties. However, the PAX-FOXO1s inhibit cell division and hence tumour growth. Genetically re-activating core cell cycle regulators rescues PAX-FOXO1 mediated cell cycle inhibition. Together, our findings bring further support to the idea that the PAX-FOXO1s are stricto sensu oncoproteins, whose oncogenicity is limited by negative effects on cell cycle.

Chromosome structural variation in tumorigenesis: mechanisms of formation and carcinogenesis

With the rapid development of next-generation sequencing technology, chromosome structural variation has gradually gained increased clinical significance in tumorigenesis. However, the molecular mechanism(s) underlying this structural variation remain poorly understood. A search of the literature shows that a three-dimensional chromatin state plays a vital role in inducing structural variation and in the gene expression profiles in tumorigenesis. Structural variants may result in changes in copy number or deletions of coding sequences, as well as the perturbation of structural chromatin features, especially topological domains, and disruption of interactions between genes and their regulatory elements. This review focuses recent work aiming at elucidating how structural variations develop and misregulate oncogenes and tumor suppressors, to provide general insights into tumor formation mechanisms and to provide potential targets for future anticancer therapies.

Chromatin Landscape During Skeletal Muscle Differentiation

Cellular commitment and differentiation involve highly coordinated mechanisms by which tissue-specific genes are activated while others are repressed. These mechanisms rely on the activity of specific transcription factors, chromatin remodeling enzymes, and higher-order chromatin organization in order to modulate transcriptional regulation on multiple cellular contexts. Tissue-specific transcription factors are key mediators of cell fate specification with the ability to reprogram cell types into different lineages. A classic example of a master transcription factor is the muscle specific factor MyoD, which belongs to the family of myogenic regulatory factors (MRFs). MRFs regulate cell fate determination and terminal differentiation of the myogenic precursors in a multistep process that eventually culminate with formation of muscle fibers. This developmental progression involves the activation and proliferation of muscle stem cells, commitment, and cell cycle exit and fusion of mononucleated myoblast to generate myotubes and myofibers. Although the epigenetics of muscle regeneration has been extensively addressed and discussed over the recent years, the influence of higher-order chromatin organization in skeletal muscle regeneration is still a field of development. In this review, we will focus on the epigenetic mechanisms modulating muscle gene expression and on the incipient work that addresses three-dimensional genome architecture and its influence in cell fate determination and differentiation to achieve skeletal myogenesis. We will visit known alterations of genome organization mediated by chromosomal fusions giving rise to novel regulatory landscapes, enhancing oncogenic activation in muscle, such as alveolar rhabdomyosarcomas (ARMS).

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.

Analysis of 4C-seq data: A comparison of methods

The circular chromosome conformation capture technique followed by sequencing (4C-seq) has been used in a number of studies to investigate chromosomal interactions between DNA fragments. Computational pipelines have been developed and published that offer various possibilities of 4C-seq data processing and statistical analysis. Here, we present an overview of four of such pipelines (fourSig, FourCSeq, 4C-ker and w4Cseq) taking into account the most important stages of computations. We provide comparisons of the methods and discuss their advantages and possible weaknesses. We illustrate the results with the use of data obtained for two different species, in a study devoted to vernalization control in Arabidopsis thaliana by the FLOWERING LOCUS C (FLC) gene and to long-range chromatin interactions in mouse embryonic stem cells.

A Conserved Noncoding Locus Regulates Random Monoallelic Xist Expression across a Topological Boundary

cis-Regulatory communication is crucial in mammalian development and is thought to be restricted by the spatial partitioning of the genome in topologically associating domains (TADs). Here, we discovered that the Xist locus is regulated by sequences in the neighboring TAD. In particular, the promoter of the noncoding RNA Linx (LinxP) acts as a long-range silencer and influences the choice of X chromosome to be inactivated. This is independent of Linx transcription and independent of any effect on Tsix, the antisense regulator of Xist that shares the same TAD as Linx. Unlike Tsix, LinxP is well conserved across mammals, suggesting an ancestral mechanism for random monoallelic Xist regulation. When introduced in the same TAD as Xist, LinxP switches from a silencer to an enhancer. Our study uncovers an unsuspected regulatory axis for X chromosome inactivation and a class of cis-regulatory effects that may exploit TAD partitioning to modulate developmental decisions.

Gene Co-Expression Networks Restructured Gene Fusion in Rhabdomyosarcoma Cancers

Rhabdomyosarcoma is subclassified by the presence or absence of a recurrent chromosome translocation that fuses the FOXO1 and PAX3 or PAX7 genes. The fusion protein (FOXO1-PAX3/7) retains both binding domains and becomes a novel and potent transcriptional regulator in rhabdomyosarcoma subtypes. Many studies have characterized and integrated genomic, transcriptomic, and epigenomic differences among rhabdomyosarcoma subtypes that contain the FOXO1-PAX3/7 gene fusion and those that do not; however, few investigations have investigated how gene co-expression networks are altered by FOXO1-PAX3/7. Although transcriptional data offer insight into one level of functional regulation, gene co-expression networks have the potential to identify biological interactions and pathways that underpin oncogenesis and tumorigenicity. Thus, we examined gene co-expression networks for rhabdomyosarcoma that were FOXO1-PAX3 positive, FOXO1-PAX7 positive, or fusion negative. Gene co-expression networks were mined using local maximum Quasi-Clique Merger (lmQCM) and analyzed for co-expression differences among rhabdomyosarcoma subtypes. This analysis observed 41 co-expression modules that were shared between fusion negative and positive samples, of which 17/41 showed significant up- or down-regulation in respect to fusion status. Fusion positive and negative rhabdomyosarcoma showed differing modularity of co-expression networks with fusion negative (n = 109) having significantly more individual modules than fusion positive (n = 53). Subsequent analysis of gene co-expression networks for PAX3 and PAX7 type fusions observed 17/53 were differentially expressed between the two subtypes. Gene list enrichment analysis found that gene ontology terms were poorly matched with biological processes and molecular function for most co-expression modules identified in this study; however, co-expressed modules were frequently localized to cytobands on chromosomes 8 and 11. Overall, we observed substantial restructuring of co-expression networks relative to fusion status and fusion type in rhabdomyosarcoma and identified previously overlooked genes and pathways that may be targeted in this pernicious disease.

An evolutionary perspective of regulatory landscape dynamics in development and disease

The organization of animal genomes into topologically associating domains (TADs) provides a structural scaffold in which cis-regulatory elements (CREs) operate on their target genes. Determining the position of CREs and genes relative to TADs has become instrumental to trace gene expression changes during evolution and in diseases. Here we will review recent studies and discuss TADs as structural units with respect to their conservation and stability during genome reorganization. Furthermore, we describe how TAD restructuring contributed to morphological novelties during evolution but also their deleterious effects associated with disease. Despite considering TADs as structural units, the nested and dynamic scaffold within TADs contributes to tissue-specific gene expression, implying that such changes can also account for gene expression differences during evolution.

Interpreting Chromosomal Rearrangements in the Context of 3-Dimentional Genome Organization: A Practical Guide for Medical Genetics

In this exciting era of "next-gen cytogenetics", the use of novel molecular methods such as comparative genome hybridization and whole genome and whole exome sequencing becomes more and more common in clinics. This results in generation of large amounts of high-resolution patient-specific data and challenges the development of new approaches for interpretation of obtained information. Usually, interpretation of chromosomal rearrangements is focused on alterations of linear genome sequence, underestimating the role of spatial chromatin organization. In this article, we describe the main features of 3-dimentional genome organization, emphasizing their role in normal and pathological development. We highlight some tips to help physicians estimating the impact of chromosomal rearrangements on the patient phenotype. A separate section describes available tools that can be used to visualize and analyze human genome architecture.

Hedgehog Pathway Drives Fusion-Negative Rhabdomyosarcoma Initiated From Non-myogenic Endothelial Progenitors

Rhabdomyosarcoma (RMS) is a pediatric soft tissue sarcoma that histologically resembles embryonic skeletal muscle. RMS occurs throughout the body and an exclusively myogenic origin does not account for RMS occurring in sites devoid of skeletal muscle. We previously described an RMS model activating a conditional constitutively active Smoothened mutant (SmoM2) with aP2-Cre. Using genetic fate mapping, we show SmoM2 expression in Cre-expressing endothelial progenitors results in myogenic transdifferentiation and RMS. We show that endothelium and skeletal muscle within the head and neck arise from Kdr-expressing progenitors, and that hedgehog pathway activation results in aberrant expression of myogenic specification factors as a potential mechanism driving RMS genesis. These findings suggest that RMS can originate from aberrant development of non-myogenic cells.