Chromosomal instability in untreated primary prostate cancer as an indicator of metastatic potential.

BACKGROUND: Metastatic prostate cancer (PC) is highly lethal. The ability to identify primary tumors capable of dissemination is an unmet need in the quest to understand lethal biology and improve patient outcomes. Previous studies have linked chromosomal instability (CIN), which generates aneuploidy following chromosomal missegregation during mitosis, to PC progression. Evidence of CIN includes broad copy number alterations (CNAs) spanning > 300 base pairs of DNA, which may also be measured via RNA expression signatures associated with CNA frequency. Signatures of CIN in metastatic PC, however, have not been interrogated or well defined. We examined a published 70-gene CIN signature (CIN70) in untreated and castration-resistant prostate cancer (CRPC) cohorts from The Cancer Genome Atlas (TCGA) and previously published reports. We also performed transcriptome and CNA analysis in a unique cohort of untreated primary tumors collected from diagnostic prostate needle biopsies (PNBX) of localized (M0) and metastatic (M1) cases to determine if CIN was linked to clinical stage and outcome. METHODS: PNBX were collected from 99 patients treated in the VA Greater Los Angeles (GLA-VA) Healthcare System between 2000 and 2016. Total RNA was extracted from high-grade cancer areas in PNBX cores, followed by RNA sequencing and/or copy number analysis using OncoScan. Multivariate logistic regression analyses permitted calculation of odds ratios for CIN status (high versus low) in an expanded GLA-VA PNBX cohort (n = 121). RESULTS: The CIN70 signature was significantly enriched in primary tumors and CRPC metastases from M1 PC cases. An intersection of gene signatures comprised of differentially expressed genes (DEGs) generated through comparison of M1 versus M0 PNBX and primary CRPC tumors versus metastases revealed a 157-gene "metastasis" signature that was further distilled to 7-genes (PC-CIN) regulating centrosomes, chromosomal segregation, and mitotic spindle assembly. High PC-CIN scores correlated with CRPC, PC-death and all-cause mortality in the expanded GLA-VA PNBX cohort. Interestingly, approximately 1/3 of M1 PNBX cases exhibited low CIN, illuminating differential pathways of lethal PC progression. CONCLUSIONS: Measuring CIN in PNBX by transcriptome profiling is feasible, and the PC-CIN signature may identify patients with a high risk of lethal progression at the time of diagnosis.

Probabilistic modeling of personalized drug combinations from integrated chemical screen and molecular data in sarcoma.

BACKGROUND: Cancer patients with advanced disease routinely exhaust available clinical regimens and lack actionable genomic medicine results, leaving a large patient population without effective treatments options when their disease inevitably progresses. To address the unmet clinical need for evidence-based therapy assignment when standard clinical approaches have failed, we have developed a probabilistic computational modeling approach which integrates molecular sequencing data with functional assay data to develop patient-specific combination cancer treatments. METHODS: Tissue taken from a murine model of alveolar rhabdomyosarcoma was used to perform single agent drug screening and DNA/RNA sequencing experiments; results integrated via our computational modeling approach identified a synergistic personalized two-drug combination. Cells derived from the primary murine tumor were allografted into mouse models and used to validate the personalized two-drug combination. Computational modeling of single agent drug screening and RNA sequencing of multiple heterogenous sites from a single patient's epithelioid sarcoma identified a personalized two-drug combination effective across all tumor regions. The heterogeneity-consensus combination was validated in a xenograft model derived from the patient's primary tumor. Cell cultures derived from human and canine undifferentiated pleomorphic sarcoma were assayed by drug screen; computational modeling identified a resistance-abrogating two-drug combination common to both cell cultures. This combination was validated in vitro via a cell regrowth assay. RESULTS: Our computational modeling approach addresses three major challenges in personalized cancer therapy: synergistic drug combination predictions (validated in vitro and in vivo in a genetically engineered murine cancer model), identification of unifying therapeutic targets to overcome intra-tumor heterogeneity (validated in vivo in a human cancer xenograft), and mitigation of cancer cell resistance and rewiring mechanisms (validated in vitro in a human and canine cancer model). CONCLUSIONS: These proof-of-concept studies support the use of an integrative functional approach to personalized combination therapy prediction for the population of high-risk cancer patients lacking viable clinical options and without actionable DNA sequencing-based therapy.

