Dissecting the Origin of Heterogeneity in Uterine and Ovarian Carcinosarcomas

Abstract

Gynecologic carcinosarcomas (CS) are biphasic neoplasms composed of carcinomatous (C) and sarcomatous (S) malignant components. Because of their rarity and histologic complexity, genetic and functional studies on CS are scarce and the mechanisms of initiation and development remain largely unknown. Whole-genome analysis of the C and S components reveals shared genomic alterations, thus emphasizing the clonal evolution of CS. Reconstructions of the evolutionary history of each tumor further reveal that C and S samples are composed of both ancestral cell populations and component-specific subclones, supporting a common origin followed by distinct evolutionary trajectories. However, while we do not find any recurrent genomic features associated with phenotypic divergence, transcriptomic and methylome analyses identify a common mechanism across the cohort, the epithelial-to-mesenchymal transition (EMT), suggesting a role for nongenetic factors in inflicting changes to cellular fate. Altogether, these data accredit the hypothesis that CS tumors are driven by both clonal evolution and transcriptomic reprogramming, essential for susceptibility to transdifferentiation upon encountering environmental cues, thus linking CS heterogeneity to genetic, transcriptomic, and epigenetic influences.

Significance: We have provided a detailed characterization of the genomic landscape of CS and identified EMT as a common mechanism associated with phenotypic divergence, linking CS heterogeneity to genetic, transcriptomic, and epigenetic influences.

FIGURE 1

Genomic landscape of CS. A, Tissue origin and histological subtypes for each tumor (from P01 to P15). B, Main component and detailed percentage content for both samples (a and b) derived from the same tumor. C, CNA analysis: tumor purity (%), ploidy, FGA, and CNA breakpoints number. D, SV number, classified as deletion, duplication, inversion, and interchromosomal translocation. E, SNV and small indel number. F, Proportion of each of four mutational signatures deciphered in the whole cohort. G, MSI phenotype: MSI score (representing the percentage of unstable microsatellite loci), MLH1 promoter methylation (beta-value), and mRNA expression level (TPM – transcripts per million). H, Kataegis number illustrating APOBEC error-prone DNA repair process. I, BRCAness phenotype: number of LST, (white stars highlight samples with LST number ≥20), BRCA1 promoter methylation level (beta-value), BRCA1/2 genes alterations and TDP score, collectively summarized by the HRD status track.

FIGURE 2

Key gene alterations in uterine and ovarian CS. Oncoprint of alterations identified in the 50 most altered genes from TCGA cancer pathways and custom lists of CRG. The type of genomic alteration (deletion, amplification, fusion, broken by SV, missense, truncated) is described in the legend. Number and proportions of alteration types are summarized by gene (horizontal bars on right) and by sample (vertical bars on top). The samples are classified into MSI and CN-high according to TCGA endometrial carcinoma classification.

FIGURE 3

Differential expression and methylation analysis (C vs. S). A, Heatmap and pathway overrepresentation analysis of differentially expressed genes between C and S (or undifferentiated) components of the 10 samples derived from the five selected CS tumors. Gene clustering method: Ward's; distance: Spearman. FDR: false discovery rate. B, Heatmaps of EMT and stemness pathways scores (left) and methylation level of EMT-related miRNAs (right). mRNA level: scaled TPM. Pathway scores: scaled ssGSEA scores. Methylation level: scaled beta-values.

FIGURE 4

Evolutionary histories of CS. Clonal lineage histories of P02 MSI phenotype tumor (A), P12 BRCAness-phenotype tumor (B), and P15 ovarian bilateral tumor (C). Left, Clonal lineage inference for the whole tumor: subclonal populations (Ci), main cancer gene alterations, genome doubling events, and representative mutational signature proportions are represented. Right, Subclonal population frequencies projection in each tumor sample. Clones colored in gray are found in both derived samples, pink (or green) clones are specific to sarcoma- (or carcinoma-) derived samples, respectively.