The patients clinical data for prostate adenocarcinoma were downloaded using TCGA assembler.73 However, downloaded data comprised of only tumor pathologic and node pathologic information. Thus, based on tumor pathologic and node pathologic data as per 'https://cancerstaging.org/references-tools/quickreferences/Documents/ProstateSmall.pdf', samples were categorized into primary and metastatic tumor. Afterwards, level3 TCGA RNA-seq data (including raw_read_count and scaled_estimate for each sample) for all primary tumor, metastatic tumor and matched normal samples were downloaded using TCGA assembler. Transcript per million values for each gene was obtained by multiplying scaled_estimate by 1 000 000. Boxplot was generated using R (https://cran.r-project.org/).