The mutational landscape of lethal castration-resistant prostate cancer.

Characterization of the prostate cancer transcriptome and genome has identified chromosomal rearrangements and copy number gains and losses, including ETS gene family fusions, PTEN loss and androgen receptor (AR) amplification, which drive prostate cancer development and progression to lethal, metastatic castration-resistant prostate cancer (CRPC). However, less is known about the role of mutations. Here we sequenced the exomes of 50 lethal, heavily pre-treated metastatic CRPCs obtained at rapid autopsy (including three different foci from the same patient) and 11 treatment-naive, high-grade localized prostate cancers. We identified low overall mutation rates even in heavily treated CRPCs (2.00 per megabase) and confirmed the monoclonal origin of lethal CRPC. Integrating exome copy number analysis identified disruptions of CHD1 that define a subtype of ETS gene family fusion-negative prostate cancer. Similarly, we demonstrate that ETS2, which is deleted in approximately one-third of CRPCs (commonly through TMPRSS2:ERG fusions), is also deregulated through mutation. Furthermore, we identified recurrent mutations in multiple chromatin- and histone-modifying genes, including MLL2 (mutated in 8.6% of prostate cancers), and demonstrate interaction of the MLL complex with the AR, which is required for AR-mediated signalling. We also identified novel recurrent mutations in the AR collaborating factor FOXA1, which is mutated in 5 of 147 (3.4%) prostate cancers (both untreated localized prostate cancer and CRPC), and showed that mutated FOXA1 represses androgen signalling and increases tumour growth. Proteins that physically interact with the AR, such as the ERG gene fusion product, FOXA1, MLL2, UTX (also known as KDM6A) and ASXL1 were found to be mutated in CRPC. In summary, we describe the mutational landscape of a heavily treated metastatic cancer, identify novel mechanisms of AR signalling deregulated in prostate cancer, and prioritize candidates for future study.

Methods and challenges in timing chromosomal abnormalities within cancer samples.

MOTIVATION: Tumors acquire many chromosomal amplifications, and those acquired early in the lifespan of the tumor may be not only important for tumor growth but also can be used for diagnostic purposes. Many methods infer the order of the accumulation of abnormalities based on their occurrence in a large cohort of patients. Recently, Durinck et al. (2011) and Greenman et al. (2012) developed methods to order a single tumor's chromosomal amplifications based on the patterns of mutations accumulated within those regions. This method offers an unprecedented opportunity to assess the etiology of a single tumor sample, but has not been widely evaluated. RESULTS: We show that the model for timing chromosomal amplifications is limited in scope, particularly for regions with high levels of amplification. We also show that the estimation of the order of events can be sensitive for events that occur early in the progression of the tumor and that the partial maximum likelihood method of Greenman et al. (2012) can give biased estimates, particularly for moderate read coverage or normal contamination. We propose a maximum-likelihood estimation procedure that fully accounts for sequencing variability and show that it outperforms the partial maximum-likelihood estimation method. We also propose a Bayesian estimation procedure that stabilizes the estimates in certain settings. We implement these methods on a small number of ovarian tumors, and the results suggest possible differences in how the tumors acquired amplifications. AVAILABILITY AND IMPLEMENTATION: We provide implementation of these methods in an R package cancerTiming, which is available from the Comprehensive R Archive Network (CRAN) at http://CRAN.R-project.org/.

PAK4 inhibition improves PD1 blockade immunotherapy in prostate cancer by increasing immune infiltration.

Antitumor immunity requires lymphocytes to localize to the tumor. Prostate cancers (PCs) are immunologically cold and tend to lack T-cell infiltration. Most advanced PCs are insensitive to PD1 blockade therapies. Using syngeneic RM1 prostate tumors, p21-activated kinase-4 (PAK4) knockdown (KD) and pharmacological inhibition was assessed in C57BL/6J mice treated with PD1 antibodies (αPD1). RNASeq was used to characterize the immune response in the tumor. Immunohistochemistry, flow cytometry, and in vivo blocking studies confirmed the role of cell surface proteins in the generation of immune responses. In The Cancer Genome Atlas, PAK4 expression was inversely correlated with immune cell infiltration. PAK4 expression was controlled by the androgen receptor and its pioneering factor, FOXA1. PAK4 KD increased CD8+ T-cell infiltration and expression of IFNγ response genes. PAK4 KD also upregulated angiogenesis and endothelial cell adhesion molecules in the tumor microenvironment, contributing to CD8+ lymphocyte recruitment. Pharmacological inhibition of PAK4 made PC more responsive to immunotherapy with αPD1. A decrease in PAK4 activity increases immune activation and vascularity, which increases CD8+ lymphocyte infiltration into the tumor. Therefore, targeting PAK4 may improve the response of human PC to immunotherapy.

