Systematic analysis of RNA-binding proteins identifies targetable therapeutic vulnerabilities in osteosarcoma.

Osteosarcoma is the most common primary malignant bone tumor with a strong tendency to metastasize, limiting the prognosis of affected patients. Genomic, epigenomic and transcriptomic analyses have demonstrated the exquisite molecular complexity of this tumor, but have not sufficiently defined the underlying mechanisms or identified promising therapeutic targets. To systematically explore RNA-protein interactions relevant to OS, we define the RNA interactomes together with the full proteome and the transcriptome of cells from five malignant bone tumors (four osteosarcomata and one malignant giant cell tumor of the bone) and from normal mesenchymal stem cells and osteoblasts. These analyses uncover both systematic changes of the RNA-binding activities of defined RNA-binding proteins common to all osteosarcomata and individual alterations that are observed in only a subset of tumors. Functional analyses reveal a particular vulnerability of these tumors to translation inhibition and a positive feedback loop involving the RBP IGF2BP3 and the transcription factor Myc which affects cellular translation and OS cell viability. Our results thus provide insight into potentially clinically relevant RNA-binding protein-dependent mechanisms of osteosarcoma.

Deciphering the role of FUS::DDIT3 expression and tumor microenvironment in myxoid liposarcoma development.

BACKGROUND: Myxoid liposarcoma (MLS) displays a distinctive tumor microenvironment and is characterized by the FUS::DDIT3 fusion oncogene, however, the precise functional contributions of these two elements remain enigmatic in tumor development. METHODS: To study the cell-free microenvironment in MLS, we developed an experimental model system based on decellularized patient-derived xenograft tumors. We characterized the cell-free scaffold using mass spectrometry. Subsequently, scaffolds were repopulated using sarcoma cells with or without FUS::DDIT3 expression that were analyzed with histology and RNA sequencing. RESULTS: Characterization of cell-free MLS scaffolds revealed intact structure and a large variation of protein types remaining after decellularization. We demonstrated an optimal culture time of 3 weeks and showed that FUS::DDIT3 expression decreased cell proliferation and scaffold invasiveness. The cell-free MLS microenvironment and FUS::DDIT3 expression both induced biological processes related to cell-to-cell and cell-to-extracellular matrix interactions, as well as chromatin remodeling, immune response, and metabolism. Data indicated that FUS::DDIT3 expression more than the microenvironment determined the pre-adipocytic phenotype that is typical for MLS. CONCLUSIONS: Our experimental approach opens new means to study the tumor microenvironment in detail and our findings suggest that FUS::DDIT3-expressing tumor cells can create their own extracellular niche.

Evolutionary mode and timing of dissemination of high-grade serous carcinomas.

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

The interplay between neoantigens and immune cells in sarcomas treated with checkpoint inhibition.

Introduction: Sarcomas are comprised of diverse bone and connective tissue tumors with few effective therapeutic options for locally advanced unresectable and/or metastatic disease. Recent advances in immunotherapy, in particular immune checkpoint inhibition (ICI), have shown promising outcomes in several cancer indications. Unfortunately, ICI therapy has provided only modest clinical responses and seems moderately effective in a subset of the diverse subtypes. Methods: To explore the immune parameters governing ICI therapy resistance or immune escape, we performed whole exome sequencing (WES) on tumors and their matched normal blood, in addition to RNA-seq from tumors of 31 sarcoma patients treated with pembrolizumab. We used advanced computational methods to investigate key immune properties, such as neoantigens and immune cell composition in the tumor microenvironment (TME). Results: A multifactorial analysis suggested that expression of high quality neoantigens in the context of specific immune cells in the TME are key prognostic markers of progression-free survival (PFS). The presence of several types of immune cells, including T cells, B cells and macrophages, in the TME were associated with improved PFS. Importantly, we also found the presence of both CD8+ T cells and neoantigens together was associated with improved survival compared to the presence of CD8+ T cells or neoantigens alone. Interestingly, this trend was not identified with the combined presence of CD8+ T cells and TMB; suggesting that a combined CD8+ T cell and neoantigen effect on PFS was important. Discussion: The outcome of this study may inform future trials that may lead to improved outcomes for sarcoma patients treated with ICI.

