The Genomic and Epigenomic Landscape of Double-Negative Metastatic Prostate Cancer.

Systemic targeted therapy in prostate cancer is primarily focused on ablating androgen signaling. Androgen deprivation therapy and second-generation androgen receptor (AR)-targeted therapy selectively favor the development of treatment-resistant subtypes of metastatic castration-resistant prostate cancer (mCRPC), defined by AR and neuroendocrine (NE) markers. Molecular drivers of double-negative (AR-/NE-) mCRPC are poorly defined. In this study, we comprehensively characterized treatment-emergent mCRPC by integrating matched RNA sequencing, whole-genome sequencing, and whole-genome bisulfite sequencing from 210 tumors. AR-/NE- tumors were clinically and molecularly distinct from other mCRPC subtypes, with the shortest survival, amplification of the chromatin remodeler CHD7, and PTEN loss. Methylation changes in CHD7 candidate enhancers were linked to elevated CHD7 expression in AR-/NE+ tumors. Genome-wide methylation analysis nominated Krüppel-like factor 5 (KLF5) as a driver of the AR-/NE- phenotype, and KLF5 activity was linked to RB1 loss. These observations reveal the aggressiveness of AR-/NE- mCRPC and could facilitate the identification of therapeutic targets in this highly aggressive disease. SIGNIFICANCE: Comprehensive characterization of the five subtypes of metastatic castration-resistant prostate cancer identified transcription factors that drive each subtype and showed that the double-negative subtype has the worst prognosis.

The 5-Hydroxymethylcytosine Landscape of Prostate Cancer.

Analysis of DNA methylation is a valuable tool to understand disease progression and is increasingly being used to create diagnostic and prognostic clinical biomarkers. While conversion of cytosine to 5-methylcytosine (5mC) commonly results in transcriptional repression, further conversion to 5-hydroxymethylcytosine (5hmC) is associated with transcriptional activation. Here we perform the first study integrating whole-genome 5hmC with DNA, 5mC, and transcriptome sequencing in clinical samples of benign, localized, and advanced prostate cancer. 5hmC is shown to mark activation of cancer drivers and downstream targets. Furthermore, 5hmC sequencing revealed profoundly altered cell states throughout the disease course, characterized by increased proliferation, oncogenic signaling, dedifferentiation, and lineage plasticity to neuroendocrine and gastrointestinal lineages. Finally, 5hmC sequencing of cell-free DNA from patients with metastatic disease proved useful as a prognostic biomarker able to identify an aggressive subtype of prostate cancer using the genes TOP2A and EZH2, previously only detectable by transcriptomic analysis of solid tumor biopsies. Overall, these findings reveal that 5hmC marks epigenomic activation in prostate cancer and identify hallmarks of prostate cancer progression with potential as biomarkers of aggressive disease. SIGNIFICANCE: In prostate cancer, 5-hydroxymethylcytosine delineates oncogene activation and stage-specific cell states and can be analyzed in liquid biopsies to detect cancer phenotypes. See related article by Wu and Attard, p. 3880.

Predicting cancer drug TARGETS - TreAtment Response Generalized Elastic-neT Signatures.

We are now in an era of molecular medicine, where specific DNA alterations can be used to identify patients who will respond to specific drugs. However, there are only a handful of clinically used predictive biomarkers in oncology. Herein, we describe an approach utilizing in vitro DNA and RNA sequencing and drug response data to create TreAtment Response Generalized Elastic-neT Signatures (TARGETS). We trained TARGETS drug response models using Elastic-Net regression in the publicly available Genomics of Drug Sensitivity in Cancer (GDSC) database. Models were then validated on additional in-vitro data from the Cancer Cell Line Encyclopedia (CCLE), and on clinical samples from The Cancer Genome Atlas (TCGA) and Stand Up to Cancer/Prostate Cancer Foundation West Coast Prostate Cancer Dream Team (WCDT). First, we demonstrated that all TARGETS models successfully predicted treatment response in the separate in-vitro CCLE treatment response dataset. Next, we evaluated all FDA-approved biomarker-based cancer drug indications in TCGA and demonstrated that TARGETS predictions were concordant with established clinical indications. Finally, we performed independent clinical validation in the WCDT and found that the TARGETS AR signaling inhibitors (ARSI) signature successfully predicted clinical treatment response in metastatic castration-resistant prostate cancer with a statistically significant interaction between the TARGETS score and PSA response (p = 0.0252). TARGETS represents a pan-cancer, platform-independent approach to predict response to oncologic therapies and could be used as a tool to better select patients for existing therapies as well as identify new indications for testing in prospective clinical trials.

