Song-Chao Li
Xiang-hui Ning- Department of Urology, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
Kidney cancer is one of the most common urological cancers worldwide, and kidney renal clear cell cancer (KIRC) is the major histologic subtype. Our previous study found that von-Hippel Lindau (VHL) gene mutation, the dominant reason for sporadic KIRC and hereditary kidney cancer-VHL syndrome, could affect VHL disease-related cancers development by inducing telomere shortening. However, the prognosis role of telomere-related genes in kidney cancer has not been well discussed. In this study, we obtained the telomere-related genes (TRGs) from TelNet. We obtained the clinical information and TRGs expression status of kidney cancer patients in The Cancer Genome Atlas (TCGA) database, The International Cancer Genome Consortium (ICGC) database, and the Clinical Proteomic Tumor Analysis Consortium (CPTAC) database. Totally 353 TRGs were differential between tumor and normal tissues in the TCGA-KIRC dataset. The total TCGA cohort was divided into discovery and validation TCGA cohorts and then using univariate cox regression, lasso regression, and multivariate cox regression method to conduct data analysis sequentially, ten TRGs (ISG15, RFC2, TRIM15, NEK6, PRKCQ, ATP1A1, ELOVL3, TUBB2B, PLCL1, NR1H3) risk model had been constructed finally. The kidney patients in the high TRGs risk group represented a worse outcome in the discovery TCGA cohort (p<0.001), and the result was validated by these four cohorts (validation TCGA cohort, total TCGA cohort, ICGC cohort, and CPTAC cohort). In addition, the TRGs risk score is an independent risk factor for kidney cancer in all these five cohorts. And the high TRGs risk group correlated with worse immune subtypes and higher tumor mutation burden in cancer tissues. In addition, the high TRGs risk group might benefit from receiving immune checkpoint inhibitors and targeted therapy agents. Moreover, the proteins NEK6, RF2, and ISG15 were upregulated in tumors both at the RNA and protein levels, while PLCL1 and PRKCQ were downregulated. The other five genes may display the contrary expression status at the RNA and protein levels. In conclusion, we have constructed a telomere-related genes risk model for predicting the outcomes of kidney cancer patients, and the model may be helpful in selecting treatment agents for kidney cancer patients.
Introduction
Kidney cancer is one of the most common urological cancers worldwide, and kidney renal clear cell cancer (KIRC) is the major histologic subtype (1, 2). The classic prognostic factors of kidney cancer are tumor features such as tumor stage, node stage, metastasis stage, and nuclear grade. Several models have been built to predict prognosis by combining the clinical features of cancer patients. These include the International Metastatic RCC Database Consortium (IMDC) risk model, UCLA Integrated Staging System (UISS), and Memorial Sloan-Kettering Cancer Center (MSKCC/Motzer) Score (3–6). Recently, owing to the development of sequencing technology, several molecular signatures have been identified to predict the overall survival of kidney cancer patients. These signatures show a robust predictive effect, and some of them could indicate the potential mechanism of cancer development (7, 8).
Telomeres are regions composed of repetitive TTAGGG DNA sequences and shelterin complex located at the end of chromosomes (9). Telomeres are essential for chromosome stability, and telomere length shortens following cell division and some disease statuses. In addition, abnormalities in the telomere might result in many diseases, such as dyskeratosis congenita, heart disease, mental health problems, and cancer (10, 11). Our previous study found that von-Hippel Lindau (VHL) gene mutation, the dominant reason for sporadic KIRC and hereditary kidney cancer- VHL syndrome, could affect VHL disease-related cancers development by inducing telomere shortening (12).
Moreover, many studies have been conducted to investigate the role of telomeres in the development and progression of cancers. Recent findings suggest that the shortening telomeres can affect the process of cancer development in two different ways. First, telomere shortening might play a tumor-suppressive role by arresting cell proliferation. On the other hand, telomere shortening could also result in extensive genome instability, which promotes cancer progression. In breast cancer, long telomere lengths correlated with a better prognosis (13). A meta-analysis indicated that short telomere length was associated with increased mortality risk and poor prognosis in cancer patients (14). In kidney cancer, telomere length in tumor cells was shorter than in normal kidney cells, but the prognostic role of telomere length remains controversial (15). However, shorter leukocyte telomere length in kidney cancer patients was an independent factor of worse outcomes, especially for stage I cancer patients (16).
Previous studies have focused on the telomere length in cancers and its role in the prognosis of cancers. There has been no study done to investigate telomere-related genes in the prognosis of cancers. Herein, we constructed a risk model using telomere-related genes to predict the prognosis of kidney cancer and then evaluated the potential role of this risk model in selecting treatment agents.
Method
Acquisition of data
Kidney cancer patients’ data files from The Cancer Genome Atlas (TCGA) dataset (TCGA-KIRC cohort), The International Cancer Genome Consortium (ICGC) database (RECA-EU cohort), and Clinical Proteomic Tumor Analysis Consortium (CPTAC) database (PDC000127) were downloaded and processed according to the operational processes of the public data provider. All the data in the TCGA-KIRC cohort were enrolled in this study to screen the differential expression genes. The analysis included the mRNA expression data and exact clinical features in all patients from these three datasets, including the survival time and tumor characteristics. The tumor mutation burden (TMB) data of the TCGA-KIRC cohorts were also downloaded and processed. In addition, the protein expression status of the CPTAC cohort and the immunohistochemical (IHC) data of renal cancer in the Human Protein Atlas (HPA) database were also acquired. In addition, the telomere-related genes were obtained from http://www.cancertelsys.org/telnet/ (17).
Screening and analysis of differential telomere-related genes
The RNA sequence profiles of 533 kidney cancer tissues and 72 adjacent normal kidney tissues were used for screening the differential genes using the limma package (18). Telomere-related genes in all the differential genes were selected and enrolled in weighted gene co-expression network analysis (WGCNA), which was conducted using the WGCNA R package (19). The most differentially expressed telomere-related genes between the tumor and the adjacent normal kidney tissues in the WCGNA results were included in the construction of the prognosis model.
Construction and verification of telomere-related genes risk model
The TCGA-KIRC cohort was randomly divided into the discovery TCGA cohort and the validation cohort in a 1:1 ratio. The discovery TCGA cohort was analyzed sequentially through univariate cox regression, Lass regression, and multivariate cox regression. Subsequently, the risk model was constructed based on the genes’ expression value and the coefficient, which were acquired in the multivariate cox regression using the Akaike Information Criterion (AIC) method. The risk score of each patient was calculated in these five cohorts: discovery TCGA cohort, validation TCGA cohort, total TCGA cohort, ICGC cohort, and CPTAC cohort. Each cohort was divided into two groups, the high-risk group, and the low-risk group, according to the median risk score, and the prognosis of each risk group was examined using the log-rank test. The risk score was also evaluated using univariate and multivariate cox regression to determine its role in predicting the overall survival of kidney cancer patients in the five cohorts. In addition, the receiver operating characteristic (ROC) curve was used to check the accuracy of the risk model in predicting the prognosis. Finally, a nomogram was constructed to predict the 1-, 3-, and 5-year survival rates using the risk score based on the total TCGA cohort.
Immune features and immune subtype analysis
The infiltration of the different types of immune cells in the kidney cancer tumor microenvironment was assessed by the specific genes’ expression value using CIBERSORT iterated 1000 times (https://cibersort.stanford.edu/) (20). The immune subtypes of pan-cancers in the TCGA database have been described in a previous study, and the immune subtypes of KIRC patients were obtained (21). Additionally, TIDE (Tumor ImmuneDysfunction and Exclusion) score was calculated online following the instructions (https://tide.dfci.harvard.edu/) (22).
IC50 prediction of the different targeted therapy agents
The targeted drugs’ half-maximal inhibitory concentrations (IC50) were predicted using the gene expression level to reflect the treatment sensitivity. This was done using the R package named “pRRophetic” (23, 24).
Cell culture, RNA Extraction and quantitative real-time PCR
Totally seven cell lines, including a normal human renal proximal tubular cell line (HK2) and six RCC cell lines (Caki-1, 769-P, OSRC-2, ACHN, 786-O, A498) were used to investigate the telomere-related gene expression status. Among these cell lines, 769-P, A498, 786-O, and OSRC-2 were cultured in RPMI 1640 medium supplemented with 10% fetal bovine serum. While HK2, Caki-1, and ACHN were cultured in DMEM medium with 10% fetal bovine serum. TRIZOL reagent (Invitrogen, USA) was used to extract the total RNA, and the reverse transcription reactions were carried out using a qPCR RT Kit (TOYOBO Life Science, Shanghai, China). Expression levels of telomere-related genes were detected by qRT-PCR). The human beta-actin gene was used as a reference gene. The primers were shown in Table 1 (25). The experiment was conducted on a Bio-Rad S1000 machine and using an SYBR Green RT-PCR Master Mix reagent (TOYOBO). TRGs’ Relative expression value was computed using the 2−ΔΔCt method and normalized with the beta-actin (26).
Table 1 The primers used in qRT-PCR of the TRG risk model genes.
Protein expression analysis and statistical analysis
The expression of proteins encoded by the genes included in the telomere-related genes model in the CPTAC cohort was compared using the limma package. In addition, the protein expression status of renal cancer and kidney tissues in the HPA database was accessed and analyzed using the “HPAanalyze” package (27). Continuous variables were compared by independent t-test, while the categorical variables were analyzed using the chi-square test. A P-value < 0.05 (two sides) was considered statistically significant.
Results
Analysis of the differential telomere-related genes
The detailed characteristics of the three cohorts are summarized in the Supplementary Table. A total of 2086 telomere-related genes were obtained from TelNet. Among these, 353 genes were differentially expressed between kidney cancer and adjacent normal kidney tissues. In particular, 234 genes were upregulated, while 119 were downregulated in the tumor tissues (Supplementary Figure 1). In addition, these differential telomere-related genes were analyzed using the WCGNA method. The results showed that these genes were clustered into six models, MEgreen, MEturquoise, ME brown, MEblue, MEyellow, and MEgrey modules. Among these modules, the ME-turquoise module, composed of 118 genes, displayed the most significant difference between normal kidney and kidney cancer tissues (R2 = 0.84, P<0.001, Figure 1).
Figure 1 WCGNA analysis of differential expressed telomere-related genes. The differential expressed telomere-related genes were clustered into six modules (A), and the ME-turquoise module shows the most significant difference (B).
Telomere-related genes risk model could predict the prognosis of kidney cancer patients
Among the 353 differential expression genes, 47 genes correlated with kidney cancer patients’ overall survival in the discovery TCGA cohort. Then 19 genes in these 47 genes were finally screened by Lasso regression. Finally, ten telomere-related genes (ISG15, RFC2, TRIM15, NEK6, PRKCQ, ATP1A1, ELOVL3, TUBB2B, PLCL1, NR1H3) in these 19 genes were identified as independent risk factors through the multivariate cox regression and were used to format the risk model, telomere related genes (TRGs) risk model. And the TRGs risk score formula was as follows: risk score=0.005134398*ISG15 +0.047733021*RFC2-0.071026696*TRIM15-0.0203981*NEK6-0.133649782*PRKCQ-0.007396692*A TP1A1+0.700714514*ELOVL3+0.206037577*TUBB2B-0.148322187*PLCL1+0.070770091*NR1H3. The risk score of patients in five cohorts was computed, and the patients in each cohort were divided into a high and low-risk group according to each cohort’s median risk score value. The patients in the high-risk group represent the worse outcomes in all these five cohorts, discovery TCGA cohort (p<0.001), validation TCGA cohort (p<0.001), total TCGA cohort (p<0.001), ICGC cohort (p=0.018) and CPTAC cohort (p=0.003) (Figures 2, 3). In addition, the TRGs risk score proved to be an independent prognostic factor for kidney cancer patients in these five cohorts (Figure 4 and Table 2). Moreover, a nomogram, which consisted of patients’ age, clinical stage, and TRGs risk score, was constructed to individually predict kidney cancer patients’ overall survival (Figure 5).
Figure 2 The distribution of survival status and TRGs risk score in these six cohorts. The patients were ordered according to the TRGs risk score, shown in the up panel, and the survival status of each patient with a different risk score was shown in the middle panel. The TRGs risk model gene expression value has presented in the lower panel. (A) discovery TCGA cohort; (B) validation TCGA cohort; (C) total TCGA cohort; (D) ICGC cohort; (E) CPTAC cohort.
Figure 3 The differential prognosis in patients with different TRGs risk groups. Log-rank methods had compared the survival status in different TRGs risk groups, the patients in the high TRGs risk group presented worse outcomes in these five cohorts. (A) discovery TCGA cohort; (B) validation TCGA cohort; (C) total TCGA cohort; (D) ICGC cohort; (E) CPTAC cohort.
Figure 4 TRGs risk model as an independent prognosis factor in these five cohorts. TRGs risk score, clinical stage, and age were the independent factors in the discovery TCGA cohort (A), validation TCGA cohort (B), and total TCGA cohort (C). Clinical stage and TRGs risk score were the independent factors in the ICGC cohort (D) and CPTAC cohort (E). The ROC curve and AUC value were also shown in each cohort. (*represents p<0.01,**represents p<0.01, *** represents p<0.001).
Table 2 Univariate and Multivariate Cox regression reveals TRGs risk score is an independent risk factor in kidney cancer.
Figure 5 A nomogram was constructed for predicting the overall survival of kidney cancer patients. TRGs risk score, clinical stage, and age composed the nomogram using the total TCGA cohort, and the risk score has the most contribution to the total points.
The patients with high-risk scores present high TMB
The patients in the high-risk group showed a higher TMB (P=0.0023) than those in the low-risk group (Figure 6A). Furthermore, the risk score correlated with the TMB (R=0.21, P=0.00019, Figure 6B).
Figure 6 The tumor mutation burden (TMB) correlated with the TRGs risk score and groups. The TMB of patients in the high TRGs risk group (red points) was higher than those in the low TRGs group (blue points). (A) The patients’ TMB positively correlated with TRGs risk score (p=0.00019). (B).
The patients in different risk groups show different immune status
Different types of immune cells exhibited different infiltration rates in the tumor microenvironment between the high and low-risk groups. Plasma cells, T cells follicular helper, T cells regulatory (Tregs), and macrophages M0 had a higher rate of infiltration in the high-risk group than in the low-risk group. T cells CD4 memory resting, monocytes, macrophages M2, dendritic cells resting, and mast cells resting showed less infiltration in the high-risk group than in the low-risk group (Figure 7A).
Figure 7 The immune cells and subtypes of the TRGs risk group. The immune cells distribution portion in the cancer tissues with different TRGs risk group patients. (*represents p<0.01, **represents p<0.01, *** represents p<0.001). (A) Four hundred of ninety-four patients in the total TCGA cohort has the immune subtype results, no C5 subtype in these patients, and C3 is the most common subtype in the two TRGs risk group. (B).
Previous studies have clustered the tumors in the TCGA database into six subtypes according to the immune status: C1 (wound healing), C2 (IFN-g dominant), C3 (inflammatory), C4 (lymphocyte depleted), C5 (immunologically quiet), and C6 (TGF-b dominant) (21). The immune subtype analysis results indicated the patients in the high-risk group had a higher proportion of C1(2%), C2(7%), C4(8%), and C6(5%) subtype than the patients in the low-risk group, and a lower proportion of C3(77%) subtype than the patients in the low-risk group (97%) (Figure 7B).
Tumor-related genes risk model could be used in choosing a treatment strategy
The TIDE score was significantly higher in the high-risk group than in the low-risk group (P<0.001, Figure 8A). There were no significant differences in the MSI and exclusion score (Figures 8B, C). The IC50 values of Axitinib, Sorafenib, Sunitinib, and Temsirolimus in the high-risk group were lower than in the low-risk group. In contrast, the IC50 value of Pazopanib was higher in the high-risk group than in the low-risk group (Figures 8D–H).
Figure 8 The TIDE score and targeted agents’ treatment sensitivity. Tumor tissues in the high TRGs risk group present the high TIDE score (A), and the MSI and Expulsion score has no significant difference in these two groups (*** represents p<0.001, ns represents no significant difference). (B, C). The treatment response of targeted agents, axitinib (D), sorafenib (E), sunitinib (F), pazopanib (G), and temsirolimus (H).
RNA and protein expression status of the tumor-related genes used in the model
Among the ten tumor-related genes used in the risk model, TRIM15, NEK6, ISG15, RFC2, and NR1H3 were upregulated. At the same time, PLCL1, ATP1A1, TUBB2B, PRKCQ, and ELOVL3 were downregulated in tumors in the TCGA-KIRC dataset at the transcription level (Figure 9A). At the protein level, seven proteins were notable in the tumors of the CPTAC cohort. ATP1A1, ISG15, NEK6, and RFC2 were significantly upregulated, while PRKCQ, PLCL1, and TRIM15 were downregulated (Figure 9B). In addition, eight proteins had IHC results in the HPA database. Though the exact differences between normal kidney tissue and tumors cannot be thoroughly evaluated, the primary data showed that ISG15 might be upregulated in tumors. At the same time, NEK6 and NR1H3 might be downregulated in tumor tissues. No significant differences were found in the expression of RFC2, ATP1A1, ELOVL3, TUBB2B, and PLCL1 (Figure 9C). In addition, these genes’ relative expression status was assessed in RCC cell lines by qRT-PCR, and the results show a significant difference in these cell lines (Figure 10).
Figure 9 The RNA and protein expression of TRGs model genes. The heatmap shows these ten TRGs model gene expression in kidney cancer and adjacent normal kidney tissues. (A) The seven proteins expression, assessed by Tandem Mass Tag (TMT) 10 method, compares the kidney cancer and adjacent normal kidney tissues in the CPTAC cohort. The thin lines represent the maximum and minimum values, and the thick line represents the median value. (B) The detection of eight proteins by IHC in the HPA database. (C).
Figure 10 The TRGs model genes’ relative expression status in RCC cell lines. These 10 genes’ relative expression was assessed by qRT-PCR in HK2 and six RCC cell lines (Caki-1, 769-P, OSRC-2, ACHN, 786-O, A498). All these genes’ expressions in RCC cell lines were compared to their expressions in the HK2 cell line. The NEK6 gene expression in 769-P was set to 2, and the true value was more than 2 (Supplementary File). The headline of this figure represents the TRGs model genes name. (*represents p<0.01, **represents p<0.01).
Discussion
Telomeres play an essential role in the development of kidney cancer. To the best of our knowledge, this study is the first to assess the role of telomere-related genes in the prognosis of renal cancer. We established a prognostic model based on telomere-related genes using public databases. We found that this model can serve as a basis for selecting therapeutic drugs for renal cancer.
Previous studies show the varied prognosis of the tumors in different immune subtypes, the C3 subtype presented with the most favorable outcome, while the C6 and C4 subtypes had the worst prognosis. Our results suggested that the telomere-related genes risk model group was significantly correlated with the immune subtypes. The high-risk group had a higher proportion of C1, C2, C4, and C6 subtypes and a lower proportion of the C3 subtype than the low-risk group. The immune subtype analysis indicated that the patients in the high-risk group might have had a poor prognosis due to the higher portion of subtypes (C1+C2+C4+C6, 23%) with a worse outcome than in the low-risk group (3%). In addition, the high-risk group showed a high TIDE score, indicating that the high-risk group patients may be more likely to experience an immune escape. Overall, the results suggest that the differences in prognosis between the high and low-risk groups might be partly due to the different immune statuses of the patients.
High TMB has been identified as a biomarker to predict the potential benefit of immune checkpoint blockade (ICB) therapy. Our study showed that the high-risk group was associated with a high TMB (28). The drug sensitivity prediction analysis results showed that the high-risk group patients might be more sensitive to Axitinib, Sorafenib, Sunitinib, and Temsirolimus treatment. In contrast, those in the low-risk group might benefit from Pazopanib treatment. Taken together, the patients with a high-risk score might be more suitable candidates to receive immune checkpoint inhibitors and targeted therapy agents. In contrast, the patients with a low-risk score had a limited selection of treatment.
The genes in the telomere-related genes risk model play varying roles in disease. ATP1A1 encodes the α-1 isoform of the familiar Na+/K+ATPase. Studies have found that ATP1A1 mutations could cause aldosterone-producing adenoma (APA). Inhibiting ATP1A1 expression in pancreatic ductal adenocarcinoma (PDAC) cells can suppress the tumor invasion. In kidney cancer, the promoter methylation rate of ATP1A1 was about 15.8% (29–32). Replication Factor C Subunit 2 (RFC2) gene encodes subunit 2 of the RFC complex, a primer recognition factor of DNA polymerase. Studies have shown that RFC2 could regulate the cell cycle and DNA replication process to promote liver cancer development and could act as an oncogene in the progression of lower-grade gliomas (33, 34). Tripartite Motif Containing 15 (TRIM15) encodes a member of the Tripartite motif family, which could exhibit E3 ubiquitin ligase activity. TRIM15 could regulate the Wnt/β-catenin signaling pathway and Keap1-Nrf2 pathway and mediate the ubiquitination of APOA1 and ERK to promote the development and progression of cancers such as non-small cell lung cancer, esophageal squamous cell carcinoma, and pancreatic cancer (35–37). NIMA Related Kinase 6 (NEK6) encodes a kinase, which plays multiple roles in the tumor, including suppression of tumor cell senescence and facilitation of breast cancer cell proliferation. Moreover, NEK6 could participate in the development of castration resistance in prostate cancer. In kidney cancer, studies revealed that miR-141-3p regulated NEK6 to influence cell proliferation, migration, and apoptosis of tumor cells (38–41). Protein Kinase C Theta (PRKCQ) encodes a serine/threonine kinase. It can regulate the immune system by controlling T cells’ activation, survival, and differentiation. PRKCQ affects different processes in cancers, including tumor cell proliferation, migration, and invasion (42–44). IFN-stimulated gene factor 15 (ISG15) encodes a protein induced by type I IFNs, and it is a ubiquitin-like protein. ISG15 exists in two forms in vivo: conjugated to other proteins or free protein. ISG15 and its conjugation affect the progression and treatment response in different cancers. In kidney cancer, ISG15 has been investigated as a novel protein adjuvant in vaccines (45–48). ELOVL Fatty Acid Elongase 3 (ELOVL3) encodes a long-chain fatty acid elongase protein expressed in the liver and brown adipose tissues. Previous studies have shown that ELOVL3 could predict the prognosis of cancers (49–51). Tubulin Beta 2B Class IIb (TUBB2B) encodes a beta isoform of tubulin. The protein could bind GTP and is the major component of microtubules. TUBB2B has been shown to participate in the construction of a prognostic model of kidney cancer (52). Phospholipase C Like 1 (PLCL1) encodes a protein that was predicted to enable phospholipase C activity. PLCL1 was downregulated in kidney cancer tissues and correlated with a poor prognosis. In addition, PLCL1 could repress the progression of kidney cancer through UCP1-mediated lipid browning (53). Nuclear Receptor Subfamily 1 Group H Member 3 (NR1H3), also known as LXR-A, encodes a protein that belongs to the nuclear receptor superfamily. Elevated LXRα expression correlates with a high tumor stage, histologic grade, and pathologic stage of ccRCC, and this could regulate ccRCC cell migration and invasion (54, 55). In our study, NEK6, RF2, and ISG15 were upregulated in tumors both at the RNA and protein levels, while PLCL1 and PRKCQ were downregulated. The other five genes may display the contrary expression status at the RNA and protein levels. The different expression status reflects the different transcriptional, and post-transcriptional mechanisms of these genes in kidney cancer cells, the detailed mechanism and role of these risk model genes should be investigated further.
Our study has some limitations. Though we used three databases to construct and validate this ten genes TRGs risk model, we cannot carry out a single trial to verify our findings due to the long translation and follow-up circle and the high cost (56, 57). In addition, as our results were based on the transcriptomics profile, which limited the clinical application and promotion of the TRGs risk model, a further simple and easy method should be developed. And the functions of the TRGs risk model genes in kidney cancer should be clarified by more basic experiments in further.
In this study, we have constructed a telomere-related genes risk model using the TCGA dataset and carefully verified the model in two datasets (CPTAC and ICGC). We have discussed the protein and RNA expression status of these genes in kidney cancer; however, there were no experiments to validate our findings. The model we have constructed may be helpful in the selection of treatment agents for kidney cancer patients.
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
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
Conceived and designed the experiment: X-HN, S-CL. Performed the experiments and analyzed the data: X-HN, S-CL, Z-KJ. Interpretation of the findings: X-HN, S-CL, Z-KJ and J-JY. 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 number 1901004), 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 manuscript preparation.
Acknowledgments
The results shown here are in whole or part, based upon data generated by the TCGA Research Network (https://www.cancer.gov/tcga), ICGC database (https://dcc.icgc.org/), and CPTAC (https://proteomic.datacommons.cancer.gov/pdc/study/PDC000127).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.975057/full#supplementary-material
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71:209–49. doi: 10.3322/caac.21660
2. Capitanio U, Montorsi F. Renal cancer. Lancet (2016) 387:894–906. doi: 10.1016/S0140-6736(15)00046-X
3. Zigeuner R, Hutterer G, Chromecki T, Imamovic A, Kampel-Kettner K, Rehak P, et al. External validation of the Mayo clinic stage, size, grade, and necrosis (SSIGN) score for clear-cell renal cell carcinoma in a single European centre applying routine pathology. Eur Urol (2010) 57:102–9. doi: 10.1016/j.eururo.2008.11.033
4. 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
5. Motzer RJ, Bacik J, Murphy BA, Russo P, Mazumdar M. Interferon-alfa as a comparative treatment for clinical trials of new therapies against advanced renal cell carcinoma. J Clin Oncol (2002) 20:289–96. doi: 10.1200/JCO.2002.20.1.289
6. Zisman A, Pantuck AJ, Wieder J, Chao DH, Dorey F, Said JW, et al. Risk group assessment and clinical outcome algorithm to predict the natural history of patients with surgically resected renal cell carcinoma. J Clin Oncol (2002) 20:4559–66. doi: 10.1200/JCO.2002.05.111
7. Ning XH, Li NY, Qi YY, Li SC, Jia ZK, Yang JJ. Identification of a hypoxia-related gene model for predicting the prognosis and formulating the treatment strategies in kidney renal clear cell carcinoma. Front Oncol (2021) 11:806264. doi: 10.3389/fonc.2021.806264
8. Sun Z, Tao W, Guo X, Jing C, Zhang M, Wang Z, et al. Construction of a lactate-related prognostic signature for predicting prognosis, tumor microenvironment, and immune response in kidney renal clear cell carcinoma. Front Immunol (2022) 13:818984. doi: 10.3389/fimmu.2022.818984
9. Shay JW, Wright WE. Telomeres and telomerase: three decades of progress. Nat Rev Genet (2019) 20:299–309. doi: 10.1038/s41576-019-0099-1
10. Savage SA. Beginning at the ends: telomeres and human disease. F1000Res (2018) 7:524. doi: 10.12688/f1000research.14068.1
11. Maciejowski J, de Lange T. Telomeres in cancer: tumour suppression and genome instability. Nat Rev Mol Cell Biol (2017) 18:175–86. doi: 10.1038/nrm.2016.171
12. Ning XH, Zhang N, Li T, Wu PJ, Wang X, Li XY, et al. Telomere shortening is associated with genetic anticipation in Chinese Von hippel-lindau disease families. Cancer Res (2014) 74:3802–9. doi: 10.1158/0008-5472.CAN-14-0024
13. Ennour-Idrissi K, Maunsell E, Diorio C. Telomere length and breast cancer prognosis: A systematic review. Cancer Epidemiol Biomarkers Prev (2017) 26:3–10. doi: 10.1158/1055-9965.EPI-16-0343
14. Zhang C, Chen X, Li L, Zhou Y, Wang C, Hou S. The association between telomere length and cancer prognosis: Evidence from a meta-analysis. PloS One (2015) 10:e0133174. doi: 10.1371/journal.pone.0133174
15. Morais M, Dias F, Teixeira AL, Medeiros R. Telomere length in renal cell carcinoma: The Jekyll and Hyde biomarker of ageing of the kidney. Cancer Manag Res (2020) 12:1669–79. doi: 10.2147/CMAR.S211225
16. Callahan CL, Schwartz K, Ruterbusch JJ, Shuch B, Graubard BI, Lan Q, et al. Leukocyte telomere length and renal cell carcinoma survival in two studies. Br J Cancer (2017) 117:752–5. doi: 10.1038/bjc.2017.237
17. Braun DM, Chung I, Kepper N, Deeg KI, Rippe K. TelNet - a database for human and yeast genes involved in telomere maintenance. BMC Genet (2018) 19:32. doi: 10.1186/s12863-018-0617-8
18. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43:e47. doi: 10.1093/nar/gkv007
19. Langfelder P, Horvath S. WGCNA: an r package for weighted correlation network analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559
20. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods (2015) 12:453–7. doi: 10.1038/nmeth.3337
21. Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of cancer. Immunity (2018) 48:812–830 e14. doi: 10.1016/j.immuni.2018.03.023
22. 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
23. Geeleher P, Cox N, Huang RS. pRRophetic: an r package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PloS One (2014) 9:e107468. doi: 10.1371/journal.pone.0107468
24. Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of drug sensitivity in cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res (2013) 41:D955–61. doi: 10.1093/nar/gks1111
25. Wang X, Spandidos A, Wang H, Seed B. PrimerBank: a PCR primer database for quantitative gene expression analysis, 2012 update. Nucleic Acids Res (2012) 40:D1144–9. doi: 10.1093/nar/gkr1013
26. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods (2001) 25:402–8. doi: 10.1006/meth.2001.1262
27. Tran AN, Dussaq AM, Kennell T Jr., Willey CD, Hjelmeland AB. HPAanalyze: an r package that facilitates the retrieval and analysis of the human protein atlas data. BMC Bioinf (2019) 20:463. doi: 10.1186/s12859-019-3059-z
28. Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol (2019) 30:44–56. doi: 10.1093/annonc/mdy495
29. Funder JW. Primary aldosteronism. Hypertension (2019) 74:458–66. doi: 10.1161/HYPERTENSIONAHA.119.12935
30. Beuschlein F, Boulkroun S, Osswald A, Wieland T, Nielsen HN, Lichtenauer UD, et al. Somatic mutations in ATP1A1 and ATP2B3 lead to aldosterone-producing adenomas and secondary hypertension. Nat Genet (2013) 45:440–4, 444e1-2. doi: 10.1038/ng.2550
31. Chen YI, Chang CC, Hsu MF, Jeng YM, Tien YW, Chang MC, et al. Homophilic ATP1A1 binding induces activin a secretion to promote EMT of tumor cells and myofibroblast activation. Nat Commun (2022) 13:2945. doi: 10.1038/s41467-022-30638-4
32. Deckers IA, van Engeland M, van den Brandt PA, Van Neste L, Soetekouw PM, Aarts MJ, et al. Promoter CpG island methylation in ion transport mechanisms and associated dietary intakes jointly influence the risk of clear-cell renal cell cancer. Int J Epidemiol (2017) 46:622–31. doi: 10.1093/ije/dyw266
33. Ji Z, Li J, Wang J. Up-regulated RFC2 predicts unfavorable progression in hepatocellular carcinoma. Hereditas (2021) 158:17. doi: 10.1186/s41065-021-00179-9
34. Zhao X, Wang Y, Li J, Qu F, Fu X, Liu S, et al. RFC2: a prognosis biomarker correlated with the immune signature in diffuse lower-grade gliomas. Sci Rep (2022) 12:3122. doi: 10.1038/s41598-022-06197-5
35. Sun Y, Ren D, Yang C, Yang W, Zhao J, Zhou Y, et al. TRIM15 promotes the invasion and metastasis of pancreatic cancer cells by mediating APOA1 ubiquitination and degradation. Biochim Biophys Acta Mol Basis Dis (2021) 1867:166213. doi: 10.1016/j.bbadis.2021.166213
36. Liang M, Wang L, Sun Z, Chen X, Wang H, Qin L, et al. E3 ligase TRIM15 facilitates non-small cell lung cancer progression through mediating Keap1-Nrf2 signaling pathway. Cell Commun Signal (2022) 20:62. doi: 10.1186/s12964-022-00875-7
37. Zhu G, Herlyn M, Yang X. TRIM15 and CYLD regulate ERK activation via lysine-63-linked polyubiquitination. Nat Cell Biol (2021) 23:978–91. doi: 10.1038/s41556-021-00732-8
38. Liu Y, Fu W, Yin F, Xia L, Zhang Y, Wang B, et al. miR-141-3p suppresses development of clear cell renal cell carcinoma by regulating NEK6. Anticancer Drugs (2022) 33:e125–33. doi: 10.1097/CAD.0000000000001158
39. He Z, Ni X, Xia L, Shao Z. Overexpression of NIMA-related kinase 6 (NEK6) contributes to malignant growth and dismal prognosis in human breast cancer. Pathol Res Pract (2018) 214:1648–54. doi: 10.1016/j.prp.2018.07.030
40. Choudhury AD, Schinzel AC, Cotter MB, Lis RT, Labella K, Lock YJ, et al. Castration resistance in prostate cancer is mediated by the kinase NEK6. Cancer Res (2017) 77:753–65. doi: 10.1158/0008-5472.CAN-16-0455
41. Jee HJ, Kim HJ, Kim AJ, Song N, Kim M, Lee HJ, et al. The inhibition of Nek6 function sensitizes human cancer cells to premature senescence upon serum reduction or anticancer drug treatment. Cancer Lett (2013) 335:175–82. doi: 10.1016/j.canlet.2013.02.012
42. Nicolle A, Zhang Y, Belguise K. The emerging function of PKCtheta in cancer. Biomolecules (2021) 11(2):221. doi: 10.3390/biom11020221
43. Hayashi K, Altman A. Protein kinase c theta (PKCtheta): a key player in T cell life and death. Pharmacol Res (2007) 55:537–44. doi: 10.1016/j.phrs.2007.04.009
44. Marrocco V, Fiore P, Madaro L, Crupi A, Lozanoska-Ochser B, Bouche M. Targeting PKCtheta in skeletal muscle and muscle diseases: good or bad? Biochem Soc Trans (2014) 42:1550–5. doi: 10.1042/BST20140207
45. Mirzalieva O, Juncker M, Schwartzenburg J, Desai S. ISG15 and ISGylation in human diseases. Cells (2022) 11(3):538. doi: 10.3390/cells11030538
46. Hermann M, Bogunovic D. ISG15: In sickness and in health. Trends Immunol (2017) 38:79–93. doi: 10.1016/j.it.2016.11.001
47. Han HG, Moon HW, Jeon YJ. ISG15 in cancer: Beyond ubiquitin-like protein. Cancer Lett (2018) 438:52–62. doi: 10.1016/j.canlet.2018.09.007
48. Wood LM, Pan ZK, Guirnalda P, Tsai P, Seavey M, Paterson Y. Targeting tumor vasculature with novel listeria-based vaccines directed against CD105. Cancer Immunol Immunother (2011) 60:931–42. doi: 10.1007/s00262-011-1002-x
49. Westerberg R, Mansson JE, Golozoubova V, Shabalina IG, Backlund EC, Tvrdik P, et al. ELOVL3 is an important component for early onset of lipid recruitment in brown adipose tissue. J Biol Chem (2006) 281:4958–68. doi: 10.1074/jbc.M511588200
50. Anzulovich A, Mir A, Brewer M, Ferreyra G, Vinson C, Baler R. Elovl3: a model gene to dissect homeostatic links between the circadian clock and nutritional status. J Lipid Res (2006) 47:2690–700. doi: 10.1194/jlr.M600230-JLR200
51. Yuan C, Yuan M, Chen M, Ouyang J, Tan W, Dai F, et al. Prognostic implication of a novel metabolism-related gene signature in hepatocellular carcinoma. Front Oncol (2021) 11:666199. doi: 10.3389/fonc.2021.666199
52. Chen L, Xiang Z, Chen X, Zhu X, Peng X. A seven-gene signature model predicts overall survival in kidney renal clear cell carcinoma. Hereditas (2020) 157:38. doi: 10.1186/s41065-020-00152-y
53. Xiong Z, Xiao W, Bao L, Xiong W, Xiao H, Qu Y, et al. Tumor cell "Slimming" regulates tumor progression through PLCL1/UCP1-mediated lipid browning. Adv Sci (Weinh) (2019) 6:1801862. doi: 10.1002/advs.201801862
54. Zhao L, Lei W, Deng C, Wu Z, Sun M, Jin Z, et al. The roles of liver X receptor alpha in inflammation and inflammation-associated diseases. J Cell Physiol (2021) 236:4807–28. doi: 10.1002/jcp.30204
55. Wang K, Xu T, Ruan H, Xiao H, Liu J, Song Z, et al. LXRalpha promotes cell metastasis by regulating the NLRP3 inflammasome in renal cell carcinoma. Cell Death Dis (2019) 10:159. doi: 10.1038/s41419-019-1345-3
56. Chang X, Lv YF, He J, Cao Y, Li CQ, Guo QN. Gene expression profile and prognostic value of m6A RNA methylation regulators in hepatocellular carcinoma. J Hepatocell Carcinoma (2021) 8:85–101. doi: 10.2147/JHC.S296438
Keywords: kidney cancer, telomere, prognosis, tumor mutation burden, immuno therapy
Citation: Li S-C, Jia Z-K, Yang J-J and Ning X-h (2022) Telomere-related gene risk model for prognosis and drug treatment efficiency prediction in kidney cancer. Front. Immunol. 13:975057. doi: 10.3389/fimmu.2022.975057
Received: 21 June 2022; Accepted: 29 August 2022;
Published: 16 September 2022.
Edited by:
Zhangqun Ye, Huazhong University of Science and Technology, ChinaReviewed by:
Xu-jie Zhou, First Hospital, Peking University, ChinaJing He, Guangzhou Medical University, China
Copyright © 2022 Li, Jia, Yang and Ning. 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=