
==== Front
Commun Med (Lond)
Commun Med (Lond)
Communications Medicine
2730-664X
Nature Publishing Group UK London

38762630
522
10.1038/s43856-024-00522-x
Article
Proteomic features of soft tissue tumours in adolescents and young adults
http://orcid.org/0000-0002-1319-5148
Tam Yuen Bun 1
http://orcid.org/0000-0003-4275-660X
Low Kaan 1
PS Hari 1
Chadha Madhumeeta 1
Burns Jessica 1
Wilding Christopher P. 1
Arthur Amani 1
Chen Tom W. 23
Thway Khin 14
Sadanandam Anguraj 1
Jones Robin L. 45
http://orcid.org/0000-0003-3972-5087
Huang Paul H. paul.huang@icr.ac.uk

1
1 https://ror.org/043jzw605 grid.18886.3f 0000 0001 1499 0189 Division of Molecular Pathology, The Institute of Cancer Research, London, United Kingdom
2 https://ror.org/03nteze27 grid.412094.a 0000 0004 0572 7815 Department of Oncology, National Taiwan University Hospital, Taipei, Taiwan
3 https://ror.org/05bqach95 grid.19188.39 0000 0004 0546 0241 Graduate Institute of Oncology, National Taiwan University College of Medicine, Taipei, Taiwan
4 https://ror.org/0008wzh48 grid.5072.0 0000 0001 0304 893X The Royal Marsden NHS Foundation Trust, London, United Kingdom
5 https://ror.org/043jzw605 grid.18886.3f 0000 0001 1499 0189 Division of Clinical Studies, The Institute of Cancer Research, London, United Kingdom
18 5 2024
18 5 2024
2024
4 9313 7 2023
7 5 2024
© The Author(s) 2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Background

Adolescents and young adult (AYA) patients with soft tissue tumours including sarcomas are an underserved group with disparities in treatment outcomes.

Methods

To define the molecular features between AYA and older adult (OA) patients, we analysed the proteomic profiles of a large cohort of soft tissue tumours across 10 histological subtypes (AYA n = 66, OA n = 243), and also analysed publicly available functional genomic data from soft tissue tumour cell lines (AYA n = 5, OA n = 8).

Results

Biological hallmarks analysis demonstrates that OA tumours are significantly enriched in MYC targets compared to AYA tumours. By comparing the patient-level proteomic data with functional genomic profiles from sarcoma cell lines, we show that the mRNA splicing pathway is an intrinsic vulnerability in cell lines from OA patients and that components of the spliceosome complex are independent prognostic factors for metastasis free survival in AYA patients.

Conclusions

Our study highlights the importance of performing age-specific molecular profiling studies to identify risk stratification tools and targeted agents tailored for the clinical management of AYA patients.

Plain language summary

Soft tissue tumours are cancers that develop in the connective and supporting tissues of the body, such as muscle or fat. These tumours arise in patients across the entire age range. However, improvements in survival outcomes in adolescent and young adult (AYA) patients have lagged behind outcomes in older adults (OA) and children. To better understand the biology of AYA patients with soft tissue tumours, we analysed protein profiles across 10 different types. We identified biological differences between AYA and OA patients and report an age-specific signature that can potentially be used to help predict which AYA patients are more likely to have aggressive cancers that will spread to other parts of the body. Our study highlights the importance of performing age-specific studies to identify new tools to predict patient outcomes and potentially find more suitable treatments.

Tam et al. compare molecular features between adolescents and young adult (AYA) and older adult (OA) patients with soft tissue tumours and identify that components of the spliceosome complex are independent prognostic factors for metastasis free survival in AYA patients.

Subject terms

Sarcoma
Paediatric cancer
https://doi.org/10.13039/501100000289 Cancer Research UK (CRUK) C56167/A29363 Huang Paul H. https://doi.org/10.13039/501100000272 DH | National Institute for Health Research (NIHR) https://doi.org/10.13039/100009756 Sarcoma UK SUK03.2019 Jones Robin L. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Soft tissue tumours are rare mesenchymal tumours that span >80 histological subtypes of distinct biology and genetics1. These include malignant cancers such as sarcomas as well as tumours with no metastatic potential such as desmoid tumours (DES). Sarcomas have a higher incidence in cancers diagnosed in the adolescent and young adult (AYA) age group (16–39 years at the time of diagnosis, 8% of all cancer diagnosis), compared to older adults (OA) (>39 years, 1% of all cancer diagnosis)2,3. Despite the increased incidence, improvements in survival in AYA patients with soft tissue sarcomas (STS) have lagged behind other age groups4. Reasons for this disparity are multi-factorial, and include under-representation in clinical trials4,5, unique psychosocial considerations6,7, inadequate age-specific services8–10, and poor knowledge of their unique biology11,12. A recent study showed that STS histological subtypes typically sensitive to chemotherapy in other age groups are instead chemoresistant in AYA patients13, which suggests that patients in this age group may have distinct biological differences compared to either OA or paediatric patients. In the case of non-rhabdomyosarcoma STS (NRSTS) where the majority of current treatment guidelines relies on drugs that have been optimised in OA patients14,15, the lack of therapies tailored to AYA patients is a major unmet need and a key barrier to improving survival rates.

Several recent pan-cancer analyses have leveraged publicly available datasets (TCGA, ICGC, GENIE) to demonstrate that there are age-associated genomic, transcriptomic and immune microenvironmental differences across multiple cancer types that include STS16–20. For instance, Lee et al., showed that sarcoma patients <50 years old had lower immune-related pathway expression compared to patients >50 years of age at the transcriptomic level17. They further determined that at the genomic level, the older sarcoma patients had higher copy number variation rates. However, the AYA age group is under-represented in all these studies, with only a small number of STS patients and histological subtypes included. Furthermore, these aggregate analyses do not consider well-established differences in the spectrum of STS histological subtypes in AYA versus OA patients21. Prior studies that have undertaken molecular profiling of AYA sarcoma specimens including a recent EORTC SPECTA-AYA study have focused exclusively on genomic and transcriptomic data20,22,23. While informative, these technologies do not provide a direct measure of proteins which are key mediators of tumour cell signalling and the largest class of targets for oncology drugs24–27, making it challenging to bridge the translational gap towards clinical applications. Given the unique tumour, microenvironmental and host differences between AYA and OA patients12, it is likely that AYA patients with STS harbour distinct molecular features which may influence clinical and treatment outcomes, although this has yet to be conclusively demonstrated. Due to the rarity and heterogeneity of STS, to date, there are few studies that have systematically evaluated the molecular differences between AYA and OA patients.

Here, we undertake a detailed analysis of the proteomic features in AYA and OA patients across 10 histological subtypes of soft tissue tumours. By interrogating clinically annotated proteomic profiles in a large cohort of AYA and OA patients and combining it with functional genomics data derived from sarcoma cell lines within the Cancer Cell Line Encyclopaedia (CCLE), we demonstrate that there are significant differences in the biological networks and intrinsic vulnerabilities between these two age groups with implications for biomarker development and therapy selection.