An integrated functional and clinical genomics approach reveals genes driving aggressive metastatic prostate cancer.

Genomic sequencing of thousands of tumors has revealed many genes associated with specific types of cancer. Similarly, large scale CRISPR functional genomics efforts have mapped genes required for cancer cell proliferation or survival in hundreds of cell lines. Despite this, for specific disease subtypes, such as metastatic prostate cancer, there are likely a number of undiscovered tumor specific driver genes that may represent potential drug targets. To identify such genetic dependencies, we performed genome-scale CRISPRi screens in metastatic prostate cancer models. We then created a pipeline in which we integrated pan-cancer functional genomics data with our metastatic prostate cancer functional and clinical genomics data to identify genes that can drive aggressive prostate cancer phenotypes. Our integrative analysis of these data reveals known prostate cancer specific driver genes, such as AR and HOXB13, as well as a number of top hits that are poorly characterized. In this study we highlight the strength of an integrated clinical and functional genomics pipeline and focus on two top hit genes, KIF4A and WDR62. We demonstrate that both KIF4A and WDR62 drive aggressive prostate cancer phenotypes in vitro and in vivo in multiple models, irrespective of AR-status, and are also associated with poor patient outcome.

Multiomics Characterization of Low-Grade Serous Ovarian Carcinoma Identifies Potential Biomarkers of MEK Inhibitor Sensitivity and Therapeutic Vulnerability.

Low-grade serous ovarian carcinoma (LGSOC) is a rare tumor subtype with high case fatality rates in patients with metastatic disease. There is a pressing need to develop effective treatments using newly available preclinical models for therapeutic discovery and drug evaluation. Here, we use multiomics integration of whole-exome sequencing, RNA sequencing, and mass spectrometry-based proteomics on 14 LGSOC cell lines to elucidate novel biomarkers and therapeutic vulnerabilities. Comparison of LGSOC cell line data with LGSOC tumor data enabled predictive biomarker identification of MEK inhibitor (MEKi) efficacy, with KRAS mutations found exclusively in MEKi-sensitive cell lines and NRAS mutations found mostly in MEKi-resistant cell lines. Distinct patterns of Catalogue of Somatic Mutations in Cancer mutational signatures were identified in MEKi-sensitive and MEKi-resistant cell lines. Deletions of CDKN2A/B and MTAP genes were more frequent in cell lines than tumor samples and possibly represent key driver events in the absence of KRAS/NRAS/BRAF mutations. These LGSOC cell lines were representative models of the molecular aberrations found in LGSOC tumors. For prediction of in vitro MEKi efficacy, proteomic data provided better discrimination than gene expression data. Condensin, minichromosome maintenance, and replication factor C protein complexes were identified as potential treatment targets in MEKi-resistant cell lines. This study suggests that CDKN2A/B or MTAP deficiency may be exploited using synthetically lethal treatment strategies, highlighting the importance of using proteomic data as a tool for molecular drug prediction. Multiomics approaches are crucial to improving our understanding of the molecular underpinnings of LGSOC and applying this information to develop new therapies. SIGNIFICANCE: These findings highlight the utility of global multiomics to characterize LGSOC cell lines as research models, to determine biomarkers of MEKi resistance, and to identify potential novel therapeutic targets.

Autoantibody Landscape in Patients with Advanced Prostate Cancer.

