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.

The long noncoding RNA landscape of neuroendocrine prostate cancer and its clinical implications.

Background: Treatment-induced neuroendocrine prostate cancer (tNEPC) is an aggressive variant of late-stage metastatic castrate-resistant prostate cancer that commonly arises through neuroendocrine transdifferentiation (NEtD). Treatment options are limited, ineffective, and, for most patients, result in death in less than a year. We previously developed a first-in-field patient-derived xenograft (PDX) model of NEtD. Longitudinal deep transcriptome profiling of this model enabled monitoring of dynamic transcriptional changes during NEtD and in the context of androgen deprivation. Long non-coding RNA (lncRNA) are implicated in cancer where they can control gene regulation. Until now, the expression of lncRNAs during NEtD and their clinical associations were unexplored. Results: We implemented a next-generation sequence analysis pipeline that can detect transcripts at low expression levels and built a genome-wide catalogue (n = 37,749) of lncRNAs. We applied this pipeline to 927 clinical samples and our high-fidelity NEtD model LTL331 and identified 821 lncRNAs in NEPC. Among these are 122 lncRNAs that robustly distinguish NEPC from prostate adenocarcinoma (AD) patient tumours. The highest expressed lncRNAs within this signature are H19, LINC00617, and SSTR5-AS1. Another 742 are associated with the NEtD process and fall into four distinct patterns of expression (NEtD lncRNA Class I, II, III, and IV) in our PDX model and clinical samples. Each class has significant (z-scores >2) and unique enrichment for transcription factor binding site (TFBS) motifs in their sequences. Enriched TFBS include (1) TP53 and BRN1 in Class I, (2) ELF5, SPIC, and HOXD1 in Class II, (3) SPDEF in Class III, (4) HSF1 and FOXA1 in Class IV, and (5) TWIST1 when merging Class III with IV. Common TFBS in all NEtD lncRNA were also identified and include E2F, REST, PAX5, PAX9, and STAF. Interrogation of the top deregulated candidates (n = 100) in radical prostatectomy adenocarcinoma samples with long-term follow-up (median 18 years) revealed significant clinicopathological associations. Specifically, we identified 25 that are associated with rapid metastasis following androgen deprivation therapy (ADT). Two of these lncRNAs (SSTR5-AS1 and LINC00514) stratified patients undergoing ADT based on patient outcome. Discussion: To date, a comprehensive characterization of the dynamic landscape of lncRNAs during the NEtD process has not been performed. A temporal analysis of the PDX-based NEtD model has for the first time provided this dynamic landscape. TFBS analysis identified NEPC-related TF motifs present within the NEtD lncRNA sequences, suggesting functional roles for these lncRNAs in NEPC pathogenesis. Furthermore, select NEtD lncRNAs appear to be associated with metastasis and patients receiving ADT. Treatment-related metastasis is a clinical consequence of NEPC tumours. Top candidate lncRNAs FENDRR, H19, LINC00514, LINC00617, and SSTR5-AS1 identified in this study are implicated in the development of NEPC. We present here for the first time a genome-wide catalogue of NEtD lncRNAs that characterize the transdifferentiation process and a robust NEPC lncRNA patient expression signature. To accomplish this, we carried out the largest integrative study that applied a PDX NEtD model to clinical samples. These NEtD and NEPC lncRNAs are strong candidates for clinical biomarkers and therapeutic targets and warrant further investigation.

Genomic Alterations in Cell-Free DNA and Enzalutamide Resistance in Castration-Resistant Prostate Cancer.

