Pan-cancer subclonal mutation analysis of 7,827 tumors predicts clinical outcome.

Intra-tumor heterogeneity is an important driver of tumor evolution and therapy response. Advances in precision cancer treatment will require understanding of mutation clonality and subclonal architecture. Currently the slow computational speed of subclonal reconstruction hinders large cohort studies. To overcome this bottleneck, we developed Clonal structure identification through Pairwise Penalization, or CliPP, which clusters subclonal mutations using a regularized likelihood model. CliPP reliably processed whole-genome and whole-exome sequencing data from over 12,000 tumor samples within 24 hours, thus enabling large-scale downstream association analyses between subclonal structures and clinical outcomes. Through a pan-cancer investigation of 7,827 tumors from 32 cancer types, we found that high subclonal mutational load (sML), a measure of latency time in tumor evolution, was significantly associated with better patient outcomes in 16 cancer types with low to moderate tumor mutation burden (TMB). In a cohort of prostate cancer patients participating in an immunotherapy clinical trial, high sML was indicative of favorable response to immune checkpoint blockade. This comprehensive study using CliPP underscores sML as a key feature of cancer. sML may be essential for linking mutation dynamics with immunotherapy response in the large population of non-high TMB cancers.

The history of chromosomal instability in genome doubled tumors.

Tumors frequently display high chromosomal instability and contain multiple copies of genomic regions. Here, we describe GRITIC, a generic method for timing genomic gains leading to complex copy number states, using single-sample bulk whole-genome sequencing data. By applying GRITIC to 6,091 tumors, we found that non-parsimonious evolution is frequent in the formation of complex copy number states in genome-doubled tumors. We measured chromosomal instability before and after genome duplication in human tumors and found that late genome doubling was followed by an increase in the rate of copy number gain. Copy number gains often accumulate as punctuated bursts, commonly after genome doubling. We infer that genome duplications typically affect the landscape of copy number losses, while only minimally impacting copy number gains. In summary, GRITIC is a novel copy number gain timing framework that permits the analysis of copy number evolution in chromosomally unstable tumors.

Age-specific induction of mutant p53 drives clonal hematopoiesis and acute myeloid leukemia in adult mice.

The investigation of the mechanisms behind p53 mutations in acute myeloid leukemia (AML) has been limited by the lack of suitable mouse models, which historically have resulted in lymphoma rather than leukemia. This study introduces two new AML mouse models. One model induces mutant p53 and Mdm2 haploinsufficiency in early development, showing the role of Mdm2 in myeloid-biased hematopoiesis and AML predisposition, independent of p53. The second model mimics clonal hematopoiesis by inducing mutant p53 in adult hematopoietic stem cells, demonstrating that the timing of p53 mutation determines AML vs. lymphoma development. In this context, age-related changes in hematopoietic stem cells (HSCs) collaborate with mutant p53 to predispose toward myeloid transformation rather than lymphoma development. Our study unveils new insights into the cooperative impact of HSC age, Trp53 mutations, and Mdm2 haploinsufficiency on clonal hematopoiesis and the development of myeloid malignancies.

Genomic profile analysis of leiomyomas with bizarre nuclei and fumarate hydratase deficient leiomyomas: Strengths, weaknesses, and limitations of array-CGH interpretation.

A close relationship has been demonstrated between genomic complexity and clinical outcome in uterine smooth muscle tumors. We studied the genomic profiles by array-CGH of 28 fumarate hydratase deficient leiomyomas and 37 leiomyomas with bizarre nuclei (LMBN) from 64 patients. Follow-up was available for 46 patients (from three to 249 months, mean 87.3 months). All patients were alive without evidence of disease. For 51 array-CGH interpretable tumors the mean Genomic Index (GI) was 16.4 (median: 9.8; from 1 to 57.8), significantly lower than the mean GI in LMS (mean GI 51.8, p < 0.001). We described three groups: (1) a group with FH deletion (24/58) with low GI (mean GI: 11 vs. 22,4, p = 0.02), (2) a group with TP53 deletion (17/58) with higher GI (22.4 vs. 11 p = 0.02), and (3) a group without genomic events on FH or TP53 genes (17/58) (mean GI:18.3; from 1 to 57.8). Because none of these tumors recurred and none showed morphological features of LMS we concluded that GI at the cut-off of 10 was not applicable in these subtypes of LM. By integration of all those findings, a GI <10 in LMBN remains a valuable argument for benignity. Conversely, in LMBN a GI >10 or alteration in tumor suppressor genes, should not alone warrant a diagnosis of malignancy. Nine tumors were tested with Nanocind CINSARC® signature and all were classified in low risk of recurrence. We propose, based on our observations, a diagnostic approach of these challenging lesions.

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring.

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