Chromosomal instability and a deregulated cell cycle are intrinsic features of high-risk gastrointestinal stromal tumours with a metastatic potential.

Patients with localised, high-risk gastrointestinal stromal tumours (GIST) benefit from adjuvant imatinib treatment. Still, approximately 40% of patients relapse within 3 years after adjuvant therapy and the clinical and histopathological features currently used for risk classification cannot precisely predict poor outcomes after standard treatment. This study aimed to identify genomic and transcriptomic profiles that could be associated with disease relapse and thus a more aggressive phenotype. Using a multi-omics approach, we analysed a cohort of primary tumours from patients with untreated, resectable high-risk GISTs. We compared patients who developed metastatic disease within 3 years after finishing adjuvant imatinib treatment and patients without disease relapse after more than 5 years of follow-up. Combining genomics and transcriptomics data, we identified somatic mutations and deregulated mRNA and miRNA genes intrinsic to each group. Our study shows that increased chromosomal instability (CIN), including chromothripsis and deregulated kinetochore and cell cycle signalling, separates high-risk samples according to metastatic potential. The increased CIN seems to be an intrinsic feature for tumours that metastasise and should be further validated as a novel prognostic biomarker for high-risk GIST.

Heritable defects in telomere and mitotic function selectively predispose to sarcomas.

Cancer genetics has to date focused on epithelial malignancies, identifying multiple histotype-specific pathways underlying cancer susceptibility. Sarcomas are rare malignancies predominantly derived from embryonic mesoderm. To identify pathways specific to mesenchymal cancers, we performed whole-genome germline sequencing on 1644 sporadic cases and 3205 matched healthy elderly controls. Using an extreme phenotype design, a combined rare-variant burden and ontologic analysis identified two sarcoma-specific pathways involved in mitotic and telomere functions. Variants in centrosome genes are linked to malignant peripheral nerve sheath and gastrointestinal stromal tumors, whereas heritable defects in the shelterin complex link susceptibility to sarcoma, melanoma, and thyroid cancers. These studies indicate a specific role for heritable defects in mitotic and telomere biology in risk of sarcomas.

Clonal evolution in primary breast cancers under sequential epirubicin and docetaxel monotherapy.

BACKGROUND: Subclonal evolution during primary breast cancer treatment is largely unexplored. We aimed to assess the dynamic changes in subclonal composition of treatment-naïve breast cancers during neoadjuvant chemotherapy. METHODS: We performed whole exome sequencing of tumor biopsies collected before, at therapy switch, and after treatment with sequential epirubicin and docetaxel monotherapy in 51 out of 109 patients with primary breast cancer, who were included in a prospectively registered, neoadjuvant single-arm phase II trial. RESULTS: There was a profound and differential redistribution of subclones during epirubicin and docetaxel treatment, regardless of therapy response. While truncal mutations and main subclones persisted, smaller subclones frequently appeared or disappeared. Reassessment of raw data, beyond formal mutation calling, indicated that the majority of subclones seemingly appearing during treatment were in fact present in pretreatment breast cancers, below conventional detection limits. Likewise, subclones which seemingly disappeared were still present, below detection limits, in most cases where tumor tissue remained. Tumor mutational burden (TMB) dropped during neoadjuvant therapy, and copy number analysis demonstrated specific genomic regions to be systematically lost or gained for each of the two chemotherapeutics. CONCLUSIONS: Sequential epirubicin and docetaxel monotherapy caused profound redistribution of smaller subclones in primary breast cancer, while early truncal mutations and major subclones generally persisted through treatment. TRIAL REGISTRATION: ClinicalTrials.gov, NCT00496795 , registered on July 4, 2007.

Clonal evolution after treatment pressure in multiple myeloma: heterogenous genomic aberrations and transcriptomic convergence.

