Toward development of a precision medicine framework for metastatic, castration-resistant prostate cancer (mCRPC), we established a multi-institutional clinical sequencing infrastructure to conduct prospective whole-exome and transcriptome sequencing of bone or soft tissue tumor biopsies from a cohort of 150 mCRPC affected individuals. Aberrations of AR, ETS genes, TP53, and PTEN were frequent (40%-60% of cases), with TP53 and AR alterations enriched in mCRPC compared to primary prostate cancer. We identified new genomic alterations in PIK3CA/B, R-spondin, BRAF/RAF1, APC, β-catenin, and ZBTB16/PLZF. Moreover, aberrations of BRCA2, BRCA1, and ATM were observed at substantially higher frequencies (19.3% overall) compared to those in primary prostate cancers. 89% of affected individuals harbored a clinically actionable aberration, including 62.7% with aberrations in AR, 65% in other cancer-related genes, and 8% with actionable pathogenic germline alterations. This cohort study provides clinically actionable information that could impact treatment decisions for these affected individuals.

Genetic variants in antioxidant pathways might decrease the efficacy of radiation therapy (RT) by suppressing the generation of reactive oxygen species. We studied the association between single nucleotide polymorphisms (SNPs) in the antioxidant gene superoxide dismutase-2 (SOD2) and cancer-specific outcomes after RT.

MicroRNAs (miRNAs) play important roles in cancer formation and progression by suppressing the production of key functional proteins at the post-transcriptional level in a sequence-specific manner. While differential expression of miRNAs is widely observed in cancers including prostate cancer (PCa), how these miRNAs are transcriptionally regulated is largely unknown. MiRNA-221 and miRNA-222 (miR-221/-222) are well-established oncogenes and overexpressed in breast, liver, pancreas, and lung cancer, but their expression and biological functions in PCa remain controversial. Both up and down regulation have been observed in patient samples. Specifically, studies have demonstrated miR-221/-222 function as oncogenes, and promote PCa cell proliferation and the development of castration-resistant prostate cancer (CRPC). However, the expression level of miR-221/-222 is downregulated in several miRNA expression profiling studies. In this study, we demonstrate miR-221/-222 are androgen receptor (AR)-repressed genes and reside in a long primary transcript (pri-miRNA). Derepression of miR-221/-222 after androgen deprivation therapy (ADT) may enhance PCa cell proliferation potential through promoting G1/S phase transition. This function is likely transient but important in the development of CRPC. Downregulation of miR-221/-222 subsequently occurs once AR activity is restored through AR overexpression in CRPC. Our findings shed light on the complexity of transcriptional regulation of miRNAs in PCa and suggest context-dependent targeting of oncogenic miRNAs.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

Current prostate cancer prognostic models are based on pre-treatment prostate specific antigen (PSA) levels, biopsy Gleason score, and clinical staging but in practice are inadequate to accurately predict disease progression. Hence, we sought to develop a molecular panel for prostate cancer progression by reasoning that molecular profiles might further improve current clinical models.

The large registry, PROVENGE Registry for the Observation, Collection, and Evaluation of Experience Data (PROCEED)(NCT01306890), evaluated sipuleucel-T immunotherapy for asymptomatic/minimally symptomatic metastatic castration-resistant prostate cancer (mCRPC).