Deep Functional and Molecular Characterization of a High-Risk Undifferentiated Pleomorphic Sarcoma.

Nonrhabdomyosarcoma soft-tissue sarcomas (STSs) are a class of 50+ cancers arising in muscle and soft tissues of children, adolescents, and adults. Rarity of each subtype often precludes subtype-specific preclinical research, leaving many STS patients with limited treatment options should frontline therapy be insufficient. When clinical options are exhausted, personalized therapy assignment approaches may help direct patient care. Here, we report the results of an adult female STS patient with relapsed undifferentiated pleomorphic sarcoma (UPS) who self-drove exploration of a wide array of personalized Clinical Laboratory Improvement Amendments (CLIAs) level and research-level diagnostics, including state of the art genomic, proteomic, ex vivo live cell chemosensitivity testing, a patient-derived xenograft model, and immunoscoring. Her therapeutic choices were also diverse, including neoadjuvant chemotherapy, radiation therapy, and surgeries. Adjuvant and recurrence strategies included off-label and natural medicines, several immunotherapies, and N-of-1 approaches. Identified treatment options, especially those validated during the in vivo study, were not introduced into the course of clinical treatment but did provide plausible treatment regimens based on FDA-approved clinical agents.

PAK4 inhibition improves PD-1 blockade immunotherapy.

Lack of tumor infiltration by immune cells is the main mechanism of primary resistance to programmed cell death protein 1 (PD-1) blockade therapies for cancer. It has been postulated that cancer cell-intrinsic mechanisms may actively exclude T cells from tumors, suggesting that the finding of actionable molecules that could be inhibited to increase T cell infiltration may synergize with checkpoint inhibitor immunotherapy. Here, we show that p21-activated kinase 4 (PAK4) is enriched in non-responding tumor biopsies with low T cell and dendritic cell infiltration. In mouse models, genetic deletion of PAK4 increased T cell infiltration and reversed resistance to PD-1 blockade in a CD8 T cell-dependent manner. Furthermore, combination of anti-PD-1 with the PAK4 inhibitor KPT-9274 improved anti-tumor response compared with anti-PD-1 alone. Therefore, high PAK4 expression is correlated with low T cell and dendritic cell infiltration and a lack of response to PD-1 blockade, which could be reversed with PAK4 inhibition.

Integrated molecular analysis of Tamoxifen-resistant invasive lobular breast cancer cells identifies MAPK and GRM/mGluR signaling as therapeutic vulnerabilities.

Invasive lobular breast cancer (ILC) is an understudied malignancy with distinct clinical, pathological, and molecular features that distinguish it from the more common invasive ductal carcinoma (IDC). Mounting evidence suggests that estrogen receptor-alpha positive (ER+) ILC has a poor response to Tamoxifen (TAM), but the mechanistic drivers of this are undefined. In the current work, we comprehensively characterize the SUM44/LCCTam ILC cell model system through integrated analysis of gene expression, copy number, and mutation, with the goal of identifying actionable alterations relevant to clinical ILC that can be co-targeted along with ER to improve treatment outcomes. We show that TAM has several distinct effects on the transcriptome of LCCTam cells, that this resistant cell model has acquired copy number alterations and mutations that impinge on MAPK and metabotropic glutamate receptor (GRM/mGluR) signaling networks, and that pharmacological inhibition of either improves or restores the growth-inhibitory actions of endocrine therapy.

IL-13 receptors as possible therapeutic targets in diffuse intrinsic pontine glioma.

