Identification of hub genes in prostate cancer using robust rank aggregation and weighted gene co-expression network analysis

The pathogenic mechanisms of prostate cancer (PCa) remain to be defined. In this study, we utilized the Robust Rank Aggregation (RRA) method to integrate 10 eligible PCa microarray datasets from the GEO and identified a set of significant differentially expressed genes (DEGs) between tumor samples and normal, matched specimens. To explore potential associations between gene sets and PCa clinical features and to identify hub genes, we utilized WGCNA to construct gene co-expression networks incorporating the DEGs screened with the use of RRA. From the key module, we selected LMNB1, TK1, ZWINT, and RACGAP1 for validation. We found that these genes were up-regulated in PCa samples, and higher expression levels were associated with higher Gleason scores and tumor grades. Moreover, ROC and K-M plots indicated these genes had good diagnostic and prognostic value for PCa. On the other hand, methylation analyses suggested that the abnormal up-regulation of these four genes likely resulted from hypomethylation, while GSEA and GSVA for single hub gene revealed they all had a close association with proliferation of PCa cells. These findings provide new insight into PCa pathogenesis, and identify LMNB1, TK1, RACGAP1 and ZWINT as candidate biomarkers for diagnosis and prognosis of PCa.