IMPORTANCE: The molecular landscape underpinning response to the androgen receptor (AR) antagonist enzalutamide in patients with metastatic castration-resistant prostate cancer (mCRPC) is undefined. Consequently, there is an urgent need for practical biomarkers to guide therapy selection and elucidate resistance. Although tissue biopsies are impractical to perform routinely in the majority of patients with mCRPC, the analysis of plasma cell-free DNA (cfDNA) has recently emerged as a minimally invasive method to explore tumor characteristics. OBJECTIVE: To reveal genomic characteristics from cfDNA associated with clinical outcomes during enzalutamide treatment. DESIGN, SETTING, AND PARTICIPANTS: Plasma samples were obtained from August 4, 2013, to July 31, 2015, at a single academic institution (British Columbia Cancer Agency) from 65 patients with mCRPC. We collected temporal plasma samples (at baseline, 12 weeks, end of treatment) for circulating cfDNA and performed array comparative genomic hybridization copy number profiling and deep AR gene sequencing. Samples collected at end of treatment were also subjected to targeted sequencing of 19 prostate cancer-associated genes. EXPOSURE: Enzalutamide, 160 mg, daily orally. MAIN OUTCOMES AND MEASURES: Prostate-specific antigen response rate (decline ≥50% from baseline confirmed ≥3 weeks later). Radiographic (as per Prostate Cancer Working Group 2 Criteria) and/or clinical progression (defined as worsening disease-related symptoms necessitating a change in anticancer therapy and/or deterioration in Eastern Cooperative Group performance status ≥2 levels). RESULTS: The 65 patients had a median (interquartile range) age of 74 (68-79) years. Prostate-specific antigen response rate to enzalutamide treatment was 38% (25 of 65), while median clinical/radiographic progression-free survival was 3.5 (95% CI, 2.1-5.0) months. Cell-free DNA was isolated from 122 of 125 plasma samples, and targeted sequencing was successful in 119 of 122. AR mutations and/or copy number alterations were robustly detected in 48% (31 of 65) and 60% (18 of 30) of baseline and progression samples, respectively. Detection of AR amplification, heavily mutated AR (≥2 mutations), and RB1 loss were associated with worse progression-free survival, with hazard ratios of 2.92 (95% CI, 1.59-5.37), 3.94 (95% CI, 1.46-10.64), and 4.46 (95% CI, 2.28-8.74), respectively. AR mutations exhibited clonal selection during treatment, including an increase in glucocorticoid-sensitive AR L702H and promiscuous AR T878A in patients with prior abiraterone treatment. At the time of progression, cfDNA sequencing revealed mutations or copy number changes in all patients tested, including clinically actionable alterations in DNA damage repair genes and PI3K pathway genes, and a high frequency (4 of 14) of activating CTNNB1 mutations. CONCLUSIONS AND RELEVANCE: Clinically informative genomic profiling of cfDNA was feasible in nearly all patients with mCRPC and can provide important insights into enzalutamide response and resistance.

Androgen Receptor Gene Aberrations in Circulating Cell-Free DNA: Biomarkers of Therapeutic Resistance in Castration-Resistant Prostate Cancer.

PURPOSE: Although novel agents targeting the androgen-androgen receptor (AR) axis have altered the treatment paradigm of metastatic castration-resistant prostate cancer (mCRPC), development of therapeutic resistance is inevitable. In this study, we examined whether AR gene aberrations detectable in circulating cell-free DNA (cfDNA) are associated with resistance to abiraterone acetate and enzalutamide in mCRPC patients. EXPERIMENTAL DESIGN: Plasma was collected from 62 mCRPC patients ceasing abiraterone acetate (n = 29), enzalutamide (n = 19), or other agents (n = 14) due to disease progression. DNA was extracted and subjected to array comparative genomic hybridization (aCGH) for chromosome copy number analysis, and Roche 454 targeted next-generation sequencing of exon 8 in the AR. RESULTS: On aCGH, AR amplification was significantly more common in patients progressing on enzalutamide than on abiraterone or other agents (53% vs. 17% vs. 21%, P = 0.02, χ(2)). Missense AR exon 8 mutations were detected in 11 of 62 patients (18%), including the first reported case of an F876L mutation in an enzalutamide-resistant patient and H874Y and T877A mutations in 7 abiraterone-resistant patients. In patients switched onto enzalutamide after cfDNA collection (n = 39), an AR gene aberration (copy number increase and/or an exon 8 mutation) in pretreatment cfDNA was associated with adverse outcomes, including lower rates of PSA decline ≥ 30% (P = 0.013, χ(2)) and shorter time to radiographic/clinical progression (P = 0.010, Cox proportional hazards regression). CONCLUSIONS: AR gene aberrations in cfDNA are associated with resistance to enzalutamide and abiraterone in mCRPC. Our data illustrate that genomic analysis of cfDNA is a minimally invasive method for interrogating mechanisms of therapeutic resistance in mCRPC.

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.

Poly-gene fusion transcripts and chromothripsis in prostate cancer.

