Comparison of MiRNA and Gene Expression Profiles Between Metastatic and Primary Prostate Cancer

The present study aimed to identify the regulatory mechanisms associated with the metastasis of prostate cancer (PC). The microRNA (miRNA/miR) microarray dataset GSE21036 and gene transcript dataset GSE21034 were downloaded from the Gene Expression Omnibus database. Following pre-processing, differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs) between samples from patients with primary prostate cancer (PPC) and metastatic prostate cancer (MPC) with |log fold change (FC)| >1 and a false discovery rate <0.05 were selected using the Linear Models for Microarray and RNA-seq Data 4 package of R. Next, a DEM-DEG regulatory network was constructed by downloading miRNA-DEG pairs from the miRNA.org database. Finally, functional annotation of each DEM-DEG module was performed using the Database for Annotation, Visualization and Integrated Discovery based on the Gene Ontology database. The upregulated miRNAs, including miR-144, miR-494 and miR-181a, exhibited a higher degree of connections compared with other nodes, including in the DEM-DEG regulatory network, and regulated a number of downregulated DEGs. According to the functional annotation of the DEM-DEG modules, miR-144 and its targeted DEGs enriched the highest number of biological process terms (36 terms), followed by miR-494 (24 terms), miR-30d (18 terms), miR-181a (15 terms), hsa-miR-196a (8 terms), miR-708 (7 terms) and miR-486-5p (2 terms). Therefore, these miRNAs may serve roles in the metastasis of PC cells via downregulation of their corresponding target DEGs.