Global gene expression data was generated using RNA isolated from SUB1 siRNA knockdown PC3 and non-target control cells in profiling analysis as well as in transcriptome sequencing analysis.6 Expression profiling was performed using the Agilent Whole Human Genome Oligo Microarray (Agilent, Santa Clara, CA, USA) and the analysis was performed according to the manufacturer's protocol. A bioconductor package ‘agilp'74 was used to extract and normalize raw data from two-channel experiment arrays. Loess normalization was applied on each array. The gene expression profiling data has been deposited at gene expression omnibus (GSE74895). The differential expression of each gene was estimated by subtracting loess normalized log-2-transformed signal intensity of control sample from that of knockdown sample. Genes with absolute fold change of 2 were selected as differentially expressed genes. The heatmap.2 function of R package ‘gplots' was used to create the heat map. To measure mRNA expression levels, total RNA was isolated from prostate cancer cell lines and prostate cancer tissue samples using the RNeasy Mini Kit (Qiagen, Valencia, CA, USA). qPCR was performed as described.37 All primers were synthesized by Integrated DNA Technologies (Coralville, IA, USA) and PCR reactions were performed in triplicate. Primer sequences are listed in Supplementary Table S2.