- Primary research
- Open access
- Published:
Epigenetic signature predicts overall survival clear cell renal cell carcinoma
Cancer Cell International volume 20, Article number: 564 (2020)
Abstract
Background
Recently, increasing study have found that DNA methylation plays an important role in tumor, including clear cell renal cell carcinoma (ccRCC).
Methods
We used the DNA methylation dataset of The Cancer Genome Atlas (TCGA) database to construct a 31-CpG-based signature which could accurately predict the overall survival of ccRCC. Meanwhile, we constructed a nomogram to predict the prognosis of patients with ccRCC.
Result
Through LASSO Cox regression analysis, we obtained the 31-CpG-based epigenetic signature which were significantly related to the prognosis of ccRCC. According to the epigenetic signature, patients were divided into two groups with high and low risk, and the predictive value of the epigenetic signature was verified by other two sets. In the training set, hazard ratio (HR) = 13.0, 95% confidence interval (CI) 8.0–21.2, P < 0.0001; testing set: HR = 4.1, CI 2.2–7.7, P < 0.0001; entire set: HR = 7.2, CI 4.9–10.6, P < 0.0001, Moreover, combined with clinical indicators, the prediction of 5-year survival of ccRCC reached an AUC of 0.871.
Conclusions
Our study constructed a 31-CpG-based epigenetic signature that could accurately predicted overall survival of ccRCC and staging progression of ccRCC. At the same time, we constructed a nomogram, which may facilitate the prediction of prognosis for patients with ccRCC.
Background
Renal cell carcinoma (RCC) is a cancer that originates in the renal epithelial cells and accounts for more than 90% of renal cancer, of which clear cell RCC (ccRCC) is most common subtype and causes the most deaths [1]. According to statistics, ccRCC caused 14,400 deaths in 2017 [2]. The TNM staging system is still the most commonly used clinical tool to stratify ccRCC patients, however, it is not enough to reflect the biological heterogeneity of tumor and accurately predict the prognosis of ccRCC patients [3], and a better marker is urgently needed to help us to accurately predict.
The aberrant methylation status of CpG islands located in the promoter region of tumor genes is becoming increasingly important in the search for new potential biomarkers for cancer [4,5,6], because these aberrant are relatively stable and potentially reversible [7]. Combining multiple markers, rather than using only one marker to construct a prognostic model, will make the results more stable and reliable, and improve its predictive value [8]. Many studies have used multiple targets to construct the prognostic model of ccRCC [9, 10]. Here, we used the methylation data of the ccRCC of TCGA database to construct a 31-CpG-based signature.
In this study, we used the methylation data of the ccRCC of TCGA database to construct a 31-CpG-based epigenetic signature, which can accurately predict the overall survival rate of ccRCC through the most recently used The Least Absolute Shrinkage And Selection Operator method (LASSO) algorithm. We then verified this epigenetic marker using two validation sets (testing set and entire set). Furthermore, we constructed a nomogram to facilitate clinicians to accurately predict overall survival in ccRCC patients. Our results identified a new 31-CpG-based epigenetic marker that may be a new target for predicting overall survival of ccRCC.
Methods
Study population and data collection
A total of 319 cases of ccRCC were randomly assigned to two groups, which had a 2:1 ratio. The former was defined as the training set (n = 213) and the latter as testing set (n = 106). The corresponding clinical follow-up information of all cases included the survival time, survival outcome, tumor staging, grading and other information of the cases. The level three of RNA-seq (Illumina RNASeqV2) and the Infinium HumanMethylation450 BeadChip array (Illumina) data were downloaded from The Cancer Genome Atlas (TCGA) database (http://cancergenome.nih.gov/), The details were listed in Table 1, and the flow chart of our entire experiment was shown in Fig. 1.
Differentially methylated positions screening
We first filtered out the probes that have too many zeros, then used the R package “ChAMP” [11] to quality control, remove the batch effects, BMIQ normalization [12] and Single Nucleotide Polymorphism (SNP) filtering (if it located in the CpG dinucleotide) [13]. After that, we used “champ.DMP” function to perform differentially methylated positions (DMP) analysis, we chosen the false discovery rate (FDR) < 0.05 and | Δβ | ≥ 0.2 as the threshold for screening DMPs.
Establishment and validation of epigenetic signature
LASSO cox regression analysis was conducted for the DMPs obtained in the previous step, based on R packet “glmnet”. The degree of LASSO regression complexity adjustment was controlled by the parameter λ, and the larger λ was, the greater the penalty was, so as to obtain a model with fewer variables [14]. We construct the model according to the protocol of official website (https://cran.r-project.org/web/packages/glmnet/glmnet.pdf). Frist, the function “glmnet” returned a sequence of models for us to choose from. Second, we used the function “cv.glmnet” to perform Cross-validation, we followed the protocol, did 100 times cross-validations, and finally got an λ average. Based on this λ value, we constructed an epigenetic signature based on 31 CpGs from training set. The accuracy of the epigenetic signature was then verified in the testing set and the entire set by using the time-dependent ROC curve and Kaplan–Meier survival curve analysis. Time-dependent ROC curve analysis was performed by R package “survivalROC” [15] while the Kaplan–Meier survival curve analysis was performed by R package “survival”.
Cox regression analysis of the epigenetic signature
To investigate whether the epigenetic signature we constructed could be independent of other factors affecting the prognosis of ccRCC, we conducted univariate and multivariate analyses by R package “survival”. First of all, the epigenetic signature score, age, gender, neoadjuvant treatment, lymph node count, histologic stage and pathologic grade were performed the univariate analysis. Then, we put the univariate analysis has significant difference factor to carry on the multivariate analysis. Finally, we used the forest map to visualize our results by R package “forestplot”.
The construction of the nomogram
Based on the results of the previous step, we used the factors with significant differences (epigenetic signature score and histologic stage) in multivariate regression analysis to construct a nomogram by using the R package “rms”. In order to verify the accuracy of the nomogram, we used the calibration curve to evaluate the nomogram. In the calibration curve, if the observed value and the actual value are more coincident, it indicates that the prediction accuracy of the nomogram is higher.
Functional annotation of the epigenetic signature
We calculated the top 500 genes with the highest correlation with epigenetic signature score, and then enriched the functions of these 500 genes, thus indirectly predicting the functional annotation of epigenetic marker score. The 500 genes were uploaded to the DAVID website (https://david.ncifcrf.gov/) for GO analysis and KEGG pathway analysis. The following, the STRING database (https://string-db.org/) was used to performed protein–protein interaction (PPI) network [16] analysis. We’ve picked the top 50 genes from the “cytoHubba” app at Cytoscape and used the Cytoscape to visualize the results. According to the scores of the epigenetic signature, patients were divided into two groups with high and low risk, and then the corresponding RNA-seq data of these two groups were analyzed by single-samples gene-set enrichment analysis. The R package “GSVA” [17] were used to perform ssGSEA for the 31-CpGs-based epigenetic signature. The C2 (c2.cp.kegg.v6.1.symbols.gmt) set was downloaded from the Molecular Signatures Database (http://software.broadinstitute.org/gsea/msigdb/index.jsp). And it was chosen as the signature to perform ssGSEA. The R package “limma” was used to identify the differentially expressed gene sets. The threshold for screening differentially expressed gene sets were the p-value < 0.05.
Results
Screening of DMPs
The R package “ChAMP” was used to screen DMPs between ccRCC and normal ccRCC samples in ccRCC dataset of the TCGA database, a total of 3858 DMPs were identified (1641 up-regulated and 576 down-regulated), under the threshold of FDR < 0.05 and | Δβ | ≥ 0.2. All DMPs were listed in Additional file 1: Table S1.
Establishment and validation of the epigenetic signature
We used the R package of “glmnet” to construct LASSO cox regression model. After 100 cross - validation, we got an optimal λ value (λ = 0.054), and finally we constructed an epigenetic signature based on 31 CpGs. The formula of epigenetic signature score were obtained by calculating the LASSO regression coefficient: epigenetic signature score = \( \sum\nolimits_{k = 1}^{n} {\beta_{ki} X_{ki} } \), here, n is the number of the CpG site; β is the LASSO regression coefficient of the CpG site; X is the methylation value of CpG k and patient i; k is the CpG site. Details of the 31 CpGs were listed in Table 2, the heatmap of the 31 CpGs were shown in Fig. 2a, and the results of ten-time cross-validation were shown in Fig. 2b.
We divided the patients in the training set into high and low risk groups based on the median epigenetic signature scores (here, the median epigenetic signature score was 2.96, Fig. 3a, left panel). Then, time-dependent ROC curve analysis revealed the epigenetic signature we constructed were able to predict the overall survival rate of patients with ccRCC with great accuracy. The AUC of time-dependent ROC for 1-year overall survival was 0.875, 3-year 0.876, 5-year 0.851 and 7-year 0.855 (Fig. 3b middle panel). Finally, Kaplan–Meier survival curve analysis showed that the high-risk group had a significantly worse prognosis than the low-risk group (HR = 13.0, 95% CI 8.0–21.2, P < 0.0001, Fig. 3a right panel). In order to avoid over-fitting effect, testing set (HR = 4.1, CI 2.2–7.7, P < 0.0001, Fig. 3b) and entire set (HR = 7.2, CI 4.9–10.6, P < 0.0001, Fig. 3c) were used for verification our results.
Univariate and multivariate regression analysis of the epigenetic signature
We first performed a univariate analysis of the epigenetic signature scores, age, gender, neoadjuvant treatment, lymph node examination, histologic grade and pathologic stage in the entire set, and the results revealed that only the epigenetic signature score (HR = 3.87, 95% CI 3.02–4.97, P < 0.0001), age (HR = 1.5, 95% CI 1.02–2.21, P = 0.0378), histologic grade (HR = 3.35, 95% CI 2.07–5.41, P < 0.0001) and pathologic stage (HR = 4.49, 95% CI 2.92–6.89, P < 0.0001) have significant prognostic value (Fig. 4a). Next, we conducted multivariate analysis of the factors with significant differences above, and found that the epigenetic signature score (HR = 3.09, 95% CI 2.35–4.08, P < 0.0001) and pathologic stage (HR = 2.42, 95% CI 1.52–3.86, P = 0.0002) could accurately predict the overall survival rate of ccRCC independently of other factors (Fig. 4b, Table 3).
Nomogram creating and calibrating
In order to more conveniently predict the overall survival rate of ccRCC patients for clinicians, we constructed a nomogram based on the R package “rms”. This nomogram could accurately predict 3- or 5-year overall survival in patients with ccRCC (Fig. 5a). Moreover, we used the calibration curve to calibrate the nomogram and found that the prediction value of nomogram was significantly correlated with the actual value (Fig. 5b, c).
Subgroup analysis of the epigenetic signature
To explore the predictive value of epigenetic signature in different subgroups, we performed subgroup analyzed patients from the entire set according to different characteristics: pathologic stage (III/IV vs I/II, Fig. 6a, b), histologic grade (III/IV vs I/II, Fig. 6c, d), lymph node examination (Positive vs negative, Fig. 6e, f) and age (>=65 vs < 65, Fig. 6g, h). The detailed information of subgroup analysis for the epigenetic signature were listed in Table 4.
Functional annotation of the epigenetic signature
Pearson correlation coefficients of all genes and the epigenetic signature scores were calculated. The 500 genes with the highest correlation were selected for subsequent analysis, and the threshold value was set as p-value less than 0.05. Functional enrichment of the 500 genes showed that this epigenetic signature was mainly enriched in the cell adhesion, notch signaling pathway, axon guidance, thyroid hormone signaling pathway and so on (Additional file 2: Figure S1). The PPI network diagram we built was shown in Additional file 3: Figure S2.
After the enrichment scores of each sample were calculated, the differences of the high and low groups were analyzed according to the epigenetic signature scores. A total of 26 differentially expressed gene sets (10 up-regulated and 16 down-regulated) were screened under the threshold FDR < 0.05 and |LogFC| > 0.015. The result of ssGSEA revealed that patients in high risk group mainly enriched in p53, nod like receptor, cytosolic DNA sensing and other signaling pathways. The patients in low risk group were mainly enriched in the fatty acid metabolism pathway, PPAR signaling pathway, renin angiotensin system and so on (Fig. 7c, Additional file 1: Table S2).
Verification of CpG island methylation levels in ccRCC tissues and renal carcinoma cell lines
We collected 4 samples diagnosed as ccRCC and 2 paracancer tissue from the Zhongnan Hospital of Wuhan University, plus human renal tubular epithelial cell line (HK2) and 2 kidney cancer cell lines (ACHN and 769-P), and sequenced them with targeted methylation (MethylTarget, Illumina Hiseq, Illumina, CA, USA). The results showed that CpG island involved in the construction of prognostic model had different expression patterns between carcinoma and paracancer (Fig. 7d).
Discussion
Patients with ccRCC exhibit a highly invasive clinical process, poor prognosis, and high recurrence rate [1]. The TNM staging system is still widely used to predict the prognosis of ccRCC, however, it relies mainly on anatomical information and has no biological characteristics [9]. Therefore, we need a new target to accurately predict the prognosis of ccRCC.
Overfitting often occurs in selecting prognostic models using higher-dimensional data [18], here, we chosen the LASSO algorithm, which could eliminate this disadvantage [19], as the method to select CpG for constructing epigenetic signature. The algorithm of LASSO has been widely used to construct predictive survival models, for instance, colon cancer [20], gastric cancer [21], T-cell lymphoblastic lymphoma [22] and ccRCC [9, 10, 23].
In recent years, epigenetics has been studied more and more in ccRCC, DNA methylation is one of the most studied patterns in epigenetic regulation [24,25,26]. And there have also been a number of studies using methylation data to construct prognostic signatures [5, 6, 27]. Here, we constructed and validated a 31-CpG-based epigenetic signature that could accurately predicted overall survival of ccRCC. Using the epigenetic signature we constructed, we divided patients into high and low risk groups, and high-risk patients had significantly lower overall survival rates in each set. Although previous studies have constructed several molecular biomarkers to predict the survival of ccRCC, few of the signature they have constructed are as accurate as our epigenetic signature. The AUC of time-dependent ROC was all greater than 0.75 in each set. After the epigenetic signature combined with staging and grading, the AUC predicted 5-year overall survival was reached 0.871. At the same time, we validated our epigenetic signature in different subgroups (pathologic stage, histologic grade, lymph node examination, and age) and found that the epigenetic signature showed excellent prognostic value.
We use ssGSEA and other methods to infer the functional annotation of epigenetic signature, and we found that the high risk group patients were mainly enriched in p53 and nod like receptor signaling pathway, which have been reported in previous studies to be associated with the pathogenesis of ccRCC [28,29,30]. The group of low risk were enriched in fatty acid metabolism and which has also been reported to be associated with ccRCC [31,32,33]. We speculate that these pathways may be responsible for the very different overall survival rates between the two groups based on the epigenetic signature. So what we’re going to do is to look at how do these pathways affect the prognosis of ccRCC.
We performed methylation sequencing on ccRCC tissue samples and kidney cancer cell lines, and the results showed that CpG islands involved in the construction of the prognosis model had significantly different expression patterns. However, our study has some limitations, we lack more data to verify it, but we have built a new biological sample database called Biological Repositories Zhongnan Hospital of Wuhan University (http://biobank.znhospital.cn), which will provide a lot of data to verify our epigenetic markers in the future.
In conclusion, we provide an effective epigenetic signature that accurately predicts overall survival of ccRCC independent of other risk factors (age, gender, neoadjuvant treatment, lymph node examination, histologic grade and pathologic stage). At the same time, we also constructed a nomogram for the convenience of clinicians to predict the prognosis of ccRCC.
Conclusion
We developed a 31-CpG-based signature using bioinformatics methods (LASSO cox regression analysis). We found that this epigenetic signature could accurately predict the overall survival rate of ccRCC, contributing to the prediction of the prognosis of ccRCC. And we constructed a nomogram based on this epigenetic signature, which could accurately and conveniently predict the prognosis of ccRCC patients for clinicians.
Availability of data and materials
The datasets analyzed during the current study are available in the Cancer Genome Atlas (TCGA) database (http://cancergenome.nih.gov/).
Abbreviations
- RCC:
-
Renal cell carcinoma
- ccRCC:
-
Clear cell renal cell carcinoma
- LASSO:
-
Least absolute shrinkage and selection operator
- TCGA:
-
The cancer genome atlas
- SNP:
-
Single nucleotide polymorphism
- DMP:
-
Differentially methylated positions
- FDR:
-
False discovery rate
- PPI:
-
Protein–protein interaction
References
Hsieh JJ, Purdue MP, Signoretti S, Swanton C, Albiges L, Schmidinger M, Heng DY, Larkin J, Ficarra V. Renal cell carcinoma. Nat Rev Dis Primers. 2017;3:17009.
Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.
Ljungberg B, Bensalah K, Canfield S, Dabestani S, Hofmann F, Hora M, Kuczyk MA, Lam T, Marconi L, Merseburger AS, et al. EAU guidelines on renal cell carcinoma: 2014 update. Eur Urol. 2015;67(5):913–24.
Esteller M. Epigenetics in cancer. N Engl J Med. 2008;358(11):1148–59.
Shen S, Wang G, Shi Q, Zhang R, Zhao Y, Wei Y, Chen F, Christiani DC. Seven-CpG-based prognostic signature coupled with gene expression predicts survival of oral squamous cell carcinoma. Clin Epigenet. 2017;9:88.
Qiu J, Peng B, Tang Y, Qian Y, Guo P, Li M, Luo J, Chen B, Tang H, Lu C, et al. CpG methylation signature predicts recurrence in early-stage hepatocellular carcinoma: results from a multicenter study. J Clin Oncol. 2017;35(7):734–42.
Garcia-Manero G, Jabbour E, Borthakur G, Faderl S, Estrov Z, Yang H, Maddipoti S, Godley LA, Gabrail N, Berdeja JG, et al. Randomized open-label phase II study of decitabine in patients with low- or intermediate-risk myelodysplastic syndromes. J Clin Oncol. 2013;31(20):2548–53.
Ng SW, Mitchell A, Kennedy JA, Chen WC, McLeod J, Ibrahimova N, Arruda A, Popescu A, Gupta V, Schimmer AD, et al. A 17-gene stemness score for rapid determination of risk in acute leukaemia. Nature. 2016;540(7633):433–7.
Qu L, Wang ZL, Chen Q, Li YM, He HW, Hsieh JJ, Xue S, Wu ZJ, Liu B, Tang H, et al. Prognostic value of a long non-coding RNA signature in localized clear cell renal cell carcinoma. Eur Urol. 2018;74(6):756–63.
Chen L, Luo Y, Wang G, Qian K, Qian G, Wu CL, Dan HC, Wang X, Xiao Y. Prognostic value of a gene signature in clear cell renal cell carcinoma. J Cell Physiol. 2019;234(7):10324–35.
Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, Beck S. ChAMP: 450k chip analysis methylation pipeline. Bioinformatics. 2014;30(3):428–30.
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29(2):189–96.
Sandoval J, Mendez-Gonzalez J, Nadal E, Chen G, Carmona FJ, Sayols S, Moran S, Heyn H, Vizoso M, Gomez A, et al. A prognostic DNA methylation signature for stage I non-small-cell lung cancer. J Clin Oncol. 2013;31(32):4140–7.
Simon N, Friedman J, Hastie T, Tibshirani R. Regularization paths for Cox’s proportional hazards model via coordinate descent. J Stat Softw. 2011;39(5):1–13.
Pencina MJ, D’Agostino RB. Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimation. Stat Med. 2004;23(13):2109–23.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447–52.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinform. 2013;14:7.
Simon R, Altman DG. Statistical aspects of prognostic factor studies in oncology. Br J Cancer. 1994;69(6):979–85.
Goeman JJ. L1 penalized estimation in the Cox proportional hazards model. Biom J. 2010;52(1):70–84.
Zhang JX, Song W, Chen ZH, Wei JH, Liao YJ, Lei J, Hu M, Chen GZ, Liao B, Lu J, et al. Prognostic and predictive value of a microRNA signature in stage II colon cancer: a microRNA expression analysis. Lancet Oncol. 2013;14(13):1295–306.
Jiang Y, Zhang Q, Hu Y, Li T, Yu J, Zhao L, Ye G, Deng H, Mou T, Cai S, et al. ImmunoScore signature: a prognostic and predictive tool in gastric cancer. Ann Surg. 2018;267(3):504–13.
Tian XP, Huang WJ, Huang HQ, Liu YH, Wang L, Zhang X, Lin TY, Rao HL, Li M, Liu F, et al. Prognostic and predictive value of a microRNA signature in adults with T-cell lymphoblastic lymphoma. Leukemia. 2019;33(10):2454–65.
Luo Y, Chen L, Wang G, Xiao Y, Ju L, Wang X. Identification of a three-miRNA signature as a novel potential prognostic biomarker in patients with clear cell renal cell carcinoma. J Cell Biochem. 2019;120(8):13751–64.
Rondinelli B, Rosano D, Antonini E, Frenquelli M, Montanini L, Huang D, Segalla S, Yoshihara K, Amin SB, Lazarevic D, et al. Histone demethylase JARID1C inactivation triggers genomic instability in sporadic renal cancer. J Clin Invest. 2015;125(12):4625–37.
Li Y, Gong Y, Ning X, Peng D, Liu L, He S, Gong K, Zhang C, Li X, Zhou L. Downregulation of CLDN7 due to promoter hypermethylation is associated with human clear cell renal cell carcinoma progression and poor prognosis. J Exp Clin Cancer Res. 2018;37(1):276.
Gooskens SL, Klasson TD, Gremmels H, Logister I, Pieters R, Perlman EJ, Giles RH, van den Heuvel-Eibrink MM. TCF21 hypermethylation regulates renal tumor cell clonogenic proliferation and migration. Mol Oncol. 2018;12(2):166–79.
Geybels MS, Wright JL, Bibikova M, Klotzle B, Fan JB, Zhao S, Feng Z, Ostrander EA, Lin DW, Nelson PS, et al. Epigenetic signature of Gleason score and prostate cancer recurrence after radical prostatectomy. Clin Epigenet. 2016;8:97.
Wang X, Cheng Y, Zhu Y, Li H, Ge W, Wu X, Zhao K, Yuan J, Li Z, Jiang S, et al. Epigenetic silencing of ASPP1 confers 5-FU resistance in clear cell renal cell carcinoma by preventing p53 activation. Int J Cancer. 2017;141(7):1422–33.
Mastropasqua F, Marzano F, Valletti A, Aiello I, Di Tullio G, Morgano A, Liuni S, Ranieri E, Guerrini L, Gasparre G, et al. TRIM8 restores p53 tumour suppressor function by blunting N-MYC activity in chemo-resistant tumours. Mol Cancer. 2017;16(1):67.
Mey L, Jung M, Roos F, Blaheta R, Hegele A, Kinscherf R, Urbschat A. NOD1 and NOD2 of the innate immune system is differently expressed in human clear cell renal cell carcinoma, corresponding healthy renal tissue, its vasculature and primary isolated renal tubular epithelial cells. J Cancer Res Clin Oncol. 2019;145(6):1405–16.
Du W, Zhang L, Brett-Morris A, Aguila B, Kerner J, Hoppel CL, Puchowicz M, Serra D, Herrero L, Rini BI, et al. HIF drives lipid deposition and cancer in ccRCC via repression of fatty acid metabolism. Nat Commun. 2017;8(1):1769.
Wettersten HI, Aboud OA, Lara PN Jr, Weiss RH. Metabolic reprogramming in clear cell renal cell carcinoma. Nat Rev Nephrol. 2017;13(7):410–9.
Theoleyre S, Mottier S, Masson D, Denis MG. HtrA3 is regulated by 15-deoxy-Delta 12,14-prostaglandin J2 independently of PPARgamma in clear cell renal cell carcinomas. Biochem Biophys Res Commun. 2010;394(3):453–8.
Acknowledgements
The excellent technical assistance of Shanshan Zhang and Danni Shan is gratefully acknowledged.
Funding
This work was supported by grants from the Health commission of Hubei Province scientific research project (Grant Number WJ2019H013) and the Fundamental Research Funds for the Central Universities (Grant Number 2042019kf0176).
Author information
Authors and Affiliations
Contributions
YW, XW, YX, and GW conceived and designed the study, YW and LC performed the analysis procedures, YW, LC, LJ, KY, GW, XL and YX analyzed the results, YW, LC and YX contributed analysis tools, YW, LC, XW and YX contributed to the writing of the manuscript. All authors reviewed the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The human tissues used in the experiment were approved by the patients and approved by the Medical Ethics Committee, Zhongnan Hospital of Wuhan University. Ethics Committee Approval (number: 2020102).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1.
Supplementary tables 1–2, Supplementary Information: Related file 1. Ethics Committee Approval (number: 2020102).
Additional file 2.
Supplementary figure S1: Functional annotation of the epigenetic signature.
Additional file 3.
Supplementary figure S2: The PPI network diagram.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, Y., Chen, L., Ju, L. et al. Epigenetic signature predicts overall survival clear cell renal cell carcinoma. Cancer Cell Int 20, 564 (2020). https://doi.org/10.1186/s12935-020-01640-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12935-020-01640-x