Methods

Patient cohort

The cohort is comprised of 309 patients from two centres (The Royal Marsden Hospital and National Taiwan University Hospital). Patients with a histopathologically confirmed diagnosis of soft tissue sarcoma (STS) or desmoid tumour, and 16 years of age or older at the time of diagnosis were included in the analysis. Soft tissue sarcoma diagnoses included angiosarcoma, alveolar soft part sarcoma, clear cell sarcoma, dedifferentiated liposarcoma, desmoplastic small round cell tumour, epithelioid sarcoma, synovial sarcoma, leiomyosarcoma, and undifferentiated pleomorphic sarcoma. Retrospective collection and analysis of associated clinical data was approved as part of the Royal Marsden Hospital (RMH) PROgnoStic and PrEdiCTive ImmUnoprofiling of Sarcomas (PROSPECTUS) study (NHS Research Ethics Committee Reference 16/EE/0213) or National Taiwan University Hospital (Research Ethics Committee Reference 201912226RINB). Written informed consent was obtained from participants, and all participants were old enough to provide informed consent according to local regulations. Baseline clinicopathological characteristics and survival data were collected by retrospective review of medical records as part of our previous study28.

Proteomic data

Proteomic data for this study was downloaded from ProteomeXchange (PXD036226)28. The SequestHT search engine in Proteome Discoverer 2.2 or 2.3 (Thermo Scientific, Waltham, MA, USA) was used to search the raw mass spectra against reviewed UniProt human protein entries (v2018_07 or later) for protein identification and quantification. Precursor mass tolerance was set at 20 ppm and fragment ion mass tolerance was 0.02 Da. Spectra were searched for fully tryptic peptides with a maximum of two missed cleavages. TMT6plex at N-terminus/lysine and carbamidomethyl at cysteine were selected as fixed modifications. Dynamic modifications were the oxidation of methionine and deamidation of asparagine/glutamine. Peptide confidence was estimated with the Percolator node. The peptide false discovery rate (FDR) was set at 0.01 and validation was based on q value and decoy database search. The reporter ion quantifier node included an integration window tolerance of 15 ppm and an integration method based on the most confident centroid peak at the MS3 level. Only unique peptides were used for quantification, considering protein groups for peptide uniqueness. Peptides with average reporter signal-to-noise >3 were used for protein quantification. Proteins with an FDR <0.01 and a minimum of two peptides were used for downstream analyses.

All data were processed using custom R scripts in R v3.5.1 or later. Proteins identified in <75% of samples were removed, and the remaining missing values were imputed using the k-nearest neighbour (k-NN) algorithm29. To normalise the data and remove batch effects, data for each patient sample was divided by the corresponding reference sample and log2 transformed, followed by median centring across samples and standardising within samples. To visualise the STS proteomic dataset, hierarchical clustering was performed using Pearson correlation distance.

Statistical methods

All statistical tests were two-sided and unless otherwise stated, p values were adjusted to false discovery rate (FDR) using the Benjamini–Hochberg (BH) procedure to account for multiple comparisons where required. Unless otherwise specified, analysis was performed using custom R scripts in R v4.1.1 or later. Two-way analysis of variance (ANOVA) and chi-square tests were implemented, with further details of statistical tests listed in the figure legends.

Differential expression analysis

To identify upregulated proteins in AYA and OA patients, a two-tailed multiple t-test was performed and corrected for multiple comparisons by the BH procedure. Logistic regression analysis was performed to adjust for confounding factors of tumour size, grade, anatomical site, performance status, histological subtype, tumour margin, and tumour depth. Univariate logistic regression first was performed to identify significantly different proteins between AYA and OA (FDR <0.05). Univariate logistic regression was then performed to identify significantly different confounding factors. Each significant protein’s expression was then combined with significant confounding factors and multiple logistic regression was performed with AYA and OA variable.

Single sample GSEA (ssGSEA)

Single sample GSEA (ssGSEA) was performed using the GenePattern online tool (www.genepattern.org, v10.1.0) to score sample-specific enrichment of hallmark gene sets (v2023.1), KEGG Spliceosome complex gene set30 and Sarcoma Proteomic Module (SPM) gene sets28 in the proteomics dataset. ssGSEA score between AYA and OA patients were analysed using a two-way analysis of variance (ANOVA) followed by Šidák correction. Univariate logistic regression was first performed to identify significantly differentially expressed hallmark gene sets between AYA and OA (FDR <0.05). Each significant gene set was then combined with significant confounding factors and multiple logistic regression was performed with AYA and OA variables.

Comparative analysis of CCLE functional genomic and proteomics data

