[Abstract] Objective: To investigate the pathogenesis of prostate cancer by analyzing the associated hub gene modules of prostate cancer and identifying key transcription factors and genes that affect these modules. Methods: WGCNA (weighted gene co-expressed network analysis) was used to identify hub gene modules associated with important clinicopathological features of prostate cancer, such as pathological staging, Gleason grading etc. The OPOSSUM online tool was used to analyze the transcription factors enriching and regulating those genes. Pathway enrichment analysis and protein-protein interaction network analysis were used to identify key genes in prostate cancer. Finally, the effects of these genes on clinical features and disease-free survival (DFS) of prostate cancer patients were analyzed. Results: Three hub modules were identified, and they were highly associated with pathologic T stage, pathologic N stage and Gleason grading of prostate cancer, respectively. Further screening revealed 13 key dysregulated transcription factors that participated in the regulation of these three hub modules. The differentially expressed genes regulated by the 13 key transcription factors were significantly enriched in Calcium signaling pathway, cGMP-PKG signaling pathway and cAMP signaling pathway. 14 key genes (PRKG1, PRKG2, CYSLTR2, GRPR, CHRM3, ADCY5, ADRA1D, EDNRA, EDNRB, CYSLTR2, AGTR1, GRPR, GRIA1 and OXT) were at important nodes in the gene network. Among them, the high expression of ADRA1A, PRKG2, CHRM3, ADRA1D and EDN3 significantly extended the DFS of patients with prostate cancer (all P<0.01). Conclusion: ADRA1A, PRKG2, CHRM3, ADRA1D and EDN3 are regulated by key dysregulated transcription factors and highly associated with clinical features of prostate cancer. Their high expressions will significantly prolong the DFS of prostate cancer patients, which may shed light to the discovery of mechanism in prostate adenocarcinoma.