PURPOSE: Autoantibody responses in cancer are of great interest, as they may be concordant with T-cell responses to cancer antigens or predictive of response to cancer immunotherapies. Thus, we sought to characterize the antibody landscape of metastatic castration-resistant prostate cancer (mCRPC). EXPERIMENTAL DESIGN: Serum antibody epitope repertoire analysis (SERA) was performed on patient serum to identify tumor-specific neoepitopes. Somatic mutation-specific neoepitopes were investigated by associating serum epitope enrichment scores with whole-genome sequencing results from paired solid tumor metastasis biopsies and germline blood samples. A protein-based immunome-wide association study (PIWAS) was performed to identify significantly enriched epitopes, and candidate serum antibodies enriched in select patients were validated by ELISA profiling. A distinct cohort of patients with melanoma was evaluated to validate the top cancer-specific epitopes. RESULTS: SERA was performed on 1,229 serum samples obtained from 72 men with mCRPC and 1,157 healthy control patients. Twenty-nine of 6,636 somatic mutations (0.44%) were associated with an antibody response specific to the mutated peptide. PIWAS analyses identified motifs in 11 proteins, including NY-ESO-1 and HERVK-113, as immunogenic in mCRPC, and ELISA confirmed serum antibody enrichment in candidate patients. Confirmatory PIWAS, Identifying Motifs Using Next-generation sequencing Experiments (IMUNE), and ELISA analyses performed on serum samples from 106 patients with melanoma similarly revealed enriched cancer-specific antibody responses to NY-ESO-1. CONCLUSIONS: We present the first large-scale profiling of autoantibodies in advanced prostate cancer, utilizing a new antibody profiling approach to reveal novel cancer-specific antigens and epitopes. Our study recovers antigens of known importance and identifies novel tumor-specific epitopes of translational interest.

Identification of conserved evolutionary trajectories in tumors.

MOTIVATION: As multi-region, time-series and single-cell sequencing data become more widely available; it is becoming clear that certain tumors share evolutionary characteristics with others. In the last few years, several computational methods have been developed with the goal of inferring the subclonal composition and evolutionary history of tumors from tumor biopsy sequencing data. However, the phylogenetic trees that they report differ significantly between tumors (even those with similar characteristics). RESULTS: In this article, we present a novel combinatorial optimization method, CONETT, for detection of recurrent tumor evolution trajectories. Our method constructs a consensus tree of conserved evolutionary trajectories based on the information about temporal order of alteration events in a set of tumors. We apply our method to previously published datasets of 100 clear-cell renal cell carcinoma and 99 non-small-cell lung cancer patients and identify both conserved trajectories that were reported in the original studies, as well as new trajectories. AVAILABILITY AND IMPLEMENTATION: CONETT is implemented in C++ and available at https://github.com/ehodzic/CONETT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Well-Differentiated Papillary Mesothelioma of the Peritoneum Is Genetically Distinct from Malignant Mesothelioma.

Well-differentiated papillary mesothelioma (WDPM) is an uncommon mesothelial proliferation that is most commonly encountered as an incidental finding in the peritoneal cavity. There is controversy in the literature about whether WDPM is a neoplasm or a reactive process and, if neoplastic, whether it is a variant or precursor of epithelial malignant mesothelioma or is a different entity. Using whole exome sequencing of five WDPMs of the peritoneum, we have identified distinct mutations in EHD1, ATM, FBXO10, SH2D2A, CDH5, MAGED1, and TP73 shared by WDPM cases but not reported in malignant mesotheliomas. Furthermore, we show that WDPM is strongly enriched with C > A transversion substitution mutations, a pattern that is also not found in malignant mesotheliomas. The WDPMs lacked the alterations involving BAP1, SETD2, NF2, CDKN2A/B, LASTS1/2, PBRM1, and SMARCC1 that are frequently found in malignant mesotheliomas. We conclude that WDPMs are neoplasms that are genetically distinct from malignant mesotheliomas and, based on observed mutations, do not appear to be precursors of malignant mesotheliomas.

Y-box binding protein-1 is crucial in acquired drug resistance development in metastatic clear-cell renal cell carcinoma.