Diffuse intrinsic pontine glioma (DIPG) is a universally fatal childhood cancer of the brain. Despite the introduction of conventional chemotherapy and radiotherapy, improvements in survival have been marginal and long-term survivorship is uncommon. Thus, new targets for therapeutics are critically needed. Early phase clinical trials exploring molecularly-targeted therapies against the epidermal growth factor receptor (EGFR) and novel immunotherapies targeting interleukin receptor-13α2 (IL-13Rα2) have demonstrated activity in this disease. To identify additional therapeutic markers for cell surface receptors, we performed exome sequencing (16 new samples, 22 previously published samples, total 38 with 26 matched normal DNA samples), RNA deep sequencing (17 new samples, 11 previously published samples, total 28 with 18 matched normal RNA samples), and immunohistochemistry (17 DIPG tissue samples) to examine the expression of the interleukin-4 (IL-4) signaling axis components (IL-4, interleukin 13 (IL-13), and their respective receptors IL-4Rα, IL-13Rα1, and IL-13Rα2). In addition, we correlated cytokine and receptor expression with expression of the oncogenes EGFR and c-MET. In DIPG tissues, transcript-level analysis found significant expression of IL-4, IL-13, and IL-13Rα1/2, with strong differential expression of IL-13Rα1/2 in tumor versus normal brain. At the protein level, immunohistochemical studies revealed high content of IL-4 and IL-13Rα1/2 but notably low expression of IL-13. Additionally, a strong positive correlation was observed between c-Met and IL-4Rα. The genomic and transcriptional landscape across all samples was also summarized. These data create a foundation for the design of potential new immunotherapies targeting IL-13 cell surface receptors in DIPG.

Rapid, ultra low coverage copy number profiling of cell-free DNA as a precision oncology screening strategy.

Current cell-free DNA (cfDNA) next generation sequencing (NGS) precision oncology workflows are typically limited to targeted and/or disease-specific applications. In advanced cancer, disease burden and cfDNA tumor content are often elevated, yielding unique precision oncology opportunities. We sought to demonstrate the utility of a pan-cancer, rapid, inexpensive, whole genome NGS of cfDNA approach (PRINCe) as a precision oncology screening strategy via ultra-low coverage (~0.01x) tumor content determination through genome-wide copy number alteration (CNA) profiling. We applied PRINCe to a retrospective cohort of 124 cfDNA samples from 100 patients with advanced cancers, including 76 men with metastatic castration-resistant prostate cancer (mCRPC), enabling cfDNA tumor content approximation and actionable focal CNA detection, while facilitating concordance analyses between cfDNA and tissue-based NGS profiles and assessment of cfDNA alteration associations with mCRPC treatment outcomes. Therapeutically relevant focal CNAs were present in 42 (34%) cfDNA samples, including 36 of 93 (39%) mCRPC patient samples harboring AR amplification. PRINCe identified pre-treatment cfDNA CNA profiles facilitating disease monitoring. Combining PRINCe with routine targeted NGS of cfDNA enabled mutation and CNA assessment with coverages tuned to cfDNA tumor content. In mCRPC, genome-wide PRINCe cfDNA and matched tissue CNA profiles showed high concordance (median Pearson correlation = 0.87), and PRINCe detectable AR amplifications predicted reduced time on therapy, independent of therapy type (Kaplan-Meier log-rank test, chi-square = 24.9, p < 0.0001). Our screening approach enables robust, broadly applicable cfDNA-based precision oncology for patients with advanced cancer through scalable identification of therapeutically relevant CNAs and pre-/post-treatment genomic profiles, enabling cfDNA- or tissue-based precision oncology workflow optimization.

Tumor evolution and progression in multifocal and paired non-invasive/invasive urothelial carcinoma.

Although multifocal tumors and non-invasive/invasive components are commonly encountered in surgical pathology, their genetic relationship is often poorly characterized. We used next-generation sequencing (NGS) to characterize somatic alterations in a patient with five spatially distinct, high-grade papillary urothelial carcinomas (UCs), with one tumor harboring an underlying invasive component. NGS of 409 cancer-related genes was performed on DNA isolated from formalin-fixed paraffin-embedded (FFPE) blocks representing each papillary tumor (n = 5), the invasive component of one tumor, and matched normal tissue. We identified nine unique non-synonymous somatic mutations across the six UC samples, including five present in each carcinoma sample, consistent with clonal origin and limited intertumoral heterogeneity. Copy number and loss of heterogeneity (LOH) profiles were similar in all six carcinomas; however, the invasive carcinoma component uniquely showed focal CDKN2A loss and chromosome 9 LOH and did not harbor gains of chromosomes 5p or X that were present in the other tumor samples. Phylogenetic analysis supported the invasive component arising from a shared progenitor prior to the outgrowth of cells in the non-invasive tumors. Results were extended to three additional cases of upper tract UC with paired non-invasive/invasive components, which identified driving alterations exclusive to both non-invasive and invasive components. Lastly, we performed targeted RNA sequencing (RNAseq) using a custom bladder cancer panel, which confirmed gene expression signature differences between paired non-invasive/invasive components. The results and approaches presented here may be useful in understanding the clonal relationships in multifocal cancers or paired non-invasive/invasive components from routine FFPE specimens.