Complex genome rearrangements are frequently observed in cancer but their impact on tumor molecular biology is largely unknown. Recent studies have identified a new phenomenon involving the simultaneous generation of tens to hundreds of genomic rearrangements, called chromothripsis. To understand the molecular consequences of these events, we sequenced the genomes and transcriptomes of two prostate tumors exhibiting evidence of chromothripsis. We identified several complex fusion transcripts, each containing sequence from three different genes, originating from different parts of the genome. One such poly-gene fusion transcript appeared to be expressed from a chain of small genomic fragments. Furthermore, we detected poly-gene fusion transcripts in the prostate cancer cell line LNCaP, suggesting they may represent a common phenomenon. Finally in one tumor with chromothripsis, we identified multiple mutations in the p53 signaling pathway, expanding on recent work associating aberrant DNA damage response mechanisms with chromothripsis. Overall, our data show that chromothripsis can manifest as massively rearranged transcriptomes. The implication that multigenic changes can give rise to poly-gene fusion transcripts is potentially of great significance to cancer genetics.

From sequence to molecular pathology, and a mechanism driving the neuroendocrine phenotype in prostate cancer.

The current paradigm of cancer care relies on predictive nomograms which integrate detailed histopathology with clinical data. However, when predictions fail, the consequences for patients are often catastrophic, especially in prostate cancer where nomograms influence the decision to therapeutically intervene. We hypothesized that the high dimensional data afforded by massively parallel sequencing (MPS) is not only capable of providing biological insights, but may aid molecular pathology of prostate tumours. We assembled a cohort of six patients with high-risk disease, and performed deep RNA and shallow DNA sequencing in primary tumours and matched metastases where available. Our analysis identified copy number abnormalities, accurately profiled gene expression levels, and detected both differential splicing and expressed fusion genes. We revealed occult and potentially dormant metastases, unambiguously supporting the patients' clinical history, and implicated the REST transcriptional complex in the development of neuroendocrine prostate cancer, validating this finding in a large independent cohort. We massively expand on the number of novel fusion genes described in prostate cancer; provide fresh evidence for the growing link between fusion gene aetiology and gene expression profiles; and show the utility of fusion genes for molecular pathology. Finally, we identified chromothripsis in a patient with chronic prostatitis. Our results provide a strong foundation for further development of MPS-based molecular pathology.

Integrated genome and transcriptome sequencing identifies a novel form of hybrid and aggressive prostate cancer.

Next-generation sequencing is making sequence-based molecular pathology and personalized oncology viable. We selected an individual initially diagnosed with conventional but aggressive prostate adenocarcinoma and sequenced the genome and transcriptome from primary and metastatic tissues collected prior to hormone therapy. The histology-pathology and copy number profiles were remarkably homogeneous, yet it was possible to propose the quadrant of the prostate tumour that likely seeded the metastatic diaspora. Despite a homogeneous cell type, our transcriptome analysis revealed signatures of both luminal and neuroendocrine cell types. Remarkably, the repertoire of expressed but apparently private gene fusions, including C15orf21:MYC, recapitulated this biology. We hypothesize that the amplification and over-expression of the stem cell gene MSI2 may have contributed to the stable hybrid cellular identity. This hybrid luminal-neuroendocrine tumour appears to represent a novel and highly aggressive case of prostate cancer with unique biological features and, conceivably, a propensity for rapid progression to castrate-resistance. Overall, this work highlights the importance of integrated analyses of genome, exome and transcriptome sequences for basic tumour biology, sequence-based molecular pathology and personalized oncology.

Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics.

We report 80,388 ESTs from 23 Atlantic salmon (Salmo salar) cDNA libraries (61,819 ESTs), 6 rainbow trout (Oncorhynchus mykiss) cDNA libraries (14,544 ESTs), 2 chinook salmon (Oncorhynchus tshawytscha) cDNA libraries (1317 ESTs), 2 sockeye salmon (Oncorhynchus nerka) cDNA libraries (1243 ESTs), and 2 lake whitefish (Coregonus clupeaformis) cDNA libraries (1465 ESTs). The majority of these are 3' sequences, allowing discrimination between paralogs arising from a recent genome duplication in the salmonid lineage. Sequence assembly reveals 28,710 different S. salar, 8981 O. mykiss, 1085 O. tshawytscha, 520 O. nerka, and 1176 C. clupeaformis putative transcripts. We annotate the submitted portion of our EST database by molecular function. Higher- and lower-molecular-weight fractions of libraries are shown to contain distinct gene sets, and higher rates of gene discovery are associated with higher-molecular weight libraries. Pyloric caecum library group annotations indicate this organ may function in redox control and as a barrier against systemic uptake of xenobiotics. A microarray is described, containing 7356 salmonid elements representing 3557 different cDNAs. Analyses of cross-species hybridizations to this cDNA microarray indicate that this resource may be used for studies involving all salmonids.