BACKGROUND: Renal cell carcinoma (RCC) is a highly vascular tumor and patients with low risk metastatic RCC of clear-cell histological sub-type (mccRCC) are treated with tyrosine-kinase inhibitors (TKIs), sunitinib, as the first-line of treatment. Unfortunately, TKI resistance eventually develops, and the underlying molecular mechanism is not well understood. METHODS: RCC cell-line with metastatic clear-cell histology (Caki-1), and patient samples were analysed to identify the role of Y-box binding protein 1 (YB-1) and ATP-binding cassette sub-family B member 1 (ABCB-1) in acquired sunitinib-resistance development. Caki-1 was conditioned with increasing sunitinib doses to recapitulate acquired resistance development in clinics. Sunitinib-conditioned and wild-type Caki-1 were subjected to cell viability assay, scratch assay, chicken embryo chorioallantoic membrane engraftment and proteomics analysis. Classical biochemical assays like flow cytometry, immunofluorescent staining, immunohistochemical staining, optical coherence tomography imaging, Western Blot and RT-PCR assays were applied to determine the possible mechanism of sunitinib-resistance development and the effect of drug treatments. Publicly available data was also used to determine the role of YB-1 upregulation in ccRCC and the patients' overall survival. RESULTS: We demonstrate that YB-1 and ABCB-1 are upregulated in sunitinib-resistant in vitro, ex vivo, in vivo and patient samples compared to the sensitive samples. This provides evidence to a mechanism of acquired sunitinib-resistance development in mccRCC. Furthermore, our results establish that inhibiting ABCB-1 with elacridar, in addition to sunitinib, has a positive impact on reverting sunitinib-resistance development in in vitro, ex vivo and in vivo models. CONCLUSION: This work proposes a targeted therapy (elacridar and sunitinib) to re-sensitize sunitinib-resistant mccRCC and, possibly, slow disease progression.

Pathway and network analysis of more than 2500 whole cancer genomes.

The catalog of cancer driver mutations in protein-coding genes has greatly expanded in the past decade. However, non-coding cancer driver mutations are less well-characterized and only a handful of recurrent non-coding mutations, most notably TERT promoter mutations, have been reported. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancer across 38 tumor types, we perform multi-faceted pathway and network analyses of non-coding mutations across 2583 whole cancer genomes from 27 tumor types compiled by the ICGC/TCGA PCAWG project that was motivated by the success of pathway and network analyses in prioritizing rare mutations in protein-coding genes. While few non-coding genomic elements are recurrently mutated in this cohort, we identify 93 genes harboring non-coding mutations that cluster into several modules of interacting proteins. Among these are promoter mutations associated with reduced mRNA expression in TP53, TLE4, and TCF4. We find that biological processes had variable proportions of coding and non-coding mutations, with chromatin remodeling and proliferation pathways altered primarily by coding mutations, while developmental pathways, including Wnt and Notch, altered by both coding and non-coding mutations. RNA splicing is primarily altered by non-coding mutations in this cohort, and samples containing non-coding mutations in well-known RNA splicing factors exhibit similar gene expression signatures as samples with coding mutations in these genes. These analyses contribute a new repertoire of possible cancer genes and mechanisms that are altered by non-coding mutations and offer insights into additional cancer vulnerabilities that can be investigated for potential therapeutic treatments.

Identification of gene signature for treatment response to guide precision oncology in clear-cell renal cell carcinoma.

Clear-cell renal cell carcinoma (ccRCC) is a common therapy resistant disease with aberrant angiogenic and immunosuppressive features. Patients with metastatic disease are treated with targeted therapies based on clinical features: low-risk patients are usually treated with anti-angiogenic drugs and intermediate/high-risk patients with immune therapy. However, there are no biomarkers available to guide treatment choice for these patients. A recently published phase II clinical trial observed a correlation between ccRCC patients' clustering and their response to targeted therapy. However, the clustering of these groups was not distinct. Here, we analyzed the gene expression profile of 469 ccRCC patients, using featured selection technique, and have developed a refined 66-gene signature for improved sub-classification of patients. Moreover, we have identified a novel comprehensive expression profile to distinguish between migratory stromal and immune cells. Furthermore, the proposed 66-gene signature was validated using a different cohort of 64 ccRCC patients. These findings are foundational for the development of reliable biomarkers that may guide treatment decision-making and improve therapy response in ccRCC patients.

