Original Research Open Access Logo

Elevated expression of SKP2 correlates with poor prognosis in osteosarcoma: a bioinformatics analysis

Tien-Manh Hoang 1 ORCID logo
Minh-Tien Nguyen 1, 2
Weisin Chen 1
Chenyang Zhuang 1
Zixiang Wang 1
Hanquan Wang 1
Juan Li 1, 3
Hong Lin 1, 3, *
  1. Department of Orthopedics, Zhongshan Hospital, Fudan University, Shanghai, China
  2. Department of Spinal Surgery, Institute of Trauma and Orthopedics, 108 Military Central Hospital, Hanoi, Viet Nam
  3. Department of Orthopedics, Shanghai Geriatrics Center, Fudan University, Shanghai, China
Correspondence to: Hong Lin, Department of Orthopedics, Zhongshan Hospital, Fudan University, Shanghai, China; Department of Orthopedics, Shanghai Geriatrics Center, Fudan University, Shanghai, China. Email: [email protected].
Volume & Issue: Vol. 8 No. 12 (2021) | Page No.: 4782-4792 | DOI: 10.15419/bmrat.v8i12.714
Published: 2021-12-31

Online metrics


Statistics from the website

  • HTML Views: 6767
  • PDF Views: 1406
  • XML Views: 0
  • PDF Views: 1472

Statistics from Dimensions

Copyright The Author(s) 2024. This article is published with open access by BioMedPress. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0) which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. 

Abstract

Introduction: Osteosarcoma (OS) is the most prevalent kind of bone cancer, but the tumorigenesis and underlying molecular drivers of OS remain unknown. The instability of the cell cycle regulation system, which leads to uncontrolled cell proliferation, is a typical characteristic of carcinogenesis. This study aimed to investigate the expression of SKP2 in human OS and assessed its prognostic value in OS patients.

Methods: Three gene expression profile datasets GSE28424, GSE42352 and GSE21257 were obtained from the GEO database. The Wilcoxon rank-sum test was performed to evaluate the differential expression of the SKP2 gene between the OS and the control groups in the GSE28424 and GSE42352 datasets, and the correlation between SKP2 expression and Huvos grade was analysed in the GSE21257 dataset. Subsequently, PPI network, GO and KEGG analyses were constructed. Furthermore, survival analysis between high SKP2 expression group and low SKP2 expression group was performed using the Kaplan–Meier method and log-rank test. In addition, the UALCAN database was used to investigate the association between SKP2 expression and sarcoma.

Results: The expression level of SKP2 in OS cells was significantly higher than that in normal bones and mesenchymal stem cell samples. Furthermore, the level of SKP2 expression was observed to decrease as Huvos grade increased. The PPI network was established, and the top ten SKP2-related genes were identified, including CDKs (1, 2, 4, and 6), Cyclin A1-2, E1-2, D1, and CDKN1A. The survival analysis showed that the elevated SKP2 expression level was significantly related to the overall survival of OS patients.

Conclusion: Our work adds to our knowledge of SKP2's function in OS and suggests that it might be used as a therapeutic target in the future.

Introduction

Osteosarcoma (OS) is the most prevalent kind of bone cancer, mostly happening in young individuals aged 10 – 25 years and comprising over one-third of all primary bone malignancies1, 2, 3. The treatment for OS patients has progressed remarkably during the last four decades, and the advances in therapy, especially the introduction of chemotherapy as well as neoadjuvant chemotherapy followed by primary tumour resection, have led to a dramatically amelioration in the long-term prognosis of individuals without distal metastasis3, 4, 5. However, the prognosis for advanced patients remains poor, and early diagnosis has been restricted due to our limited understanding of the aetiology. Besides, the tumorigenesis and underlying molecular drivers of OS remain unknown. Thus, identifying novel effective biomarkers and drug targets is the key to understanding the development and progression of OS and to improve patient survival rates.

The instability of the cell cycle regulation system, which leads to uncontrolled cell proliferation, is a typical characteristic of carcinogenesis, and a variety of regulatory factors in the cell cycle are degraded through the ubiquitin-proteasome pathway. SKP2 (S-Phase Kinase-Associated Protein 2), also named FBXL1 (F-Box/LRR-Repeat Protein 1) or p45 (CDK2/cyclin A-associated protein p45), encodes for a 45 kDa protein that belongs to the F-box family, it could ubiquitinate and degrade numerous tumour suppressor proteins such as p27, and so plays a crucial role in regulating the cell cycle. SKP2 was initially identified as a pro-tumour factor in 1995, when it was discovered to be a crucial part for the S-phase of the cyclin A — CDK2 kinase in several transformed cells6. It was found that SKP2 was highly expressed and played a major role in carcinogenesis in numerous kinds of malignancies, including soft tissue sarcomas (STS)7, breast cancer8, hepatocellular carcinoma (HCC)9, pancreatic ductal adenocarcinoma (PDAC)10, astrocytic gliomas11. To our knowledge, no bioinformatics study has been conducted to investigate the involvement of SKP2 in OS. Hence, in the present study, we aimed to investigate the expression of SKP2 in human OS cell lines and samples and assessed its prognostic value in OS patients.