Although there are molecularly distinct subtypes of prostate cancer, no molecular classification system is used clinically. The ribonucleotide reductase small subunit M2 (RRM2) gene plays an oncogenic role in many cancers. Our previous study elucidated comprehensive molecular mechanisms of RRM2 in prostate cancer (PC). Given the potent functions of RRM2, we set out to determine whether the RRM2 signature can be used to identify aggressive subtypes of PC. We applied gene ontology and pathway analysis in RNA-seq datasets from PC cells overexpressing RRM2. We refined the RRM2 signature by integrating it with two molecular classification systems (PCS and PAM50 subtypes) that define aggressive PC subtypes (PCS1 and luminal B) and correlated signatures with clinical outcomes in six published cohorts comprising 4000 cases of PC. Increased expression of genes in the RRM2 signature was significantly correlated with recurrence, high Gleason score, and lethality of PC. Patients with high RRM2 levels showed higher PCS1 score, suggesting the aggressive PC feature. Consistently, RRM2-regulated genes were highly enriched in the PCS1 signature from multiple PC cohorts. A simplified RRM2 signature (12 genes) was identified by intersecting the RRM2 signature, PCS1 signature, and the PAM50 classifier. Intriguingly, inhibition of RRM2 specifically targets PCS1 and luminal B genes. Furthermore, 11 genes in the RRM2 signature were correlated with enzalutamide resistance by using a single-cell RNA-seq dataset from PC circulating tumor cells. Finally, high expression of RRM2 was associated with an immunosuppressive tumor-immune microenvironment in both primary prostate cancer and metastatic prostate cancer using CIBERSORT analysis and LM22, a validated leukocyte gene signature matrix. These data demonstrate that RRM2 is a driver of aggressive prostate cancer subtypes and contributes to immune escape, suggesting that RRM2 inhibition may be of clinical benefit for patients with PC.

Phosphatase and tensin homologue deleted on chromosome 10 (PTEN) is a well-characterized tumour-suppressor gene that is lost or mutated in about half of metastatic castration-resistant prostate cancers and in many other human cancers. The restoration of functional PTEN as a treatment for prostate cancer has, however, proven difficult. Here, we show that PTEN messenger RNA (mRNA) can be reintroduced into PTEN-null prostate cancer cells in vitro and in vivo via its encapsulation in polymer-lipid hybrid nanoparticles coated with a polyethylene glycol shell. The nanoparticles are stable in serum, elicit low toxicity and enable high PTEN mRNA transfection in prostate cancer cells. Moreover, significant inhibition of tumour growth is achieved when delivered systemically in multiple mouse models of prostate cancer. We also show that the restoration of PTEN function in PTEN-null prostate cancer cells inhibits the phosphatidylinositol 3-kinase (PI3K)-AKT pathway and enhances apoptosis. Our findings provide proof-of-principle evidence of the restoration of mRNA-based tumour suppression in vivo.

Heterogeneity in the genomic landscape of metastatic prostate cancer has become apparent through several comprehensive profiling efforts, but little is known about the impact of this heterogeneity on clinical outcome. Here, we report comprehensive genomic and transcriptomic analysis of 429 patients with metastatic castration-resistant prostate cancer (mCRPC) linked with longitudinal clinical outcomes, integrating findings from whole-exome, transcriptome, and histologic analysis. For 128 patients treated with a first-line next-generation androgen receptor signaling inhibitor (ARSI; abiraterone or enzalutamide), we examined the association of 18 recurrent DNA- and RNA-based genomic alterations, including androgen receptor (AR) variant expression, AR transcriptional output, and neuroendocrine expression signatures, with clinical outcomes. Of these, only RB1 alteration was significantly associated with poor survival, whereas alterations in RB1, AR, and TP53 were associated with shorter time on treatment with an ARSI. This large analysis integrating mCRPC genomics with histology and clinical outcomes identifies RB1 genomic alteration as a potent predictor of poor outcome, and is a community resource for further interrogation of clinical and molecular associations.

Most pancreatic neuroendocrine tumors (PNETs) do not produce excess hormones and are therefore considered 'non-functional'1-3. As clinical behaviors vary widely and distant metastases are eventually lethal2,4, biological classifications might guide treatment. Using enhancer maps to infer gene regulatory programs, we find that non-functional PNETs fall into two major subtypes, with epigenomes and transcriptomes that partially resemble islet α- and β-cells. Transcription factors ARX and PDX1 specify these normal cells, respectively5,6, and 84% of 142 non-functional PNETs expressed one or the other factor, occasionally both. Among 103 cases, distant relapses occurred almost exclusively in patients with ARX+PDX1- tumors and, within this subtype, in cases with alternative lengthening of telomeres. These markedly different outcomes belied similar clinical presentations and histology and, in one cohort, occurred irrespective of MEN1 mutation. This robust molecular stratification provides insight into cell lineage correlates of non-functional PNETs, accurately predicts disease course and can inform postoperative clinical decisions.

