Prostate cancer, a complex disease, can be relatively harmless or extremely aggressive. To identify candidate genes involved in causal pathways of aggressive prostate cancer, we implemented a systems biology approach by combining differential expression analysis and coexpression network analysis to evaluate transcriptional profiles using lymphoblastoid cell lines from 62 prostate cancer patients with aggressive phenotype (Gleason grade >or= 8) and 63 prostate cancer patients with nonaggressive phenotype (Gleason grade <or= 5). From 13,935 mRNA genes and 273 microRNAs (miRNA) tested, we identified significant differences in 1,100 mRNAs and 7 miRNAs with a false discovery rate (FDR) of <0.01. We also identified a coexpression module demonstrating significant association with the aggressive phenotype of prostate cancer (P = 3.67 x 10(-11)). The module of interest was characterized by overrepresentation of cell cycle-related genes (FDR = 3.50 x 10(-50)). From this module, we further defined 20 hub genes that were highly connected to other genes. Interestingly, 5 of the 7 differentially expressed miRNAs have been implicated in cell cycle regulation and 2 (miR-145 and miR-331-3p) are predicted to target 3 of the 20 hub genes. Ectopic expression of these two miRNAs reduced expression of target hub genes and subsequently resulted in cell growth inhibition and apoptosis. These results suggest that cell cycle is likely to be a molecular pathway causing aggressive phenotype of prostate cancer. Further characterization of cell cycle-related genes (particularly, the hub genes) and miRNAs that regulate these hub genes could facilitate identification of candidate genes responsible for the aggressive phenotype and lead to a better understanding of prostate cancer etiology and progression.