Abstract
Clear cell renal cell carcinoma (ccRCC) is a major histological subtype of renal cell carcinoma and can be clinically divided into four stages according to the TNM criteria. Identifying clinical stage-related genes is beneficial for improving the early diagnosis and prognosis of ccRCC. By using bioinformatics analysis, we aim to identify clinical stage-relevant genes that are significantly associated with the development of ccRCC. First, we analyzed the gene expression microarray data sets: GSE53757 and GSE73731. We divided these data into five groups by staging information-normal tissue and ccRCC stages I, II, III, and IV-and eventually identified 500 differentially expressed genes (DEGs). To obtain precise stage-relevant genes, we subsequently applied weighted gene coexpression network analysis (WGCNA) to the GSE73731 dataset and KIRC data from The Cancer Genome Atlas (TCGA). Two modules from each dataset were identified to be related to the tumor TNM stage. Several genes with high inner connection inside the modules were considered hub genes. The intersection results between hub genes of key modules and 500 DEGs revealed UBE2C, BUB1B, RRM2, and TPX2 as highly associated with the stage of ccRCC. In addition, the candidate genes were validated at both the RNA expression level and the protein level. Survival analysis also showed that 4 genes were significantly correlated with overall survival. In conclusion, our study affords a deeper understanding of the molecular mechanisms associated with the development of ccRCC and provides potential biomarkers for early diagnosis and individualized treatment for patients at different stages of ccRCC.