ABSTRACT Extension of prostate cancer beyond the primary site into the surrounding organs by local invasion or nodal metastasis is associated with poor prognosis. The emergence and evolution of cancer clones at this early stage of expansion and spread has not been studied in detail. We performed whole genome sequencing on 42 prostate cancer samples from the prostate, seminal vesicles and regional lymph nodes of five treatment-naive patients with locally advanced disease who underwent radical prostatectomy. Using cancer cell fractions computed from single nucleotide variants and copy number alterations, we reconstructed the tumour phylogenies, which in turn allowed us to infer key molecular steps in the progression of prostate cancer in these individuals. We mapped the clonal composition of cancer sampled across the prostate in each individual and inferred the routes of spread of cancer cells within the prostate and to seminal vesicles and lymph nodes. Based on these data, we delineated the route of tumour progression and metastasis following the transformation of adenocarcinoma to amphicrine morphology, the molecular events leading to whole genome duplication associated with a single clonal expansion and identified putative driver events associated with local invasion and lymph node metastasis. We also correlated genomic changes associated with differences in morphology and identified putative driver events associated with spread to seminal vesicle invasion and lymph node metastasis. Taken together, these findings have implications for diagnosis and risk stratification, in addition to providing a rationale for further studies to characterise the genetic changes associated with morphological transformation. Our results demonstrate the value of integrating multi-region sequencing with histopathological data to study tumour evolution and identify mechanisms of prostate cancer spread.

Purpose: It is important to identify molecular features that improve prostate cancer (PCa) risk stratification before radical treatment with curative intent. Molecular analysis of historical diagnostic formalin-fixed paraffin-embedded (FFPE) prostate biopsies from cohorts with post-radiotherapy (RT) long-term clinical follow-up has been limited. Utilizing parallel sequencing modalities, we performed a proof-of-principle sequencing analysis of historical diagnostic FFPE prostate biopsies. We compared patients with i) stable PCa post-primary or salvage RT (sPCa), ii) progressing PCa post-RT (pPCa), and iii) de novo metastatic PCa (mPCa). Experimental Design: A cohort of 19 patients with diagnostic prostate biopsies (n=6 sPCa, n=5 pPCa, n=8 mPCa) and mean 4 years 10 months follow-up (diagnosed 2009-2016) underwent nucleic acid extraction from demarcated malignancy. Samples underwent 3"RNA sequencing (3"RNAseq) (n=19), nanoString analysis (n=12) and Illumina 850k methylation (n=8) sequencing. Bioinformatic analysis was performed to coherently identify differentially expressed genes (DEGs) and methylated genomic regions (MGRs). Results: 18 of 19 samples provided useable 3"RNAseq data. Principal Component Analysis (PCA) demonstrated similar expression profiles between pPCa and mPCa cases, versus sPCa. Coherently differentially methylated probes between these groups identified ~600 differentially MGRs. The top 50 genes with increased expression in pPCa patients were associated with reduced progression-free survival post-RT (p<0.0001) in an external cohort. Conclusions: 3"RNAseq, nanoString and 850K-methylation analyses are each achievable from historical FFPE diagnostic pre-treatment prostate biopsies, unlocking the potential to utilize large cohorts of historic clinical samples. Profiling similarities between individuals with pPCa and mPCa suggests biological similarities and historical radiological staging limitations, which warrant further investigation.