Methods

Microarray data

The Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) was used in this study. GEO, which is established and managed by the National Center for Biotechnology Information, National Library of Medicine, is an invaluable international resource that preserves and freely accesses comprehensive sets of microarray data submitted by researchers globally12, 13. The database includes data from a wide variety of technologies, such as Serial Analysis of Gene Expression, DNA microarray, protein array, high-throughput nucleic acid sequencing, and Real-Time PCR14, 15.

The search strategy we employed to identify OS-related datasets in the GEO database was as follows: “(osteosarcoma OR OS) AND Expression profiling by array[DataSet Type] AND Homo sapiens[Organism]”. Datasets GSE2842416, GSE4235217 and GSE2125718 were eventually included and extracted after a thorough search.

The platform for GSE28424 is GPL13376 — Illumina HumanWG-6 v2.0 expression beadchip, which contains nineteen OS cell lines and four normal bones. The platform for GSE42352 is GPL10295 — Illumina human-6 v2.0 expression beadchip, which comprises nineteen OS samples and twelve mesenchymal stem cell (MSC) samples used as controls. The GSE21257 dataset, which is also based on the GPL10295 platform, has fifty-three OS samples with complete clinicopathological characteristics.

Data processing

The R software (version 4.0.5) was used to acquire and process the expression profiling data. First, the raw data of all datasets was obtained using the GEOquery package (v2.58.0) and normalized and standardized using the limma package (v3.46.0). Then, the Wilcoxon rank-sum test was performed to evaluate the differential expression of SKP2 gene between the OS and the control groups in the GSE28424 and GSE42352 datasets. In the GSE21257 dataset, the samples were divided into a high expression group and low expression group based on the median value of SKP2 expression, and survival analysis was performed using the Kaplan–Meier method and log-rank test.

UALCAN database

