The probability of tumor progression in intermediate/high-risk clear cell renal cell carcinoma (ccRCC) is highly variable, underlining the lack of predictive accuracy of the current clinicopathological factors. To develop an accurate prognostic classifier for these patients, we analyzed global gene expression patterns in 13 tissue samples from progressive and non-progressive ccRCC using Illumina Hi-seq 4000. Expression levels of 22 selected differentially expressed genes (DEG) were assessed by nCounter analysis in an independent series of 71 ccRCCs. A clinicopathological-molecular model for predicting tumor progression was developed and in silico validated in a total of 202 ccRCC patients using the TCGA cohort. A total of 1202 DEGs were found between progressive and non-progressive intermediate/high-risk ccRCC in RNAseq analysis, and seven of the 22 DEGs selected were validated by nCounter. Expression of HS6ST2, pT stage, tumor size, and ISUP grade were found to be independent prognostic factors for tumor progression. A risk score generated using these variables was able to distinguish patients at higher risk of tumor progression (HR 7.27; p < 0.001), consistent with the results obtained from the TCGA cohort (HR 2.74; p < 0.002). In summary, a combined prognostic algorithm was successfully developed and validated. This model may aid physicians to select high-risk patients for adjuvant therapy.
Keywords: RNA sequencing; biomarkers; clear cell renal cell carcinoma; disease progression; gene expression; prognostic factors.