Prognostic value of p16, p53, and pcna in sarcoma and an evaluation of immune infiltration
Journal of Orthopaedic Surgery and Research volume 17, Article number: 305 (2022)
p16, p53, and proliferating cell nuclear antigen (pcna) genes play significant roles in many chromatin modifications and have been found to be highly expressed in a variety of tumor tissues. Therefore, they have been used as target genes for some tumor therapies. However, the differential expressions of the p16, p53, and pcna genes in human sarcomas and their effects on prognosis have not been widely reported.
The Oncomine dataset was used to analyze the transcription levels of p16, p53, and pcna genes, and the gene expression profile interactive analysis (GEPIA) dataset was used to analyze the differential expressions of p16, p53, and pcna. The expression levels of p16, p53, and pcna were further analyzed by Western Blotting. GEPIA and Kaplan–Meier analyses were used to analyze the prognostic value of p16, p53, and pcna. Furthermore, p16, p53, and pcna gene mutations and their association with overall survival (OS) and disease-free survival (DFS) were analyzed using cBioPortal datasets. In addition, genes co-expressed with p16, p53, and pcna were analyzed using Oncomine. The DAVID dataset was used to analyze the functional enrichment of p16, p53, pcna, and their co-expressed genes by Gene Ontology (GO) and Metascape were used to construct a network map. Finally, the immune cell infiltration of p16, p53, and pcna in patients with sarcoma was reported by Tumor Immune Estimation Resource (TIMER).
p16, p53, and pcna were up-regulated in human sarcoma tissues and almost all sarcoma cell lines. Western Blotting showed that the expression of p16, p53, and pcna was elevated in osteosarcoma cell lines. The expression of pcna was correlated with OS, the expression of p16, p53, and pcna was correlated with relapse-free survival, and the genetic mutation of p16 was negatively correlated with OS and DFS. We also found that p16, p53, and pcna genes were positively/negatively correlated with immune cell infiltration in sarcoma.
The results of this study showed that p16, p53, and pcna can significantly affect the survival and immune status of sarcoma patients. Therefore, p16, p53, and pcna could be used as potential biomarkers of prognosis and immune infiltration in human sarcoma and provide a possible therapeutic target for sarcoma.
Sarcoma is a rare group of heterogeneous tumors, mainly from the bone and soft tissue, and is highly invasive malignant tumors [1, 2]. Sarcomas can occur in almost all age groups; however, they are more common in adolescents and children and account for 10% of the malignant tumors in children and adolescents . The exact etiopathogenesis of sarcomas is still unclear and needs further research. Osteosarcoma is one of the leading causes of cancer deaths in young people [4, 5]. Childhood and adolescents osteosarcoma has a high incidence [6, 7] and is reported at about 5% . Osteosarcoma has a high rate of metastasis, and the effect of treatment and prognosis are poor . Although surgery combined with chemotherapy slightly improves the survival rate of patients, there has been no substantial improvement in the survival rate in the past 40 years. In addition, the prognosis of patients with metastatic osteosarcoma is very poor, and the overall survival (OS) rate is 30% [10, 11].
Of the tumor suppressor genes, p16 and p53 are very important. Once inactivated, there is a proliferation of malignant cells. With a tumor suppressor gene protein product cdkn2a, p16 mutations seriously affect the progression and prognosis of many tumors. p53 is an important tumor suppressor gene [12, 13]. The p53 gene mutation rate is very high. Due to changes in its spatial conformation, the p53 gene loses the ability to regulate cell growth, cell apoptosis, and DNA repair, thus progressing from a tumor suppressor gene to a cancer gene. Proliferating cell nucleus antigen (pcna) plays an important role in DNA replication and repair. It forms a DNA containing homology trimer ring, anchoring along its side DNA polymerase and other editing enzymes that regulate the DNA and protein sequence motif of pcna—interacting protein tape (PIP—box), thereby influencing the development of tumor .
According to previous research, p16, p53, and pcna genes affect the development of various tumors in many ways. Ishida et al. reported, for example, that the lack of p16 gene was associated with a 91% prevalence of esophageal squamous cell carcinoma and showed that the overexpression of the p16 was caused by an imbalance in the RB1 pathway . Wang and others reported the highest p53 gene mutation frequency in liver cancer patients in China . Ho et al. show that pcna was highly expressed in colon adenocarcinoma, with the expression of pcna in the distant metastasis of tumor as well . However, the mechanisms of p16, p53, and pcna genes in sarcoma have not been widely reported.
The aim of this work was to study the expression of p16, p53, and pcna gene in sarcomas and evaluate its prognostic significance and association with immune cell infiltration.
Transcription of p16, p53, and pcna in patients with sarcoma
Based on Oncomine database data, we found high p16, p53, and pcna expressions in sarcoma (Fig. 1). In addition, the differential transcriptional expression levels of p16, p53, and pcna in sarcoma subtypes and corresponding normal tissues were analyzed using the Detwiller and Barretina sarcoma databases. The expression level of p16 in the Detwiller's sarcoma database was up-regulated, and the multiple changes in these genes in leiomyosarcoma, pleomorphic liposarcoma, round cell liposarcoma, dedifferentiated liposarcoma, and malignant fibrous histiocytoma were 21.899, 25.731, 11.632, 4.897, and 4.789, respectively. Analysis of the Barretina sarcoma database showed that compared with normal tissue, multiple changes of p16 expression levels in dedifferentiated liposarcoma, leiomyosarcoma, pleomorphic liposarcoma, myxoid fibrosarcoma, and myxoid/round cell liposarcoma were 6.870, 5.963, 6.395, 4.967, and 1.300, respectively.
According to the Detwiller sarcoma database analysis, multiple changes in p53 expression in synovial sarcoma and fibrosarcoma compared with normal tissues were 2.660 and 1.885, respectively. In the Barretina sarcoma database, multiple changes in p53 expression in myxoid/round cell liposarcoma and dedifferentiated liposarcoma compared with normal tissues were 4.235 and 1.755, respectively.
Based on the Detwiller sarcoma database analysis, the multiple changes in pcna expression in pleomorphic liposarcoma, leiomyosarcoma, fibrosarcoma, synovial sarcoma, and malignant fibrous histiocytoma were 3.836, 3.155, 3.750, 2.385, and 3.444, respectively. In the Barretina sarcoma database, the multiple changes of pcna expression in pleomorphic liposarcoma, leiomyosarcoma, myxofibrosarcoma, dedifferentiated liposarcoma, and myxoid/round cell liposarcoma were 4.323, 5.521, 4.905, 2.874, and 2.492, respectively.
mRNA expression of p16, p53, and pcna in sarcoma
Using the GEPIA tool, we found that the p16 and pcna mRNA levels were significantly raised in sarcoma compared to normal tissue (P > 0.05), although the level of mRNA expression of p53 showed no significant increase (Fig. 2A–G).
To further evaluate the expression of p16, p53, and pcna in osteosarcoma, we selected two osteosarcoma cell lines to detect their protein expression. Normal osteoblasts constituted the control group. Western Blotting results showed that the expression of p16, p53, and pcna in various osteosarcoma cell lines was significantly higher than that in normal osteoblasts (Fig. 2H), consistent with the results of the Oncomine and GEPIA database analysis.
Prognostic value of p16, p53, and pcna gene in sarcoma patients
Kaplan–Meier plots and GEPIA evaluation demonstrated the prognostic effect of p16, p53, and pcna gene expression in sarcoma patients. The expression of pcna was related to a poor 5- and 10-year OS. With the increase in p53 expression, the 5-year recurrence-free survival (RFS) was up-regulated. The increase in pcna expression was associated with a lowering of the 5- and 10-year RFS in sarcoma, and in recurrent sarcomas, pcna was found to be differentially expressed (Fig. 3). p16 tended to affect the RFS in sarcoma, but this was not statistically significant. We found an inverse correlation between pcna expression and the OS in sarcoma patients. The expression of p16, p53, and pcna was not related to DFS in patients with sarcoma (Fig. 4).
Mutations in p16, p53, and pcna genes and their effects on OS and DFS in sarcoma patients
Next, we analyzed the mutations of p16, p53, and pcna genes and its influence on OS and DFS in sarcoma. High mutation rates of p16, p53, and pcna were observed in sarcomas; 154 of the 254 sarcomas sequenced showed gene mutations with a mutation rate of 61%. Of the three, p16 and p53 gene mutation rates were relatively high, at 17 and 50%, respectively. pcna gene mutation rate was 6% (Fig. 5A). In addition, the results of Kaplan–Meier survival analysis showed that genetic mutations in p16 represented a decrease in OS (Fig. 5B, P = 8.164E − 4) and DFS (Fig. 5C, P = 0.0130) in sarcoma patients. The high frequency of p53 gene mutations tends to be negatively correlated with DFS in sarcoma, but it is not significant (P = 0.0563). The pcna gene mutations were not shown to significantly affect the OS and DFS in sarcoma patients. These results suggest that p16 and p53 gene mutations may significantly influence the outcome of patients with sarcoma.
Identification of genes co-expressed with p16, p53, and pcna genes in osteosarcoma
We used the Oncomine dataset to detect genes co-expressed with p16, p53, and pcna in osteosarcoma. The genes co-expressed with p16 were rnf17, gcnt4, sp140, dio1, acan, clc, tyrobp, grik3, a4gnt, and cacna1c. The genes co-expressed with p53 were chaf1a, gdf15, ap1g2, abcb9, tk1, tle2, mapkapk3, fli1, hip1r, and rfc2. The genes co-expressed with pcna were hist2h4a, ccl2, s100a2, cbs, pdgfrl, nr4a1, nck2, trim27, atp6v1b2, and tk2 (Fig. 6A).
Then, we analyzed the co-expression between p16, p53, and the pcna gene. The results showed that p16 and pcna were positively correlated (R = 0.29, P < 0.05), p53, and p16 showed negative correlation (R = 0.17, P < 0.05). There was no significant correlation between p53 and pcna (Fig. 6B).
Next, the enrichment of David analysis showed that p16, p53, and pcna and co-expressed genes mainly involved were: GO:0,046,104 (thymidine metabolic process), GO:0,007,265 (Ras protein signal transduction), GO:0,048,010 (vascular endothelial growth factor receptor signaling pathway), GO:0,009,157 (deoxyribonucleoside monophosphate biosynthetic process), and GO:0,031,497 (chromatin assembly) of biological processes. Cellular components (CCs), GO:0,005,663 (DNA replication factor C complex), GO:0,016,605 (promyelocytic leukemia [PML] body), GO:0,005,654 (nucleoplasm), GO:0,005,657 (replication fork), and GO:0,032,839 (dendrite cytoplasm) were enriched, and molecular functions (MFs): GO:0,042,802 (the same protein binding): GO:0,004,797 (thymidine kinase activity): GO:0,005,515 (protein binding): GO:0,019,206 (nucleoside kinase activity), GO: 0,019,899 (binding) were enriched (Fig. 7 and Table 1).
To intuitively show the relationship between the enrichment of each gene (p16, p53, and pcna and their co-expressed genes), we used Metascape to construct a network diagram, the cytoscape visualization network with each node representing a rich group (Fig. 8).
Expression of p16, p53, and pcna genes and degree of immune infiltration in sarcoma
In this study, the TIMER was used to explore the correlation between the expression of p16, p53, and pcna genes and immune cell infiltration within the sarcoma. TIMER analysis showed that p16 and infiltrating neutrophils showed a significant correlation (R = 0.102, P < 0.05). Similarly, the expression of p53 and B lymphocyte infiltration also showed a significant correlation (R = 0.103, P < 0.05), and so did CD8 + T cells (R = 0.117, P < 0.05), macrophages (R = 0.184, P < 0.05), neutrophils (R = 0.122, P < 0.05), and dendritic cells (R = 0.197, P < 0.05). The expression of pcna and B lymphocyte infiltration showed a significant correlation (R = 0.1, P < 0.05), so did CD4 + T cells (R = 0.17, P < 0.05), and macrophages (R = 0.122, P < 0.05) (Fig. 9).
p16, p53, and pcna genes have been shown to play an irreplaceable role in the occurrence and development of a variety of malignant tumors. However, gene differential expression of p16, p53, and pcna in sarcoma and its influence on prognosis have not been widely reported. At present, tumor immunotherapy has become a research hotspot, and reports have said that effective immunotherapy has been carried out for many different types of malignant tumors [18,19,20,21,22]. In addition, clinical studies have made some progress in detecting personalized new vaccine antigens in patients with tumors . There is growing attention toward immune-related research in sarcomas and evaluation of immunotherapy combined with radiotherapy, chemotherapy, and some targeted treatment. We conducted this study to reveal the relationship between p16, p53, and pcna genes in human sarcomas with patient prognosis and immune infiltration.
The p16 gene is a very important anticancer gene. Once inactivated, malignant proliferation of cells occurs. At present, the lack of homozygous, nonsense, missense, and p16 frameshift mutations was found in many types of tumor tissue, showing that p16 gene deletions and mutations are widely involved in tumor formation. Zeng et al.  reported a method to monitor the direct effects of the genetic damage of p16 on low-channel primary human melanocytes using precision engineering and digital holography. They observed that the deletion of p16 promotes melanocyte movement through transcriptional activation of brn2, as well as increases the invading and metastasizing ability of the melanoma cells, thereby showing that p16 deletion leads to melanoma progression. Sarun et al.  reported that homozygous deletion of cdkn2a is one of the most common gene mutations in malignant pleural mesothelioma. Tirode et al.  reported that p16 gene mutation frequency of Ewing's sarcoma was 12%. Campbell et al.  reported the important role of cdkn2a (p16) gene in the etiology of pancreatic cancer and thought it to be one of the most common mutations in pancreatic cancer. They also found that cdkn2a promoter methylation plays a key role in the pathogenesis and prognosis of pancreatic cancer . Previous research has shown that changes in the p16 gene can significantly affect the prognosis of patients with many types of tumors. Botana-Rial et al.  reported that patients with p16 ink4a/methylation have a shorter survival time than those without methylation. Pessôan et al.  reported that p16 has a high mutation rate in gliomas, and the higher the degree of malignancy of gliomas, the higher the frequency of p16 mutation (16.7%) and the worse the prognosis. Trietsch and others reported that p16 was often mutated in human papilloma virus-negative vulvar squamous cell carcinoma, and the 5-year survival rate of patients with obvious mutations was worse than patients without mutations . El-mokadem et al.  found that when patients with renal clear cell carcinoma developed p16 gene mutation, the degree of malignancy was usually high, and p16 gene mutation increased the risk of tumor metastasis, leading to higher tumor recurrence and metastasis, which was associated with poor prognosis. These results are similar to the results of our study. The genetic variation rate of p16 is more common in sarcoma patients (17%), and genetic mutations of p16 are associated with poor OS and DFS. It is reported that overexpression of p16 in traditional esophageal squamous cell carcinoma is associated with a better clinical prognosis [32, 33]. Zhou et al.  found that p16 positive lung adenocarcinoma had a good prognosis. Kommoss et al.  found that p16-negative status was an indicator of poor prognosis in clear cell ovarian cancer and mucinous ovarian cancer subgroups. Barber et al.  found that in oropharyngeal squamous cell carcinoma, up-regulation of p16 expression implied better RFS. In the current study, the up-regulation of p16 expression was found to increase the RFS of sarcoma patients, but this was not statistically significant (P = 0.052). There was no significant correlation between the expression of p16 and OS in patients.
As a transcription factor, p53 suppresses cancer mainly through selective transcriptional regulation of multiple target genes, including regulation of apoptosis, cell cycle arrest, senescence, DNA repair, and metabolism [37,38,39]. The existing research is the same as the results of this work; the mutation rate of the p53 gene is so high that it occurs in more than 50% of human malignant tumors [40,41,42]. In many human tumors, p53 plays an important role in tumorigenesis and progression through gene mutation and other mechanisms, including the amplification and/or overexpression of p53-negative regulators, such as mdm2 and mdm4 . It has been reported that mutation or downregulation of p53 (encoded by tp53) accelerates the occurrence and malignant progression of esophageal squamous cell carcinoma . Hou et al.  found that the Claudin 7 closed protein gene (CLDN7) located downstream of p53 on the short arm of chromosome 17 is regulated by WTp53 by binding to the promoter region of colorectal cancer. Once p53 mutation or deletion occurs, the tumor inhibitory function of CLDN7 is lost, indicating that the tumor inhibitory effect of CLDN7 in colorectal cancer is closely related to the status of p53. Rodrigues et al.  have shown that 80% of patients with mantle cell lymphoma with tp53 mutations die within the first 5 years of diagnosis. Marcus and Ladds proposed the use of dihydrowhey acid dehydrogenase and p53 activation as tumor treatment targets, thus effectively killing tumor cells . Takamatsu et al. proposed that inhibiting p53 aggregation through various methods, that is, reducing the evolutionary ability of p53 to nuclear stress, may enhance the effectiveness of cancer treatment and reduce cancer recurrence . George et al.  reported that the mutation of the tumor suppressor gene tp53 is associated with poor survival of patients with acute myeloid leukemia and proposed a new molecular targeted therapy with the main goal of degrading or inactivating mutant p53 or restoring WTp53 to restore normal tp53 function. Saleh et al.  reported that p53 gene mutations are very frequent in high-grade serous ovarian cancer (HGSC), while the frequency of p53 mutations in normal organs is relatively low, making p53 an attractive target for specific therapy of HGSC. These results were similar to ours; the high-frequency genetic changes in p53 were associated with poor DFS in sarcoma patients, but this was not statistically significant (P = 0.0563).
As a hub protein, pcna is considered to be a key regulator of DNA and cell cycle regulation. The expression of pcna is found to be up-regulated in many tumor types, and its overexpression is thought to be related to cancer virulence . Tumor cells have strong proliferative activity and can be used as an index to evaluate the state of cell proliferation. Therefore, pcna has been studied in many tumors, involving the relationship between pcna and tumor occurrence, development, prognosis, recurrence and metastasis, tumor markers, and so on. It is known that high expressions of pcna can promote the proliferation of lung cancer cells and the ability to invade adjacent tissues; therefore, it can be used as a new molecular targeted marker for the diagnosis and treatment of lung cancer [52,53,54]. Jin et al.  have shown that excessive up-regulation of pcna expression is one of the factors that can directly affect the prognosis of rectal cancer, which significantly reduces the survival rate of patients. Wang et al.  found that the transforming growth factor (TGF)-β1 signaling pathway may affect cell growth, cell cycle distribution, and apoptosis by regulating the expression of pcna and other molecules in hepatocellular carcinoma. Hu et al.  reported that pcna expression in gastric cancer tissues was significantly up-regulated compared with normal tissues, and the positive expression of pcna is a risk factor for the prognosis of patients with gastric cancer. Qin et al.  reported that in glioblastoma (GBM), a type of microRNA (miR-1258), inhibits pcna transcription by directly targeting e2f1, which provides a new potential target for GBM therapy and other e2f1-driven cancers. According to the above reports, the up-regulated expression of pcna is related to the poor prognosis of patients with tumors, which is similar to the results we analyzed using Kaplan–Meier graphs, indicating that the increased expression of pcna is related to poor 5- and 10-year OS and RFS in sarcoma patients. Thus, the up-regulation of pcna expression means that patients with sarcoma are more likely to relapse.
A growing body of evidence indicates that the tumor microenvironment plays an important role in the regulation of the proliferation and metastasis [59, 60]. A variety of immune cells in the tumor microenvironment can promote or inhibit the activity of tumor cells. The tumor microenvironment is thought to be an important determinant of the clinical prognosis and response to immunotherapy [61, 62]. Adib et al.  found that patients with cdkn2a -altered urothelial cancer had decreased pd-l1 expression in tumor-infiltrating immune cells, as well as decreased T cell receptors, antigen processing, and activation of interferon-gamma pathways. In urothelial carcinoma, the high expression of pd-l1 in tumor-infiltrating immune cells is associated with an enhanced efficacy of immune checkpoint inhibitors and a good prognosis. After systematically exploring the differentially expressed genes related to immunity in hepatocellular carcinoma, Luo et al.  concluded that the high expression of cdkn2a is associated with a poor prognosis of hepatocellular carcinoma and a decrease in immune infiltration. The up-regulation of cdkn2a expression in hepatocellular carcinoma may be related to the involvement of cdkn2a in the MAPK signaling pathway and the diversity of hepatocellular carcinoma. The degree of tumor immune infiltration will lead to the destruction of the immune microenvironment and immune escape. Tan et al. have proved that the new Golgi programming caused by p53 mutation is an important driver of cancer over secretion. The mechanism is as follows: p53 deletion increases the expression of the Golgi scaffold protein, progesterone, and fat receptor 11 (PAQR11). The secretion of PAQR11-dependent protease PLAU can stimulate the autocrine of PLAU receptor/signal transduction and transcriptional 3-dependent pathway activator, which up-regulates the expression of PAQR11, thus completing a feedforward cycle that amplifies the secretion of premetastatic proteins. Blocking the secretion of PAQR11 can improve the immunosuppressive process in the tumor microenvironment . Kurie et al.  and Tan et al. found that tp53 deletion alleviates G55 (stacking protein 55 kDa) and myosin IIA activated G55-dependent secretion from miR-34a-dependent silencing. The G55-dependent secretory protein enhances the proliferation and invasion of tp53-deficient LUAD cells and the angiogenesis and CD8 + T cell depletion in the tumor microenvironment. Wang et al.  found that inhibiting the Y211 phosphorylation of pcna causes the failure of the replication fork, which drives the biosynthesis of cellular solute ssDNA (single-stranded DNA) into cellular solute. It is suggested that the progress of pcna Y211 phosphorylation in nuclear DNA replication is related to the DNA sensing cascade initiated by cGAS (cyclic GMP-AMP synthase) in the cytoplasm, thus regulating immune surveillance to inhibit tumor metastasis. Understanding the tumor microenvironment may contribute to a better understanding of tumor immune cell interactions to predict response to immunotherapy . In this work, we found that p16, p53, and pcna can be used as surrogate immune markers for sarcoma.
This work had the following limitations: First of all, the data were analyzed from an online database, and there was a lack of adequate in vivo and in vitro testing. Second, the sample size was small. Therefore, a follow-up study with larger sample size is required.
In short, our research suggests that p16, p53, and pcna are overexpressed in human sarcomas. The expression of pcna was correlated with OS, the expression of p16, p53, and pcna was correlated with RFS, and genetic mutations in p16 were negatively correlated with OS and DFS. The p16, p53, and pcna genes were positively/negatively correlated with immune cell infiltration in sarcoma. These results suggest that the p16, p53 and pcna genes significantly affect the prognosis of human sarcomas and are of vital value in immunotherapy. Our work may provide a reference for the choice of new prognostic biomarkers and tumor immunity, and the goal of treatment of sarcomas, providing directions for further research.
Materials and methods
Oncomine database (https://www.oncomine.org/resource/login.html) is usually used to analyze the DNA or RNA sequence of cancer, and genome-wide expression data from malignant tumors . In this work, Oncomine transcription was used to explore the differential expression of p16, p53, and pcna in various cancer tissues.
Gene expression patterns of interaction analysis (GEPIA) assessment of datasets
GEPIA datasets (http://gepia.cancer-pku.cn/) are often used to analyze data generated by the cancer genome atlas project . In this work, we use the GEPIA to explore the differential expression of p16, p53, and pcna genes.
Reagents and cell lines
Human osteosarcoma cell lines MG63 and 143B were purchased from the Shanghai Institute of Cell Biology, Chinese Academy of Sciences. Human primary osteoblasts (Cat.No.GN-H109) were obtained from the Gaining Biological Company (Shanghai, China). The cells were cultured in DMEM medium supplemented with 10% FBS (fetal bovine serum), streptomycin, and penicillin at 37 °C and 5% CO2. Anti-p16 antibody (item no. AF5484), anti-p53 antibody (item no. AF0865), anti-PCNA antibody (item no. AF0239), and anti-β-actin antibody (item no. AF7018) were purchased from Affinity Biosciences (Jiangsu, China).
Protein extraction and Western Blotting test
Total proteins were isolated from cells using the RIPA buffer (Solarbio). The Bradford assay (Bio-Rad Laboratories) was used to quantify protein. The protein was separated on 10% sodium dodecyl sulfate–polyacrylamide gel (SDS-PAGE) and then transferred to a polyvinylidene fluoride (PVDF) membrane (Millipore, USA). The resulting PVDF film was immersed in 5% skim milk and sealed for about 1 h. Next, the membrane was incubated overnight with anti-p16 (1VO1000, Affinity, China), anti-p53 (1VO1000, Affinity, China), and anti-PCNA (1REO1000, Affinity, China) antibodies at 4 °C and incubated with the corresponding secondary antibody (rabbit) at room temperature for 2 h. Finally, autoradiography was performed with the ECL kit (Thermo Fisher Scientific).
We use the Kaplan–Meier plotter database (https://kmplot.com/analysis/index.php? P = service&cancer = pancancer_rna seq) to investigate the effects of p16, p53, and pcna expression on the survival in sarcoma patients [71, 72]. This was compared with the survival curve obtained by the GEPIA dataset analysis.
CBioPortal data analysis
Enrichment analysis of co-expressed genes
David database (http://david.ncifcrf.gov) provides a comprehensive set of functional annotation tool for researchers to understand the biological significance behind a string of genes . It was used in this study to perform ontology (go) analysis of enrichment of p16, p53, and pcna, and the genes. Metascape (https://metascape.org) has been used to provide a comprehensive analysis of the gene annotation and resources, integrating a variety of functions, such as enrichment, interactive group analysis, gene function annotation and members of the search, and providing outputs as unique visual graphics . In this work, we used Metascape to build a network graph of related genes.
Tumor immune to assess resource (TIMER) analysis
Timer web server (https://cistrome.shinyapps.io/timer/) is a comprehensive resource system that analyzes immune cell infiltration across different types of cancer [76, 77]. In this work, the TIMER was used to study the relationship between the expression of p16, p53, and pcna and immune cells penetrating the sarcoma.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the TCGA repository, [https://www.tcga.org/].
Clear cell ovarian cancer
Database for annotation, visualization, and integrated discovery
Gene expression profile interactive analysis
High-grade serous ovarian cancer
Mucinous ovarian cancer
- pcna :
Proliferating cell nuclear antigen
The cancer genome atlas
Tumor immune estimation resource
Fetal bovine serum
Sodium dodecyl sulfate–polyacrylamide gel
Oza J, Doshi SD, Hao L, et al. Homologous recombination repair deficiency as a therapeutic target in sarcoma. Semin Oncol. 2020;47:380–9.
Koelsche C, Schrimpf D, Stichel D, et al. Sarcoma classification by DNA methylation profiling. Nat Commun. 2021;12:498.
Reed DR, Metts J, Pressley M, et al. An evolutionary framework for treating pediatric sarcomas. Cancer. 2020;126:2577–87.
Samsa WE, Mamidi MK, Bashur LA, et al. The crucial p53-dependent oncogenic role of JAB1 in osteosarcoma in vivo. Oncogene. 2020;39:4581–91.
Zhou Y, Yang D, Yang Q, et al. Single-cell RNA landscape of intratumoral heterogeneity and immunosuppressive microenvironment in advanced osteosarcoma. Nat Commun. 2020;11:6322.
Gianferante DM, Mirabello L, Savage SA. Germline and somatic genetics of osteosarcoma - connecting aetiology, biology and therapy. Nat Rev Endocrinol. 2017;13:480–91.
Zhang C, Hansen HM, Semmes EC, et al. Common genetic variation and risk of osteosarcoma in a multi-ethnic pediatric and adolescent population. Bone. 2020;130: 115070.
Anderson ME. Update on survival in osteosarcoma. Orthop Clin North Am. 2016;47:283–92.
Zhang W, Zhuang Y, Zhang Y, et al. Uev1A facilitates osteosarcoma differentiation by promoting Smurf1-mediated Smad1 ubiquitination and degradation. Cell Death Dis. 2017;8: e2974.
Chen C, Xie L, Ren T, et al. Immunotherapy for osteosarcoma: fundamental mechanism, rationale, and recent breakthroughs. Cancer Lett. 2021;500:1–10.
Angulo P, Kaushik G, Subramaniam D, et al. Natural compounds targeting major cell signaling pathways: a novel paradigm for osteosarcoma therapy. J Hematol Oncol. 2017;10:10.
Liu J, Zhang C, Hu W, et al. Tumor suppressor p53 and metabolism. J Mol Cell Biol. 2019;11:284–92.
Levine AJ. The many faces of p53: something for everyone. J Mol Cell Biol. 2019;11:524–30.
González-Magaña A, Blanco FJ. Human PCNA structure function and interactions. Biomolecules. 2020;10(4):570.
Ishida H, Kasajima A, Fujishima F, et al. p16 in highly malignant esophageal carcinomas: the correlation with clinicopathological factors and human papillomavirus infection. Virchows Arch. 2021;478:219–29.
Wang S, Shi H, Liu T, et al. Mutation profile and its correlation with clinicopathology in Chinese hepatocellular carcinoma patients. Hepatobiliary Surg Nutr. 2021;10:172–9.
Ho YJ, Lin YM, Huang YC, et al. Prognostic significance of high YY1AP1 and PCNA expression in colon adenocarcinoma. Biochem Biophys Res Commun. 2017;494:173–80.
Sivori S, Pende D, Quatrini L, et al. NK cells and ILCs in tumor immunotherapy. Mol Aspects Med. 2021;80: 100870.
Anderson NR, Minutolo NG, Gill S, et al. Macrophage-based approaches for cancer immunotherapy. Cancer Res. 2021;81:1201–8.
Morales E, Olson M, Iglesias F, et al. Role of immunotherapy in Ewing sarcoma. J Immunother Cancer. 2020. https://doi.org/10.1136/jitc-2020-000653.
Wang H, Najibi AJ, Sobral MC, et al. Biomaterial-based scaffold for in situ chemo-immunotherapy to treat poorly immunogenic tumors. Nat Commun. 2020;11:5696.
Guo W, Ma J, Guo S, et al. A20 regulates the therapeutic effect of anti-PD-1 immunotherapy in melanoma. J Immunother Cancer. 2020. https://doi.org/10.1136/jitc-2020-001866.
Blass E, Ott PA. Advances in the development of personalized neoantigen-based therapeutic cancer vaccines. Nat Rev Clin Oncol. 2021;18:215–29.
Zeng H, Jorapur A, Shain AH, et al. Bi-allelic loss of CDKN2A initiates melanoma invasion via BRN2 activation. Cancer Cell. 2018;34:56-68.e59.
Sarun KH, Lee K, Williams M, et al. Genomic deletion of BAP1 and CDKN2A are useful markers for quality control of malignant pleural mesothelioma (MPM) primary cultures. Int J Mol Sci. 2018. https://doi.org/10.3390/ijms19103056.
Tirode F, Surdez D, Ma X, et al. Genomic landscape of Ewing sarcoma defines an aggressive subtype with co-association of STAG2 and TP53 mutations. Cancer Discov. 2014;4:1342–53.
Campa D, Pastore M, Gentiluomo M, et al. Functional single nucleotide polymorphisms within the cyclin-dependent kinase inhibitor 2A/2B region affect pancreatic cancer risk. Oncotarget. 2016;7:57011–20.
Botana-Rial M, De Chiara L, Valverde D, et al. Prognostic value of aberrant hypermethylation in pleural effusion of lung adenocarcinoma. Cancer Biol Ther. 2012;13:1436–42.
Pessôa IA, Amorim CK, Ferreira WAS, et al. Detection and correlation of single and concomitant TP53, PTEN, and CDKN2A alterations in gliomas. Int J Mol Sci. 2019. https://doi.org/10.3390/ijms20112658.
Trietsch MD, Spaans VM, ter Haar NT, et al. CDKN2A(p16) and HRAS are frequently mutated in vulvar squamous cell carcinoma. Gynecol Oncol. 2014;135:149–55.
El-Mokadem I, Kidd T, Pratt N, et al. Tumour suppressor gene (CDKNA2) status on chromosome 9p in resected renal tissue improves prognosis of localised kidney cancer. Oncotarget. 2016;7:73045–54.
Cao F, Zhang W, Zhang F, et al. Prognostic significance of high-risk human papillomavirus and p16(INK4A) in patients with esophageal squamous cell carcinoma. Int J Clin Exp Med. 2014;7:3430–8.
Fujiwara S, Noguchi T, Takeno S, et al. Hypermethylation of p16 gene promoter correlates with loss of p16 expression that results in poorer prognosis in esophageal squamous cell carcinomas. Dis Esophagus. 2008;21:125–31.
Zhou Y, Höti N, Ao M, et al. Expression of p16 and p53 in non-small-cell lung cancer: clinicopathological correlation and potential prognostic impact. Biomark Med. 2019;13:761–71.
Kommoss S, Gilks CB, du Bois A, et al. Ovarian carcinoma diagnosis: the clinical impact of 15 years of change. Br J Cancer. 2016;115:993–9.
Barber BR, Biron VL, Klimowicz AC, et al. Molecular predictors of locoregional and distant metastases in oropharyngeal squamous cell carcinoma. J Otolaryngol Head Neck Surg. 2013;42:53.
Demir Ö, Barros EP, Offutt TL, et al. An integrated view of p53 dynamics, function, and reactivation. Curr Opin Struct Biol. 2021;67:187–94.
D’Orazi G. Recent advances in p53. Biomolecules. 2021. https://doi.org/10.3390/biom11020211.
Valente LJ, Tarangelo A, Li AM, et al. p53 deficiency triggers dysregulation of diverse cellular processes in physiological oxygen. J Cell Biol. 2020. https://doi.org/10.1083/jcb.201908212.
Redman-Rivera LN, Shaver TM, Jin H, et al. Acquisition of aneuploidy drives mutant p53-associated gain-of-function phenotypes. Nat Commun. 2021;12:5184.
Zhu G, Pan C, Bei JX, et al. Mutant p53 in cancer progression and targeted therapies. Front Oncol. 2020;10: 595187.
Tiwari B, Jones AE, Caillet CJ, et al. p53 directly represses human LINE1 transposons. Genes Dev. 2020;34:1439–51.
Liu J, Zhang C, Wang J, et al. The Regulation of ferroptosis by tumor suppressor p53 and its pathway. Int J Mol Sci. 2020. https://doi.org/10.3390/ijms21218387.
Wang H, Song X, Wang Y, et al. CircCNTNAP3-TP53-positive feedback loop suppresses malignant progression of esophageal squamous cell carcinoma. Cell Death Dis. 2020;11:1010.
Hou Y, Hou L, Liang Y, et al. The p53-inducible CLDN7 regulates colorectal tumorigenesis and has prognostic significance. Neoplasia. 2020;22:590–603.
Rodrigues JM, Hassan M, Freiburghaus C, et al. p53 is associated with high-risk and pinpoints TP53 missense mutations in mantle cell lymphoma. Br J Haematol. 2020;191:796–805.
Ladds M, Popova G, Pastor-Fernández A, et al. Exploitation of dihydroorotate dehydrogenase (DHODH) and p53 activation as therapeutic targets: a case study in polypharmacology. J Biol Chem. 2020;295:17935–49.
Takamatsu Y, Ho G, Wada R, et al. Adiponectin paradox as a therapeutic target of the cancer evolvability in aging. Neoplasia. 2021;23:112–7.
George B, Kantarjian H, Baran N, et al. TP53 in acute myeloid leukemia: molecular aspects and patterns of mutation. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms221910782.
Saleh A, Perets R. Mutated p53 in HGSC-from a common mutation to a target for therapy. Cancers. 2021. https://doi.org/10.3390/cancers13143465.
Shemesh A, Kundu K, Peleg R, et al. NKp44-derived peptide binds proliferating cell nuclear antigen and mediates tumor cell death. Front Immunol. 2018;9:1114.
Chen X, Sun J, Wang Y. Expressions of CD44, PCNA and MRP1 in lung cancer tissues and their effects on proliferation and invasion abilities of lung cancer cell line 95D. J buon. 2021;26:72–8.
Wu C, Zhu XT, Xia L, et al. High expression of long noncoding RNA PCNA-AS1 promotes non-small-cell lung cancer cell proliferation and oncogenic activity via upregulating CCND1. J Cancer. 2020;11:1959–67.
Wang L, Kong W, Liu B, et al. Proliferating cell nuclear antigen promotes cell proliferation and tumorigenesis by up-regulating STAT3 in non-small cell lung cancer. Biomed Pharmacother. 2018;104:595–602.
Jin LJ, Chen WB, Zhang XY, et al. Analysis of factors potentially predicting prognosis of colorectal cancer. World J Gastrointest Oncol. 2019;11:1206–17.
Wang XH, Chen ZG, Xu RL, et al. TGF-β1 signaling pathway serves a role in HepG2 cell regulation by affecting the protein expression of PCNA, gankyrin, p115, XIAP and survivin. Oncol Lett. 2017;13:3239–46.
Hu L, Li HL, Li WF, et al. Clinical significance of expression of proliferating cell nuclear antigen and E-cadherin in gastric carcinoma. World J Gastroenterol. 2017;23:3721–9.
Qin H, Gui Y, Ma R, et al. miR-1258 attenuates tumorigenesis through targeting E2F1 to inhibit PCNA and MMP2 transcription in glioblastoma. Front Oncol. 2021;11: 671144.
Sun Y, Li Y, Shi S, et al. Exploiting a new approach to destroy the barrier of tumor microenvironment: nano-architecture delivery systems. Molecules. 2021. https://doi.org/10.3390/molecules26092703.
Pires A, Greenshields-Watson A, Jones E, et al. Immune remodeling of the extracellular matrix drives loss of cancer stem cells and tumor rejection. Cancer Immunol Res. 2020;8:1520–31.
Jang H, Kim EH, Chi SG, et al. Nanoparticles targeting innate immune cells in tumor microenvironment. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms221810009.
Ling B, Huang Z, Huang S, et al. Microenvironment analysis of prognosis and molecular signature of immune-related genes in lung adenocarcinoma. Oncol Res. 2021;28:561–78.
Adib E, Nassar AH, Akl EW, et al. CDKN2A alterations and response to immunotherapy in solid tumors. Clin Cancer Res. 2021;27:4025–35.
Luo JP, Wang J, Huang JH. CDKN2A is a prognostic biomarker and correlated with immune infiltrates in hepatocellular carcinoma. 2021. Biosci Rep. https://doi.org/10.1042/BSR20211103.
Tan X, Banerjee P, Shi L, et al. p53 loss activates prometastatic secretory vesicle biogenesis in the Golgi. Sci Adv. 2021. https://doi.org/10.1126/sciadv.abf4885.
Tan X, Shi L, Banerjee P, et al. A protumorigenic secretory pathway activated by p53 deficiency in lung adenocarcinoma. J Clin Invest. 2021. https://doi.org/10.1172/JCI137186.
Wang YL, Lee CC, Shen YC, et al. Evading immune surveillance via tyrosine phosphorylation of nuclear pcna. Cell Rep. 2021;36: 109537.
Sadeghi Rad H, Monkman J, Warkiani ME, et al. Understanding the tumor microenvironment for effective immunotherapy. Med Res Rev. 2021;41:1474–98.
Rhodes DR, Yu J, Shanker K, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6:1–6.
Tang Z, Li C, Kang B, et al. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98-w102.
Nagy Á, Lánczky A, Menyhárt O, et al. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep. 2018;8:9227.
Győrffy B. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput Struct Biotechnol J. 2021;19:4101–9.
Gao J, Aksoy BA, Dogrusoz U, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013. https://doi.org/10.1126/scisignal.20040886:pl1.
Huang DW, Sherman BT, Tan Q, et al. DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35:W169-175.
Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Li T, Fu J, Zeng Z, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48:W509-w514.
Li T, Fan J, Wang B, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77:e108–10.
This work was supported by Anhui Science and Technology Department [S202007d0705001].
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cai, D., Ma, X., Guo, H. et al. Prognostic value of p16, p53, and pcna in sarcoma and an evaluation of immune infiltration. J Orthop Surg Res 17, 305 (2022). https://doi.org/10.1186/s13018-022-03193-3
- p16 (cdkn2a) gene
- Tp53 gene
- Pcna gene
- Immune infiltration