Combinatorial Detection of Conserved Alteration Patterns for Identifying Cancer Subnetworks.

BACKGROUND: Advances in large-scale tumor sequencing have led to an understanding that there are combinations of genomic and transcriptomic alterations specific to tumor types, shared across many patients. Unfortunately, computational identification of functionally meaningful and recurrent alteration patterns within gene/protein interaction networks has proven to be challenging. FINDINGS: We introduce a novel combinatorial method, cd-CAP (combinatorial detection of conserved alteration patterns), for simultaneous detection of connected subnetworks of an interaction network where genes exhibit conserved alteration patterns across tumor samples. Our method differentiates distinct alteration types associated with each gene (rather than relying on binary information of a gene being altered or not) and simultaneously detects multiple alteration profile conserved subnetworks. CONCLUSIONS: In a number of The Cancer Genome Atlas datasets, cd-CAP identified large biologically significant subnetworks with conserved alteration patterns, shared across many tumor samples.

BAP1 haploinsufficiency predicts a distinct immunogenic class of malignant peritoneal mesothelioma.

BACKGROUND: Malignant peritoneal mesothelioma (PeM) is a rare and fatal cancer that originates from the peritoneal lining of the abdomen. Standard treatment of PeM is limited to cytoreductive surgery and/or chemotherapy, and no effective targeted therapies for PeM exist. Some immune checkpoint inhibitor studies of mesothelioma have found positivity to be associated with a worse prognosis. METHODS: To search for novel therapeutic targets for PeM, we performed a comprehensive integrative multi-omics analysis of the genome, transcriptome, and proteome of 19 treatment-naïve PeM, and in particular, we examined BAP1 mutation and copy number status and its relationship to immune checkpoint inhibitor activation. RESULTS: We found that PeM could be divided into tumors with an inflammatory tumor microenvironment and those without and that this distinction correlated with haploinsufficiency of BAP1. To further investigate the role of BAP1, we used our recently developed cancer driver gene prioritization algorithm, HIT'nDRIVE, and observed that PeM with BAP1 haploinsufficiency form a distinct molecular subtype characterized by distinct gene expression patterns of chromatin remodeling, DNA repair pathways, and immune checkpoint receptor activation. We demonstrate that this subtype is correlated with an inflammatory tumor microenvironment and thus is a candidate for immune checkpoint blockade therapies. CONCLUSIONS: Our findings reveal BAP1 to be a potential, easily trackable prognostic and predictive biomarker for PeM immunotherapy that refines PeM disease classification. BAP1 stratification may improve drug response rates in ongoing phases I and II clinical trials exploring the use of immune checkpoint blockade therapies in PeM in which BAP1 status is not considered. This integrated molecular characterization provides a comprehensive foundation for improved management of a subset of PeM patients.

High-throughput detection of RNA processing in bacteria.

