Development and validation of the prognostic value of the immune-related genes in clear cell renal cell carcinoma

Background: Clear cell renal cell carcinoma (ccRCC) is a highly heterogeneous tumor, resulting a challenge of developing target therapeutics. Not long ago, immune checkpoint blockade regimens combine with tyrosin kinase inhibitors have evolved frontline options in metastatic RCC, which implies arrival of the era of tumor immunotherapy. Studies have demonstrated immune-related genes (IRGs) could characterize tumor milieu and related to patient survival. Nevertheless, the clinical significance of classifier depending on IRGs in ccRCC has not been well established. Methods: The R package limma, univariate and LASSO cox regression analysis were used to screen the prognostic related IRGs from TCGA database. Multivariate cox regression was utilized to establish a risk prediction model for candidate genes. Quantitative real-time PCR was used to confirm the expression of candidates in clinical samples from our institution. CIBERSORT algorithm and correlation analysis were applied to explore tumor-infiltrating immune cells signature between different risk groups. A clinical nomogram was also developed to predict OS by using the rms R package based on the risk prediction model and other independent risk factors. The ICGC data was used for external validation of either gene risk model or nomogram. Results: We identified 382 differentially expressed immune related genes. Four unique prognostic IRGs (CRABP2, LTB4R, PTGER1 and TEK) were finally affirmed to associate with tumor survival independently and utilized to establish the risk score model. All candidates' expression was successfully laboratory confirmed by q-PCR. CIBERSORT analysis implied patients in unfavorable-risk group with high CD8 T cell, regulatory T cell and NK cell infiltration, as well as high expression of PD- 1, CTLA4, TNFRSF9, TIGIT and LAG3. A nomogram combined IRGs risk score with age, gender, TNM stage, Fuhrman grade, necrosis was further generated to predict of 3- and 5-year OS, which exhibited superior discriminative power (AUCs were 0.811 and 0.795). Conclusions: Our study established and validated a survival prognostic model system based on 4 unique immune related genes in ccRCC, which expands knowledge in tumor immune status and provide a potent prediction tool in future.