We investigated genomic and transcriptomic changes in paired tumor samples of 29 in-house multiple myeloma (MM) patients and 28 patients from the MMRF CoMMpass study before and after treatment. A change in clonal composition was found in 46/57 (82%) of patients, and single-nucleotide variants (SNVs) increased from median 67 to 86. The highest increase in prevalence of genetic aberrations was found in RAS genes (60% to 72%), amp1q21 (18% to 35%), and TP53 (9% to 18%). The SBS-MM1 mutation signature was detected both in patients receiving high and low dose melphalan. A total of 2589 genes were differentially expressed between early and late samples (FDR < 0.05). Gene set enrichment analysis (GSEA) showed increased expression of E2F, MYC, and glycolysis pathways and a decreased expression in TNF-NFkB and TGFbeta pathways in late compared to early stage. Single sample GSEA (ssGSEA) scores of differentially expressed pathways revealed that these changes were most evident in end-stage disease. Increased expression of several potentially targetable genes was found at late disease stages, including cancer-testis antigens, XPO1 and ABC transporters. Our study demonstrates a transcriptomic convergence of pathways supporting increased proliferation and metabolism during disease progression in MM.

miR-486-5p expression is regulated by DNA methylation in osteosarcoma.

BACKGROUND: Osteosarcoma is the most common primary malignant tumour of bone occurring in children and young adolescents and is characterised by complex genetic and epigenetic changes. The miRNA miR-486-5p has been shown to be downregulated in osteosarcoma and in cancer in general. RESULTS: To investigate if the mir-486 locus is epigenetically regulated, we integrated DNA methylation and miR-486-5p expression data using cohorts of osteosarcoma cell lines and patient samples. A CpG island in the promoter of the ANK1 host gene of mir-486 was shown to be highly methylated in osteosarcoma cell lines as determined by methylation-specific PCR and direct bisulfite sequencing. High methylation levels were seen for osteosarcoma patient samples, xenografts and cell lines based on quantitative methylation-specific PCR. 5-Aza-2'-deoxycytidine treatment of osteosarcoma cell lines caused induction of miR-486-5p and ANK1, indicating common epigenetic regulation in osteosarcoma cell lines. When overexpressed, miR-486-5p affected cell morphology. CONCLUSIONS: miR-486-5p represents a highly cancer relevant, epigenetically regulated miRNA in osteosarcoma, and this knowledge contributes to the understanding of osteosarcoma biology.

The expressed mutational landscape of microsatellite stable colorectal cancers.

BACKGROUND: Colorectal cancer is the 2nd leading cause of cancer-related deaths with few patients benefiting from biomarker-guided therapy. Mutation expression is essential for accurate interpretation of mutations as biomarkers, but surprisingly, little has been done to analyze somatic cancer mutations on the expression level. We report a large-scale analysis of allele-specific mutation expression. METHODS: Whole-exome and total RNA sequencing was performed on 137 samples from 121 microsatellite stable colorectal cancers, including multiregional samples of primary and metastatic tumors from 4 patients. Data were integrated with allele-specific resolution. Results were validated in an independent set of 241 colon cancers. Therapeutic associations were explored by pharmacogenomic profiling of 15 cell lines or patient-derived organoids. RESULTS: The median proportion of expressed mutations per tumor was 34%. Cancer-critical mutations had the highest expression frequency (gene-wise mean of 58%), independent of frequent allelic imbalance. Systematic deviation from the general pattern of expression levels according to allelic frequencies was detected, including preferential expression of mutated alleles dependent on the mutation type and target gene. Translational relevance was suggested by correlations of KRAS/NRAS or TP53 mutation expression levels with downstream oncogenic signatures (p < 0.03), overall survival among patients with stage II and III cancer (KRAS/NRAS: hazard ratio 6.1, p = 0.0070), and targeted drug sensitivity. The latter was demonstrated for EGFR and MDM2 inhibition in pre-clinical models. CONCLUSIONS: Only a subset of mutations in microsatellite stable colorectal cancers were expressed, and the "expressed mutation dose" may provide an opportunity for more fine-tuned biomarker interpretations.

Cancer Predisposition Sequencing Reporter (CPSR): A flexible variant report engine for high-throughput germline screening in cancer.