BACKGROUND: Understanding the RNA processing of an organism's transcriptome is an essential but challenging step in understanding its biology. Here we investigate with unprecedented detail the transcriptome of Pseudomonas aeruginosa PAO1, a medically important and innately multi-drug resistant bacterium. We systematically mapped RNA cleavage and dephosphorylation sites that result in 5'-monophosphate terminated RNA (pRNA) using monophosphate RNA-Seq (pRNA-Seq). Transcriptional start sites (TSS) were also mapped using differential RNA-Seq (dRNA-Seq) and both datasets were compared to conventional RNA-Seq performed in a variety of growth conditions. RESULTS: The pRNA-Seq library revealed known tRNA, rRNA and transfer-messenger RNA (tmRNA) processing sites, together with previously uncharacterized RNA cleavage events that were found disproportionately near the 5' ends of transcripts associated with basic bacterial functions such as oxidative phosphorylation and purine metabolism. The majority (97%) of the processed mRNAs were cleaved at precise codon positions within defined sequence motifs indicative of distinct endonucleolytic activities. The most abundant of these motifs corresponded closely to an E. coli RNase E site previously established in vitro. Using the dRNA-Seq library, we performed an operon analysis and predicted 3159 potential TSS. A correlation analysis uncovered 105 antiparallel pairs of TSS that were separated by 18 bp from each other and were centered on single palindromic TAT(A/T)ATA motifs (likely - 10 promoter elements), suggesting that, consistent with previous in vitro experimentation, these sites can initiate transcription bi-directionally and may thus provide a novel form of transcriptional regulation. TSS and RNA-Seq analysis allowed us to confirm expression of small non-coding RNAs (ncRNAs), many of which are differentially expressed in swarming and biofilm formation conditions. CONCLUSIONS: This study uses pRNA-Seq, a method that provides a genome-wide survey of RNA processing, to study the bacterium Pseudomonas aeruginosa and discover extensive transcript processing not previously appreciated. We have also gained novel insight into RNA maturation and turnover as well as a potential novel form of transcription regulation. NOTE: All sequence data has been submitted to the NCBI sequence read archive. Accession numbers are as follows: [NCBI sequence read archive: SRX156386, SRX157659, SRX157660, SRX157661, SRX157683 and SRX158075]. The sequence data is viewable using Jbrowse on www.pseudomonas.com .

HIT'nDRIVE: patient-specific multidriver gene prioritization for precision oncology.

Prioritizing molecular alterations that act as drivers of cancer remains a crucial bottleneck in therapeutic development. Here we introduce HIT'nDRIVE, a computational method that integrates genomic and transcriptomic data to identify a set of patient-specific, sequence-altered genes, with sufficient collective influence over dysregulated transcripts. HIT'nDRIVE aims to solve the "random walk facility location" (RWFL) problem in a gene (or protein) interaction network, which differs from the standard facility location problem by its use of an alternative distance measure: "multihitting time," the expected length of the shortest random walk from any one of the set of sequence-altered genes to an expression-altered target gene. When applied to 2200 tumors from four major cancer types, HIT'nDRIVE revealed many potentially clinically actionable driver genes. We also demonstrated that it is possible to perform accurate phenotype prediction for tumor samples by only using HIT'nDRIVE-seeded driver gene modules from gene interaction networks. In addition, we identified a number of breast cancer subtype-specific driver modules that are associated with patients' survival outcome. Furthermore, HIT'nDRIVE, when applied to a large panel of pan-cancer cell lines, accurately predicted drug efficacy using the driver genes and their seeded gene modules. Overall, HIT'nDRIVE may help clinicians contextualize massive multiomics data in therapeutic decision making, enabling widespread implementation of precision oncology.

BIRC6 Targeting as Potential Therapy for Advanced, Enzalutamide-Resistant Prostate Cancer.

Purpose: Enzalutamide resistance has emerged as a major problem in the management of castration-resistant prostate cancer (CRPC). Research on therapy resistance of CRPCs has primarily focused on the androgen receptor pathway. In contrast, there is limited information on antiapoptotic mechanisms that may facilitate the treatment resistance. The inhibitor of apoptosis proteins (IAP) family is well recognized for its role in promoting treatment resistance of cancers by inhibiting drug-induced apoptosis. Here, we examined whether BIRC6, an IAP family member, has a role in enzalutamide resistance of CRPCs and could provide a therapeutic target for enzalutamide-resistant CRPC.Experimental Design: Use of enzalutamide-resistant CRPC models: (i) the transplantable, first high-fidelity LTL-313BR patient-derived enzalutamide-resistant CRPC tissue xenograft line showing primary enzalutamide resistance, (ii) MR42D and MR49F CRPC cells/xenografts showing acquired enzalutamide resistance. Specific BIRC6 downregulation in these models was produced using a BIRC6-targeting antisense oligonucleotide (ASO-6w2). Gene expression was determined by qRT-PCR and gene expression profiling. Molecular pathways associated with growth inhibition were assessed via gene enrichment analysis.Results: Of eight IAPs examined, BIRC6 was the only one showing elevated expression in both enzalutamide-resistant CRPC models. Treatment with ASO-6w2 markedly suppressed growth of LTL-313BR xenografts and increased tumor apoptosis without inducing major host toxicity. Pathway enrichment analysis indicated that GPCR and matrisome signaling were the most significantly altered pathways. Furthermore, ASO-6w2 inhibited expression of prosurvival genes that were upregulated in the LTL-313BR line.Conclusions:BIRC6 targeting inhibited the growth of enzalutamide-resistant CRPC models and may represent a new option for clinical treatment of advanced, enzalutamide-resistant prostate cancer. Clin Cancer Res; 23(6); 1542-51. ©2016 AACR.

