Molecular Feature-Based Classification of Retroperitoneal Liposarcoma: A Prospective Cohort Study ================================================================================================= * Mengmeng Xiao * Xiangji Li * Fanqin Bu * Shixiang Ma * Xiaohan Yang * Jun Chen * Yu Zhao * Ferdinando Cananzi * Chenghua Luo * Li Min ## ABSTRACT **Background** Retroperitoneal liposarcoma (RPLS) is a critical malignant disease with various clinical outcomes. However, the molecular heterogeneity of RPLS was poorly elucidated, and few biomarkers were proposed to monitor its progression. **Methods** RNA sequencing was performed on a training cohort of 88 RPLS patients to identify dysregulated genes and pathways using clusterprofiler. The GSVA algorithm was utilized to assess signaling pathways levels in each sample, and unsupervised clustering was employed to distinguish RPLS subtypes. Differentially expressed genes (DEGs) between RPLS subtypes were identified to construct a simplified dichotomous clustering via nonnegative matrix factorization. The feasibility of this classification was validated in a separate validation cohort (n=241) using immunohistochemistry (IHC) from the Retroperitoneal SArcoma Registry (RESAR). The study is registered with [ClinicalTrials.gov](https://ClinicalTrials.gov) under number [NCT03838718](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03838718&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom). **Results** Cell cycle, DNA damage & repair, and Metabolism were identified as the most aberrant biological processes in RPLS, enabling the division of RPLS patients into two distinct subtypes with unique molecular signatures, tumor microenvironment, clinical features and outcomes (overall survival, OS and disease-free survival, DFS). A simplified RPLS classification based on representative biomarkers (LEP and PTTG1) demonstrated high accuracy (AUC>0.99), with patients classified as LEP+ and PTTG1-showing lower aggressive pathological composition ratio and fewer surgery times, along with better OS (HR=0.41, *P*<0.001) and DFS (HR=0.60, *P*=0.005). **Conclusions** Our study provided an ever-largest gene expression landscape of RPLS and established an IHC-based molecular classification that was clinically relevant and cost-effective for guiding treatment decisions. Keywords * LEP * PTTG1 * Retroperitoneal liposarcoma (RPLS) * Molecular classification ## INTRODUCTION Retroperitoneal liposarcoma (RPLS) is a soft tissue sarcoma (STS) originating in the retroperitoneum with an insidious onset. Traditional surgical resection has been regarded as a primary and curable treatment strategy of RPLS for the past fifty years (***Ecker et al., 2016***). However, the anatomical complexity and biological properties of sarcoma brought great difficulty to achieve microscopically margin-negative resection, leading to a high postoperative recurrence rate in RPLS patients. During the past decade, scientists tried to improve the postoperative survival of RPLS patients by personalized surgical resection and neoadjuvant/adjuvant therapies, but the effect was not satisfactory (***Littau et al., 2020; Gronchi et al., 2015; Gronchi et al., 2009; Gronchi et al., 2013; Pisters et al., 2009***). Recently, precision medicine greatly enriched the therapeutic approaches and reformed the clinical decision-making chain of tumor diagnosis and treatment, prolonging the median survival of main tumor types 2-10 times (***Kam et al., 2021; Zeng et al., 2022; Alifrangis et al., 2019; Frese et al., 2021***). Biomarker-based patient stratification and targeted therapy together make up the kernel of precision medicine, which is intrinsically based on the molecular profiling of cancers. However, our knowledge of the molecular features of RPLS is limited, and few clinically applicable molecular biomarkers and targeted drugs are available for RPLS treatment. Only sporadic molecules such as CDK4 (***Pilotti et al., 2000***), MDM2 (***Binh et al., 2005***), AURK4 (***Yen et al., 2019***) and CCNDBP1 (***Yang et al., 2021***) have been reported as prognostic and diagnostic biomarkers, but these biomarkers were poorly represented and verified. Therefore, it is crucial to reveal the molecular landscape of RPLS and explore a feasible classification for its diagnosis and treatment. Here, we conducted a comprehensive inbestigation into the molecular characteristics of RPLS through the delineation of the largest gene expression landscape ever assembled for this rare disease entity. By identifying both RPLS-specific genes and prognostic biomarkers, we unveiled their intricate relationships with clinical parameters. Our findings revealed the existence of two distinct molecular subtypes within all RPLS patients, characterized by diverse pathological compositions, enriched signaling pathways, and varying clinical outcomes. This highlights the limitations of relying solely on traditional pathological classification for surgical decision-making in certain cases where patients exhibit favorable histological features but poor prognoses. Emphasizing the pivotal role of molecular subtyping in guiding individualized treatment strategies and enhancing patient management. To facilitate practical application in clinical settings, we developed a simplified RPLS classification system based on key biomarkers (LEP and PTTG1) representative of each subtype. Notably, this classification scheme was validated in a larger cohort of RPLS patients through immunohistochemistry assays (Figure 1), laying the groundwork for precise surgical interventions guided by molecular insights in the realm of RPLS treatment. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F1) Figure 1. Flow diagram of exploring RPLS dichotomous classification ## MATERIALS AND METHODS ### Patients and tissue specimens Patients who diagnosed with RPLS amenable to surgical resection were eligible for the study. The RPLS histology was confirmed according to the WHO criteria done on biopsy or surgical specimen by dedicated sarcoma pathologist. The exclusion criteria included the age<18 years; serious psychiatric disease that precludes informed consent or limits compliance; impossibility to ensure adequate follow-up. Tumor specimens from 88 RPLS patients (Training cohort 1, Table S1; Training cohort 2, Table S2) and another cohort of 241 RPLS patients (Validation cohort, Table S3) were obtained from our local Hospital. These cohorts are sourced from Retroperitoneal SArcoma Registry (RESAR, [NCT03838718](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03838718&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom)). All the patients underwent curative resection from January 2015 to May 2019. RPLS tissue specimens were snap-frozen in liquid nitrogen within 1 h and then stored in a -80 ℃ refrigerator before use. Clinical information was collected from the medical records, and no patient had undergone previous chemotherapy or radiotherapy. Overall survival (OS) was defined as the interval between the latest surgery and death from tumors or between the latest surgery and the last observation taken for surviving patients. Disease-free survival (DFS) was defined as the interval between the latest surgery and diagnosis of relapse or death. Informed consent for surgical procedures and specimen collection were obtained from each patient. This study has been reported in line with the REMARK criteria (***McShane et al., 2005***). ### RNA sequencing, primary data processing, and analysis Total RNA was extracted from Training cohort 1 (Table S4) and Training cohort 2 (Table S2) using TRIzol Reagent (Invitrogen). RNA degradation and contamination were monitored with 1% agarose gel. RNA purity was checked by the NanoPhotometer spectrophotometer (IMPLEN, Los Angeles, CA, USA). RNA concentration was measured using the Qubit RNA Assay Kit with the Qubit 2.0 Fluorometer (Life Technologies, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 System (Agilent Technologies, CA, USA). A total amount of 3-5 ug RNA per sample was used as input material for the RNA library. Sequencing libraries were generated using NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (NEB, USA) following the manufacturer’s recommendations and index codes were added to attribute sequences to each sample. The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq SR Cluster Kit v3-cBot-HS (Illumia). After cluster generation, the strand-specific cDNA were sequenced on an Illumina NovaSeq 6000 platform, and single-end reads were generated (Novogene Bioinformatic Technology, Beijing, China). FPKMs of mRNAs and non-coding RNAs in each sample were calculated by Cuffdiff (v2.1.1). FPKMs were calculated based on the length of the fragments and read counts mapped to this fragment. These sequencing data have been deposited at the Open Archive for Miscellaneous Data (OMIX) database of China National Center for Bioinformation (CNCB) under the accession number OMIX002786. ### Identification of differential genes Gene difference analysis was performed to determine the differential genes (DEGs). An adjusted *FDR*<0.05 and |log2FC|>0.585 was considered significant. This process was conducted with the R package “limma”. ### Identification of prognostic genes Cox univariate regression analysis was used to screen the prognostic genes of RPLS. Results of *P*<0.05 was considered significant. This process was conducted with the R package “survival”. ### Gene set enrichment analysis (GSEA) and immune infiltrate analysis GSEA was performed in the tumor and normal groups to explore the biological signaling pathways. Pathway annotation files were downloaded from the msigdb ([www.gsea-msigdb.org](https://www.gsea-msigdb.org)) platform. This process was conducted by the GSEA R package to elucidate the representative HALLMARK and REACTOME pathways enriched in RPLS patients. Immunocyte infiltration (immune score and stromal score) was measured by the Estimation of STromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) algorithm. This process was completed via the “estimate” R package. ### Functional annotation Functional enrichment analyses were performed to elucidate the possible biological processes and signaling pathways of the prognostic genes. Gene ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) analyses were conducted by R package “clusterprofiler”, and the false discovery rate<0.05 was considered significantly enriched. ### Consensus clustering with t-distributed stochastic neighbor embedding (t-SNE) After evaluated the relative abundance level of related pathways, the Euclidean distance was calculated between any two samples and condensed into two-dimensional points using t-distributed stochastic neighbor embedding (t-SNE) (***Guo et al., 2019***) and subsequently visualized automatically with the density-based spatial clustering of applications with noise (DBSCAN) algorithm. This consensus clustering was conducted with the R packages “Rtsne” and “dbscan”. ### Consensus clustering with nonnegative matrix factorization (NMF) Nonnegative matrix factorization (NMF) was used to perform RPLS subtyping. Specifically, NMF was applied to gene expression matrix *A* which contained gene sets of major signaling pathways and prognostic genes. Matrix *A* was factorized into 2 nonnegative matrices *W* and *H*. Repeated factorization of matrix *A* was performed and its outputs were aggregated to obtain consensus clustering of RPLS samples. The optimal number of subtypes was selected according to cophenetic, dispersion, and silhouette coefficients. This consensus clustering was conducted with the R package “NMF”. ### Construction of machine learning models Machine learning models based on biomarkers were constructed by logistic regression (LR), support vector machine (SVM), and random forest (RF). These models were specifically tailored to analyze biomarker data in order to predict clinical outcomes in surgical patients. LR is a statistical method that establishes a relationship between asset of independent variables and a binary outcome. It calculates the probability of an event occurring based on the input features derived from biomarkers relevant to surgical patient. SVM is a supervised learning algorithm that categorizes data points by identifying the optimal hyperplane that separates distinct classes within a high-dimensional space. This approach effectively maps biomarker data into a multidimensional space to facilitate accurate classification of patient outcomes. RF is an ensemble learning technique that generated multiple decision trees during training and aggregates the results to make predictions. By leveraging this method, we can enhance predictive accuracy by mitigating overfitting and increasing model robustness when analyzing biomarker-driven patient data. The performance of these machine learning models was assessed using the area under the curve (AUC) metric. A higher AUC value indicates superior discriminatory power of model in distinguishing different clinical outcomes. An AUC value closer to 1.0 signifies strong predictive capability, while 0.5 indicates no discriminatory ability at all. By evaluating the AUC values generated by LR, SVM, and RF models, clinicians can identify which algorithm yields the most reliable predictions based on biomarker profiles for surgical patients. ### Immunohistochemistry The protocol was performed as previously described (***Li et al., 2022***). In brief, the LEP and PTTG1 antibodies for immunohistochemistry were purchased from Proteintech (Cat No: bs-0409R and bs-1881R). With deparaffinization for 15min × 3 in dimethylbenzene and routine hydration, the tissues were soaked in phosphate buffer saline (PBS) for 10min and then performed high-pressure antigen retrieval (Tris-EDTA, PH=9.0) for 2.5min. After being treated with a 3% endogenous catalase blocker (ZSBIO, PV-6000) for 10min, the tissues were incubated in goat serum (ZSBIO, ZLI-9022) for the blocking of nonspecific reaction and then incubated with primary antibody (LEP=1:300 and PTTG1=1:300) at 4℃ overnight. The next day, tissues were washed and incubated with goat anti-rabbit secondary antibody (ZSBIO, PV-9000) for 1h at room temperature, then washed and stained with DAB reagents (ZSBIO, ZLI-9018). Then hematoxylin staining, 1% hydrochloric acid alcohol differentiation, ammonia water anti-blue, and neutral gum sealing. The IHC results were evaluated by pathologists, the staining extent was scored as 0-100%. The intensity score was defined as negative, low-expression, medium-expression, and high-expression, which were documented as 0, 1, 2, and 3 respectively. The final scores were calculated by the formula: *IHC score = Staining extent score* × *Staining intensity score*. ### Statistical Methods R software (version 4.1.3) was used in this study. For quantitative variables, differences between the two groups and among multiple groups were analyzed by Wilcoxon’s test and One-way analysis of variance (ANOVA), respectively. For categoric variables, groups were compared by use of Chi-square test. Survival curves were determined by Log-rank test. The clinicopathological features and levels of immune infiltration were conducted by Wilcoxon’s test. A difference of *P*<0.05 indicated statistical significance unless specified otherwise. ## RESULTS Baseline characteristics were shown in Table 1. Of 329 RPLS patients, 88 in training cohort and 241 in validation cohort. No statistically significant differences were found in the age, sex, pathology, surgery times, tumor size, and multilocation between the two cohorts (*P*>0.05). The IHC score of LEP and PTTG1 were 1.62 (0.820) and 0.830 (0.75) in validation cohort. View this table: [Table 1.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/T1) Table 1. Baseline characteristics of training cohort and validation cohort ### Cell cycle, DNA damage and repair, and Metabolism are dysregulated in RPLS To reveal the general molecular features of RPLS compared to noncancerous adipose tissues, we first recruited 8 RPLS patients and collected paired tumor and normal tissues for differentially expressed gene (DEG) analysis. A total of 1354 DEGs, 554 upregulated and 800 downregulated, were identified (Figure 2A-B). To assess the underlying pathways of RPLS, GSEA analyses were performed for those DEGs. We found that proliferation-associated pathways, such as Mitotic spindle, E2F target, G2/M checkpoint and Separation of sister chromatids, were mainly enriched in tumors; while metabolism-related pathways, such as Bile acid metabolism, Heme and fatty acid metabolism and Integration of energy metabolism, were enriched in normal controls (Figure 2C). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F2) Figure 2. Cell cycle, DNA damage and repair, and metabolism are dysregulated in RPLS Volcano plot of the DEGs in 8 normal vs 8 RPLS tissues **(A)**. Venn diagram showed shared genes between DEGs and prognostic genes **(B)**. GSEA analysis of RPLS tumors, including HALLMARK gene sets and REACTOME gene sets **(C)**. Circular plots of the prognostic genes in GO, KEGG, and enrichWP **(D)**. Then, we collected another 80 samples to investigate the molecular heterogeneity of RPLS. Gene expression profiles showed 918 and 3244 genes associated with overall survival (OS) and disease-free survival (DFS), respectively. Among 497 candidate genes associated with both OS and DFS, 83 of them also overlapped with DEGs (Figure 2B). Functional annotation (GO, KEGG, and enrichWP) demonstrated that Cell cycle, DNA damage and repair, and Metabolism-related pathways were significantly enriched (Figure 2D), suggesting these signaling pathways were dysregulated in RPLS. ### RPLS subgroups based on molecular features show different clinical outcomes To evaluate heterogeneous molecular clustering characteristics in RPLS, ssGSEA emerged as a widely adopted method for computing the enrichment level of specific biological signaling pathways for each sample based on gene expression data. This aids in gaining insights into the overall activity level of signaling pathways. Here, we scored each sample on the dysregulated pathways by ssGSEA and divided RPLS patients into two subgroups (Figure 3A). Subgroup 1 (G1) showed better OS and DFS compared to subgroup 2 (G2) (Figure 3B-C), G1 displayed elevated ssGSEA scores associated with Metabolism, whereas G2 exhibited heightened ssGSEA scores linked on Cell cycle and DNA damage and repair (Figure 3D). These features suggested that effective monitoring of the prognosis of RPLS patients can be achieved based on the activation status of specific pathways. We also evaluated the clinical features and immune infiltration levels of those samples. The results showed that G1 had lower aggressive pathological composition ratio (Figure 3E), MDM2 (Figure 3G) and Ki67 (Figure S1A) expression, larger tumor size (Figure S1B), and higher tumor microenvironment (TME) level compared to G2 (Figure 3H). Surgery times for G1 were also tended to decrease (Figure 3F). Taken together, the above restuls indicated that RPLS subgroups based on molecular features showed distinct clinical features and clinical outcomes. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F3) Figure 3. RPLS subgroups (G1 and G2) based on cell cycle, DNA damage and repair, and metabolism tSNE exhibited the subgroups (G1 and G2) of RPLS **(A)**. Survival cures of OS **(B)** and DFS **(C)** in G1 and G2. The hierarchical clustering heatmap of dysregulated pathways in G1 and G2 **(D)**. Histograms revealed the difference of pathological composition ratio **(E)**, surgery times **(F)**, and MDM2 **(G)** in G1 and G2. Violin plot of the microenvironmental scores in G1 and G2 **(H)**. ### A simplified RPLS classification strategy derived from RPLS G1/G2 subgroups To explore representative biomarkers for different RPLS subgroups, we performed a DEG analysis between G1 and G2. There were 1258 genes downregulated, among which 112 of them indicated good prognosis (protective genes). Correspondingly, 754 genes were upregulated, and 28 of them indicated poor prognosis (aggressive genes) (Figure S1C-D). Enrichment analysis suggested that those DEGs were also associated with cell cycle regulation and metabolism (Figure S1E), which was consistent with previous results (Figure 2C). To develop a simplified RPLS clustering based on DEGs, we adopted NMF and tSNE for a re-classification of those patients. The results showed RPLS patients were also divided into two clusters (Figure 4A and Figure S1F-G). We then annotated the samples of two clusters by ssGSEA and found Cluster1 (C1) was related to metabolic processes, and Cluster2 (C2) was mainly related to the processes of Cell cycle and DNA damage and repair (Figure S2A). Also, C1 showed better OS and DFS, lower pathological composition ratio and MDM2 expression, and fewer surgery times (Figure 4B-F). Lower Ki67 expression and larger tumor size were observed in C2 (Figure S2B-C). Interestingly, the biological annotations of the C1/C2 classification were greatly consistent with G1/G2. Therefore, a simplified RPLS classification strategy derived from RPLS subgroups was provisionally established. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F4) Figure 4. RPLS classification strategy (C1 and C2) derived from RPLS subgroups NMF for a re-classification of training cohort 1 (C1 and C2) **(A)**. Survival cures of OS **(B)** and DFS **(C)** in C1 and C2. Histograms revealed the difference in pathological composition ratio **(D)**, MDM2 **(E)**, and surgery times **(F)** in C1 and C2. ### Development of a dichotomous RPLS classification model For NMF classification of RPLS patients, LEP and PTTG1 were identified as representative biomarkers of C1 and C2, respectively (Figure 5A). We aimed to replicate the RPLS classification of C1 and C2 by integrating these two biomarkers with the assistance of machine learning algorithms, and this two-gene panel achieved promising results (Logistic, AUC=0.995; SVM, AUC=0.997; RF, AUC=1.000; Figure 5B). Also, a linear negative correlation between LEP and PTTG1 expression was detected (Figure 5C). Considering the enhanced interpretability and generalization of linear models, we adopted the results of Logistic regression for subsequent analysis (Risk values=2.182×PTTG1-2.204×LEP). The patients marked as high-risk (Cluster\_H) exhibited worse OS and DFS than those marked as low-risk (Cluster\_L) (Figure 5D-E). Dysregulated pathways, such as DNA repair and Cell cycle regulation, were enriched in Cluster_H (Figure 5F and Figure S2D), and Cluster\_H presented more aggressive pathological composition ratio, higher MDM2 levels and marginally increased in surgery times than Cluster_L (Figure 5G-I). Similarly, the Cluster_H showed higher Ki67 level and smaller tumor sizes (Figure S2E-F). Moreover, a Sankey diagram was drawn to show the correlation among G1/G2, C1/C2, and Cluster_L/H. Cluster_L/H were well-matched to C1/C2 and G1/G2, suggesting LEP and PTTG1 were promising biomarkers for a dichotomous RPLS classification (Figure 5J). ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F5) Figure 5. RPLS dichotomous classification (Cluster_C1 and Cluster_C2) derived from RPLS clusters Heatmap of biomarkers identified (LEP and PTTG1) in C1 and C2 **(A)**. ROC curves of the machine learning models to identify C1 and C2 **(B)**. Correlation between LEP and PTTG1 expression **(C)**. Survival curves of OS **(D)** and DFS **(E)** in Cluster\_C1 (low risk) and Cluster\_C2 (high risk) groups. GSEA of HALLMARK gene sets in Cluster\_C1 and Cluster_C2 **(F)**. Histograms revealed the difference of pathological composition ratio **(G)**, MDM2 level **(H)**, and surgery times **(I)** in Cluster_C1 and Cluster_C2. Sankey diagram indicated the correlation among G1/G2, C1/C2, and Cluster_C1/C2 **(J)**. ### Validation of the dichotomous RPLS classification in another 241 RPLS patients To validate LEP and PTTG1 as biomarkers for a dichotomous RPLS classification, we performed IHC staining of two biomarkers in validation cohort. The representative images of LEP and PTTG1 with different expression levels were shown in Figure 6A-B. The IHC scores were integrated with the previously fitted coefficients to evaluate the prognosis of RPLS patients (Risk values=2.182×PTTG1IHC-2.204×LEPIHC). The cutoff value of validation cohort is the median of risk value. The high-risk group had worse OS and DFS (Figure 6C-D), along with more surgery times and more aggressive pathological composition ratio (Figure 7A-B), but the difference of tumor size was not observed between the two groups (Figure 7C). Then we constructed visual nomograms for a precise survival prediction of RPLS patients by combining the risk score with clinical features. The predictive abilities of the 1-, 2-, and 3-year OS (Figure 7D-F) and DFS (Figure S3A-C) were 0.743-0.788. Together, we proposed a simple and clinically applicable molecular classification strategy for RPLS patients. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F6) Figure 6. Validation of the RPLS dichotomous classification in another 241 RPLS cohort Representative IHC staining images of LEP **(A)** and PTTG1 **(B)**. Survival curves of OS **(C)** and DFS **(D)** in high-risk and low-risk groups. ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/09/2024.08.09.24311657/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2024/08/09/2024.08.09.24311657/F7) Figure 7. Survival nomogram of LEP+PTTG1 model in validation cohort The difference of surgery times **(A)**, pathological composition **(B)**, and surgery times **(C)** in high-risk and low-risk groups. Nomograms for OS was developed in REASR cohort with four factors: sex, age, risk score, and differentiation **(D)**. ROC curves of 1-, 2-, and 3-year OS in validation cohort **(E)**. Calibration curves of predicting 1-, 2-, and 3-year OS in validation cohort **(F)**. ## DISCUSSION Here we divided RPLS patients into two subgroups based on Cell cycle, DNA damage & repair, and Metabolism-related pathways. G1 was annotated as Metabolism-active, which exhibited high ssGSEA scores on Metabolism-associated pathways, while G2 showed high ssGSEA scores on Cell cycle and DNA damage and repair with a high Ki67 and MDM2 level. G2 had more aggressive molecular features and worse clinical outcomes compared to G1, in accordance with previously reported tumor classification (***Lindskrog et al., 2021; Yu et al., 2021; Zhang et al., 2022***). In fact, ***Demicco et al., 2017*** has integrated SCNA and DNA methylation divide dedifferentiated liposarcoma into two subtypes (S1 and S2), the unfavorable cluster was characterized as JUN amplified (An oncogene that promotes proliferation and metastasis) and lower inferred fraction of immature dendritic cells. However, the patients in ***Demicco et al., 2017*** were of complex origin (mixed limbs, trunk and retroperitoneum), providing limited guidance for RPLS molecular classification. Here we reported the first clinically applicable RPLS molecular classification based on RNA sequencing and IHC validation cohorts. To facilitate the clinical application, we constructed a simplified RPLS molecular classification derived from the original Cell cycle/Metabolism subgroups. By NMF algorithm, we identified LEP and PTTG1 as representative biomarkers for each subtype. A model based on IHC staining of LEP and PTTG1 successfully approximated the original dichotomous RPLS classification in biological features and survival outcomes. LEP is an important regulator of basal metabolism and food intake, which is considered a linkage between metabolism and the immune system (***Jiménez-Cortegana et al., 2021***). Although LEP-based targeting therapies have not yet been fully applied, LEP has already been identified as a potent metabolic reprogramming agent to support antitumor responses in aggressive melanomas (***Waldman et al., 2020; De la et al., 2019; Rivadeneira et al., 2019***). In addition, LEP improves the immunotherapeutic effects by regulating innate and adaptive immune responses via increasing the cytotoxicity of NK cells (***Vera et al., 2018***), stimulating the proliferation of T/B cells (***Vera et al., 2018; Bernotiene et al., 2006***), and activating DC cells (***Hu et al., 2019***). Those reported roles of LEP provided a good mechanism explanation on the features of metabolism pathway-enriched, better prognosis, higher TME level of Metabolism subgroup (LEP+). In contrast, PTTG1 acts as a regulator of sister chromatid separation during cell division under physiological conditions (Zou et al., 1999), which is closely linked to genetic instability, aneuploidy, tumor progression, invasion, and metastasis (***Heaney et al., 2000; Ramaswamy et al., 2003; Kim et al., 2005; Yu et al., 2003; Teveroni et al., 2021; Romero et al., 2001***). PTTG1 also regulates the cell cycle and the transactivation of growth factors as an initiator and promoter of tumorigenesis (***Zou et al., 1999; Mora-Santos et al., 2013; McCabe et al., 2002; Ishikawa et al 2001; Hamid et al., 2005***). Overexpressing PTTG1 was correlated with worse prognosis in tumors, such as ovarian cancer (***Parte et al., 2019***), cervical cancer (***Guo et al., 2019***), renal cell carcinoma (***Tian et al., 2022***), and colorectal cancer (***Heaney et al., 2000***). Therefore, the biological functions of PTTG1 provided a good mechanism explanation of the pathway-enriched of cell cycle/DNA damage & repair-associated, worse prognosis, and more aggressive pathological composition ratio the Cell cycle subgroup (PTTG1+). ## CONCLUSION Our study presented a comprehensive gene expression landscape of RPLS, revealing distinct molecular features. Through categorizing RPLS into Metabolism and Cell Cycle subtypes and identifying key biomarkers LEP and PTTG1, we established a dichotomous classification system verified by IHC assays. This innovative approach enables precise guidance for surgeons in adjusting treatment strategies for patients with histologically favorable but prognostically challenging RPLS cases, thereby advancing the implementation of precision medicine in guiding surgical interventions for RPLS. ## Authors’ contributors Mengmeng Xiao: Methodology, Formal analysis, Data curation, Writing-original draft, Writing-review & editing. Xiang Ji Li: Methodology, Formal analysis, Investigation, Data curation, Writing-original draft. Fanqin Bu: Conceptualization, Methodology, Formal analysis, Investigation, Data curation, Writing-review & editing. Shixiang Ma: Methodology, Formal analysis. Xiaohan Yang: Methodology, Investigation. Jun Chen: Formal analysis. Yu Zhao: Investigation. Ferdinando Cananzi: Investigation. Chenghua Luo: Conceptualization, Methodology, Writing-review & editing, Supervision. Li Min: Conceptualization, Resources, Writing-review & editing, Supervision. ## Ethics approval and consent to participate Specimens of RPLS were obtained from Peking University International Hospital. The study protocol was approved by the Ethics Committee of Peking University International Hospital, Peking University Health Science Center (WA2020RW29) and conducted in accordance with Helsinki Declaration. All patients signed the informed consent. ## Finding information This work was supported by grants from the Beijing Municipal Science and Technology Project (Z191100006619081), National Natural Science Foundation of China (82073390), and Young Elite Scientists Sponsorship Program (2023QNRC001). The study sponsors had no role in the design and preparation of this manuscript. ## Research registration unique identifying number (UIN) 1. Name of the registry: ClinicalTrials.gov 2. Unique identifying number or registration ID: [NCT03838718](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03838718&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 3. Hyperlink to your specific registration: [https://clinicaltrials.gov](https://clinicaltrials.gov) ## Consent for publication All authors have read and approved the manuscript and agree with submission to ***eLife***. ## Data availability Data supporting the conclusions of this article are presented within the article and its supplementary files. ## Competing interests The authors declare that they have no competing interests. ## Supplementary Information **Table S1.** The detailed clinicopathological characteristics of the training cohort 1 **Table S2.** The RNA-seq data of the training cohort 2 **Table S3.** The detailed clinicopathological characteristics of the validation cohort **Table S4.** The RNA-seq data of the training cohort 1 **Figure S1. Clinical features and re-classification of RPLS subgroups (G1 and G2)** The difference of Ki67 **(A)** and tumor size **(B)** in G1 and G2. Volcano plot of the DEGs (G1 vs G2) **(C)**. Venn diagram showed shared genes between DEGs and prognostic genes **(D)**. Bubble plot of the DEGs enrichment **(E)**. NMF for a re-classification of training cohort 1 **(F)**. tSNE exhibits the RPLS clusters (C1 and C2) **(G)**. **Figure S2. Dysregulated pathways and clinical features of RPLS clusters and high-/low-risk groups** The hierarchical clustering heatmap of dysregulated pathways in C1 and C2 **(A)**. The difference of Ki67 **(B)** and tumor size **(C)** in C1 and C2. The hierarchical clustering heatmap of dysregulated pathways in high- and low-risk groups **(D)**. The difference of tumor size **(E)** and Ki67 **(F)** in high- and low-risk groups. **Figure S3. Survival nomogram of LEP+PTTG1 model in validation cohort** Nomograms for DFS was developed in REASR cohort with four factors: sex, age, risk score, and differentiation **(A)**. ROC curves of 1-, 2-, and 3-year DFS in validation cohort **(B)**. Calibration curves of predicting 1-, 2-, and 3-year DFS **(C)**. ## Abbreviations RPLS : Retroperitoneal liposarcoma DEGs : Differentially expressed genes IHC : Immunohistochemistry STS : Soft tissue sarcoma OS : Overall survival DFS : Disease-free survival OMIX : Open archive for miscellaneous data CNCB : China national center for bioinformation GSEA : Gene set enrichment analysis GO : Gene ontology KEGG : Kyoto encyclopedia of gene and genomes t-SNE : t-distributed stochastic neighbor embedding DBSCAN : Density-based spatial clustering of applications with noise NMF : Nonnegative matrix factorization LR : Logistic regression SVM : Support vector machine RF : Random forest AUC : Area under curve * Received August 9, 2024. * Revision received August 9, 2024. * Accepted August 9, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## REFERENCES 1. Ecker BL, Peters MG, McMillan MT, Sinnamon AJ, Zhang PJ, Fraker DL, Levin WP, Roses RE, Karakousis GC. Preoperative radiotherapy in the management of retroperitoneal liposarcoma. Br J Surg. 2016; 103(13): 1839–1846. doi: 10.1002/bjs.10305 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/bjs.10305&link_type=DOI) 2. Littau MJ, Kulshrestha S, Bunn C, Agnew S, Sweigert P, Luchette FA, Baker MS. The importance of the margin of resection and radiotherapy in retroperitoneal liposarcoma. Am J Surg. 2021;221(3):554–560. doi:10.1016/j.amjsurg.2020.11.041 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amjsurg.2020.11.041&link_type=DOI) 3. Gronchi A, Miceli R, Allard MA, Callegaro D, Le Péchoux C, Fiore M, Honoré C, Sanfilippo R, Coppola S, Stacchiotti S, Terrier P, Casali PG, Cesne AL, Mariani L, Colombo C, Bonvalot S. Personalizing the approach to retroperitoneal soft tissue sarcoma: histology-specific patterns of failure and postrelapse outcome after primary extended resection. Ann Surg Oncol. 2015;22(5):1447–1454. doi:10.1245/s10434-014-4130-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1245/s10434-014-4130-7&link_type=DOI) 4. Gronchi A, Lo Vullo S, Fiore M, Mussi C, Stacchiotti S, Collini P, Lozza L, Pennacchioli E, Mariani L, Casali PG. Aggressive surgical policies in a retrospectively reviewed single-institution case series of retroperitoneal soft tissue sarcoma patients. J Clin Oncol. 2009;27(1):24–30. doi:10.1200/JCO.2008.17.8871 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjc6IjI3LzEvMjQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8wOS8yMDI0LjA4LjA5LjI0MzExNjU3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. Gronchi A, Pollock RE. Quality of local treatment or biology of the tumor: which are the trump cards for loco-regional control of retroperitoneal sarcoma?. Ann Surg Oncol. 2013;20(7):2111–2113. doi:10.1245/s10434-013-2971-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1245/s10434-013-2971-0&link_type=DOI) 6. Pisters PW. Resection of some -- but not all -- clinically uninvolved adjacent viscera as part of surgery for retroperitoneal soft tissue sarcomas. J Clin Oncol. 2009;27(1):6–8. doi:10.1200/JCO.2008.18.7138 [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjY6IjI3LzEvNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzA5LzIwMjQuMDguMDkuMjQzMTE2NTcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 7. Kam AE, Masood A, Shroff RT. Current and emerging therapies for advanced biliary tract cancers. Lancet Gastroenterol Hepatol. 2021;6(11):956–969. doi:10.1016/S2468-1253(21)00171-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-1253(21)00171-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34626563&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 8. Zeng Y, Jin RU. Molecular pathogenesis, targeted therapies, and future perspectives for gastric cancer. Semin Cancer Biol. 2022;86(Pt 3):566–582. doi:10.1016/j.semcancer.2021.12.004 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.semcancer.2021.12.004&link_type=DOI) 9. Alifrangis C, McGovern U, Freeman A, Powles T, Linch M. Molecular and histopathology directed therapy for advanced bladder cancer. Nat Rev Urol. 2019;16(8):465–483. doi:10.1038/s41585-019-0208-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41585-019-0208-0&link_type=DOI) 10. Frese KK, Simpson KL, Dive C. Small cell lung cancer enters the era of precision medicine. Cancer Cell. 2021;39(3):297–299. doi:10.1016/j.ccell.2021.02.002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ccell.2021.02.002&link_type=DOI) 11. Pilotti S, Della Torre G, Mezzelani A, Tamborini E, Azzarelli A, Sozzi G, Pierotti MA. The expression of MDM2/CDK4 gene product in the differential diagnosis of well differentiated liposarcoma and large deep-seated lipoma. Br J Cancer. 2000;82(7):1271–1275. doi:10.1054/bjoc.1999.1090 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1054/bjoc.1999.1090&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10755400&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085917800006&link_type=ISI) 12. Binh MB, Sastre-Garau X, Guillou L, Pinieux G, Terrier P, Lagacé R, Aurias A, Hostein I, Coindre JM. MDM2 and CDK4 immunostainings are useful adjuncts in diagnosing well-differentiated and dedifferentiated liposarcoma subtypes: a comparative analysis of 559 soft tissue neoplasms with genetic data. Am J Surg Pathol. 2005;29(10):1340–1347. doi:10.1097/01.pas.0000170343.09562.39 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/01.pas.0000170343.09562.39&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16160477&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000232120800011&link_type=ISI) 13. Yen CC, Chen SC, Hung GY, Wu PK, Chua WY, Lin YC, Yen CH, Chen YC, Wang JY, Yang MH, Chao Y, Chang MC, Chen WM. Expression profile-driven discovery of AURKA as a treatment target for liposarcoma. Int J Oncol. 2019;55(4):938–948. doi:10.3892/ijo.2019.4861 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3892/ijo.2019.4861&link_type=DOI) 14. Yang L, Wu Z, Sun W, Luo P, Chen S, Chen Y, Yan W, Li Y, Wang C. CCNDBP1, a Prognostic Marker Regulated by DNA Methylation, Inhibits Aggressive Behavior in Dedifferentiated Liposarcoma *via* Repressing Epithelial Mesenchymal Transition. Front Oncol. 2021;11:687012. doi:10.3389/fonc.2021.687012 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fonc.2021.687012&link_type=DOI) 15. McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM. Reporting recommendations for tumour MARKer prognostic studies (REMARK). Br J Cancer. 2005;93(4):387–391. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.bjc.6602678&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16106245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231260400002&link_type=ISI) 16. Guo L, Chen G, Zhang W, Zhou L, Xiao T, Di X, Wang Y, Feng L, Zhang K. A high-risk luminal A dominant breast cancer subtype with increased mobility. Breast Cancer Res Treat. 2019;175(2):459–472. doi:10.1007/s10549-019-05135-w [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10549-019-05135-w&link_type=DOI) 17. Li X, Bu F, Ma S, Cananzi F, Zhao Y, Xiao M, Min L, Luo C. The Janus-faced role of TRPM2-S in retroperitoneal liposarcoma via increasing ROS levels. Cell Commun Signal. 2022;20(1):128. doi:10.1186/s12964-022-00873-9 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12964-022-00873-9&link_type=DOI) 18. Lindskrog SV, Prip F, Lamy P, Taber A, Groeneveld CS, Birkenkamp-Demtröder K, Jensen JB, Strandgaard T, Nordentoft I, Christensen E, Sokac M, Birkbak NJ, Maretty L, Hermann GG, Petersen AC, Weyerer V, Grimm MO, Horstmann M, Sjödahl G, Höglund M, Steiniche T, Mogensen K, Reyniès A, Nawroth R, Jordan B, Lin X, Dragicevic D, Ward DG, Goel A, Hurst C, Kessel K, Maurer T, Meeks JJ, DeGraff DJ, Bryan RT, Knowles MA, Simic T, Hartmann A, Zwarthoff EC, Malmström PU, Malats N, Real FX, Dyrskjøt L. An integrated multi-omics analysis identifies prognostic molecular subtypes of non-muscle-invasive bladder cancer. Nat Commun. 2021;12(1):2301. 19. Yu Z, Deng P, Chen Y, Liu S, Chen J, Yang Z, Chen J, Fan X, Wang P, Cai Z, Wang Y, Hu P, Lin D, Xiao R, Zou Y, Huang Y, Yu Q, Lan P, Tan J, Wu X. Inhibition of the PLK1-Coupled Cell Cycle Machinery Overcomes Resistance to Oxaliplatin in Colorectal Cancer. Adv Sci (Weinh). 2021;8(23):e2100759. 20. Zhang F, Zhang Q, Zhu J, Yao B, Ma C, Qiao N, He S, Ye Z, Wang Y, Han R, Feng J, Wang Y, Qin Z, Ma Z, Li K, Zhang Y, Tian S, Chen Z, Tan S, Wu Y, Ran P, Wang Y, Ding C, Zhao Y. Integrated proteogenomic characterization across major histological types of pituitary neuroendocrine tumors. Cell Res. 2022;32(12):1047–1067. 21. Cancer Genome Atlas Research Network. Electronic address: elizabeth.demicco{at}sinaihealthsystem.ca; Cancer Genome Atlas Research Network. Comprehensive and Integrated Genomic Characterization of Adult Soft Tissue Sarcomas. Cell. 2017;171(4):950–965.e28. doi:10.1016/j.cell.2017.10.014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2017.10.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 22. Jiménez-Cortegana C, López-Saavedra A, Sánchez-Jiménez F, Pérez-Pérez A, Castiñeiras J, Virizuela-Echaburu JA, Cruz-Merino L, Sánchez-Margalet V. Leptin, Both Bad and Good Actor in Cancer. Biomolecules. 2021;11(6):913. doi:10.3390/biom11060913 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/biom11060913&link_type=DOI) 23. Waldman AD, Fritz JM, Lenardo MJ. A guide to cancer immunotherapy: from T cell basic science to clinical practice. Nat Rev Immunol. 2020;20(11):651–668. doi:10.1038/s41577-020-0306-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41577-020-0306-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32433532&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 24. De la Cruz-Merino L, Palazón-Carrión N, Henao-Carrasco F, Nogales-Fernández E, Álamo-de Ia Gala A, Vallejo-Benítez A, Chiesa M, Sánchez-Margalet V. New horizons in breast cancer: the promise of immunotherapy. Clin Transl Oncol. 2019;21(2):117–125. doi:10.1007/s12094-018-1907-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12094-018-1907-3&link_type=DOI) 25. Rivadeneira DB, DePeaux K, Wang Y, Kulkarni A, Tabib T, Menk AV, Sampath P, Lafyatis R, Ferris RL, Sarkar SN, Thorne SH, Delgoffe GM. Oncolytic Viruses Engineered to Enforce Leptin Expression Reprogram Tumor-Infiltrating T Cell Metabolism and Promote Tumor Clearance. Immunity. 2019;51(3):548–560.e4. doi:10.1016/j.immuni.2019.07.003 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2019.07.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 26. Vera F, Pino J, Campos-Cabaleiro V, Ruiz-Fernández C, Mera A, Gonzalez-Gay MA, Gómez R, Gualillo O. Obesity, Fat Mass and Immune System: Role for Leptin. Front Physiol. 2018;9:640. doi: 10.3389/fphys.2018.00640 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fphys.2018.00640&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29910742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 27. Bernotiene E, Palmer G, Gabay C. The role of leptin in innate and adaptive immune responses. Arthritis Res Ther. 2006;8(5):217. doi:10.1186/ar2004 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/ar2004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16879738&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 28. Hu W, Wang G, Huang D, Sui M, Xu Y. Cancer Immunotherapy Based on Natural Killer Cells: Current Progress and New Opportunities. Front Immunol. 2019;10:1205. doi:10.3389/fimmu.2019.01205 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2019.01205&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31214177&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 29. Zou H, McGarry TJ, Bernal T, Kirschner MW. Identification of a vertebrate sister-chromatid separation inhibitor involved in transformation and tumorigenesis. Science. 1999;285(5426):418–422. doi:10.1126/science.285.5426.418 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIyODUvNTQyNi80MTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8wOS8yMDI0LjA4LjA5LjI0MzExNjU3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 30. Heaney AP, Singson R, McCabe CJ, Nelson V, Nakashima M, Melmed S. Expression of pituitary-tumour transforming gene in colorectal tumours. Lancet. 2000;355(9205):716–719. doi:10.1016/S0140-6736(99)10238-1 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(99)10238-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10703804&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085615600014&link_type=ISI) 31. Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003;33(1):49–54. doi:10.1038/ng1060 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng1060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12469122&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180136100016&link_type=ISI) 32. Kim D, Pemberton H, Stratford AL, Buelaert K, Watkinson JC, Lopes V, Franklyn JA, McCabe CJ. Pituitary tumour transforming gene (PTTG) induces genetic instability in thyroid cells. Oncogene. 2005;24(30):4861–4866. doi:10.1038/sj.onc.1208659 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.onc.1208659&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15897900&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230477900013&link_type=ISI) 33. Yu R, Lu W, Chen J, McCabe CJ, Melmed S. Overexpressed pituitary tumor-transforming gene causes aneuploidy in live human cells. Endocrinology. 2003;144(11):4991–4998. doi:10.1210/en.2003-0305 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.2003-0305&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12960092&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186178300043&link_type=ISI) 34. Teveroni E, Di Nicuolo F, Bianchetti G, Epstein AL, Grande G, Maulucci G, Spirito MD, Pontecorvi A, Milardi D, Mancini F. Nuclear Localization of PTTG1 Promotes Migration and Invasion of Seminoma Tumor through Activation of MMP-2. Cancers (Basel). 2021;13(2):212. doi:10.3390/cancers13020212 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/cancers13020212&link_type=DOI) 35. Romero F, Multon MC, Ramos-Morales F, Domínguez A, Bernal JA, Pintor-Toro JA, Tortolero M. Human securin, hPTTG, is associated with Ku heterodimer, the regulatory subunit of the DNA-dependent protein kinase. Nucleic Acids Res. 2001;29(6):1300–1307. doi:10.1093/nar/29.6.1300 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/29.6.1300&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11238996&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167529200006&link_type=ISI) 36. Mora-Santos M, Castilla C, Herrero-Ruiz J, Giráldez S, Limón-Mortés MC, Sáez C, Japón MÁ, Tortolero M, Romero F. A single mutation in Securin induces chromosomal instability and enhances cell invasion. Eur J Cancer. 2013;49(2):500–510. doi:10.1016/j.ejca.2012.06.024 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejca.2012.06.024&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22819078&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 37. McCabe CJ, Boelaert K, Tannahill LA, Heaney AP, Stratford AL, Khaira JS, Hussain S, Sheppard MC, Franklyn JA, Gittoes NJL. Vascular endothelial growth factor, its receptor KDR/Flk-1, and pituitary tumor transforming gene in pituitary tumors. J Clin Endocrinol Metab. 2002;87(9):4238–4244. doi:10.1210/jc.2002-020309 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2002-020309&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12213878&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178201200037&link_type=ISI) 38. Ishikawa H, Heaney AP, Yu R, Horwitz GA, Melmed S. Human pituitary tumor-transforming gene induces angiogenesis. J Clin Endocrinol Metab. 2001;86(2):867–874. doi:10.1210/jcem.86.2.7184 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.86.2.867&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11158059&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167004900069&link_type=ISI) 39. Hamid T, Malik MT, Kakar SS. Ectopic expression of PTTG1/securin promotes tumorigenesis in human embryonic kidney cells. Mol Cancer. 2005;4(1):3. Published 2005 Jan 13. doi:10.1186/1476-4598-4-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1476-4598-4-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15649325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F09%2F2024.08.09.24311657.atom) 40. Parte S, Virant-Klun I, Patankar M, Batra SK, Straughn A, Kakar SS. PTTG1: a Unique Regulator of Stem/Cancer Stem Cells in the Ovary and Ovarian Cancer. Stem Cell Rev Rep. 2019;15(6):866–879. doi:10.1007/s12015-019-09911-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12015-019-09911-5&link_type=DOI) 41. Guo XC, Li L, Gao ZH, Zhou HW, Li J, Wang QQ. The long non-coding RNA PTTG3P promotes growth and metastasis of cervical cancer through PTTG1. Aging (Albany NY). 2019;11(5):1333–1341. doi:10.18632/aging.101830 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18632/aging.101830&link_type=DOI) 42. Tian X, Xu WH, Xu FJ, Li H, Anwaier A, Wang HK, Wan FN, Zhu Y, Cao DL, Zhu YP, Shi GH, Qu YY, Zhang HL, Ye DW. Identification of prognostic biomarkers in papillary renal cell carcinoma and PTTG1 may serve as a biomarker for predicting immunotherapy response. Ann Med. 2022;54(1):211–226.