- 1Department of Urology, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 2Department of Nephrology, the First Affiliated Hospital, Zhengzhou University, Zhengzhou, China
Purpose: The present study aimed to establish a hypoxia related genes model to predict the prognosis of kidney clear cell carcinoma (KIRC) patients using data accessed from The Cancer Genome Atlas (TCGA) database and International Cancer Genome Consortium (ICGC) database.
Methods: Patients’ data were downloaded from the TCGA and ICGC databases, and hypoxia related genes were accessed from the Molecular Signatures Database. The differentially expressed genes were evaluated and then the differential expressions hypoxia genes were screened. The TCGA cohort was randomly divided into a discovery TCGA cohort and a validation TCGA cohort. The discovery TCGA cohort was used for constructing the hypoxia genes risk model through Lasso regression, univariate and multivariate Cox regression analysis. Receiver operating characteristic (ROC) curves were used to assess the reliability and sensitivity of our model. Then, we established a nomogram to predict the probable one-, three-, and five-year overall survival rates. Lastly, the Tumor Immune Dysfunction and Exclusion (TIDE) score of patients was calculated.
Results: We established a six hypoxia-related gene prognostic model of KIRC patients in the TCGA database and validated in the ICGC database. The patients with high riskscore present poorer prognosis than those with low riskscore in the three TCGA cohorts and ICGC cohort. ROC curves show our six-gene model with a robust predictive capability in these four cohorts. In addition, we constructed a nomogram for KIRC patients in the TCGA database. Finally, the high risk-group had a high TIDE score than the patients with low riskscore.
Conclusions: We established a six hypoxia-related gene risk model for independent prediction of the prognosis of KIRC patients was established and constructed a robust nomogram. The different riskscores might be a biomarker for immunotherapy strategy.
Introduction
Renal cell carcinoma is one of the most common malignant tumors in the urinary system (1). Clear cell renal cell carcinoma (ccRCC) is the most common type of RCC, and accounts for approximately 80% of cases, with a 10-year cancer specific-survival rate of 62% (2). Radical nephrectomy is the standard surgical treatment for localized and locally advanced ccRCC. However, there is a lack of sufficient effective treatment for metastatic ccRCC and its 5-year survival rate of 12% (1). Now the immune checkpoint inhibition agents have been recommended as the first line agents for metastatic RCC instead of targeted therapy (3). Constructing a prognostic model that can accurately predict the prognosis has become a key issue in the diagnosis and formulation of treatment strategies for renal cancer.
Nowadays, a few prognostic systems and nomograms have been developed based on the clinical characters (such as TNM stage and Fuhrman grade) and laboratory factors (such as hemoglobin, neutrophil count, and platelet count) (2). With the development of high-throughput sequencing technology, the nomogram based on the RNA-seq data or methylation data have been constructed (4–9). These nomograms have not only improved the capability of predicting the prognosis of ccRCC but have also indicated the pathogenic mechanism of ccRCC which may help to develop drug candidates.
VHL gene mutation plays a central role in the initiation and progression of ccRCC (10, 11). As VHL gene is a key regulator of hypoxia, and hypoxia related genes can facilitate tumorigenesis (12), whether the hypoxia pathway could predict the prognosis of ccRCC was interesting. Now, we used the TCGA and ICGC data to investigate the role of hypoxia pathway in predicting the prognosis of ccRCC, and our results constructed a hypoxia related genes model and a nomogram to improve the prognostic value of ccRCC based on bioinformatics approaches. We also investigate the role of our prognosis model in the immunotherapy strategy.
Materials and Methods
Patients and Public Datasets
RNA-seq data of 533 KIRC and 78 adjacent normal kidney tissues were downloaded from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) in October 2018. These original data were processed according to the guidelines and the data were merged into a single expression file for further study. All these expression data were used for screening the differential expression genes between normal and KIRC tissues. The clinical data of KIRC patients were also obtained, and 512 patients who had identified information of age, gender, clinical stage, TNM stage, overall survival status, and survival time were finally enrolled in the prognosis study (Supplementary File; Supplementary Table). The International Cancer Genome Consortium (ICGC) cohort (RECA-EU project, 91 patients) clinical information and RNA expression data were downloaded from the ICGC data portal in April 2020 (Supplementary File, Supplementary Table). The TCGA was used as a public open database, and the relevant information retrieved from it did not require further ethical approval.
Screening of Differential Expression Hypoxia Genes (DEHGs)
RNA-seq data were standardized using the “limma” package of R software. Then differentially expressed genes (DEGs) were evaluated using the “DESeq” package. Genes with adjusted p-value <0.05 and absolute log2 fold change (FC) >2 in normal and tumors tissues were considered as DEGs. The hypoxia related genes were downloaded from the Molecular Signatures Database (https://www.gsea-msigdb.org/gsea/msigdb/cards/HAL LMARK_HYPOXIA.html). The hypoxia related genes in DEGs and the ICGC genes list were considered as differential expression hypoxia genes (DEHGs).
Prognosis Related Genes Filtering and Gene Risk Model Construction
The whole 512 patients were named Total TCGA cohort. Total TCGA cohort was then divided into a discovery TCGA cohort (256 patients) and a validation TCGA cohort (256 patients) by a ratio of 1 to 1 randomly (Supplementary File). The discovery TCGA cohort was used for constructing the hypoxia genes risk model. Firstly, univariate Cox expression was used to assess the prognosis role of these DEHGs. After identifying the prognosis related DEHGs, LASSO (least absolute shrinkage and selection operator) regression was used to select a panel of genes, and the analysis was performed using “glmnet” package in R software. Lastly, the genes which were extracted by lasso regression, were enrolled in multivariate Cox regression to further screen the prognosis related DEHGs and then construct a risk model. Riskscore was calculated according to the coefficient and expression value of each DEHG which had the significant meaning in multivariate Cox regression.
Genes Risk Model Validation, Nomogram Construction and TIDE Score Calculation
Risk-groups were deemed as two categorical variables (high risk-group and low risk-group) according to the median value of riskscore. Kaplan–Meier curve and log-rank test were used to assess the prognosis role of risk-group. Next, univariate and multivariate Cox regression were used to identify whether riskscore was an independent factor for the prognosis of KIRC patients. In the Cox regression analysis, gender, T stage, N stage, M stage, clinical stage, and risk-group were considered as categorical variables. Meanwhile, age and riskscore were considered as continuous variables. Validation TCGA cohort, total TCGA cohort, and ICGC cohort underwent the same analysis procedure to validate the genes risk model. Moreover, Receiver operating characteristic curve (ROC) analysis was performed to assess the predictive accuracy of the riskscore in all these four cohorts. In addition, a nomogram was established based on the total TCGA cohort for clinical application. Statistical analyses were performed using R software. Tumor Immune Dysfunction and Exclusion (TIDE) score was calculated online (HTTP://tide.dfci.harvard.edu/) (13). All tests were two-tailed, and a P-value of <0.05 was considered statistically significant.
Results
Data Processing and Screening of Differentially Expressed Genes
A total of 5,825 DEGs were screened (3,370 upregulated and 2,455 downregulated) from a total of 34,827 genes using the adjust p-value <0.05 and log FC (fold change) >1 threshold (Supplementary File). Two hundred hypoxia-related genes were acquired from the Molecular Signatures Database. Finally, 72 genes DEHGs were identified both in total TCGA and ICGC cohorts (Supplementary File; Supplementary Figure).
Screening of Prognostic DEHGs
The processed survival data of each tumor sample in the training set were subjected to univariate Cox proportional hazards regression analysis, in which the significant threshold was set at value <0.05. Therefore, 15 prognosis-related DEHGs containing 10 risky genes and 5 protective genes were identified (Supplementary File).
Establishment and Validation of the 6-DEHGs Prognosis Model
LASSO regression with tenfold cross validation was performed to further screen the DEHGs that significantly correlated with the prognosis of ccRCC patients. The optimal lambda value was obtained from the minimum partial likelihood deviance (Figures 1A, B). Then, the optimal 6-DEHG model was obtained which contained GPC3, KIF5A, PLAUR, ANKZF1, ETS1, and SELENBP1. The Cox coefficients of the DEHGs were obtained from the multivariate Cox proportional hazards regression analysis. Then, the riskscore was constructed based on the coefficients and categorical values of expression level as the following: Riskscore = (0.1611 ∗ GPC3) + (0.0781 ∗ KIF5A) + (0.2598 ∗ PLAUR) + (0.4296 ∗ ANKZF1) − (0.2765 ∗ ETS1) + (0.4578 ∗ SELENBP1). The riskscore of each tumor patient in the discovery TCGA cohort was calculated and the patients were divided into a high risk-group and a low risk-group based on the median riskscore (Figure 2). Then, the other three cohorts were also divided into a high risk-group and a low risk-group based on the median riskscore of each cohort individually. Patients with the high risk-group presented with poorer survival (Figure 3A) and it had been validated in the other three cohorts (Figures 3B–D).
Figure 1 LASSO regression analysis of 15 DEHGs. Ten-fold cross-validation was applied to calculate the best lambda, which leads to a minimum mean cross-validated error (A). A total of 6 DEHGs were adopted for the LASSO model (B).
Figure 2 Correlation between the prognostic model and the overall survival (OS) of patients in the Discovery TCGA Cohort (A), the Validation TCGA Cohort (B), the Total TCGA Cohort (C) and the ICGC Cohort (D). The distribution of riskscore (upper), survival time (middle) and hypoxia genes expression levels (below). Patients were classified into low-risk and high-risk groups by using the median score as a cut-off value. The red dots and lines represent the patients in high-risk groups. The green dots and lines represent the patients in low-risk groups.
Figure 3 Kaplan–Meier survival of the 6-gene prognostic model in the Discovery TCGA Cohort (A), the Validation TCGA Cohort (B), the Total TCGA Cohort (C) and the ICGC Cohort (D).
Construction of a 6-DEHG Prognosis Model-Based Nomogram
In addition, Cox regression analysis remained that riskscore was an independent prognostic factor influencing patients with ccRCC in these four cohorts (Figure 4, Table 1). The AUC value for the DEHG Prognosis model was 0.711 in the 1-year ROC curve, 0.708 in the 3-year ROC curve, and 0.779 in the 5-year ROC curve of discovery TCGA cohort (Figure 5A). The AUC value of the other three validation cohorts is shown in Figures 5B–D. The riskscore, patients’ age, and clinical stage could be an independent prognostic factor, respectively, for OS prediction of ccRCC patients in the training set after the univariate and multivariate Cox proportional hazards regression analyses. Then, these independent prognostic factors were integrated together into this nomogram to predict the 1-, 3-, and 5-year OS of ccRCC patients (Figure 6).
Figure 4 Multivariate Cox regression analyses. Discovery TCGA Cohort (A), Validation TCGA Cohort (B), Total TCGA Cohort (C) and ICGC Cohort (D). *p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 5 ROC curves with calculated AUCs were established to evaluate the prognostic value of the 6-gene model in 5-, 3-, 1-years in the Discovery TCGA Cohort (A), the Validation TCGA Cohort (B), the Total TCGA Cohort (C), and the ICGC Cohort (D).
Figure 6 Construction of a nomogram for overall survival prediction in KIRC. The nomogram consists of age, clinical stage and the riskscore based on the six-hypoxia gene model.
The Probably Benefit of Patients for Immune Checkpoint Inhibitor (ICI) Therapy in Different Risk Groups
The TIDE-score was used to assess the potential clinical efficacy of immunotherapy in different risk groups. In our results, the high-risk group had a higher TIDE score than the low risk-group (Figure 7A). Also, we found that the high risk-group had a lower microsatellite instability (MSI) score, while the high risk-group had a higher T cell dysfunction score, but there was no difference in T cell exclusion between the two groups (Figures 7B–D).
Figure 7 (A) TIDE score, (B) MSI, (C) T cell Dysfunction and (D) Exclusion in different risk-groups. The score between the two risk-groups were compared through Wilcoxon test (ns, not significant, ***p < 0.001).
Discussion
To achieve the precision medicine of kidney cancer, accurate diagnosis, individualized therapeutic strategies are needed. Providing the detail prognosis information for patients make the precision medicine more meaningful. Hypoxia pathway plays the central function in ccRCC pathogenesis and development. In our study, we screened the hypoxia-related genes for predicting the prognosis of ccRCC and then constructed a six hypoxia-related gene prognosis model using the LASSO regression method. The model was validated by the validation TCGA cohort, the total TCGA cohort, and the ICGC cohort. A nomogram, based on the clinical features and our model, was built for predicting the overall survival probability of ccRCC patients. In addition, we found the different immune function statuses in the different risk groups.
Recently, prognostic model integrating hypoxia-related gene expression and clinical features emerged (14–16). In our study, the hypoxia-related genes prognosis model was made up of six genes: GPC3, KIF5A, PLAUR, ANKZF1, ETS1, and SELENBP1. In the calculation formula of riskscore, the coefficient of ETS1 was a negative numbers, while the coefficient of the other five genes was a positive number. Therefore, there was a negative relationship between riskscore and ETS1, while there was a positive relationship between riskscore and GPC3, KIF5A, PLAUR, ANKZF1, and SELENBP1. Glypican 3 (GPC3) is an oncofetal proteoglycan anchored on the cell membrane (17). High GPC3 expression was also extensively associated with worse tumor differentiation, later tumor stage, presence of vascular invasion, and hepatitis B virus (HBV) infection (18). This protein is expressed in the liver and the kidney of healthy fetuses but is hardly expressed in adults, except in the placenta. Contrarily, GPC3 is specifically expressed in hepatocellular carcinoma (HCC), ovarian clear cell carcinoma, melanoma, squamous cell carcinoma of the lung, hepatoblastoma, nephroblastoma (Wilms tumor), yolk sac tumor, and some pediatric cancers (19). In addition, GPC3 has a low expression in ccRCC tissues than normal kidney tissues, and it can reduce the proliferation of RCC cell lines (20). Kinesin family member 5A (KIF5A) is a member of the kinesin superfamily which can modulate many cell physiological behavior such as the proliferation of cell cycle (21). KIF5A mutation causes familial Amyotrophic lateral sclerosis (ALS; OMIM: 05400) (22). Tian et al. have proved that KIF5A can regulate the bladder cancer development and progression (21). Plasminogen activator, urokinase receptor (PLAUR) encodes the receptor for urokinase plasminogen activator and is a key molecule in regulating of cell-surface plasminogen activation (23). Previous studies have shown that PLAUR was involved in cancer cell migration, invasion and metastasis processing and could predict the prognosis of many cancers, such as glioma and oral squamous cell carcinoma (24–27). Ankyrin repeat and zinc finger peptidyl tRNA hydrolase 1 (ANKZF1) is essential for the oxidative stress and the maintenance of mitochondrial integrity (28). Some bioinformatics analyses have proved that ANKZF1 based genes signature can predict the prognosis of patients with colon cancer, renal cell carcinoma, and prostate cancer (14, 29–33). ETS proto-oncogene 1, transcription factor (ETS1) is a transcription factor belong to the ETS domain family. Ets1 expression is linked to poor survival of some cancers and contributes to the cancer cell invasiveness, EMT, and drug resistance. In addition, as a major MAPK downstream molecule, ETS1 might be a prospective therapeutic target for cancer (34, 35). Moreover, high expression levels of ETS1 were associated with poor survival in RCC tissues (36). Selenium binding protein 1 (SELENBP1) encodes protein and plays a selenium-dependent role in many physiological functions, such as protein degradation, cell differentiation and redox modulation. SELENBP1 is downregulated in many cancer types and its low expression levels are associated with poorer prognosis (37). In RCC, SELENBP1 has been proved as a tumor suppressor gene and low SELENBP1 mRNA expression predicts a worse cancer-specific survival (38).
Although tumor stage, lymph node metastasis and grade are significantly associated with the prognosis of KIRC (39), prognostic model integrating gene expression and clinical features are gained increasing attention along with the development of sequencing technology. Here, we developed a hypoxia-related genes prognostic model and a significant difference in overall survival was observed between high-and low-risk subgroups. The 1- (TCGA 0.68, ICGC 0.72), 3- (TCGA 0.67, ICGC 0.70), and 5- (TCGA 0.71, ICGC 0.66) year AUC showed a good accuracy of ROC curves in the TCGA cohort and the ICGC cohort, respectively. Hypoxia-related genes can be promising prognostic biomarkers for KIRC which gain evidence for hypoxia-related genes targeted therapy.
TIDE score could predict patient response to immunotherapy as it can reflect the potential capacity for the tumor's immune evasion (13).In our results, the high risk-group had a higher TIDE score than the low risk-group, indicating that patients with low riskscore could benefit more from ICI therapy than patients with high riskscore (Figure 7A). In addition, high TIDE score predicted a worse outcome, which was accordance with our results and it might interpret the high riskscore group presented a worse prognosis .
In conclusion, a six hypoxia-related gene risk model for independent prediction of the prognosis of KIRC patients was established. The different riskscores might be a biomarker for immunotherapy strategy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by the TCGA Research Network: http://cancergenome.nih.gov/. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
X-hN and J-jY conceived and designed the experiment. X-hN, Y-yQ, N-yL, S-cL, and Z-kJ performed the experiments and analyzed the data. Interpretation of the findings were done by X-hN, Y-yQ, and Z-kJ. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Postdoctoral Research Grant in Henan Province (grant numbers 1901004, 1902005), the Henan Science and Technology Research Program (grant number 2018020142) and The Natural Science Foundation of Henan Province (212300410265). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The results shown here are in whole or part based upon data generated by the TCGA Research Network: http://cancergenome.nih.gov/.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.806264/full#supplementary-material
References
1. Padala SA, Barsouk A, Thandra KC, Saginala K, Mohammed A, Vakiti A, et al. Epidemiology of Renal Cell Carcinoma. World J Oncol (2020) 11:79–87. doi: 10.14740/wjon1279
2. Heng DY, Xie W, Regan MM, Harshman LC, Bjarnason GA, Vaishampayan UN, et al. External Validation and Comparison With Other Models of the International Metastatic Renal-Cell Carcinoma Database Consortium Prognostic Model: A Population-Based Study. Lancet Oncol (2013) 14:141–8. doi: 10.1016/S1470-2045(12)70559-4
3. Bedke J, Albiges L, Capitanio U, Giles RH, Hora M, Lam TB, et al. Updated European Association of Urology Guidelines on Renal Cell Carcinoma: Nivolumab Plus Cabozantinib Joins Immune Checkpoint Inhibition Combination Therapies for Treatment-Naive Metastatic Clear-Cell Renal Cell Carcinoma. Eur Urol (2021) 79:339–42. doi: 10.1016/j.eururo.2020.12.005
4. Gao R, Qin H, Lin P, Ma C, Li C, Wen R, et al. Development and Validation of a Radiomic Nomogram for Predicting the Prognosis of Kidney Renal Clear Cell Carcinoma. Front Oncol (2021) 11:613668. doi: 10.3389/fonc.2021.613668
5. Li H, Hu J, Yu A, Othmane B, Guo T, Liu J, et al. RNA Modification of N6-Methyladenosine Predicts Immune Phenotypes and Therapeutic Opportunities in Kidney Renal Clear Cell Carcinoma. Front Oncol (2021) 11:642159. doi: 10.3389/fonc.2021.642159
6. Yoo SH, Yun J, Keam B, Hong SP, Ock CY, Koh J, et al. Discovery of Acquired Molecular Signature on Immune Checkpoint Inhibitors in Paired Tumor Tissues. Cancer Immunol Immunother (2021) 70:1755–69. doi: 10.1007/s00262-020-02799-y
7. Yu J, Mao W, Sun S, Hu Q, Wang C, Xu Z, et al. Identification of an M6a-Related lncRNA Signature for Predicting the Prognosis in Patients With Kidney Renal Clear Cell Carcinoma. Front Oncol (2021) 11:663263. doi: 10.3389/fonc.2021.663263
8. Zhang W, Li C, Wu F, Li N, Wang Y, Hu Y, et al. Analyzing and Validating the Prognostic Value of a TNF-Related Signature in Kidney Renal Clear Cell Carcinoma. Front Mol Biosci (2021) 8:689037. doi: 10.3389/fmolb.2021.689037
9. Zheng W, Zhang S, Guo H, Chen X, Huang Z, Jiang S, et al. Multi-Omics Analysis of Tumor Angiogenesis Characteristics and Potential Epigenetic Regulation Mechanisms in Renal Clear Cell Carcinoma. Cell Commun Signal (2021) 19:39. doi: 10.1186/s12964-021-00728-9
10. Clark PE. The Role of VHL in Clear-Cell Renal Cell Carcinoma and its Relation to Targeted Therapy. Kidney Int (2009) 76:939–45. doi: 10.1038/ki.2009.296
11. Gossage L, Eisen T, Maher ER. VHL, the Story of a Tumour Suppressor Gene. Nat Rev Cancer (2014) 15:55–64. doi: 10.1038/nrc3844
12. Schodel J, Grampp S, Maher ER, Moch H, Ratcliffe PJ, Russo P, et al. Hypoxia, Hypoxia-Inducible Transcription Factors, and Renal Cancer. Eur Urol (2016) 69:646–57. doi: 10.1016/j.eururo.2015.08.007
13. Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T Cell Dysfunction and Exclusion Predict Cancer Immunotherapy Response. Nat Med (2018) 24:1550–8. doi: 10.1038/s41591-018-0136-1
14. Li Z, Du G, Zhao R, Yang W, Li C, Huang J, et al. Identification and Validation of a Hypoxia-Related Prognostic Signature in Clear Cell Renal Cell Carcinoma Patients. Med (Baltimore) (2021) 100:e27374. doi: 10.1097/MD.0000000000027374
15. Zhang Z, Li Q, Wang F, Ma B, Meng Y, Zhang Q. Identifying Hypoxia Characteristics to Stratify Prognosis and Assess the Tumor Immune Microenvironment in Renal Cell Carcinoma. Front Genet (2021) 12:606816. doi: 10.3389/fgene.2021.606816
16. Zhong W, Zhong H, Zhang F, Huang C, Lin Y, Huang J. Characterization of Hypoxia-Related Molecular Subtypes in Clear Cell Renal Cell Carcinoma to Aid Immunotherapy and Targeted Therapy via Multi-Omics Analysis. Front Mol Biosci (2021) 8:684050. doi: 10.3389/fmolb.2021.684050
17. Mu W, Jiang D, Mu S, Liang S, Liu Y, Zhang N. Promoting Early Diagnosis and Precise Therapy of Hepatocellular Carcinoma by Glypican-3-Targeted Synergistic Chemo-Photothermal Theranostics. ACS Appl Mater Interf (2019) 11:23591–604. doi: 10.1021/acsami.9b05526
18. Zhang J, Zhang M, Ma H, Song X, He L, Ye X, et al. Overexpression of Glypican-3 is a Predictor of Poor Prognosis in Hepatocellular Carcinoma: An Updated Meta-Analysis. Med (Baltimore) (2018) 97:e11130. doi: 10.1097/MD.0000000000011130
19. Shimizu Y, Suzuki T, Yoshikawa T, Endo I, Nakatsura T. Next-Generation Cancer Immunotherapy Targeting Glypican-3. Front Oncol (2019) 9:248. doi: 10.3389/fonc.2019.00248
20. Valsechi MC, Oliveira AB, Conceicao AL, Stuqui B, Candido NM, Provazzi PJ, et al. GPC3 Reduces Cell Proliferation in Renal Carcinoma Cell Lines. BMC Cancer (2014) 14:631. doi: 10.1186/1471-2407-14-631
21. Tian DW, Wu ZL, Jiang LM, Gao J, Wu CL, Hu HL. KIF5A Promotes Bladder Cancer Proliferation In Vitro and In Vivo. Dis Markers (2019) 2019:4824902. doi: 10.1155/2019/4824902
22. Brenner D, Yilmaz R, Muller K, Grehl T, Petri S, Meyer T, et al. Hot-Spot KIF5A Mutations Cause Familial ALS. Brain (2018) 141:688–97. doi: 10.1093/brain/awx370
23. Borglum AD, Byskov A, Ragno P, Roldan AL, Tripputi P, Cassani G, et al. Assignment of the Urokinase-Type Plasminogen Activator Receptor Gene (PLAUR) to Chromosome 19q13.1-Q13.2. Am J Hum Genet (1992) 50:492–7. doi: 0002-9297/92/5003-0006
24. Gonias SL, Hu J. Urokinase Receptor and Resistance to Targeted Anticancer Agents. Front Pharmacol (2015) 6:154. doi: 10.3389/fphar.2015.00154
25. Gilder AS, Natali L, Van Dyk DM, Zalfa C, Banki MA, Pizzo DP, et al. The Urokinase Receptor Induces a Mesenchymal Gene Expression Signature in Glioblastoma Cells and Promotes Tumor Cell Survival in Neurospheres. Sci Rep (2018) 8:2982. doi: 10.1038/s41598-018-21358-1
26. Tan YQ, Li YT, Yan TF, Xu Y, Liu BH, Yang JA, et al. Six Immune Associated Genes Construct Prognostic Model Evaluate Low-Grade Glioma. Front Immunol (2020) 11:606164. doi: 10.3389/fimmu.2020.606164
27. Wang J, Wang Y, Kong F, Han R, Song W, Chen D, et al. Identification of a Six-Gene Prognostic Signature for Oral Squamous Cell Carcinoma. J Cell Physiol (2020) 235:3056–68. doi: 10.1002/jcp.29210
28. van Haaften-Visser DY, Harakalova M, Mocholi E, van Montfrans JM, Elkadri A, Rieter E, et al. Ankyrin Repeat and Zinc-Finger Domain-Containing 1 Mutations are Associated With Infantile-Onset Inflammatory Bowel Disease. J Biol Chem (2017) 292:7904–20. doi: 10.1074/jbc.M116.772038
29. Guo K, Lai C, Shi J, Tang Z, Liu C, Li K, et al. A Novel Risk Factor Model Based on Glycolysis-Associated Genes for Predicting the Prognosis of Patients With Prostate Cancer. Front Oncol (2021) 11:605810. doi: 10.3389/fonc.2021.605810
30. Xing Q, Zeng T, Liu S, Cheng H, Ma L, Wang Y. A Novel 10 Glycolysis-Related Genes Signature Could Predict Overall Survival for Clear Cell Renal Cell Carcinoma. BMC Cancer (2021) 21:381. doi: 10.1186/s12885-021-08111-0
31. Chen S, Cao G, Wu W, Lu Y, He X, Yang L, et al. Mining Novel Cell Glycolysis Related Gene Markers That Can Predict the Survival of Colon Adenocarcinoma Patients. Biosci Rep (2020) 40(8):BSR20201427. doi: 10.1042/BSR20201427
32. Xu F, Guan Y, Xue L, Huang S, Gao K, Yang Z, et al. The Effect of a Novel Glycolysis-Related Gene Signature on Progression, Prognosis and Immune Microenvironment of Renal Cell Carcinoma. BMC Cancer (2020) 20:1207. doi: 10.1186/s12885-020-07702-7
33. Zhou X, Shang YN, Lu R, Fan CW, Mo XM. High ANKZF1 Expression Is Associated With Poor Overall Survival and Recurrence-Free Survival in Colon Cancer. Future Oncol (2019) 15:2093–106. doi: 10.2217/fon-2018-0920
34. Dittmer J. The Role of the Transcription Factor Ets1 in Carcinoma. Semin Cancer Biol (2015) 35:20–38. doi: 10.1016/j.semcancer.2015.09.010
35. Tetsu O, McCormick F. ETS-Targeted Therapy: Can it Substitute for MEK Inhibitors? Clin Transl Med (2017) 6:16. doi: 10.1186/s40169-017-0147-4
36. Zhai W, Ma J, Zhu R, Xu C, Zhang J, Chen Y, et al. MiR-532-5p Suppresses Renal Cancer Cell Proliferation by Disrupting the ETS1-Mediated Positive Feedback Loop With the KRAS-Nap1l1/P-ERK Axis. Br J Cancer (2018) 119:591–604. doi: 10.1038/s41416-018-0196-5
37. Elhodaky M, Diamond AM. Selenium-Binding Protein 1 in Human Health and Disease. Int J Mol Sci (2018) 19(11):3437. doi: 10.3390/ijms19113437
38. Ha YS, Lee GT, Kim YH, Kwon SY, Choi SH, Kim TH, et al. Decreased Selenium-Binding Protein 1 mRNA Expression Is Associated With Poor Prognosis in Renal Cell Carcinoma. World J Surg Oncol (2014) 12:288. doi: 10.1186/1477-7819-12-288
Keywords: kidney cancer, hypoxia gene, prognosis, immunotherapy, nomogram
Citation: Ning X-h, Li N-y, Qi Y-y, Li S-c, Jia Z-k and Yang J-j (2022) Identification of a Hypoxia-Related Gene Model for Predicting the Prognosis and Formulating the Treatment Strategies in Kidney Renal Clear Cell Carcinoma. Front. Oncol. 11:806264. doi: 10.3389/fonc.2021.806264
Received: 31 October 2021; Accepted: 23 December 2021;
Published: 24 January 2022.
Edited by:
Jonathan Olivier, Université de Lille, FranceReviewed by:
Di Gu, First Affiliated Hospital of Guangzhou Medical University, ChinaPasquale Ditonno, University of Bari, Italy
Copyright © 2022 Ning, Li, Qi, Li, Jia and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiang-hui Ning, bmluZ3hpYW5naHVpQDEyNi5jb20=; Jin-jian Yang, eWFuZ2ppbmppYW4yMDExQDE2My5jb20=