Gene Set Enrichment Analysis (GSEA)30 was performed using the GenePattern online tool (www.genepattern.org, v20.4.0). To identify MSigDB Reactome gene sets (v2023.1) enriched in the two phenotypic classes (in this case, the OA or AYA age groups) within the functional genomic CCLE dataset, only cell lines with a NRSTS subtype and a recorded patient age of >16 years were included for analysis. Cell lines derived from patients between 16–39 years old were grouped as AYA and above 39 years old as OA. The full list of cell lines included is provided in Supplementary Data 1. Genome-scale CRISPR-Cas9 screening data of cell lines was downloaded from the CCLE portal (https://sites.broadinstitute.org/ccle)31. The full list of genes from the CRISPRGeneEffect dataset (DepMap Public 22Q4) was used for GSEA with an FDR <0.05 cut-off for significance. In parallel, GSEA was separately performed on the proteomic dataset to analyse Reactome gene sets that were enriched in soft tissue sarcoma specimens from OA versus AYA patients (FDR <0.05 cut-off). The top enriched gene sets from the proteomic and functional genomic datasets were compared, and shared significant hits were reported.

SPM analysis and protein–protein interaction (PPI) networks

PPI networks were built in Cytoscape v3.9.1 or later32. Previously described SPMs were utilised to identify proteomic signatures enriched in AYA and OA patients28 (Supplementary Data 2). The full SPM network was visualised using protein co-occurrence scores and the group attributes layout using SPM membership. For the full SPM network, a co-occurrence score threshold of >0.05 was applied. To inspect the network of individual spliceosome components, protein networks were constructed using the STRING scores obtained from the STRING database v11.033, with a confidence cut-off score of 0.7 and a grid layout used.

Survival analyses

Patients were split into -high or -low expressing groups for SPM6 or spliceosome components based on the median protein expression level. The association of patient groups with survival outcome were evaluated based on Kaplan–Meier survival estimates and univariable Cox analysis with two-sided Wald test. Multivariable Cox analysis was used to adjust for clinicopathological variables. Three survival outcome endpoints were used. Overall survival (OS) is defined as the time from primary disease surgery to death from any cause. Metastasis free survival (MFS) is defined as the time from primary disease surgery to radiologically confirmed metastatic disease or death. Local recurrence free survival (LRFS) is defined as time from primary disease surgery to radiologically confirmed local recurrence or death. Patients who did not have an event were censored at their last follow-up time, up to 60 months.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Results

Cohort and clinicopathological data

The cohort comprises primary tumour specimens from 309 patients (AYA = 66, OA = 243), for which comprehensive proteomic profiles by mass spectrometry (MS) have previously been generated by our laboratory28. Nine sarcoma subtypes are represented, including alveolar soft part sarcoma (ASPS), angiosarcoma (AS), clear cell sarcoma (CCS), dedifferentiated liposarcoma (DDLPS), desmoplastic small round cell tumour (DSRCT), epithelioid sarcoma (EPS), leiomyosarcoma (LMS), synovial sarcoma (SS) and undifferentiated pleomorphic sarcoma (UPS) (full clinicopathological information provided in Supplementary Table 1). Additionally, DES, a locally infiltrative soft tissue tumour with no metastatic potential and a relatively high incidence in the AYA age group, was included in the cohort34. When broken down by age groups, AYA patients are enriched for ASPS (100%) and DSRCT (75%), while OA patients are enriched for UPS (98%), DDLPS (95%), LMS (91%), AS (90%) and CCS (67%) (Fig. 1a). There are almost equal number of AYA and OA patients in SS (44% AYA), DES (49% AYA) and EPS (50% AYA). The proportion of different histological subtypes in the two age groups is reflective of real-world incidence35. The cohort has a female predominance (male [37%], female [63%]), with a broad distribution of different anatomical sites in each age group (Fig. 1A). Consistent with a previous study of a large cohort of ~5000 STS patients by the Scandinavian Sarcoma Group (SSG)3, our cohort has a higher proportion of patients with grade 3 (OA: 53%, AYA: 15%) and large tumours ≥15 cm (OA: 26%, AYA: 12%) in the OA age group (Supplementary Table 1). Similar to the SSG study, univariate Cox regression analysis showed that patients belonging to the AYA age group were at statistically significant lower risk of death compared to patients in the OA group (Hazard Ratio (HR) = 0.431, 95% Confidence Interval (CI) = 0.232-0.8, p = 0.0077) (Fig. 1b). There was no statistically significant difference in the two age groups for metastasis free survival (MFS) and local relapse-free survival (LRFS) (Supplementary Fig. 1). Note that DES was not included in any of the survival analysis undertaken in this study because they are not malignant.Fig. 1 Overview of the adolescent and young adult (AYA) and older adult (OA) patients in the cohort.

a Distribution of patient sex, anatomical site, and histological subtype within the AYA and OA cohorts. b Kaplan–Meier plot of overall survival (OS) for AYA and OA patients. Hazard ratio (HR), 95% confidence intervals (CI) and p value determined by univariable Cox regression. AS angiosarcoma, ASPS alveolar soft part sarcoma, CCS clear cell sarcoma, DDLPS dedifferentiated liposarcoma, DSRCT desmoplastic small round cell tumour, DE desmoid tumour, EPS epithelioid sarcoma, LMS leiomyosarcoma, SS synovial sarcoma, UPS undifferentiated pleomorphic sarcoma.

Analysis of the AYA and OA proteomic landscape

A total of 8148 proteins were identified, with 3299 proteins quantified across all samples (Fig. 2a and Supplementary Data 3). We defined the proteins that are significantly upregulated in AYA versus OA patients. Following multiple testing correction, 32 and 35 proteins were identified to be significantly upregulated in AYA or OA respectively (FDR <0.05, fold change >2) (Fig. 2b and Supplementary Data 4). Our analysis finds that OA patients harboured a significant upregulation of proteins involved in DNA replication (MCM complex), cell cycle regulation (CDK1 and CDKN2A) and immune regulation (CD163, B2M, IL4I1), while AYA patients displayed an upregulation of proteins involved in mitochondrial metabolism (NDUFA9, SUCLA2, FDXR and ACADVL) and skeletal and cardiac myosin chains (MYL1, MYL2, MYLPF and MYH7). Multivariable logistic regression analysis was performed to adjust for potential confounding factors (tumour size, grade, anatomical site, performance status and histological subtype), which led to five proteins remaining significant between the two age groups (AYA: NDUFA9, SUCLA2, TUBB2B, MACROH2A2 and OA: CDK1).Fig. 2 Analysis of the adolescent and young adult (AYA) and older adult (OA) proteomic landscape.

a Annotated heatmap showing the unsupervised clustering (Pearson’s distance) of 3299 proteins across the study cohort. Patients are ordered from youngest (left) to oldest (right). From top to bottom, panels indicate tumour size, anatomical site, tumour grade, patient sex, histological subtype, and patient age group. b Volcano plot showing significantly upregulated proteins in AYA and OA patient tumours. Significant proteins (false discovery rate <0.05, fold change >2) determined by multiple t-test followed by Benjamini–Hochberg procedure are shown in yellow. Proteins that remained significant after multivariable logistic regression analysis to account for confounding factors are shown in purple. c Significantly enriched hallmark gene sets (q < 0.05) in OA patients using single sample gene set enrichment analysis (ssGSEA) scores, as determined by two-way analysis of variance (ANOVA) followed by Šidák correction. Following multivariable logistic regression to account for confounding factors, the MYC targets hallmark gene set remained significant (p = 0.047). Boxplots show ssGSEA scores, with boxes indicating the 25th and 75th percentile and the black dot indicating the 50th percentile. Whiskers extend from the 25th percentile − (1.5* interquartile range) to the 75th percentile + (1.5* interquartile range), and outliers are plotted as grey points. AS angiosarcoma, ASPS alveolar soft part sarcoma, CCS clear cell sarcoma, DDLPS dedifferentiated liposarcoma, DSRCT desmoplastic small round cell tumour, DES desmoid tumour, EPS epithelioid sarcoma, LMS leiomyosarcoma, SS synovial sarcoma, UPS undifferentiated pleomorphic sarcoma.

By undertaking single sample gene set enrichment analysis (ssGSEA), we show that compared to AYA patients, OA patients are significantly enriched (q < 0.01) for distinct biological hallmark features (Fig. 2c), including gene sets involved in cell cycle regulation (E2F targets, G2M checkpoint), oncogenic signalling (KRAS signalling, MYC targets, TNFα signalling, mTORC1 and PI3K signalling) and inflammatory pathways (inflammatory response, INFα response) (Fig. 2b). AYA patients are significantly enriched for oxidative phosphorylation (q = 0.004) and coagulation (q = 0.006) hallmarks (Supplementary Fig. 2). Multivariable logistic regression analysis was performed to adjust for potential confounding factors (tumour size, grade, anatomical site, performance status and histological subtype) which led to only the MYC targets hallmark remaining significant between the two age groups (p = 0.047).

Sarcoma proteomic modules highlight distinct biological pathways in the two age groups

We have previously identified 14 protein signatures based on sarcoma protein co-expression patterns termed Sarcoma Proteomic Modules (SPMs) which capture a broad spectrum of STS biology and transcend histological subtype (Fig. 3a and Supplementary Data 2)28. We first compared the enrichment of SPMs in the AYA and OA patients in the full cohort. Calculating patient-specific ssGSEA scores for each SPM showed significant differences between the two age groups for the majority of SPMs (Fig. 3b). To account for the possibility that histological subtype imbalances in the two age groups within the proteomic dataset may impact the SPM results, we repeated the ssGSEA in a sub-cohort with a balanced number of cases in both age groups for each histological subtype (n = 120, 60 AYA, 60 OA) (Fig. 3c). We find that SPM4 (splicing proteins) was significantly upregulated in OA patients compared to AYA patients in both the full cohort and the balanced sub-cohort. Conversely, SPM7 (immune proteins) was significantly upregulated in the AYA age group versus the OA age group in both the full cohort and the balanced sub-cohort.Fig. 3 Differential expression of sarcoma proteomic modules (SPM) in adolescent and young adult (AYA) and older adult (OA) patients.

a Protein co-expression network showing the 14 previously described SPM identified in the full proteomic cohort. Nodes indicate proteins and are coloured based on SPM membership. Edges indicate a correlation between protein expression and thickness of edges are scaled to the correlation score. b Plot showing single sample gene set enrichment analysis (ssGSEA) scores of the 14 SPMs for AYA and OA patients in the full cohort (n = 309, 66 AYA, 243 OA). c Plot showing single sample gene set enrichment analysis (ssGSEA) scores of the 14 SPMs for AYA and OA patients in a balanced number of cases in both age groups for each histological subtype (n = 120, 60 AYA, 60 OA). Error bars indicate the mean ± 1 standard deviation. Statistical significance was determined by two-way analysis of variance (ANOVA) using an uncorrected Fisher’s least significant difference (LSD) test. ****p < 0.0001, ***p < 0.001, **p < 0.01, *p < 0.05. d Sankey plot showing the distribution of each patient age group and histological subtype that falls into the SPM6-high and low expression groups. e Kaplan–Meier plot of metastasis free survival (MFS) for AYA and OA patients with high and low median expression levels of SPM6 proteins. Hazard ratio (HR), 95% confidence intervals (CI) and p value determined by univariable Cox regression. AS angiosarcoma, ASPS alveolar soft part sarcoma, CCS clear cell sarcoma, DDLPS dedifferentiated liposarcoma, DSRCT desmoplastic small round cell tumour, ECM extracellular matrix, EPS epithelioid sarcoma, LMS leiomyosarcoma, SS synovial sarcoma, UPS undifferentiated pleomorphic sarcoma.

SPM6 is a DNA replication module which we have previously shown to be prognostic for metastasis free survival (MFS) across the whole age range28. To further refine this candidate biomarker signature, we evaluated the histological subtype distribution of cases classified into SPM6-high and SPM6-low subgroups based on the median protein expression levels of the 41 proteins that make up SPM6 (Supplementary Data 2) in each of the two age groups. The Sankey plot shows that there is broad representation of histotypes in each SPM6 group, with the exception of ASPS which is only found in the SPM6-low group (Fig. 3d). In addition, there was a similar distribution of AYA or OA patients in both the SPM6-high and SPM6-low groups (Fig. 3d). Combining SPM6 and age stratified patients into four subgroups (OA-SPM6-high, OA-SPM6-low, AYA-SPM6-high, and AYA-SPM6-low) (Fig. 3e). While univariate Cox regression analysis found that the SPM6 module was able to stratify the OA patients into two groups with significantly different MFS outcomes (HR = 0.365, 95% CI 0.230–0.58, p = 1.97 × 10−5), there was no significant difference in the AYA patients (Fig. 3e). Multivariable Cox proportional hazards analysis showed that that the prognostic value of SPM6 in OA patients was independent of known prognostic factors of tumour size, grade, performance status, histological subtype and anatomical location36,37 (HR = 0.381, 95% CI = 0.219–0.661, p = 6.02 × 10−4) (Supplementary Table 2). This analysis demonstrates that SPM6 can be used as an independent risk stratification tool to identify a subgroup of OA patients (SPM6-high) with a high risk of distant relapse but has limited utility in AYA patients.

Comparative functional genomic and proteomic analyses reveal prognostic significance of the spliceosome complex

We reasoned that biological pathways that are enriched in both CRISPR-based functional genomics and MS-based proteomics datasets may yield useful candidate drug targets and biomarkers for AYA and OA patients. Here, we sought to compare the pathway information gained from in vitro functional genomic data in a panel of sarcoma cell lines with the proteomic dataset generated from soft tissue tumour patients (Fig. 4a). We first undertook an analysis of the genome-scale CRISPR-Cas9 loss-of-function screen data focusing on the NRSTS panel of cell lines within the CCLE database31,38. We identified 13 NRSTS cell lines where clinical information of patient age was available (AYA:n = 5, OA:n = 8, Supplementary Data 1) and undertook gene set enrichment analysis (GSEA) using the Reactome pathway gene sets to determine the biological pathways with selective dependencies in OA compared to AYA lines (and vice versa). Separately, we performed GSEA on the proteomic dataset of 309 soft tissue tumour patients to define differential pathways that are enriched in the two age groups. Comparing the top hits that were significantly enriched in the OA group in both the CCLE functional genomics and patient-derived proteomic datasets identified multiple overlapping shared gene sets related to mRNA splicing (Fig. 4b), while gene sets comprising the respiratory electron transport were significantly enriched in both datasets in the AYA patients (Supplementary Fig. 3). To account for the possibility that histological subtype imbalances in the two age groups within the proteomic dataset may impact the GSEA results, we repeated the GSEA of the proteomic data in two sub-cohorts: (A) SS cases only (n = 43, 19 AYA, 24 OA) and (B) Balanced number of cases in both age groups for each histological subtype (n = 120, 60 AYA, 60 OA). When GSEA was performed in each of these sub-cohorts, the mRNA splicing gene sets remained in the top two hits that were enriched in the OA population compared to the AYA group (Supplementary Fig. 4), confirming that the enrichment of this biological pathway was age-specific and not dependent on histological subtype differences between the two age groups.Fig. 4 Comparative functional genomic and proteomic analysis identifies the prognostic value of the spliceosome.

a Overview schematic of the comparative functional genomic and proteomic analysis. b Gene set enrichment analysis (GSEA) results showing the top five enriched reactome gene sets based on normalised enrichment score (NES) in the older adult (OA) group in the cell line functional genomics data from The Cancer Cell Line Encyclopaedia (CCLE) (left) and patient proteomics data (right). Overlapping gene sets enriched in both datasets are highlighted in blue. c Small nuclear ribonucleoprotein (snRNP) subunit genes that comprise the spliceosome complex found within the mRNA splicing Reactome gene set, which are identified in the CCLE cell line functional genomic and patient proteomic datasets. Genes that are only identified in the cell line functional genomic dataset are indicated with a dotted border. Genes that are core enriched as determined by GSEA leading-edge analysis in OA versus adolescent and young adult (AYA) patients in the CCLE cell line functional genomic and patient proteomic datasets are shown in green and yellow, respectively. Linkages indicate protein–protein interaction scores obtained from the STRING database, with a darker line indicating a higher score (cut-off >0.7). d Overview of multivariable Cox results for AYA and OA patients with high and low expression of each of the 21 subunits of the spliceosome complex as defined by Hegele et al. (based on median protein expression), and metastasis free survival (MFS). e Kaplan–Meier plot of MFS for AYA and OA patients with high and low expression of U2 snRNP proteins. Hazard ratio (HR), 95% confidence intervals (CI) and p value determined by univariable Cox regression. FDR false discovery rate.

Given that the OA and AYA cell lines harboured distinct dependencies within the mRNA splicing gene sets, we hypothesised that components in this pathway may serve as prognostic signatures for patients in the two age groups. In particular, we focused on the spliceosome complex, which regulates the removal of introns from precursor mRNA during the splicing process (Fig. 4c). The spliceosome is a large macromolecular complex that is comprised of >200 splicing factors that vary in their composition in a spatiotemporal manner39. We first demonstrated that ssGSEA scores for the spliceosome complex were significantly upregulated in OA patients compared to AYA patients in the full proteomic cohort (Supplementary Fig. 5A). We then undertook the same analysis in a cohort of SS cases only or a cohort comprising a balanced number of cases in both age groups for each histological subtype. In both analyses, ssGSEA scores for the spliceosome complex remain significantly higher in OA patients compared to AYA patients (Supplementary Fig. 5B, C), demonstrating that the upregulation of this complex is age-specific.

We and others have previously shown that co-regulation of splicing factors is important in the pathology of mesenchymal and epithelial tumours27,40,41. As an exemplar, Fig. 4c shows the proteins in the U1, U2, U4/5/6 small nuclear ribonucleoprotein (snRNP) subunits of the spliceosome complex that are found in OA and AYA patients in the GSEA of either the CCLE functional genomics dataset or the patient proteomic dataset. We systematically assessed the prognostic value of proteins that make up each of the spliceosome functional component subunits as defined by ref. 42. Out of a total of 21 spliceosome subunits found in our dataset, only the U2 snRNP and a non-core miscellaneous (MISC) group of splicing factors were identified to be prognostic for MFS in AYA patients in multivariable Cox analysis (Fig. 4d). Whereas there was no significant difference in MFS between the two age groups (Supplementary Fig. 1), when categorised by median expression levels of U2 snRNP proteins (n = 12 proteins, Fig. 4c), AYA patients with high U2 snRNP expression (U2-high) had superior MFS outcomes compared to those with low U2 snRNP levels (U2-low) (multivariable: HR = 4.03, 95% CI = 1.13–14.4, p = 0.0319) (Fig. 4e and Supplementary Table 3). No spliceosome subunit protein signatures were identified to be prognostic for OA in multivariable Cox analysis. These findings highlight the utility of comparative analysis of functional genomics data with proteomic profiling as a means of defining new prognostic factors in AYA patients.

Discussion

AYA patients with soft tissue tumours are an understudied age group with disparities in treatment and survival outcomes, where improvements in 5-year survival rates over the past two decades have lagged behind other age groups4. Furthermore, the use of intensive multi-modal therapy often leads to chronic health conditions and secondary malignancies in AYA patients43,44. Rather than the current “one size fits all” approach where AYA patients, in particular those with NRSTS, are offered treatments which have been optimised in OA, tailored strategies using targeted agents and risk stratification tools could have a substantial impact on survivorship and management of late effects. As a result of the under-representation of AYA patients in most molecular and biological studies in all cancer types, including STS11, there is a poor knowledge of the biological pathways that are unique to the AYA age group, which is an obstacle to developing precision medicine approaches for these patients. Notably, each of the large-scale proteomic profiling studies published thus far by The Clinical Proteomic Tumour Analysis Consortium (CPTAC) include less than 10 AYA patients45–52. Here, we present a large-scale analysis of the proteomic features of AYA and OA patients with soft tissue tumours. We show that there are inherent biological and pathway differences in the two age groups, which are maintained even when confounding variables such as tumour grade, size and histological subtypes are considered. We further demonstrate that comparative analysis of in vitro functional genomic data in a panel of NRSTS cell lines with the patient-level proteomic data leads to the prioritisation of age-specific vulnerabilities and independent prognostic factors, which provide new avenues for personalised treatment of AYA patients.

Several age-associated pan-cancer genomic analyses have shown that aging is associated with chronic inflammation and reprogramming of the immune cell landscape53. Our study finds that OA patients with soft tissue tumours are enriched in proteins involved in inflammatory response and INFα signalling hallmarks. This is consistent with a previous report by ref. 17, who demonstrated using GSEA and immune cell deconvolution of transcriptomic data that sarcoma patients that are <50 years of age have lower interferon responses and lymphocyte infiltration than those >50 years. We also determined that OA patients are enriched in proteins involved in cell cycle regulation, including the E2F targets and G2M checkpoint hallmarks. In agreement with our study, Chatsirisupachai et al., has shown in a pan-cancer analysis that mutations and somatic copy number alterations of genes within the cell cycle pathway are strongly enriched in tumours from older patients16. Interestingly, our data indicates that tumours from AYA patients harbour elevated levels of proteins involved in mitochondrial metabolism and the oxidative phosphorylation pathway, which could be indicative of metabolic rewiring in younger patients. Future investigation on the functional role of metabolic rewiring in STT subtypes in this age group is warranted54. Given several confounding factors contribute to biological differences between the AYA and OA age groups (for instance, tumour size, grade and histological subtype), following multivariable logistic regression to account for these factors, only the MYC targets hallmark remained significant. Our findings are in line with a recent study of age-specific proteomic features in colorectal cancer, which similarly showed enrichment of MYC targets in colorectal patients above 50 years of age55.

A previous analysis of genome-scale CRISPR-Cas9-based loss of function screening data in a panel of paediatric cancer cell lines identified vulnerabilities that were distinct from cell lines derived from adult patients56, suggesting that oncology drugs that are used in adult patients may not always be applicable to childhood cancers. In this study, we focused on the pathway vulnerabilities that are specific to sarcoma cell lines in either the AYA or OA age groups and undertook analysis to compare the GSEA outputs from the genome-scale CRISPR screening data with the patient-level proteomics data. Our analysis finds that in vitro pathway dependencies observed in sarcoma cell lines derived from patients of different age groups correspond to significantly higher expression levels of pathway proteins in either AYA (respiratory electron transport) or OA (mRNA splicing) patients. Sarcomas are a group of diseases of unmet need with a lack of effective therapies and novel agents. Investigational drugs that target components for each of these pathways are available for repurposing57–59 and should be evaluated in prospective studies in the different age groups. Our data further suggest that high protein expression levels of these proteins in sarcoma tissue specimens may facilitate the selection of patients who are most likely to benefit from these investigational agents and therefore should be incorporated as candidate biomarkers in clinical trial design. It should be noted that this analysis relies on the functional genomics dataset from CCLE, which suffers from several limitations. CCLE comprise immortalised cell lines which have previously been shown to include sarcoma lines which may not be fully representative of the tumours from which they were derived due to adaptation to long-term culture conditions60,61. Furthermore, there are multiple reports of inconsistencies and discordance between different functional genomics screening efforts using the same cell lines (for example, CCLE versus the Genomics of Drug Sensitivity in Cancer (GDSC) databases)62–64. It is, therefore, important to acknowledge that while functional experiments are beyond the scope of this study, there is a need for future functional validation of some of our findings in an independent set of patient-derived AYA NRSTS cell line models.

Despite optimal clinical management, a substantial proportion of NRSTS patients (up to 50%) with localised disease experience distant relapse following surgery15. Stratification of these high-risk patients has been limited to the use of nomograms, which consider known prognostic factors including tumour grade, size, histological subtype and age, amongst other variables65–67. There are currently very few molecular prognostic signatures for NRSTS and none which are optimised for AYA patients68. Here, we show that specific subunits of the spliceosome complex are independent prognostic factors in AYA patients. In particular, AYA patients with low tumour protein expression levels of the U2 snRNP spliceosome subunit are at higher risk of developing metastasis compared to those with high expression levels. These protein signatures have potential utility as precision medicine tools to tailor more aggressive treatment strategies such as peri-operative chemo/radiotherapy in AYA patients that are predicted to have higher risk of distant relapse. Conversely, low-risk AYA patients may be spared potential overtreatment, thereby reducing the risk of chronic health conditions and late effects. Mechanistically, it is not clear why AYA tumours with reduced spliceosome levels appear to have more aggressive features and future functional experiments are required to dissect the role of individual splicesome protein components in AYA sarcoma cell lines. Our study further highlights the importance of performing age-specific studies to delineate biomarkers tailored for the clinical management of AYA patients.

We acknowledge several limitations of this study. This is a retrospective cohort, which is prone to selection bias therefore, the study is hypothesis generating and our findings need to be validated in independent cohorts. Soft tissue tumours comprise a broad range of histological subtypes, and our study is limited to 10 and 7 histologies in the proteomic and functional genomics datasets, respectively. Future studies which include wider histological subtype representation is needed to determine if our findings are generalisable to all AYA patients, although this may be challenging given the limited number of publicly available AYA NRSTS cell line models available for functional studies69. There is a notable imbalance of subtypes between the two age groups in our cohort where tumours comprising complex karyotypes are prevalent in OA while those with simple genomes are enriched in the AYA group. It is important to highlight that this is reflective of the real-world incidence, with some sarcoma subtypes more likely to arise in younger patients or vice versa35. This is a key study limitation, and where possible, we have accounted for these differences using multivariable regression analysis to identify age-specific effects. Since our study does not include a comparative analysis of proteomic profiles of normal tissue from AYA and OA individuals, we cannot exclude the possibility that some of the enriched protein signatures and pathways identified in this study are the result of physiological aging rather than being tumour-specific. Despite this limitation, our data identify age-specific protein signatures with prognostic value in MFS in both AYA and OA patients which is indicative of pathological disease relevance.

In summary, we have undertaken a deep analysis of the biological differences in the proteomic profiles of soft tissue tumours from patients in the AYA and OA age groups. We highlight important protein-specific pathways and genetic vulnerabilities that are enriched in AYA patients and identify age-specific prognostic signatures to facilitate tailored clinical management of this underserved patient group.

Supplementary information

Supplementary Information

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Reporting Summary

Supplementary information

The online version contains supplementary material available at 10.1038/s43856-024-00522-x.

Acknowledgements

This study is funded by grants from the Sarah Burkeman Trust, Cancer Research UK (C56167/A29363), Children’s Cancer and Leukaemia Group and Little Princess Trust (CCLG 2023 09), Royal Marsden Cancer Charity, The Institute of Cancer Research, and the National Institute for Health Research (NIHR) Biomedical Research Centre at The Royal Marsden NHS Foundation Trust and The Institute of Cancer Research, and a charitable donation from Geoff Crocker and Bristol Care Homes to P.H.H.; Sarcoma UK (SUK03.2019) to R.L.J.; Ministry of Technology and Science of Taiwan grant (109-BOT-I-002-502) to T.W.C.

Author contributions

Conceptualisation: P.H.H.; Formal analysis: Y.B.T., K.L., H.P., and A.S.; Validation: M.C.; Data curation: J.B., C.P.W., and A.A.; Resources: T.W.C., K.T., and R.L.J.; Writing—original draft: Y.B.T. and P.H.H.; Writing—review and editing: all authors.; Supervision: A.S. and P.H.H.; and Funding acquisition: T.W.C, R.L.J., and P.H.H.

Peer review

Peer review information

Communications Medicine thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Dr Amand Schmidt and Dr Ben Abbott.

Data availability

The raw proteomic data generated in this study have been deposited in the ProteomeXchange Consortium via the PRIDE partner repository70,71 with the dataset identifier PXD036226. The normalised proteomic dataset used for the study is provided as Supplementary Data 3. The clinical data is available under restricted access due to data privacy legislation, access can be obtained by contacting the corresponding author (P.H.H.) and will require researchers to sign a data access agreement with the Institute of Cancer Research after approval by the Data Access Committee (DAC) of the Institute of Cancer Research. The DAC will determine the length of permitted access, with an expected response time frame of 2 weeks for access requests. The genome-scale CRISPR-Cas9 screening data of cell lines (CRISPRGeneEffect DepMap Public 22Q4) is available from the Cancer Cell Line Encyclopaedia (CCLE) portal (https://sites.broadinstitute.org/ccle). Additional data used for figures are provided with this paper in Supplementary Data 5.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. WHO Classification of Tumours Editorial Board. Soft Tissue and Bone Tumours (International Agency for Research on Cancer, 2020).
2. Trama A Botta L Steliarova-Foucher E Cancer burden in adolescents and young adults: a review of epidemiological evidence Cancer J. 2018 24 256 266 10.1097/PPO.0000000000000346 30480570
3. Papworth KE Soft-tissue sarcoma in adolescents and young adults compared with older adults: A report among 5000 patients from the Scandinavian Sarcoma Group Central Register Cancer 2019 125 3595 3602 10.1002/cncr.32367 31287163
4. Bleyer A Montello M Budd T Saxman S National survival trends of young adults with sarcoma: lack of progress is associated with lack of clinical trial participation Cancer 2005 103 1891 1897 10.1002/cncr.20995 15795902
5. Shaw PH Hayes-Lattin B Johnson R Bleyer A Improving enrollment in clinical trials for adolescents with cancer Pediatrics 2014 133 S109 S113 10.1542/peds.2014-0122F 24918208
6. Smith AW Health-related quality of life of adolescent and young adult patients with cancer in the United States: the Adolescent and Young Adult Health Outcomes and Patient Experience study J. Clin. Oncol. 2013 31 2136 2145 10.1200/JCO.2012.47.3173 23650427
7. Zebrack BJ Psychological, social, and behavioral issues for young adults with cancer Cancer 2011 117 2289 2294 10.1002/cncr.26056 21523748
8. van der Graaf WTA Orbach D Judson IR Ferrari A Soft tissue sarcomas in adolescents and young adults: a comparison with their paediatric and adult counterparts Lancet Oncol. 2017 18 e166 e175 10.1016/S1470-2045(17)30099-2 28271871
9. Younger E Age-related sarcoma patient experience: results from a national survey in England BMC Cancer 2018 18 991 10.1186/s12885-018-4866-8 30333006
10. Taylor RM Description of the BRIGHTLIGHT cohort: the evaluation of teenage and young adult cancer services in England BMJ Open 2019 9 e027797 10.1136/bmjopen-2018-027797 31005941
11. Ferrari A Adolescents and young adults (AYA) with cancer: a position paper from the AYA Working Group of the European Society for Medical Oncology (ESMO) and the European Society for Paediatric Oncology (SIOPE) ESMO Open 2021 6 100096 10.1016/j.esmoop.2021.100096 33926710
12. Bleyer A The distinctive biology of cancer in adolescents and young adults Nat. Rev. Cancer 2008 8 288 298 10.1038/nrc2349 18354417
13. Younger E Metastatic soft tissue sarcomas in adolescents and young adults: a specialist center experience J. Adolesc. Young Adult Oncol. 2020 9 628 638 10.1089/jayao.2020.0010 32379517
14. Reed DR Naghavi A Binitie O Sarcoma as a model for adolescent and young adult care J. Oncol. Pract. 2019 15 239 247 10.1200/JOP.18.00684 31075215
15. Gronchi A Soft tissue and visceral sarcomas: ESMO-EURACAN-GENTURIS Clinical Practice Guidelines for diagnosis, treatment and follow-up(✰) Ann. Oncol. 2021 32 1348 1365 10.1016/j.annonc.2021.07.006 34303806
16. Chatsirisupachai K Lesluyes T Paraoan L Van Loo P de Magalhaes JP An integrative analysis of the age-associated multi-omic landscape across cancers Nat. Commun. 2021 12 2345 10.1038/s41467-021-22560-y 33879792
17. Lee W Wang Z Saffern M Jun T Huang KL Genomic and molecular features distinguish young adult cancer from later-onset cancer Cell Rep. 2021 37 110005 10.1016/j.celrep.2021.110005 34788626
18. Li CH Haider S Boutros PC Age influences on the molecular presentation of tumours Nat. Commun. 2022 13 208 10.1038/s41467-021-27889-y 35017538
19. Shah Y Pan-cancer analysis reveals molecular patterns associated with age Cell Rep. 2021 37 110100 10.1016/j.celrep.2021.110100 34879281
20. Wang X Langevin AM Houghton PJ Zheng S Genomic disparities between cancers in adolescent and young adults and in older adults Nat. Commun. 2022 13 7223 10.1038/s41467-022-34959-2 36433963
21. Steppan DA Pratilas CA Loeb DM Targeted therapy for soft tissue sarcomas in adolescents and young adults Adolesc. Health Med. Ther. 2017 8 41 55 28408855
22. Morfouace M Comprehensive molecular profiling of sarcomas in adolescent and young adult patients: results of the EORTC SPECTA-AYA international proof-of-concept study Eur. J. Cancer 2023 178 216 226 10.1016/j.ejca.2022.10.020 36470093
23. Gounder MM Clinical genomic profiling in the management of patients with soft tissue and bone sarcoma Nat. Commun. 2022 13 3406 10.1038/s41467-022-30496-0 35705558
24. Burns J Wilding CP Jones RL Huang PH Proteomic research in sarcomas - current status and future opportunities Semin. Cancer Biol. 2020 61 56 70 10.1016/j.semcancer.2019.11.003 31722230
25. Chadha M Huang PH Proteomic and metabolomic profiling in soft tissue sarcomas Curr. Treat. Options Oncol. 2022 23 78 88 10.1007/s11864-022-00947-3 35171456
26. Vyse S Thway K Huang PH Jones RL Next-generation sequencing for the management of sarcomas with no known driver mutations Curr. Opin. Oncol. 2021 33 315 322 10.1097/CCO.0000000000000741 33927108
27. Milighetti M Proteomic profiling of soft tissue sarcomas with SWATH mass spectrometry J. Proteomics 2021 241 104236 10.1016/j.jprot.2021.104236 33895336
28. Burns J The proteomic landscape of soft tissue sarcomas Nat. Commun. 2023 14 3834 10.1038/s41467-023-39486-2 37386008
29. Troyanskaya O Missing value estimation methods for DNA microarrays Bioinformatics 2001 17 520 525 10.1093/bioinformatics/17.6.520 11395428
30. Subramanian A Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles Proc. Natl Acad. Sci. USA 2005 102 15545 15550 10.1073/pnas.0506580102 16199517
31. Barretina J The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity Nature 2012 483 603 607 10.1038/nature11003 22460905
32. Shannon P Cytoscape: a software environment for integrated models of biomolecular interaction networks Genome Res. 2003 13 2498 2504 10.1101/gr.1239303 14597658
33. Szklarczyk D STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets Nucleic Acids Res. 2019 47 D607 D613 10.1093/nar/gky1131 30476243
34. Gounder MM Thomas DM Tap WD Locally aggressive connective tissue tumors J. Clin. Oncol. 2018 36 202 209 10.1200/JCO.2017.75.8482 29220303
35. Weiss, A. R. & Harrison, D. J. Soft tissue sarcomas in adolescents and young adults. J. Clin. Oncol. 42, 675–685 (2023).
36. Pisters PW Leung DH Woodruff J Shi W Brennan MF Analysis of prognostic factors in 1,041 patients with localized soft tissue sarcomas of the extremities J. Clin. Oncol. 1996 14 1679 1689 10.1200/JCO.1996.14.5.1679 8622088
37. Coindre JM Prognostic factors in adult patients with locally controlled soft tissue sarcoma. A study of 546 patients from the French Federation of Cancer Centers Sarcoma Group J. Clin. Oncol. 1996 14 869 877 10.1200/JCO.1996.14.3.869 8622035
38. Ghandi M Next-generation characterization of the Cancer Cell Line Encyclopedia Nature 2019 569 503 508 10.1038/s41586-019-1186-3 31068700
39. Wahl MC Will CL Luhrmann R The spliceosome: design principles of a dynamic RNP machine Cell 2009 136 701 718 10.1016/j.cell.2009.02.009 19239890
40. Koedoot E Co-regulated gene expression of splicing factors as drivers of cancer progression Sci. Rep. 2019 9 5484 10.1038/s41598-019-40759-4 30940821
41. Krasny L Proteomic profiling identifies co-regulated expression of splicing factors as a characteristic feature of intravenous leiomyomatosis Cancers 2022 14 2907 10.3390/cancers14122907 35740573
42. Hegele A Dynamic protein-protein interaction wiring of the human spliceosome Mol. Cell 2012 45 567 580 10.1016/j.molcel.2011.12.034 22365833
43. Livingston JA Improving outcomes for adolescents and young adults with sarcoma: a focus on cancer care delivery J. Oncol. Pract. 2019 15 253 254 10.1200/JOP.19.00170 31075216
44. Youn P Milano MT Constine LS Travis LB Long-term cause-specific mortality in survivors of adolescent and young adult bone and soft tissue sarcoma: a population-based study of 28,844 patients Cancer 2014 120 2334 2342 10.1002/cncr.28733 24752471
45. Cao L Proteogenomic characterization of pancreatic ductal adenocarcinoma Cell 2021 184 5031 5052 e5026 10.1016/j.cell.2021.08.023 34534465
46. Clark DJ Integrated proteogenomic characterization of clear cell renal cell carcinoma Cell 2020 180 207 10.1016/j.cell.2019.12.026 31923397
47. Dou Y Proteogenomic characterization of endometrial carcinoma Cell 2020 180 729 748 e726 10.1016/j.cell.2020.01.026 32059776
48. Gillette MA Proteogenomic characterization reveals therapeutic vulnerabilities in lung adenocarcinoma Cell 2020 182 200 225 e235 10.1016/j.cell.2020.06.013 32649874
49. Krug K Proteogenomic landscape of breast cancer tumorigenesis and targeted therapy Cell 2020 183 1436 1456 e1431 10.1016/j.cell.2020.10.036 33212010
50. Satpathy S A proteogenomic portrait of lung squamous cell carcinoma Cell 2021 184 4348 4371 e4340 10.1016/j.cell.2021.07.016 34358469
51. Wang LB Proteogenomic and metabolomic characterization of human glioblastoma Cancer Cell 2021 39 509 528 e520 10.1016/j.ccell.2021.01.006 33577785
52. Zhang H Integrated proteogenomic characterization of human high-grade serous ovarian cancer Cell 2016 166 755 765 10.1016/j.cell.2016.05.069 27372738
53. Chatsirisupachai K Lagger C de Magalhaes JP Age-associated differences in the cancer molecular landscape Trends Cancer 2022 8 962 971 10.1016/j.trecan.2022.06.007 35811230
54. Miallot R Galland F Millet V Blay JY Naquet P Metabolic landscapes in sarcomas J. Hematol. Oncol. 2021 14 114 10.1186/s13045-021-01125-y 34294128
55. Gong Y Age-associated proteomic signatures and potential clinically actionable targets of colorectal cancer Mol. Cell Proteomics 2021 20 100115 10.1016/j.mcpro.2021.100115 34129943
56. Dharia NV A first-generation pediatric cancer dependency map Nat. Genet. 2021 53 529 538 10.1038/s41588-021-00819-w 33753930
57. Murphy AJ Li AH Li P Sun H Therapeutic targeting of alternative splicing: a new frontier in cancer treatment Front. Oncol. 2022 12 868664 10.3389/fonc.2022.868664 35463320
58. Lee SC Abdel-Wahab O Therapeutic targeting of splicing in cancer Nat. Med. 2016 22 976 986 10.1038/nm.4165 27603132
59. Greene J Segaran A Lord S Targeting OXPHOS and the electron transport chain in cancer; Molecular and therapeutic implications Semin. Cancer Biol. 2022 86 851 859 10.1016/j.semcancer.2022.02.002 35122973
60. Hemming ML Oncogenic gene-expression programs in leiomyosarcoma and characterization of conventional, inflammatory, and uterogenic subtypes Mol. Cancer Res. 2020 18 1302 1314 10.1158/1541-7786.MCR-20-0197 32518213
61. Sandberg R Ernberg I Assessment of tumor characteristic gene expression in cell lines using a tissue similarity index (TSI) Proc. Natl Acad. Sci. USA 2005 102 2052 2057 10.1073/pnas.0408105102 15671165
62. Haibe-Kains B Inconsistency in large pharmacogenomic studies Nature 2013 504 389 393 10.1038/nature12831 24284626
63. Haverty PM Reproducible pharmacogenomic profiling of cancer cell line panels Nature 2016 533 333 337 10.1038/nature17987 27193678
64. Mpindi JP Consistency in drug response profiling Nature 2016 540 E5 E6 10.1038/nature20171 27905421
65. Pasquali S Neoadjuvant chemotherapy in high-risk soft tissue sarcomas: a sarculator-based risk stratification analysis of the ISG-STS 1001 randomized trial Cancer 2022 128 85 93 10.1002/cncr.33895 34643947
66. van Praag VM A prediction model for treatment decisions in high-grade extremity soft-tissue sarcomas: personalised sarcoma care (PERSARC) Eur. J. Cancer 2017 83 313 323 10.1016/j.ejca.2017.06.032 28797949
67. Callegaro D Development and external validation of two nomograms to predict overall survival and occurrence of distant metastases in adults after surgical resection of localised soft-tissue sarcomas of the extremities: a retrospective analysis Lancet Oncol 2016 17 671 680 10.1016/S1470-2045(16)00010-3 27068860
68. Merry E Thway K Jones RL Huang PH Predictive and prognostic transcriptomic biomarkers in soft tissue sarcomas npj Precis. Oncol. 2021 5 1 8 33479506
69. Hattori E Oyama R Kondo T Systematic review of the current status of human sarcoma cell lines Cells 2019 8 157 10.3390/cells8020157 30781855
70. Deutsch EW The ProteomeXchange consortium in 2020: enabling ‘big data’ approaches in proteomics Nucleic Acids Res. 2020 48 D1145 D1152 31686107
71. Perez-Riverol Y The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences Nucleic Acids Res. 2022 50 D543 D552 10.1093/nar/gkab1038 34723319