Heterogeneity in the inter-tumor transcriptome of high risk prostate cancer.

BACKGROUND: Genomic analyses of hundreds of prostate tumors have defined a diverse landscape of mutations and genome rearrangements, but the transcriptomic effect of this complexity is less well understood, particularly at the individual tumor level. We selected a cohort of 25 high-risk prostate tumors, representing the lethal phenotype, and applied deep RNA-sequencing and matched whole genome sequencing, followed by detailed molecular characterization. RESULTS: Ten tumors were exposed to neo-adjuvant hormone therapy and expressed marked evidence of therapy response in all except one extreme case, which demonstrated early resistance via apparent neuroendocrine transdifferentiation. We observe high inter-tumor heterogeneity, including unique sets of outlier transcripts in each tumor. Interestingly, outlier expression converged on druggable cellular pathways associated with cell cycle progression, translational control or immune regulation, suggesting distinct contemporary pathway affinity and a mechanism of tumor stratification. We characterize hundreds of novel fusion transcripts, including a high frequency of ETS fusions associated with complex genome rearrangements and the disruption of tumor suppressors. Remarkably, several tumors express unique but potentially-oncogenic non-ETS fusions, which may contribute to the phenotype of individual tumors, and have significance for disease progression. Finally, one ETS-negative tumor has a striking tandem duplication genotype which appears to be highly aggressive and present at low recurrence in ETS-negative prostate cancer, suggestive of a novel molecular subtype. CONCLUSIONS: The multitude of rare genomic and transcriptomic events detected in a high-risk tumor cohort offer novel opportunities for personalized oncology and their convergence on key pathways and functions has broad implications for precision medicine.

A meta-analysis approach for characterizing pan-cancer mechanisms of drug sensitivity in cell lines.

Understanding the heterogeneous drug response of cancer patients is essential to precision oncology. Pioneering genomic analyses of individual cancer subtypes have begun to identify key determinants of resistance, including up-regulation of multi-drug resistance (MDR) genes and mutational alterations of drug targets. However, these alterations are sufficient to explain only a minority of the population, and additional mechanisms of drug resistance or sensitivity are required to explain the remaining spectrum of patient responses to ultimately achieve the goal of precision oncology. We hypothesized that a pan-cancer analysis of in vitro drug sensitivities across numerous cancer lineages will improve the detection of statistical associations and yield more robust and, importantly, recurrent determinants of response. In this study, we developed a statistical framework based on the meta-analysis of expression profiles to identify pan-cancer markers and mechanisms of drug response. Using the Cancer Cell Line Encyclopaedia (CCLE), a large panel of several hundred cancer cell lines from numerous distinct lineages, we characterized both known and novel mechanisms of response to cytotoxic drugs including inhibitors of Topoisomerase 1 (TOP1; Topotecan, Irinotecan) and targeted therapies including inhibitors of histone deacetylases (HDAC; Panobinostat) and MAP/ERK kinases (MEK; PD-0325901, AZD6244). Notably, our analysis implicated reduced replication and transcriptional rates, as well as deficiency in DNA damage repair genes in resistance to TOP1 inhibitors. The constitutive activation of several signaling pathways including the interferon/STAT-1 pathway was implicated in resistance to the pan-HDAC inhibitor. Finally, a number of dysregulations upstream of MEK were identified as compensatory mechanisms of resistance to the MEK inhibitors. In comparison to alternative pan-cancer analysis strategies, our approach can better elucidate relevant drug response mechanisms. Moreover, the compendium of putative markers and mechanisms identified through our analysis can serve as a foundation for future studies into these drugs.