UALCAN (http://ualcan.path.uab.edu) is a web-portal with user-friendly interactive features for conducting in-depth analysis of gene expression data from The Cancer Genome Atlas (TCGA)19. The mRNA expression of SKP2 in sarcoma tissues and adjacent normal tissues were explored. In addition, we also investigated the effect of SKP2 expression level on sarcoma patient survival through the UALCAN database by comparing the overall survival between the SKP2 high expression group and the SKP2 low/medium expression group.

Protein-Protein Interaction (PPI) network construction

The STRING database, which is a tool for retrieving and displaying the genes that a query gene appears within clusters on the genome20, was employed to establish the PPI network of SKP2 protein and other proteins. The PPI network was then visualized using the Cytoscape program21. Furthermore, Cytohubba, a Cytoscape plug-in, was performed to detect the hub genes that are significantly related to SKP2.

Functional and pathway enrichment analysis

Gene Ontology (GO) analysis aims to standardize the representation of gene and gene product properties across species, comprising cell components (CC), biological process (BP) and molecular function (MF)22. Kyoto Encyclopaedia of Genes and Genomes (KEGG) analysis serves as a combined database resource for biologically interpreting totally sequenced genomes, the process to map genes in the genome to generated pathway maps23.

GO and KEGG were utilized to probe the molecular mechanisms and pathways related to SKP2 for the development and prognosis of OS, and the visualization of the enrichment analyses was conducted through the cluster-Profiler package (v3.18.1), the GOplot package (v1.0.2), and the enrichplot package (v1.10.2) in the R language software.

Figure 1

Expression of SKP2 mRNA in OS and control cells. A. GSE28424 dataset, OS cell lines and normal bone cells; B. GSE42352 dataset, OS samples and MSC samples; C. TCGA database, sarcoma tissues and the corresponding normal tissues.

Results

Increased SKP2 expression in OS

The expression of SKP2 in OS cells and control cells was compared. In GSE28424 (Figure 1A), the expression level of SKP2 in OS cell lines was significantly higher than that in normal bones (P < 0.001). We also found that SKP2 expressed highly in OS samples compared to MSC samples in GSE42352 (P = 0.0014), as shown in Figure 1B. Besides, the expression difference of SKP2 between the tumour tissues of sarcoma and the corresponding normal tissues in the TCGA project was additionally evaluated through the UALCAN web portal (Figure 1C), and it was statistically significant (P < 0.001). These findings indicate that SKP2 is expressed at an obviously high level in OS cells.

Figure 2

SKP2 expression in osteosarcoma with different Huvos pathological grades in the GSE21257 dataset.

Relationship between SKP2 expression and Huvos grade

The Huvos grading system is widely used in the chemotherapy assessment of OS, and it is an important parameter for predicting long-term prognosis. Grades I and II are considered to be poor in response to chemotherapy, while grades III and IV are regarded to have a favourable response. We used the GSE21257 dataset to perform a correlation analysis between SKP2 expression and Huvos grade. The level of SKP2 expression was observed to decrease as Huvos grade increased (P < 0.05) (Figure 2), indicating that the expression level of SKP2 falls as the rate of tumour necrosis tissue increases the following chemotherapy. This suggests that SKP2 might be used to predict the effectiveness of OS chemotherapy.

Figure 3

Protein-ProteinInteraction (PPI) network. A. PPI network of SKP2 and the other related proteins; B. the top ten SKP2-related genes using cyto Hubba plug-in.

Figure 4

Functional enrichment analysis. A. Biological processes (BP) analysis; B. Cell composition (CC) analysis; C. Molecular function (MF)analysis; D. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis.

Table 1

Top 5 results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses

GO/KEGG

ID

Description

BgRatio

P value

Adjusted P value

q value

Count

BP

GO:0044843

cell cycle G1/S phase transition

310/18866

1.98E-20

1.11E-17

5.05E-18

11

BP

GO:0000082

G1/S transition of mitotic cell cycle

287/18866

6.16E-18

1.73E-15

7.88E-16

10

BP

GO:0000079

regulation of cyclin-dependent protein serine/threonine kinase activity

102/18866

3.55E-14

6.56E-12

2.99E-12

7

BP

GO:1904029

regulation of cyclin-dependent protein kinase activity

106/18866

4.68E-14

6.56E-12

2.99E-12

7

BP

GO:0000086

G2/M transition of mitotic cell cycle

254/18866

2.33E-11

2.61E-09

1.19E-09

7

CC

GO:0000307

cyclin-dependent protein kinase holoenzyme complex

43/19559

9.35E-27

2.99E-25

1.08E-25

10

CC

GO:1902554

serine/threonine protein kinase complex

89/19559

2.47E-23

3.96E-22

1.43E-22

10

CC

GO:1902911

protein kinase complex

104/19559

1.27E-22

1.35E-21

4.90E-22

10

CC

GO:0061695

transferase complex, transferring phosphorus-containing groups

253/19559

1.19E-18

9.55E-18

3.46E-18

10

CC

GO:0016592

mediator complex

41/19559

0.000233

0.001491

0.00054

2

MF

GO:0016538

cyclin-dependent protein serine/threonine kinase regulator activity

50/18352

2.35E-16

7.29E-15

3.22E-15

7

MF

GO:0019887

protein kinase regulator activity

185/18352

3.01E-12

3.90E-11

1.72E-11

7

MF

GO:0030332

cyclin binding

30/18352

3.77E-12

3.90E-11

1.72E-11

5

MF

GO:0019207

kinase regulator activity

216/18352

9.00E-12

6.98E-11

3.08E-11

7

MF

GO:0004693

cyclin-dependent protein serine/threonine kinase activity

29/18352

1.65E-09

8.51E-09

3.75E-09

4

KEGG

hsa04110

Cell cycle

126/8095

8.35E-21

5.18E-19

1.76E-19

11

KEGG

hsa05203

Viral carcinogenesis

204/8095

1.99E-18

6.18E-17

2.10E-17

11

KEGG

hsa04218

Cellular senescence

156/8095

5.73E-17

1.18E-15

4.02E-16

10

KEGG

hsa05169

Epstein-Barr virus infection

202/8095

8.08E-16

1.25E-14

4.25E-15

10

KEGG

hsa04115

p53 signaling pathway

73/8095

4.76E-15

5.91E-14

2.01E-14

8

Figure 5

Kaplan–Meier analysis for overall survival based on SKP2 expression. A. Between high expression group and low expression group in OS. B. Between high expression group and low/medium expression group insarcoma.

PPI network construction

The PPI network created by STRING, which consisted of 51 nodes and 709 edges, was visualized using Cytoscape to explore the functional connections between SKP2 and the other proteins related to it (Figure 3A). Subsequently, as shown in Figure 3B, the top ten related genes identified using the Maximal Clique Centrality method in the Cytohubba plug-in were CDKs (1, 2, 4, and 6), Cyclin A1-2, E1-2, D1, and CDKN1A.

GO and KEGG enrichment analyses

Enrichment analyses of GO and KEGG pathway for SKP2 and its top ten associated genes were performed and displayed using the R software packages clusterProfiler, GOplot, and enrichplot (Figure 4, Table 1 and Table S1). According to the GO enrichment results, regarding BP, the genes were significantly enriched in cell cycle G1/S phase transition, G1/S transition of mitotic cell cycle, regulation of cyclin−dependent protein serine/threonine kinase activity, regulation of cyclin−dependent protein kinase activity, and G2/M transition of mitotic cell cycle (Figure 4A). For CC, the genes were significantly enriched in cyclin−dependent protein kinase holoenzyme complex, serine/threonine protein kinase complex, protein kinase complex, and transferase complex, transferring phosphorus−containing groups (Figure 4B). Regarding MF, the genes were significantly enriched in cyclin−dependent protein serine/threonine kinase regulator activity, protein kinase regulator activity, cyclin binding, and kinase regulator activity (Figure 4C). The findings of the KEGG pathway analysis, as shown in Figure 4D, revealed that the genes were mainly engaged in the cell cycle, viral carcinogenesis, and cellular senescence.

Relationship between SKP2 expression and prognosis in OS patients

We further studied the role of SKP2 in the prognosis of OS patients. In the GSE21257 dataset, based on the cut-off value, patients were split into two groups: low expression group (26 patients) and high expression group (27 patients). The Kaplan-Meier curve and log-rank test shown that the elevated SKP2 expression level was significantly related to the overall survival (P < 0.01) of OS patients (Figure 5A). Moreover, the dataset in UALCAN also indicated that there was a significant difference in survival time between the SKP2 high expression group (65 patients) and the SKP2 medium/low expression group (194 patients) (P = 0.017), as shown in Figure 5B. Taken together, OS patients with a high level of SKP2 expression were predicted to have a poor prognosis.

Discussion

OS is one of the most prevalent bone cancers, affecting primarily adolescents and children and has a 5-year overall survival rate of approximately 50% — 70%1. The histology of OS is very variable, and its genetic structure is unstable, resulting in poor survival for the patients24. Hence, new prognostic biomarkers are required to assist patients with OS to improve their outcomes.

There have been several bioinformatics studies on OS in the past. HLA-DRA, EGR1, CXCL10, MYC, and CXCR4 were discovered to be hub normal-primary differentially expressed genes (DEGs) in OS25, 26. Besides, Zhang .27measured the amounts of tumour-infiltrating immune cells in the OS tumour microenvironment via the ESTIMATE algorithm-based immune score, and they detected 191 immune-related DEGs, with PPARG, IGHG3, PDK, CD209, and CCL8 ranking first through fifth. Another research studied the differential expression of circRNA in OS, they determined that circ_20403 and circ_2137 were highly upregulated and circ_24831 and circ_32279 were highly downregulated in OS samples, these circRNAs were significantly associated with BP and MF28.

SKP2 encodes for a type of Leucine-rich repeat protein that belongs to the F-box family29. SKP2 was initially recognized as an important cell cycle regulator because the ubiquitin-mediated degradation of several regulators of cell cycle required the promotion of SKP230, 31. Owing to the fact that the majority of its ubiquitinated and degraded substrates are tumour suppressors, SKP2 is classified as an oncoprotein. Besides, many studies also suggested that SKP2 functions as an oncogene in human malignancies.

In a large-scale analysis, elevated SKP2 expression in patients with STS was demonstrated as an independent risk factor for disease-specific survival in females and those who did not undergo chemo- or radiotherapy32. Another research showed that medium-to-high SKP2 expression in the nucleus in ≥ 10% of the cells was recognized in 37% of the total cases of STS and was more frequent in the high-grade group than in low-grade group7. Regarding breast cancer, cytoplasmic over-expression of SKP2 is presented in invasive ductal carcinoma, positively correlated with tumour grade and tumour size and inversely correlated with patient outcomes33. In HCC, Lee . revealed that LKB1 polyubiquitination mediated by Skp2 was crucial for tumour development in vivo34; in addition, enhanced SKP2 expression was shown to be negatively associated with apoptotic process and prognosis and positively associated with the increased level of tumour cells and microvascularization9. About PDAC, SKP2 tended to overexpress in patients with higher tumour grade, more lymph node metastases, greater amount of lymphatic permeation, and worse survival10.

In this study, the role of SKP2 in OS was investigated using a bioinformatics approach. According to our findings, all three datasets GSE28424, GSE42352 and TCGA determined that SKP2 was upregulated in OS, this result was consistent with prior research. When compared to normal human osteoblasts or human MSC-derived osteoblasts, Zhang .35 found that the levels of Skp2 mRNA were significantly higher in both conventional and patient-derived OS cell lines. In OS cells, elevated expression of SKP2 was determined to stimulate cell proliferation and cell cycle progression, reduce apoptosis, and accelerate cell invasion and migration36. Notably, Nagao .37 revealed that overexpression of SKP2 induced overexpression of GLI2 in OS cells, resulting in enhancing cell proliferation and promoting cell cycle progression, showing that Skp2 was an important link in the control of cell growth in OS cells. Another study discovered that SKP2 interacts with its substrate p27 to promote OS development and stemness38. Furthermore, SKP2 might be used as a potential prognostic and pro-metastatic biomarker. The findings of our survival study determined that SKP2 was tightly related to prognosis in OS patients, with greater levels of SKP2 expression having a worse prognosis. Similar to this result, in a study of metastasis-free survival analysis in 88 untreated high-grade OS patients, the high Skp2 group had a poorer prognosis than the low Skp2 group35. In addition, Liu .39 also demonstrated that OS patients with lung metastasis with high SKP2 TPM (Transcripts Per Million) had a worse survival than those with low SKP2 TPM.

Moreover, on the basis of established criteria, the top ten SKP2-related genes were selected as hub genes for further investigation, which were CDK1-2, 4, and 6, Cyclin A1-2, D1, E1-2, and CDKN1A. These genes are all closely related to the cell cycle, which is a critical process in malignant transformation. Disturbances in the cell cycle are frequently documented in tumorigenesis40. A vast body of literature has shown CDK (cyclin-dependent kinase) family members exert as oncoproteins and are deregulated in various malignancies, changes in CDK activity are frequently involved in tumour-associated cell cycle abnormalities41. CDK6 overexpression may promote OS cell proliferation and migration42, whereas SKP2 overexpression could increase the level of CDK6 through its substrate p27. Regarding CDK2, suppressing or downregulating it in OS cell lines could stop them from proliferating43, 44. Cyclin family regulates the cell progression through triggering CDK enzymes or synthesis-related enzymes for cell cycle45. Various prior reports have demonstrated that many members of the cyclin family, including cyclin A, D, and E, are significantly expressed and/or involved in cell proliferation in OS46, 47, 48, 49, 50, 51. Li . discovered that hsa_circ_0003732 may enhance OS cell proliferation by increasing cyclin A2 expression through miR-54552. Besides, numerous biomarkers, such as URG4, JARID1B, lncRNA LINC01296, promote OS cell proliferation, metastasis and cell cycle via cyclin D147, 48, 53. Additionally, cyclin E1 overexpression was identified as associated with clinicopathology and a potential marker predicting outcome for OS patients50.

The Huvos system was used to observe the effect of chemotherapy on high-grade OS in the GSE21257 dataset and in the present study. The Huvos grading system is a histologic evaluation system for bone sarcomas in which the level of necrosis in relation to the percentage of remaining viable tumour is determined for histologic assessment purposes using a semiquantitative method54. Grades 1 – 4 indicate necrosis of < 50%, 50% — 90%, 90% — 99%, and 100%, respectively, and higher grades are regarded to have a favourable response to chemotherapy. Our finding showed that individuals with high grades had a significantly low expression of SKP2 than those with low grades, suggesting that chemotherapy may exert its effect through inhibition of SKP2 expression. Given the crucial participation in regulating cancer progression, we believed that SKP2 might be a new target for the therapeutic treatment of OS. A prior report showed that downregulating the expression of SKP2 by PPARG could effectively induce impaired cell proliferation and apoptosis, suggesting that targeting SKP2 expression with PPARG agonists could be a viable treatment approach55. Hirotsu .56 discovered that suppressing SMO with cyclopamine promoted OS cell cycle arrest as well as prevented cell development by decreasing the transcription of Skp2. Moreover, upregulation of miR-506 induces Skp2 to be reduced, which stimulates apoptosis and slows motility in OS cells57. Interestingly, both SKP2 inhibitors (compound C1, pevonedistat and FKA) and SKP1-SKP2 pocket inhibitor (compound C25) have anti-proliferative effects on pRB/p53 double-deficient OS cells38. Conclusively, developing SKP2 inhibitors is a feasible and promising method for treating OS.

Conclusions

In conclusion, the expression and prognostic value of SKP2 in OS were formulated in the present study. SKP2 was overexpressed in OS cell lines and samples and correlated with poor outcomes in OS patients. Given its importance in the cell cycle, our research suggests that SKP2 might be an important biomarker in OS tumorigenesis as well as a potential target for drug development and treatment.

Abbreviations

BP: Biological process, CC: Cell components, CDK: Cyclin-dependent kinase, DEG: Differentially expressed gene, GEO: Gene Expression Omnibus, GO: Gene Ontology, HCC: Hepatocellular carcinoma, KEGG: Kyoto Encyclopedia of Genes and Genomes, MF: Molecular function, MSC: Mesenchymal stem cell, OS: Osteosarcoma, PDAC: Pancreatic ductal adenocarcinoma, PPI: Protein-Protein Interaction, SKP2: S-Phase Kinase-Associated Protein 2, STS: Soft tissue sarcomas, TCGA: The Cancer Genome Atlas, TPM: Transcripts Per Million

Acknowledgments

None.

Author’s contributions

Juan Li and Hong Lin, conception and design; Hong Lin, administrative support; Tien-Manh Hoang, Minh-Tien Nguyen and Weisin Chen, collection and assembly of data; Tien-Manh Hoang, Minh-Tien Nguyen, Chenyang Zhuang, Hanquan Wang, data analysis and interpretation. All authors wrote the article. All authors read and approved the final manuscript.

Funding

This work was supported by the Non-profit Central research Institute Fund of Chinese academy of Medical sciences [2020-JKCS-013]; and Three-Year Initiative Plan for Strengthening Public Health System Construction in Shanghai [GWV-10.2-XD11].

Availability of data and materials

All data were obtained from the GEO database (http://www.ncbi.nlm.nih.gov/geo/). The authors will also provide researchers with the full dataset upon request.

Ethics approval and consent to participate

The study was exempted from institutional review board approval since the clinical data was obtained from a publicly accessible database.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References

  1. L. Mirabello, R.J. Troisi, S.A. Savage. Osteosarcoma incidence and survival rates from 1973 to 2004: data from the Surveillance, Epidemiology, and End Results Program. Cancer 2009; 115 (7) : 1531-1543.
  2. M.U. Jawad, M.C. Cheung, J. Clarke, L.G. Koniaris, S.P. Scully. Osteosarcoma: improvement in survival limited to high-grade patients only. Journal of Cancer Research and Clinical Oncology (ISSN: 1432-1335). 2011; 137 (4) : 597-607.
  3. M.S. Isakoff, S.S. Bielack, P. Meltzer, R. Gorlick. Osteosarcoma: Current Treatment and a Collaborative Pathway to Success. J Clin Oncol 2015; 33 (27) : 3029-3035.
  4. J. Whelan, B. Seddon, M. Perisoglou. Management of osteosarcoma. Current Treatment Options in Oncology (ISSN: 1527-2729). 2006; 7 (6) : 444-55.
  5. W. Zhu, L. Zhu, Y. Bao, X. Zhong, Y. Chen, Q. Wu. Clinical evaluation of neoadjuvant chemotherapy for osteosarcoma. Journal of the Balkan Union of Oncology (ISSN: 2241-6293). 2019; 24 (3) : 1181-5.
  6. H. Zhang, R. Kobayashi, K. Galaktionov, D. Beach. p19Skp1 and p45Skp2 are essential elements of the cyclin A-CDK2 S phase kinase. Cell (ISSN: 0092-8674). 1995; 82 (6) : 915-25.
  7. A.M. Oliveira, S.H. Okuno, A.G. Nascimento, R.V. Lloyd. Skp2 protein expression in soft tissue sarcomas. Journal of Clinical Oncology (ISSN: 0732-183X). 2003; 21 (4) : 722-7.
  8. C. Li, L. Du, Y. Ren, X. Liu, Q. Jiao, D. Cui. SKP2 promotes breast cancer tumorigenesis and radiation tolerance through PDCD4 ubiquitination. J Exp Clin Cancer Res 2019; 38 (1) : 76-
  9. D.F. Calvisi, S. Ladu, F. Pinna, M. Frau, M.L. Tomasi, M. Sini. SKP2 and CKS1 promote degradation of cell cycle regulators and are associated with hepatocellular carcinoma prognosis. Gastroenterology 2009; 137 (5) : 1816-1826.e1811-1810.
  10. T. Einama, Y. Kagata, H. Tsuda, D. Morita, S. Ogata, S. Ueda. High-level Skp2 expression in pancreatic ductal adenocarcinoma: correlation with the extent of lymph node metastasis, higher histological grade, and poorer patient outcome. Pancreas (ISSN: 1536-4828). 2006; 32 (4) : 376-81.
  11. D. Schiffer, P. Cavalla, V. Fiano, C. Ghimenti, R. Piva. Inverse relationship between p27/Kip.1 and the F-box protein Skp2 in human astrocytic gliomas by immunohistochemistry and Western blot. Neuroscience Letters (ISSN: 0304-3940). 2002; 328 (2) : 125-8.
  12. T. Barrett, S.E. Wilhite, P. Ledoux, C. Evangelista, I.F. Kim, M. Tomashevsky. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res 2013; 41 (Database issue) : D991-995.
  13. R. Edgar, M. Domrachev, A.E. Lash. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 2002; 30 (1) : 207-210.
  14. E. Clough, T. Barrett. The Gene Expression Omnibus Database. Methods Mol Biol 2016; (1418) : 93-110.
  15. D. Toro-Domínguez, J. Martorell-Marugán, R. López-Domínguez, A. García-Moreno, V. González-Rumayor, M.E. Alarcón-Riquelme. ImaGEO: integrative gene expression meta-analysis from GEO database. Bioinformatics (Oxford, England) (ISSN: 1367-4811). 2019; 35 (5) : 880-2.
  16. H.M. Naml∅s, L.A. Meza-Zepeda, T. Bar∅y, I.H. Ostensen, S.H. Kresse, M.L. Kuijjer. Modulation of the osteosarcoma expression phenotype by microRNAs. PLoS One 2012; 7 (10) : e48086-
  17. M.L. Kuijjer, E.F. Peterse, B.E. van den Akker, I.H. Briaire-de Bruijn, M. Serra, L.A. Meza-Zepeda. IR/IGF1R signaling as potential target for treatment of high-grade osteosarcoma. BMC Cancer 2013; 13245-
  18. E.P. Buddingh, M.L. Kuijjer, R.A. Duim, H. Bürger, K. Agelopoulos, O. Myklebost. Tumor-infiltrating macrophages are associated with metastasis suppression in high-grade osteosarcoma: a rationale for treatment with macrophage activating agents. Clinical Cancer Research (ISSN: 1557-3265). 2011; 17 (8) : 2110-9.
  19. D.S. Chandrashekar, B. Bashel, S.A.H. Balasubramanya, C.J. Creighton, I. Ponce-Rodriguez, B. Chakravarthi. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017; 19 (8) : 649-658.
  20. B. Snel, G. Lehmann, P. Bork, M.A. Huynen. STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res 2000; 28 (18) : 3442-3444.
  21. S.M. Menon, A. Elengoe. Evaluation of the role of kras gene in colon cancer pathway using string and Cytoscape software. Biomedical Research and Therapy (ISSN: 2198-4093). 2020; 7 (6) : 3835-42.
  22. . Gene Ontology Consortium. The Gene Ontology in 2010: extensions and refinements. Nucleic Acids Res 2010; 38 (suppl_1) : D331-335.
  23. M. Kanehisa, M. Furumichi, M. Tanabe, Y. Sato, K. Morishima. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res 2017; 45 (D1) : D353-d361.
  24. A.J. Mutsaers, C.R. Walkley. Cells of origin in osteosarcoma: mesenchymal stem cells or osteoblast committed cells?. Bone (ISSN: 1873-2763). 2014; 6256-63.
  25. J. Liu, S. Wu, X. Xie, Z. Wang, Q. Lei. Identification of potential crucial genes and key pathways in osteosarcoma. Hereditas 2020; 157 (1) : 29-
  26. J. Niu, T. Yan, W. Guo, W. Wang, Z. Zhao, T. Ren. Identification of Potential Therapeutic Targets and Immune Cell Infiltration Characteristics in Osteosarcoma Using Bioinformatics Strategy. Front Oncol 2020; 101628-
  27. C. Zhang, J.H. Zheng, Z.H. Lin, H.Y. Lv, Z.M. Ye, Y.P. Chen. Profiles of immune cell infiltration and immune-related genes in the tumor microenvironment of osteosarcoma. Aging (Albany NY) 2020; 12 (4) : 3486-3501.
  28. Y. Xi, M. Fowdur, Y. Liu, H. Wu, M. He, J. Zhao. Differential expression and bioinformatics analysis of circRNA in osteosarcoma. Biosci Rep 2019; 39 (5) : BSR20181514-
  29. Z. Wang, P. Liu, H. Inuzuka, W. Wei. Roles of F-box proteins in cancer. Nat Rev Cancer 2014; 14 (4) : 233-247.
  30. Z.K. Yu, J.L. Gervais, H. Zhang. Human CUL-1 associates with the SKP1/SKP2 complex and regulates p21(CIP1/WAF1) and cyclin D proteins. Proc Natl Acad Sci U S A 1998; 95 (19) : 11324-11329.
  31. A.C. Carrano, E. Eytan, A. Hershko, M. Pagano. SKP2 is required for ubiquitin-mediated degradation of the CDK inhibitor p27. Nature Cell Biology (ISSN: 1465-7392). 1999; 1 (4) : 193-9.
  32. S.W. Sorbye, T.K. Kilvaer, A. Valkov, T. Donnem, E. Smeland, K. Al-Shibli. Prognostic impact of Jab1, p16, p21, p62, Ki67 and Skp2 in soft tissue sarcomas. PLoS One 2012; 7 (10) : e47068-
  33. J. Liu, X.L. Wei, W.H. Huang, C.F. Chen, J.W. Bai, G.J. Zhang. Cytoplasmic Skp2 expression is associated with p-Akt1 and predicts poor prognosis in human breast carcinomas. PLoS One 2012; 7 (12) : e52675-
  34. S.W. Lee, C.F. Li, G. Jin, Z. Cai, F. Han, C.H. Chan. Skp2-dependent ubiquitination and activation of LKB1 is essential for cancer cell survival under energy stress. Mol Cell 2015; 57 (6) : 1022-1033.
  35. Y. Zhang, Y.S. Zvi, B. Batko, N. Zaphiros, E.F. ODonnell, J. Wang. Down-regulation of Skp2 expression inhibits invasion and lung metastasis in osteosarcoma. Sci Rep 2018; 8 (1) : 14294-
  36. L. Ding, R. Li, R. Sun, Y. Zhou, Y. Zhou, X. Han. S-phase kinase-associated protein 2 promotes cell growth and motility in osteosarcoma cells. Cell Cycle 2017; 16 (16) : 1547-1555.
  37. H. Nagao, K. Ijiri, M. Hirotsu, Y. Ishidou, T. Yamamoto, S. Nagano. Role of GLI2 in the growth of human osteosarcoma. The Journal of Pathology (ISSN: 1096-9896). 2011; 224 (2) : 169-79.
  38. J. Wang, O. Aldahamsheh, A. Ferrena, H. Borjihan, A. Singla, S. Yaguare. The interaction of SKP2 with p27 enhances the progression and stemness of osteosarcoma. Annals of the New York Academy of Sciences (ISSN: 1749-6632). 2021; 1490 (1) : 90-104.
  39. D. Liu, R. Zhou, A. Zhou. Identification of key biomarkers and functional pathways in osteosarcomas with lung metastasis: Evidence from bioinformatics analysis. Medicine (Baltimore) 2021; 100 (6) : e24471-
  40. G.K. Schwartz, M.A. Shah. Targeting the cell cycle: a new approach to cancer therapy. Journal of Clinical Oncology (ISSN: 0732-183X). 2005; 23 (36) : 9408-21.
  41. M. Malumbres, M. Barbacid. Cell cycle, CDKs and cancer: a changing paradigm. Nature Reviews. Cancer (ISSN: 1474-1768). 2009; 9 (3) : 153-66.
  42. K. Zhu, L. Liu, J. Zhang, Y. Wang, H. Liang, G. Fan. MiR-29b suppresses the proliferation and migration of osteosarcoma cells by targeting CDK6. Protein Cell 2016; 7 (6) : 434-444.
  43. O. Tetsu, F. McCormick. Proliferation of cancer cells despite CDK2 inhibition. Cancer Cell (ISSN: 1535-6108). 2003; 3 (3) : 233-45.
  44. S. van den Heuvel, E. Harlow. Distinct roles for cyclin-dependent kinases in cell cycle control. Science (ISSN: 0036-8075). 1993; 262 (5142) : 2050-4.
  45. U. Galderisi, F.P. Jori, A. Giordano. Cell cycle regulation and neural differentiation. Oncogene (ISSN: 0950-9232). 2003; 22 (33) : 5208-19.
  46. R. Shekhar, P. Priyanka, P. Kumar, T. Ghosh, M.M. Khan, P. Nagarajan. The microRNAs miR-449a and miR-424 suppress osteosarcoma by targeting cyclin A2 expression. J Biol Chem 2019; 294 (12) : 4381-4400.
  47. Y. Liu, Y. Xi, G. Chen, X. Wu, M. He. URG4 mediates cell proliferation and cell cycle in osteosarcoma via GSK3β/β-catenin/cyclin D1 signaling pathway. J Orthop Surg Res 2020; 15 (1) : 226-
  48. W. Wang, K. Zheng, Y. Pei, X. Zhang. Histone Demethylase JARID1B Is Overexpressed in Osteosarcoma and Upregulates Cyclin D1 Expression via Demethylation of H3K27me3. Oncol Res 2018; 26 (3) : 373-384.
  49. P. Yang, K. Yin, D. Zhong, Q. Liao, K. Li. Inhibition of osteosarcoma cell progression by MacroH2A via the downregulation of cyclin D and cyclin dependent kinase genes. Molecular Medicine Reports (ISSN: 1791-3004). 2015; 11 (3) : 1905-10.
  50. R. Wei, P. Thanindratarn, D.C. Dean, F.J. Hornicek, W. Guo, Z. Duan. Cyclin E1 is a prognostic biomarker and potential therapeutic target in osteosarcoma. Journal of Orthopaedic Research (ISSN: 1554-527X). 2020; 38 (9) : 1952-64.
  51. A. Kumar, S. Kaur, K. Pandit, V. Kaur, S. Thakur, S. Kaur. Onosma bracteata Wall. induces G0/G1 arrest and apoptosis in MG-63 human osteosarcoma cells via ROS generation and AKT/GSK3β/cyclin E pathway. Environmental Science and Pollution Research International (ISSN: 1614-7499). 2021; 28 (12) : 14983-5004.
  52. L. Li, X.A. Kong, M. Zang, J. Dong, Y. Feng, B. Gui. Hsa_circ_0003732 promotes osteosarcoma cells proliferation via miR-545/CCNA2 axis. Biosci Rep 2020; 40 (6) : BSR20191552-
  53. X. Yu, L. Pang, T. Yang, P. Liu. lncRNA LINC01296 regulates the proliferation, metastasis and cell cycle of osteosarcoma through cyclin D1. Oncol Rep 2018; 40 (5) : 2507-2514.
  54. J.S. Wunder, G. Paulian, A.G. Huvos, G. Heller, P.A. Meyers, J.H. Healey. The histological response to chemotherapy as a predictor of the oncological outcome of operative treatment of Ewing sarcoma. The Journal of Bone and Joint Surgery. American Volume (ISSN: 0021-9355). 1998; 80 (7) : 1020-33.
  55. J. Meng, Y. Ding, A. Shen, M. Yan, F. He, H. Ji. Overexpression of PPARγ can down-regulate Skp2 expression in MDA-MB-231 breast tumor cells. Molecular and Cellular Biochemistry (ISSN: 1573-4919). 2010; 345 (1-2) : 171-80.
  56. M. Hirotsu, T. Setoguchi, H. Sasaki, Y. Matsunoshita, H. Gao, H. Nagao. Smoothened as a new therapeutic target for human osteosarcoma. Mol Cancer 2010; 95-
  57. L. Ding, R. Sun, Q. Yan, C. Wang, X. Han, Y. Cui. MiR-506 exerts antineoplastic effects on osteosarcoma cells via inhibition of the Skp2 oncoprotein. Aging (Albany NY) 2021; 13 (5) : 6724-6739.

Comments