Next-Gen Sequencing Exposes Frequent MED12 Mutations and Actionable Therapeutic Targets in Phyllodes Tumors.

UNLABELLED: Phyllodes tumors are rare fibroepithelial tumors with variable clinical behavior accounting for a small subset of all breast neoplasms, yet little is known about the genetic alterations that drive tumor initiation and/or progression. Here, targeted next-generation sequencing (NGS) was used to identify somatic alterations in formalin-fixed paraffin-embedded (FFPE) patient specimens from malignant, borderline, and benign cases. NGS revealed mutations in mediator complex subunit 12 (MED12) affecting the G44 hotspot residue in the majority (67%) of cases spanning all three histologic grades. In addition, loss-of-function mutations in p53 (TP53) as well as deleterious mutations in the tumor suppressors retinoblastoma (RB1) and neurofibromin 1 (NF1) were identified exclusively in malignant tumors. High-level copy-number alterations (CNA) were nearly exclusively confined to malignant tumors, including potentially clinically actionable gene amplifications in IGF1R and EGFR. Taken together, this study defines the genomic landscape underlying phyllodes tumor development, suggests potential molecular correlates to histologic grade, expands the spectrum of human tumors with frequent recurrent MED12 mutations, and identifies IGF1R and EGFR as potential therapeutic targets in malignant cases. IMPLICATIONS: Integrated genomic sequencing and mutational profiling provides insight into the molecular origin of phyllodes tumors and indicates potential druggable targets in malignant disease. Visual Overview: http://mcr.aacrjournals.org/content/early/2015/04/02/1541-7786.MCR-14-0578/F1.large.jpg.

Functionally defined therapeutic targets in diffuse intrinsic pontine glioma.

Diffuse intrinsic pontine glioma (DIPG) is a fatal childhood cancer. We performed a chemical screen in patient-derived DIPG cultures along with RNA-seq analyses and integrated computational modeling to identify potentially effective therapeutic strategies. The multi-histone deacetylase inhibitor panobinostat demonstrated therapeutic efficacy both in vitro and in DIPG orthotopic xenograft models. Combination testing of panobinostat and the histone demethylase inhibitor GSK-J4 revealed that the two had synergistic effects. Together, these data suggest a promising therapeutic strategy for DIPG.

Development and validation of a scalable next-generation sequencing system for assessing relevant somatic variants in solid tumors.

Next-generation sequencing (NGS) has enabled genome-wide personalized oncology efforts at centers and companies with the specialty expertise and infrastructure required to identify and prioritize actionable variants. Such approaches are not scalable, preventing widespread adoption. Likewise, most targeted NGS approaches fail to assess key relevant genomic alteration classes. To address these challenges, we predefined the catalog of relevant solid tumor somatic genome variants (gain-of-function or loss-of-function mutations, high-level copy number alterations, and gene fusions) through comprehensive bioinformatics analysis of >700,000 samples. To detect these variants, we developed the Oncomine Comprehensive Panel (OCP), an integrative NGS-based assay [compatible with <20 ng of DNA/RNA from formalin-fixed paraffin-embedded (FFPE) tissues], coupled with an informatics pipeline to specifically identify relevant predefined variants and created a knowledge base of related potential treatments, current practice guidelines, and open clinical trials. We validated OCP using molecular standards and more than 300 FFPE tumor samples, achieving >95% accuracy for KRAS, epidermal growth factor receptor, and BRAF mutation detection as well as for ALK and TMPRSS2:ERG gene fusions. Associating positive variants with potential targeted treatments demonstrated that 6% to 42% of profiled samples (depending on cancer type) harbored alterations beyond routine molecular testing that were associated with approved or guideline-referenced therapies. As a translational research tool, OCP identified adaptive CTNNB1 amplifications/mutations in treated prostate cancers. Through predefining somatic variants in solid tumors and compiling associated potential treatment strategies, OCP represents a simplified, broadly applicable targeted NGS system with the potential to advance precision oncology efforts.