Prostate cancer is the second most common cause of male cancer deaths in the United States. However, the full range of prostate cancer genomic alterations is incompletely characterized. Here we present the complete sequence of seven primary human prostate cancers and their paired normal counterparts. Several tumours contained complex chains of balanced (that is, 'copy-neutral') rearrangements that occurred within or adjacent to known cancer genes. Rearrangement breakpoints were enriched near open chromatin, androgen receptor and ERG DNA binding sites in the setting of the ETS gene fusion TMPRSS2-ERG, but inversely correlated with these regions in tumours lacking ETS fusions. This observation suggests a link between chromatin or transcriptional regulation and the genesis of genomic aberrations. Three tumours contained rearrangements that disrupted CADM2, and four harboured events disrupting either PTEN (unbalanced events), a prostate tumour suppressor, or MAGI2 (balanced events), a PTEN interacting protein not previously implicated in prostate tumorigenesis. Thus, genomic rearrangements may arise from transcriptional or chromatin aberrancies and engage prostate tumorigenic mechanisms.

Comprehensive analyses of cancer genomes promise to inform prognoses and precise cancer treatments. A major barrier, however, is inaccessibility of metastatic tissue. A potential solution is to characterize circulating tumor cells (CTCs), but this requires overcoming the challenges of isolating rare cells and sequencing low-input material. Here we report an integrated process to isolate, qualify and sequence whole exomes of CTCs with high fidelity using a census-based sequencing strategy. Power calculations suggest that mapping of >99.995% of the standard exome is possible in CTCs. We validated our process in two patients with prostate cancer, including one for whom we sequenced CTCs, a lymph node metastasis and nine cores of the primary tumor. Fifty-one of 73 CTC mutations (70%) were present in matched tissue. Moreover, we identified 10 early trunk and 56 metastatic trunk mutations in the non-CTC tumor samples and found 90% and 73% of these mutations, respectively, in CTC exomes. This study establishes a foundation for CTC genomics in the clinic.

Tumour expression of selected microRNAs (miRs) correlates with cisplatin efficacy in multiple cancers. We investigated the role of selected miRs in patients receiving cisplatin-based therapy for advanced urothelial carcinoma (UC).

Absence of pathologic complete response (pCR) to neoadjuvant chemotherapy (NACT) correlates with poor long-term survival in patients with triple negative breast cancer (TNBC). These incomplete treatment responses are likely determined by mechanisms that enable cancer cells to resist being killed. However, the detailed characterization of a drug-resistant cancer cell state in residual TNBC tissue after NACT has remained elusive. AKT1low quiescent cancer cells (QCCs) are a quiescent, epigenetically plastic, and chemotherapy-resistant subpopulation initially identified in experimental cancer models. Here, we asked whether QCCs exist in primary tumors from patients with TNBC and persist after treatment with NACT.

While progression from normal prostatic epithelium to invasive cancer is driven by molecular alterations, tumor cells and cells in the cancer microenvironment are co-dependent and co-evolve. Few human studies to date have focused on stroma. Here, we performed gene expression profiling of laser capture microdissected normal non-neoplastic prostate epithelial tissue and compared it to non-transformed and neoplastic low-grade and high-grade prostate epithelial tissue from radical prostatectomies, each with its immediately surrounding stroma. Whereas benign epithelium in prostates with and without tumor were similar in gene expression space, stroma away from tumor was significantly different from that in prostates without cancer. A stromal gene signature reflecting bone remodeling and immune-related pathways was upregulated in high compared to low-Gleason grade cases. In validation data, the signature discriminated cases that developed metastasis from those that did not. These data suggest that the microenvironment may influence prostate cancer initiation, maintenance, and metastatic progression.Stromal cells contribute to tumor development but the mechanisms regulating this process are still unclear. Here the authors analyze gene expression profiles in the prostate and show that stromal gene signature changes ahead of the epithelial gene signature as prostate cancer initiates and progresses.