The value of high-throughput germline genetic testing is increasingly recognized in clinical cancer care. Disease-associated germline variants in cancer patients are important for risk management and surveillance, surgical decisions and can also have major implications for treatment strategies since many are in DNA repair genes. With the increasing availability of high-throughput DNA sequencing in cancer clinics and research, there is thus a need to provide clinically oriented sequencing reports for germline variants and their potential therapeutic relevance on a per-patient basis. To meet this need, we have developed the Cancer Predisposition Sequencing Reporter (CPSR), an open-source computational workflow that generates a structured report of germline variants identified in known cancer predisposition genes, highlighting markers of therapeutic, prognostic and diagnostic relevance. A fully automated variant classification procedure based on more than 30 refined American College of Medical Genetics and Genomics (ACMG) criteria represents an integral part of the workflow. Importantly, the set of cancer predisposition genes profiled in the report can be flexibly chosen from more than 40 virtual gene panels established by scientific experts, enabling customization of the report for different screening purposes and clinical contexts. The report can be configured to also list actionable secondary variant findings, as recommended by ACMG. CPSR demonstrates comparable sensitivity and specificity for the detection of pathogenic variants when compared to other algorithms in the field. Technically, the tool is implemented in Python/R, and is freely available through Docker technology. Source code, documentation, example reports and installation instructions are accessible via the project GitHub page: https://github.com/sigven/cpsr.

Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden.

Sequencing of cell-free DNA in the blood of cancer patients (liquid biopsy) provides attractive opportunities for early diagnosis, assessment of treatment response, and minimally invasive disease monitoring. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we introduce an integrated genetic/epigenetic analysis method and demonstrate its utility on 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. Our method achieves sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. Moreover, we benchmark different metrics for cell-free DNA fragmentation analysis, and we introduce the LIQUORICE algorithm for detecting circulating tumor DNA based on cancer-specific chromatin signatures. Finally, we combine several fragmentation-based metrics into an integrated machine learning classifier for liquid biopsy analysis that exploits widespread epigenetic deregulation and is tailored to cancers with low mutation rates. Clinical associations highlight the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent genetic aberrations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

Methylation-dependent SUMOylation of the architectural transcription factor HMGA2.

High mobility group A2 (HMGA2) is a chromatin-associated protein involved in the regulation of stem cell function, embryogenesis and cancer development. Although the protein does not contain a consensus SUMOylation site, it is shown to be SUMOylated. In this study, we demonstrate that the first lysine residue in the reported K66KAE SUMOylation motif in HMGA2 can be methylated in vitro and in vivo by the Set7/9 methyltransferase. By editing the lysine, the increased hydrophobicity of the resulting 6-N-methyl-lysine transforms the sequence into a consensus SUMO motif. This post-translational editing dramatically increases the subsequent SUMOylation of this site. Furthermore, similar putative methylation-dependent SUMO motifs are found in a number of other chromatin factors, and we confirm methylation-dependent SUMOylation of a site in one such protein, the Polyhomeotic complex 1 homolog (PHC1). Together, these results suggest that crosstalk between methylation and SUMOylation is a general mode for regulation of chromatin function.

Clinical and molecular implications of NAB2-STAT6 fusion variants in solitary fibrous tumour.

Solitary fibrous tumour (SFT) is a mesenchymal neoplasm characterised by pathognomonic NAB2-STAT6 gene fusions. The clinical implications and prognostic value of different fusion variants has not been clarified. In the current study, we explore the clinicopathological, prognostic and molecular differences between tumours with different fusions. Thirty-nine patients with localised, extrameningeal SFT were included, of whom 20 developed distant recurrence and 19 were without recurrence after long term follow-up. Capture-based RNA sequencing identified 12 breakpoint variants, which were categorised into two groups based on the STAT6 domain composition in the predicted chimeric proteins. Twenty-one of 34 (62%) sequenced tumours had fusions with most of the STAT6 domains intact and were classified as STAT6-Full. Thirteen tumours (38%) contained only the transactivation domain of STAT6 and were classified as STAT6-TAD. Tumours with STAT6-TAD fusions had a higher mitotic count (p=0.016) and were associated with inferior recurrence-free interval (p=0.004) and overall survival (p=0.012). Estimated 10-year recurrence-free survival was 25% for patients with STAT6-TAD tumours compared to 78% for the STAT6-Full group. Distinct transcriptional signatures between the fusion groups were identified, including higher expression of FGF2 in the STAT6-TAD group and IGF2, EGR2, PDGFRB, STAT6 and several extracellular matrix genes in STAT6-Full tumours. In summary, we demonstrate that NAB2-STAT6 fusion variants are associated with distinct clinicopathological and molecular characteristics and have prognostic significance in extrameningeal SFT.