HRAS mutations are frequent in inverted urothelial neoplasms.

Inverted urothelial papilloma (IUP) is an uncommon neoplasm of the urinary bladder with distinct morphologic features. Studies regarding the role of human papillomavirus (HPV) in the etiology of IUP have provided conflicting evidence of HPV infection. In addition, little is known regarding the molecular alterations present in IUP or other urothelial neoplasms, which might demonstrate inverted growth pattern like low-grade or high-grade urothelial carcinoma (UCA). Here, we evaluated for the presence of common driving somatic mutations and HPV within a cohort of IUPs, (n = 7) noninvasive low-grade papillary UCAs with inverted growth pattern (n = 5), and noninvasive high-grade papillary UCAs with inverted growth pattern (n = 8). HPV was not detected in any case of IUP or inverted UCA by either in situ hybridization or by polymerase chain reaction. Next-generation sequencing identified recurrent mutations in HRAS (Q61R) in 3 of 5 IUPs, described for the first time in this neoplasm. Additional mutations of Ras pathway members were detected including HRAS, KRAS, and BRAF. The presence of Ras pathway member mutations at a relatively high rate suggests this pathway may contribute to pathogenesis of inverted urothelial neoplasms. In addition, we did not find any evidence supporting a role for HPV in the etiology of IUP.

Detection of somatic copy number alterations in cancer using targeted exome capture sequencing.

The research community at large is expending considerable resources to sequence the coding region of the genomes of tumors and other human diseases using targeted exome capture (i.e., "whole exome sequencing"). The primary goal of targeted exome sequencing is to identify nonsynonymous mutations that potentially have functional consequences. Here, we demonstrate that whole-exome sequencing data can also be analyzed for comprehensively monitoring somatic copy number alterations (CNAs) by benchmarking the technique against conventional array CGH. A series of 17 matched tumor and normal tissues from patients with metastatic castrate-resistant prostate cancer was used for this assessment. We show that targeted exome sequencing reliably identifies CNAs that are common in advanced prostate cancer, such as androgen receptor (AR) gain and PTEN loss. Taken together, these data suggest that targeted exome sequencing data can be effectively leveraged for the detection of somatic CNAs in cancer.

Personalized oncology through integrative high-throughput sequencing: a pilot study.

Individual cancers harbor a set of genetic aberrations that can be informative for identifying rational therapies currently available or in clinical trials. We implemented a pilot study to explore the practical challenges of applying high-throughput sequencing in clinical oncology. We enrolled patients with advanced or refractory cancer who were eligible for clinical trials. For each patient, we performed whole-genome sequencing of the tumor, targeted whole-exome sequencing of tumor and normal DNA, and transcriptome sequencing (RNA-Seq) of the tumor to identify potentially informative mutations in a clinically relevant time frame of 3 to 4 weeks. With this approach, we detected several classes of cancer mutations including structural rearrangements, copy number alterations, point mutations, and gene expression alterations. A multidisciplinary Sequencing Tumor Board (STB) deliberated on the clinical interpretation of the sequencing results obtained. We tested our sequencing strategy on human prostate cancer xenografts. Next, we enrolled two patients into the clinical protocol and were able to review the results at our STB within 24 days of biopsy. The first patient had metastatic colorectal cancer in which we identified somatic point mutations in NRAS, TP53, AURKA, FAS, and MYH11, plus amplification and overexpression of cyclin-dependent kinase 8 (CDK8). The second patient had malignant melanoma, in which we identified a somatic point mutation in HRAS and a structural rearrangement affecting CDKN2C. The STB identified the CDK8 amplification and Ras mutation as providing a rationale for clinical trials with CDK inhibitors or MEK (mitogen-activated or extracellular signal-regulated protein kinase kinase) and PI3K (phosphatidylinositol 3-kinase) inhibitors, respectively. Integrative high-throughput sequencing of patients with advanced cancer generates a comprehensive, individual mutational landscape to facilitate biomarker-driven clinical trials in oncology.