Spatiotemporal genomic profiling of intestinal metaplasia reveals clonal dynamics of gastric cancer progression.

Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. Analyzing 1,256 gastric samples (1,152 IMs) across 692 subjects from a prospective 10-year study, we identify 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9). Single-cell and spatial profiles highlight changes in tissue ecology and IM lineage heterogeneity, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling reveals expression-based molecular subtypes of IM associated with incomplete histology, antral/intestinal cell types, ARID1A mutations, inflammation, and microbial communities normally associated with the healthy oral tract. We demonstrate that combined clinical-genomic models outperform clinical-only models in predicting IMs likely to transform to GC. By highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression, our results raise opportunities for GC precision prevention and interception.

Cancer origin tracing and timing in two high-risk prostate cancers using multisample whole genome analysis: prospects for personalized medicine.

BACKGROUND: Prostate cancer (PrCa) genomic heterogeneity causes resistance to therapies such as androgen deprivation. Such heterogeneity can be deciphered in the context of evolutionary principles, but current clinical trials do not include evolution as an essential feature. Whether or not analysis of genomic data in an evolutionary context in primary prostate cancer can provide unique added value in the research and clinical domains remains an open question. METHODS: We used novel processing techniques to obtain whole genome data together with 3D anatomic and histomorphologic analysis in two men (GP5 and GP12) with high-risk PrCa undergoing radical prostatectomy. A total of 22 whole genome-sequenced sites (16 primary cancer foci and 6 lymph node metastatic) were analyzed using evolutionary reconstruction tools and spatio-evolutionary models. Probability models were used to trace spatial and chronological origins of the primary tumor and metastases, chart their genetic drivers, and distinguish metastatic and non-metastatic subclones. RESULTS: In patient GP5, CDK12 inactivation was among the first mutations, leading to a PrCa tandem duplicator phenotype and initiating the cancer around age 50, followed by rapid cancer evolution after age 57, and metastasis around age 59, 5 years prior to prostatectomy. In patient GP12, accelerated cancer progression was detected after age 54, and metastasis occurred around age 56, 3 years prior to prostatectomy. Multiple metastasis-originating events were identified in each patient and tracked anatomically. Metastasis from prostate to lymph nodes occurred strictly ipsilaterally in all 12 detected events. In this pilot, metastatic subclone content analysis appears to substantially enhance the identification of key drivers. Evolutionary analysis' potential impact on therapy selection appears positive in these pilot cases. CONCLUSIONS: PrCa evolutionary analysis allows tracking of anatomic site of origin, timing of cancer origin and spread, and distinction of metastatic-capable from non-metastatic subclones. This enables better identification of actionable targets for therapy. If extended to larger cohorts, it appears likely that similar analyses could add substantial biological insight and clinically relevant value.

Analysis of the Peritumoral Tissue Unveils Cellular Changes Associated with a High Risk of Recurrence.

BACKGROUND: The management of soft-tissue sarcoma (STS) relies on a multidisciplinary approach involving specialized oncological surgery combined with other adjuvant therapies to achieve optimal local disease control. Purpose and Results: Genomic and transcriptomic pseudocapsules of 20 prospective sarcomas were analyzed and revealed to be correlated with a higher risk of recurrence after surgery. CONCLUSIONS: A peritumoral environment that has been remodeled and infiltrated by M2 macrophages, and is less expressive of healthy tissue, would pose a significant risk of relapse and require more aggressive treatment strategies.

CINSARC in high-risk soft tissue sarcoma patients treated with neoadjuvant chemotherapy: Results from the ISG-STS 1001 study.