Discovery of novel candidates for anti-liposarcoma therapies by medium-scale high-throughput drug screening.

Sarcomas are a heterogeneous group of mesenchymal orphan cancers and new treatment alternatives beyond traditional chemotherapeutic regimes are much needed. So far, tumor mutation analysis has not led to significant treatment advances, and we have attempted to bypass this limitation by performing direct drug testing of a library of 353 anti-cancer compounds that are either FDA-approved, in clinical trial, or in advanced stages of preclinical development on a panel of 13 liposarcoma cell lines. We identified and validated six drugs, targeting different mechanisms and with good efficiency across the cell lines: MLN2238 -a proteasome inhibitor, GSK2126458 -a PI3K/mTOR inhibitor, JNJ-26481585 -a histone deacetylase inhibitor, triptolide-a multi-target drug, YM155 -a survivin inhibitor, and APO866 (FK866)-a nicotinamide phosphoribosyl transferase inhibitor. GR50s for those drugs were mostly in the nanomolar range, and in many cases below 10 nM. These drugs had long-lasting effect upon drug withdrawal, limited toxicity to normal cells and good efficacy also against tumor explants. Finally, we identified potential genomic biomarkers of their efficacy. Being approved or in clinical trials, these drugs are promising candidates for liposarcoma treatment.

Cell Fusion of Mesenchymal Stem/Stromal Cells and Breast Cancer Cells Leads to the Formation of Hybrid Cells Exhibiting Diverse and Individual (Stem Cell) Characteristics.

Cancer is one of the most common diseases worldwide, and treatment bears many challenges such as drug and radioresistance and formation of metastases. These difficulties are due to tumor heterogeneity, which has many origins. One may be cell fusion, a process that is relevant in both physiological (e.g., wound healing) and pathophysiological (cancer and viral infection) processes. In this study, we examined if cell fusion between mesenchymal stem/stromal cells (MSCs) and breast cancer (BC) cells occurs and if newly generated hybrid cells may exhibit cancer stem/initiating cell (CS/IC) characteristics. Therefore, several methods such as mammosphere assay, AldeRed assay, flow cytometry (CD24, CD44, CD104) and Western blot analysis (of epithelial to mesenchymal transition markers such as SNAIL, SLUG and Twist) were applied. In short, four different hybrid clones, verified by short tandem repeat (STR) analysis, were analyzed; each expressed an individual phenotype that seemed not to be explicitly related to either a more stem cell or cancer cell phenotype. These results show that cancer cells and MSCs are able to fuse spontaneously in vitro, thereby giving rise to hybrid cells with new properties, which likely indicate that cell fusion may be a trigger for tumor heterogeneity.

Accurate 3-gene-signature for early diagnosis of liposarcoma progression.

BACKGROUND: Well- and dedifferentiated liposarcoma (WD/DDLPS) are rare mesenchymal malignant tumors that account for 20% of all sarcomas in adults. The WD form is a low-grade malignancy with a favourable prognosis which may progress to DDLPS, a high-grade aggressive counterpart. WDLPS is referred to as atypical lipomatous tumour (ALT) when localised in extremities, due to its better prognosis. Currently the final differential diagnosis to distinguish between more aggressive and less aggressive form is based on post-surgical histological examination and no molecular biomarkers for early detection are available. METHODS: Quantitative polymerase chain reaction (qPCR) analysis of 11 metabolic genes involved in general and adipose tissue-specific metabolism, was performed on ALT (= 8), WDLPS (= 9) and DDLPS (= 20) samples. Subsequent statistical analysis was carried out to determine genes that most accurately can predict DDLPS differential diagnosis. Selected genes were further validated in a separate cohort by qPCR and the data statistically analysed. Deep sequencing was performed on DDLPS specimen from the metastatic patient and on five random WDLPS specimens. RESULTS: We established a three-gene signature based on PNPLA2, LIPE and PLIN1, which identified DDLPS with 100% sensitivity and 90% specificity, even in specimens from the WD component of DDLPS tumors. Interestingly, the PNPLA2 gene is deleted in 45% of DDLPS samples analyzed under TCGA project, and the deletion is associated with significantly lower PNPLA2 expression level. However, other mechanisms causing loss or downregulation of the expression of these three genes may be involved. Moreover, the significantly lower level of PNPLA2 is associated with R1 surgical margins, compare to R0 margins, which suggests the more invasive tumor phenotype in the absence of PNPLA2. CONCLUSIONS: The identified metabolic signature allows highly accurate differential diagnosis between WD- and DDLPS even in samples containing lipid droplets, a marker of differentiation, which makes it very suitable for the use on biopsies. In respect to the pathogenesis of the disease, our results give a new insight into possible molecular mechanisms involved and support the recent observation that deletion of PNPLA2 is a novel factor in liposarcoma progression.

