The proteomic landscape of soft tissue sarcomas

Figures

• 

  Fig. 1. Schematic overview of the study. a, b Pie charts showing the count and percentage breakdown of histological subtypes (a) and anatomical sites (b) within the study cohort. c Overview of the proteomic analysis workflow. AS angiosarcoma, ASPS alveolar soft part sarcoma, CCS clear cell sarcoma, DDLPS dedifferentiated liposarcoma, DES desmoid tumour, DSRCT desmoplastic small round cell tumour, EPS epithelioid sarcoma, LMS leiomyosarcoma, RT rhabdoid tumour, SS synovial sarcoma, UPS undifferentiated pleomorphic sarcoma, FFPE formalin-fixed paraffin-embedded, LC liquid chromatography, MS mass spectrometry.

• 

  Fig. 2. The proteomic landscape of soft tissue sarcoma. a Annotated heatmap showing the unsupervised clustering (Pearson’s distance) of 3290 proteins across the study cohort. From top to bottom, panels indicate histological subtype, anatomical site, tumour grade, patient sex, patient age, and tumour size. b Uniform manifold approximation and projection (UMAP) of the proteomic data with individual cases coloured by histological subtype. c Heatmap showing the proteins (n = 1362) uniquely upregulated in histological subtypes with greater than 20 cases in the cohort (FDR < 1%, fold change ≥1.5), sorted by histology. Annotations indicate key proteins (DDLPS & SS) identified by significant analysis of microarray (SAM) and gene sets (AS, DES, LMS, UPS) identified by overrepresentation analysis in each histological subtype (Supplementary Data 4). AS angiosarcoma, ASPS alveolar soft part sarcoma, CCS clear cell sarcoma, DDLPS dedifferentiated liposarcoma, DES desmoid tumour, DSRCT desmoplastic small round cell tumour, EPS epithelioid sarcoma, LMS leiomyosarcoma, RT rhabdoid tumour, SS synovial sarcoma, UPS undifferentiated pleomorphic sarcoma.

• 

  Fig. 3. Leiomyosarcoma (LMS) is comprised of three proteomic subtypes. a Annotated heatmap showing the unsupervised clustering (Spearman distance) of 3262 proteins across LMS cases (n = 80), arranged by proteomic subtype (top annotation). Bottom annotations indicate key tumour and patient characteristics, and significant (one-way ANOVA; FDR < 0.001) biological features obtained from single sample Gene Set Enrichment Analysis (ssGSEA) of the MSigDB Hallmark gene sets (Supplemental Data 6A, B). b Pie charts depicting the breakdown of LMS proteomic subtypes at different anatomical sites. c Kaplan–Meier plot of local recurrence-free survival (LRFS) across the LMS proteomic subtypes stratified by P3 and P1/P2 combined. Hazard ratio (HR), 95% confidence intervals (CI) and p-value determined by univariable Cox regression. d Multivariable Cox regression assessing local recurrence-free survival (LRFS) in patients categorised by leiomyosarcoma (LMS) proteomic subtype. I-A intra-abdominal, RP retroperitoneal.

• 

  Fig. 4. Characterisation of the immune profiles of dedifferentiated liposarcoma (DDLPS) and undifferentiated pleomorphic sarcoma (UPS). a Representative images of high and low CD3+ tumour infiltrating lymphocyte (TIL) staining from an exemplar in DDLPS (green) and UPS (purple) cases in the cohort. Samples were stratified as high and low based on median TIL counts (107 cells/mm2). b Kaplan–Meier plot of overall survival (OS) in CD3+ TIL-high and -low patients (n = 82). Hazard ratio (HR), 95% confidence intervals (CI) and p-value determined by univariable Cox regression. c Boxplots comparing expression of 21 immune-related genes in CD3+ TIL-high and -low cases. Boxplots indicate 25th, 50th, and 75th percentile, with whiskers extending from 25th percentile−(1.5*IQR) to 75th percentile+(1.5*IQR), and outliers plotted as points. p values determined by Kruskal–Wallis tests and adjusted to false discovery rate (FDR). d Gene set enrichment analysis (GSEA) results showing the top 15 gene sets enriched in CD3+ TIL-high and and-low cases based on normalised enrichment score (NES) with gene sets related to complement activity (blue) and coagulation processes (orange) highlighted. e To inspect the…

• 

  Fig. 5. Comparative analysis of transcriptomic and proteomic profiles of angiosarcomas (AS). a Volcano plot showing Spearman’s correlation and -log10 transformed p-values for the 3383 genes/proteins. Negatively correlated genes/proteins with p-value < 0.005 and positively correlated genes/proteins with FDR < 0.001 are annotated on the plot. b Annotated heatmap of proteomic data (3383 proteins) for 25 AS cases. The samples were clustered using M3C method with K-means. From top to bottom, panels indicate age, sex, size, performance status, tumour grade, depth, margin, size and aetiology/subtype. The corresponding RNA-seq clusters (Fig. S7A) are shown. c Scatter plots of log2-transformed hazard ratios from univariate Cox regression models fitted using OS (left panel), LRFS (middle panel) and MFS (right panel) using gene/protein expression. Blue dots are proteins with a p-value < 0.05, green dots are genes with a p-value < 0.05 and red dots are the gene/proteins where both datasets returned a p-value < 0.05. d Venn diagram showing the overlap of the genes and proteins that are significantly associated with all the survival endpoint measures (OS, LRFS, and MFS). e Likelihood ratios…

• 

  Fig. 6. Sarcoma proteomic modules (SPM) are associated with patient survival outcomes. a Co-expression heatmap showing the correlation of protein expression based on topological overlap matrix (TOM) dissimilarity (1−TOM)7. Cluster dendrogram height indicates 1−Pearson’s correlation. b Protein co-expression network comprising 3290 nodes and 168,574 edges. Nodes indicate proteins and are coloured based on SPM membership. Edges show a correlation between protein expression, where a thicker line indicates a stronger correlation. Representative biological features are annotated for each module. c Overview of univariable Cox regression results for each SPM and local recurrence-free survival (LRFS), metastasis-free survival (MFS), and overall survival (OS). d Protein–protein interaction (PPI) network of SPM 10 comprising 94 nodes and 233 edges. Nodes are proteins and edges represent the StringDB database score between proteins, where a thicker line indicates a higher score (range = 0.401– 0.999). e Sankey diagram illustrating the distribution of histological subtype (excluding DES and RT) across three SPM10 subgroups. Subgroups identified by tertile stratification based on median SPM 10…