BACKGROUND: The Complexity INdex in SARComas (CINSARC) is a transcriptional signature derived from the expression of 67 genes involved in mitosis control and chromosome integrity. This study aims to assess CINSARC value of in an independent series of high-risk patients with localized soft tissue sarcoma (STS) treated with preoperative chemotherapy within a prospective, randomized, phase III study (ISG-STS 1001). PATIENTS AND METHODS: Patients with available pre-treatment samples, treated with 3 cycles of either standard (ST) preoperative or histotype-tailored (HT) chemotherapy, were scored according to CINSARC (low-risk, C1; high-risk, C2). The 10-year overall survival probability (pr-OS) according to SARCULATOR was calculated, and patients were classified accordingly (low-risk, Sarc-LR, 10-year pr-OS>60%; high-risk, Sarc-HR, 10-year pr-OS<60%). Survival functions were estimated using the Kaplan-Meier method and compared using log-rank test. RESULTS: Eighty-six patients were included, 30 C1 and 56 C2, 49 Sarc-LR and 37 Sarc-HR. A low level of agreement between CINSARC and SARCULATOR was observed (Cohen's Kappa = 0.174). The 5-year relapse-free survival in C1 and C2 were 0.57 and 0.55 (p = 0.481); 5-year metastases-free survival 0.63 and 0.64 (p = 0.740); 5-year OS 0.80 and 0.72 (p = 0.460). The 5-year OS in C1 treated with ST and HT chemotherapy was 0.84 and 0.76 (p = 0.251) respectively; in C2 treated it was 0.72 and 0.70 (p = 0.349). The 5-year OS in Sarc-LR treated with S and HT chemotherapy was 0.80 and 0.82 (p = 0.502) respectively; in Sarc-HR it was 0.70 and 0.61 (p = 0.233). CONCLUSIONS: Our results, although constrained by the small size of the series, suggest that CINSARC has weak prognostic power in high-risk, localized STS treated with neoadjuvant chemotherapy.

A genetic model for central chondrosarcoma evolution correlates with patient outcome.

BACKGROUND: Central conventional chondrosarcoma (CS) is the most common subtype of primary malignant bone tumour in adults. Treatment options are usually limited to surgery, and prognosis is challenging. These tumours are characterised by the presence and absence of IDH1 and IDH2 mutations, and recently, TERT promoter alterations have been reported in around 20% of cases. The effect of these mutations on clinical outcome remains unclear. The purpose of this study was to determine if prognostic accuracy can be improved by the addition of genomic data, and specifically by examination of IDH1, IDH2, and TERT mutations. METHODS: In this study, we combined both archival samples and data sourced from the Genomics England 100,000 Genomes Project (n = 356). Mutations in IDH1, IDH2, and TERT were profiled using digital droplet PCR (n = 346), whole genome sequencing (n=68), or both (n = 64). Complex events and other genetic features were also examined, along with methylation array data (n = 84). We correlated clinical features and patient outcomes with our genetic findings. RESULTS: IDH2-mutant tumours occur in older patients and commonly present with high-grade or dedifferentiated disease. Notably, TERT mutations occur most frequently in IDH2-mutant tumours, although have no effect on survival in this group. In contrast, TERT mutations are rarer in IDH1-mutant tumours, yet they are associated with a less favourable outcome in this group. We also found that methylation profiles distinguish IDH1- from IDH2-mutant tumours. IDH wild-type tumours rarely exhibit TERT mutations and tend to be diagnosed in a younger population than those with tumours harbouring IDH1 and IDH2 mutations. A major genetic feature of this group is haploidisation and subsequent genome doubling. These tumours evolve less frequently to dedifferentiated disease and therefore constitute a lower risk group. CONCLUSIONS: Tumours with IDH1 or IDH2 mutations or those that are IDHwt have significantly different genetic pathways and outcomes in relation to TERT mutation. Diagnostic testing for IDH1, IDH2, and TERT mutations could therefore help to guide clinical monitoring and prognostication.

Signatures of copy number alterations in human cancer.

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

A pan-cancer compendium of chromosomal instability.

Chromosomal instability (CIN) results in the accumulation of large-scale losses, gains and rearrangements of DNA1. The broad genomic complexity caused by CIN is a hallmark of cancer2; however, there is no systematic framework to measure different types of CIN and their effect on clinical phenotypes pan-cancer. Here we evaluate the extent, diversity and origin of CIN across 7,880 tumours representing 33 cancer types. We present a compendium of 17 copy number signatures that characterize specific types of CIN, with putative aetiologies supported by multiple independent data sources. The signatures predict drug response and identify new drug targets. Our framework refines the understanding of impaired homologous recombination, which is one of the most therapeutically targetable types of CIN. Our results illuminate a fundamental structure underlying genomic complexity in human cancers and provide a resource to guide future CIN research.

E3 ubiquitin ligase HECTD2 mediates melanoma progression and immune evasion.

The ubiquitin-proteasome system maintains protein homoeostasis, underpins the cell cycle, and is dysregulated in cancer. However, the role of individual E3 ubiquitin ligases, which mediate the final step in ubiquitin-mediated proteolysis, remains incompletely understood. Identified through screening for cancer-specific endogenous retroviral transcripts, we show that the little-studied E3 ubiquitin ligase HECTD2 exerts dominant control of tumour progression in melanoma. HECTD2 cell autonomously drives the proliferation of human and murine melanoma cells by accelerating the cell cycle. HECTD2 additionally regulates cancer cell production of immune mediators, initiating multiple immune suppressive pathways, which include the cyclooxygenase 2 (COX2) pathway. Accordingly, higher HECTD2 expression is associated with weaker anti-tumour immunity and unfavourable outcome of PD-1 blockade in human melanoma and counteracts immunity against a model tumour antigen in murine melanoma. This central, multifaceted role of HECTD2 in cancer cell-autonomous proliferation and in immune evasion may provide a single target for a multipronged therapy of melanoma.

ATRX Alteration Contributes to Tumor Growth and Immune Escape in Pleomorphic Sarcomas.

Whole genome and transcriptome sequencing of a cohort of 67 leiomyosarcomas has been revealed ATRX to be one of the most frequently mutated genes in leiomyosarcomas after TP53 and RB1. While its function is well described in the alternative lengthening of telomeres mechanism, we wondered whether its alteration could have complementary effects on sarcoma oncogenesis. ATRX alteration is associated with the down-expression of genes linked to differentiation in leiomyosarcomas, and to immunity in an additional cohort of 60 poorly differentiated pleomorphic sarcomas. In vitro and in vivo models showed that ATRX down-expression increases tumor growth rate and immune escape by decreasing the immunity load of active mast cells in sarcoma tumors. These data indicate that an alternative to unsuccessful targeting of the adaptive immune system in sarcoma could target the innate system. This might lead to a better outcome for sarcoma patients in terms of ATRX status.

An integrative analysis of the age-associated multi-omic landscape across cancers.

Age is the most important risk factor for cancer, as cancer incidence and mortality increase with age. However, how molecular alterations in tumours differ among patients of different age remains largely unexplored. Here, using data from The Cancer Genome Atlas, we comprehensively characterise genomic, transcriptomic and epigenetic alterations in relation to patients' age across cancer types. We show that tumours from older patients present an overall increase in genomic instability, somatic copy-number alterations (SCNAs) and somatic mutations. Age-associated SCNAs and mutations are identified in several cancer-driver genes across different cancer types. The largest age-related genomic differences are found in gliomas and endometrial cancer. We identify age-related global transcriptomic changes and demonstrate that these genes are in part regulated by age-associated DNA methylation changes. This study provides a comprehensive, multi-omics view of age-associated alterations in cancer and underscores age as an important factor to consider in cancer research and clinical practice.

Cell-cell fusion of mesenchymal cells with distinct differentiations triggers genomic and transcriptomic remodelling toward tumour aggressiveness.

Cell-cell fusion is a physiological process that is hijacked during oncogenesis and promotes tumour evolution. The main known impact of cell fusion is to promote the formation of metastatic hybrid cells following fusion between mobile leucocytes and proliferating tumour cells. We show here that cell fusion between immortalized myoblasts and transformed fibroblasts, through genomic instability and expression of a specific transcriptomic profile, leads to emergence of hybrid cells acquiring dissemination properties. This is associated with acquisition of clonogenic ability by fused cells. In addition, by inheriting parental properties, hybrid tumours were found to mimic the histological characteristics of a specific histotype of sarcomas: undifferentiated pleomorphic sarcomas with incomplete muscular differentiation. This finding suggests that cell fusion, as macroevolution event, favours specific sarcoma development according to the differentiation lineage of parent cells.

A Global and Integrated Analysis of CINSARC-Associated Genetic Defects.

The Complexity Index in Sarcomas (CINSARC) signature is a transcriptomic marker that identifies high-risk soft-tissue sarcomas and is associated with high metastatic potential. During the last decade, CINSARC has been successfully developed and validated and is currently being assessed in two prospective phase III clinical trials for stratification of therapy. Although the link between CINSARC expression and tumor aggressiveness is well established, questions remain about how CINSARC genes are regulated. In this study, we leveraged a The Cancer Genome Atlas multiomics study on sarcomas with complex genetics to appraise the association between CINSARC profile, genomic features, and two potential regulation mechanisms, DNA methylation and miRNA expression. CINSARC expression was associated with an increase of ploidy, intratumor heterogeneity, copy-number alteration, altered expression of 37 miRNAs, and a decrease of DNA methylation. These genetic changes are not independent, but rather act together to promote or repress CINSARC expression. These findings depict new insights into CINSARC regulation. SIGNIFICANCE: These findings demonstrate that CINSARC is associated with a variety of genomic aberrations that contribute to higher risk for metastasis and may serve as a prognostic factor in sarcomas and beyond.

Value of peri-operative chemotherapy in patients with CINSARC high-risk localized grade 1 or 2 soft tissue sarcoma: study protocol of the target selection phase III CHIC-STS trial.

BACKGROUND: The value of chemotherapy in soft tissue sarcoma (STS) remains controversial. Several expert teams consider that chemotherapy provides a survival advantage and should be proposed in high-risk (HR) patients. However, the lack of accuracy in identifying HR patients with conventional risk factors (large, deep, FNCLCC grade 3, extremity STS) is an issue that cannot be neglected. For example, while the FNCLCC grading system is a powerful tool, it has several limitations. CINSARC, a 67-gene signature, has proved to be an additional independent factor for predicting metastatic spread and outperforms histological grade. Regardless of FNCLCC grade, CINSARC stratifies patients into two separate prognostic groups: one with an excellent prognosis (low-risk (LR) CINSARC) and the other with a worse outcome (HR-CINSARC) in terms of metastatic relapse. Here we evaluate the role of chemotherapy in grade 1-2 STS patients with HR-CINSARC and assess the prognostic value of CINSARC in patients treated with standard of care. METHODS: CHIC is a parallel, randomized, open-label, multicenter study evaluating the effect on metastasis-free survival of adding perioperative chemotherapy to standard of care in patients with grade ½ STS sarcoma defined as HR by CINSARC. In this target selection design, 600 patients will be screened with CINSARC to randomize 250 HR-CINSARC patients between standard of care and standard of care plus chemotherapy (4 cycles of 3 weeks of intravenous chemotherapy with doxorubicin in combination with dacarbazine or ifosfamide according to histologic subtype). LR-CINSARC patients will be treated by standard of care according to the investigator. The primary endpoint is metastasis-free survival. Secondary endpoints include overall survival, disease-free survival and safety. Furthermore, the prognostic value of CINSARC will be evaluated by comparing LR-CINSARC patients to HR-CINSARC patients randomized in standard of care. DISCUSSION: CHIC is a prospective randomized phase III trial designed to comprehensively evaluate the benefit of chemotherapy in HR-CINSARC patients and to prospectively validate the prognostic value of CINSARC in grade ½ STS sarcoma patients. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT04307277 Date of registration: 13 March 2020.

LIX1 regulates YAP activity and controls gastrointestinal cancer cell plasticity.

Gastrointestinal stromal tumours (GISTs), the most common mesenchymal neoplasm of the gastrointestinal tract, result from deregulated proliferation of transformed KIT-positive interstitial cells of Cajal that share mesenchymal progenitors with smooth muscle cells. Despite the identification of selective KIT inhibitors, primary resistance and relapse remain a major concern. Moreover, most patients develop resistance partly through reactivation of KIT and its downstream signalling pathways. We previously identified the Limb Expression 1 (LIX1) gene as a unique marker of digestive mesenchyme immaturity. We also demonstrated that LIX1 regulates mesenchymal progenitor proliferation and differentiation by controlling the Hippo effector YAP1, which is constitutively activated in many sarcomas. Therefore, we wanted to determine LIX1 role in GIST development. We found that LIX1 is strongly up-regulated in GIST samples and this is associated with unfavourable prognosis. Moreover, LIX1 controls GIST cell proliferation in vitro and in vivo. Upon LIX1 inactivation in GIST cells, YAP1/TAZ activity is reduced, KIT (the GIST signature) is down-regulated, and cells acquire smooth muscle lineage features. Our data highlight LIX1 role in digestive mesenchyme-derived cell-fate decisions and identify this novel regulator as a target for drug design for GIST treatment by influencing its differentiation status.

Tetraploidization of Immortalized Myoblasts Induced by Cell Fusion Drives Myogenic Sarcoma Development with DMD Deletion.

Whole-genome doubling is the second most frequent genomic event, after TP53 alterations, in advanced solid tumors and is associated with poor prognosis. Tetraploidization step will lead to aneuploidy and chromosomic rearrangements. The mechanism leading to tetraploid cells is important since endoreplication, abortive cytokinesis and cell fusion could have distinct consequences. Unlike processes based on duplication, cell fusion involves the merging of two different genomes, epigenomes and cellular states. Since it is involved in muscle differentiation, we hypothesized that it could play a role in the oncogenesis of myogenic cancers. Spontaneous hybrids, but not their non-fused immortalized myoblast counterparts they are generated from, induced tumors in mice. Unstable upon fusion, the hybrid genome evolved from initial mitosis to tumors with a highly rearranged genome. This genome remodeling finally produced targeted DMD deletions associated with replicative stress, isoform relocalization and metastatic spreading, exactly as observed in human myogenic sarcomas. In conclusion, these results draw a model of myogenic oncogenesis in which cell fusion and oncogene activation combine to produce pleomorphic aggressive sarcomas.