Preclinical Evaluation of the Pan-FGFR Inhibitor LY2874455 in FRS2-Amplified Liposarcoma.

Background: FGFR inhibition has been proposed as treatment for dedifferentiated liposarcoma (DDLPS) with amplified FRS2, but we previously only demonstrated transient cytostatic effects when treating FRS2-amplified DDLPS cells with NVP-BGJ398. Methods: Effects of the more potent FGFR inhibitor LY2874455 were investigated in three DDLPS cell lines by measuring effects on cell growth and apoptosis in vitro and also testing efficacy in vivo. Genome, transcriptome and protein analyses were performed to characterize the signaling components in the FGFR pathway. Results: LY2874455 induced a stronger, longer-lasting growth inhibitory effect and moderate level of apoptosis for two cell lines. The third cell line, did not respond to FGFR inhibition, suggesting that FRS2 amplification alone is not sufficient to predict response. Importantly, efficacy of LY2874455 was confirmed in vivo, using an independent FRS2-amplified DDLPS xenograft model. Expression of FRS2 was similar in the responding and non-responding cell lines and we could not find any major difference in downstream FGFR signaling. The only FGF expressed by unstimulated non-responding cells was the intracellular ligand FGF11, whereas the responding cell lines expressed extracellular ligand FGF2. Conclusion: Our study supports LY2874455 as a better therapy than NVP-BGJ398 for FRS2-amplified liposarcoma, and a clinical trial is warranted.

Multifocal Primary Prostate Cancer Exhibits High Degree of Genomic Heterogeneity.

BACKGROUND: Most primary prostate cancers are multifocal with individual tumors harboring different aggressiveness; however, the genomic heterogeneity among these tumors is poorly understood. OBJECTIVE: To better understand the biological basis for clinical variability among different lesions, we sought to comprehensively characterize the heterogeneity of somatic gene mutations in multifocal prostate cancer. DESIGN, SETTING, AND PARTICIPANTS: High-coverage whole-exome sequencing of 153 frozen tissue samples, taken from two to three distinct tumor foci and one non-cancerous area from each of 41 patients, covering a total of 89 tumor foci. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: State-of-the-art bioinformatics tools for mutation calling and copy number determination from whole-exome sequencing data. RESULTS AND LIMITATIONS: We found a very high degree of interfocal heterogeneity among tumors, that is, 76% of pairwise-compared tumor foci from the same prostatectomy specimen had no point mutations in common and DNA copy number changes were rarely shared across cancer foci. The few point mutations shared across tumor foci were seldom in cancer-critical genes. CONCLUSIONS: In this first large genomic heterogeneity study of primary prostate cancer, we observe that different tumor foci within the same patient are genetically distinct, only rarely sharing any somatic gene mutations, including those in cancer driver genes. This heterogeneity affects how genomics-based management of prostate cancer can be implemented, as information from all tumor foci is necessary to draw valid conclusions about the cancer's genomic alterations. PATIENT SUMMARY: Most primary prostate cancers consist of multiple tumors within the same organ, but little is known about their relationships. We have compared the sets of gene mutations among such tumors and found that they only exceptionally have any in common. This will influence treatment decisions in the future as each tumor's mutations will render it unique and have to be considered to gain the best treatment results.