Drugs that block the activity of the methyltransferase EZH2 are in clinical development for the treatment of non-Hodgkin lymphomas harboring EZH2 gain-of-function mutations that enhance its polycomb repressive function. We have previously reported that EZH2 can act as a transcriptional activator in castration-resistant prostate cancer (CRPC). Now we show that EZH2 inhibitors can also block the transactivation activity of EZH2 and inhibit the growth of CRPC cells. Gene expression and epigenomics profiling of cells treated with EZH2 inhibitors demonstrated that in addition to derepressing gene expression, these compounds also robustly down-regulate a set of DNA damage repair (DDR) genes, especially those involved in the base excision repair (BER) pathway. Methylation of the pioneer factor FOXA1 by EZH2 contributes to the activation of these genes, and interaction with the transcriptional coactivator P300 via the transactivation domain on EZH2 directly turns on the transcription. In addition, CRISPR-Cas9-mediated knockout screens in the presence of EZH2 inhibitors identified these BER genes as the determinants that underlie the growth-inhibitory effect of EZH2 inhibitors. Interrogation of public data from diverse types of solid tumors expressing wild-type EZH2 demonstrated that expression of DDR genes is significantly correlated with EZH2 dependency and cellular sensitivity to EZH2 inhibitors. Consistent with these findings, treatment of CRPC cells with EZH2 inhibitors dramatically enhances their sensitivity to genotoxic stress. These studies reveal a previously unappreciated mechanism of action of EZH2 inhibitors and provide a mechanistic basis for potential combination cancer therapies.

Understanding the biologic mechanisms underlying the development of lethal prostate cancer is critical for improved therapeutic and prevention strategies. In this study we explored the role of tumor metabolism in prostate cancer progression using mRNA expression profiling of seven metabolic pathways; fatty acid metabolism, glycolysis/gluconeogenesis, oxidative phosphorylation, pentose phosphate, purine metabolism, pyrimidine metabolism and the tricarboxylic acid cycle.

The primary objective of this study is to detect biomarkers and develop models that enable the identification of clinically significant prostate cancer and to understand the biologic implications of the genes involved. Peripheral blood samples (1018 patients) were split chronologically into independent training (n = 713) and validation (n = 305) sets. Whole transcriptome RNA sequencing was performed on isolated phagocytic CD14+ and non-phagocytic CD2+ cells and their gene expression levels were used to develop predictive models that correlate to adverse pathologic features. The immune-transcriptomic model with the highest performance for predicting adverse pathology, based on a subtraction of the log-transformed expression signals of the two cell types, displayed an area under the curve (AUC) of the receiver operating characteristic of 0.70. The addition of biomarkers in combination with traditional clinical risk factors (age, serum prostate-specific antigen (PSA), PSA density, race, digital rectal examination (DRE), and family history) enhanced the AUC to 0.91 and 0.83 for the training and validation sets, respectively. The markers identified by this approach uncovered specific pathway associations relevant to (prostate) cancer biology. Increased phagocytic activity in conjunction with cancer-associated (mis-)regulation is also represented by these markers. Differential gene expression of circulating immune cells gives insight into the cellular immune response to early tumor development and immune surveillance.

To date, there have been no reports characterizing the genome-wide somatic DNA chromosomal copy-number alteration landscape in metastatic urothelial carcinoma. We sought to characterize the DNA copy-number profile in a cohort of metastatic samples and compare them to a cohort of primary urothelial carcinoma samples in order to identify changes that are associated with progression from primary to metastatic disease.

Mounting evidence suggests that long noncoding RNAs (lncRNAs) can function as microRNA sponges and compete for microRNA binding to protein-coding transcripts. However, the prevalence, functional significance and targets of lncRNA-mediated sponge regulation of cancer are mostly unknown. Here we identify a lncRNA-mediated sponge regulatory network that affects the expression of many protein-coding prostate cancer driver genes, by integrating analysis of sequence features and gene expression profiles of both lncRNAs and protein-coding genes in tumours. We confirm the tumour-suppressive function of two lncRNAs (TUG1 and CTB-89H12.4) and their regulation of PTEN expression in prostate cancer. Surprisingly, one of the two lncRNAs, TUG1, was previously known for its function in polycomb repressive complex 2 (PRC2)-mediated transcriptional regulation, suggesting its sub-cellular localization-dependent function. Our findings not only suggest an important role of lncRNA-mediated sponge regulation in cancer, but also underscore the critical influence of cytoplasmic localization on the efficacy of a sponge lncRNA.