Search

Background & Aims Continuing recalcitrance to therapy cements pancreatic cancer (PC) as the most lethal malignancy, which is set to become the second leading cause of cancer death in our society. The study aim was to investigate the association between DNA damage response (DDR), replication stress, and novel therapeutic response in PC to develop a biomarker-driven therapeutic strategy targeting DDR and replication stress in PC. Methods We interrogated the transcriptome, genome, proteome, and functional characteristics of 61 novel PC patient–derived cell lines to define novel therapeutic strategies targeting DDR and replication stress. Validation was done in patient-derived xenografts and human PC organoids. Results Patient-derived cell lines faithfully recapitulate the epithelial component of pancreatic tumors, including previously described molecular subtypes. Biomarkers of DDR deficiency, including a novel signature of homologous recombination deficiency, cosegregates with response to platinum (P < .001) and PARP inhibitor therapy (P < .001) in vitro and in vivo. We generated a novel signature of replication stress that predicts response to ATR (P < .018) and WEE1 inhibitor (P < .029) treatment in both cell lines and human PC organoids. Replication stress was enriched in the squamous subtype of PC (P < .001) but was not associated with DDR deficiency. Conclusions Replication stress and DDR deficiency are independent of each other, creating opportunities for therapy in DDR-proficient PC and after platinum therapy., Graphical abstract

Pancreatic ductal adenocarcinoma (PDAC) can be divided into transcriptomic subtypes with two broad lineages referred to as classical (pancreatic) and squamous. We find that these two subtypes are driven by distinct metabolic phenotypes. Loss of genes that drive endodermal lineage specification, HNF4A and GATA6, switch metabolic profiles from classical (pancreatic) to predominantly squamous, with glycogen synthase kinase 3 beta (GSK3β) a key regulator of glycolysis. Pharmacological inhibition of GSK3β results in selective sensitivity in the squamous subtype; however, a subset of these squamous patient-derived cell lines (PDCLs) acquires rapid drug tolerance. Using chromatin accessibility maps, we demonstrate that the squamous subtype can be further classified using chromatin accessibility to predict responsiveness and tolerance to GSK3β inhibitors. Our findings demonstrate that distinct patterns of chromatin accessibility can be used to identify patient subgroups that are indistinguishable by gene expression profiles, highlighting the utility of chromatin-based biomarkers for patient selection in the treatment of PDAC.

Purpose: Glioblastoma is the most frequent and lethal primary brain tumor. Development of novel therapies relies on the availability of relevant preclinical models. We have established a panel of 96 glioblastoma patient-derived xenografts (PDX) and undertaken its genomic and phenotypic characterization. Experimental Design: PDXs were established from glioblastoma, IDH-wildtype (n = 93), glioblastoma, IDH-mutant (n = 2), diffuse midline glioma, H3 K27M-mutant (n = 1), and both primary (n = 60) and recurrent (n = 34) tumors. Tumor growth rates, histopathology, and treatment response were characterized. Integrated molecular profiling was performed by whole-exome sequencing (WES, n = 83), RNA-sequencing (n = 68), and genome-wide methylation profiling (n = 76). WES data from 24 patient tumors was compared with derivative models. Results: PDXs recapitulate many key phenotypic and molecular features of patient tumors. Orthotopic PDXs show characteristic tumor morphology and invasion patterns, but largely lack microvascular proliferation and necrosis. PDXs capture common and rare molecular drivers, including alterations of TERT, EGFR, PTEN, TP53, BRAF, and IDH1, most at frequencies comparable with human glioblastoma. However, PDGFRA amplification was absent. RNA-sequencing and genome-wide methylation profiling demonstrated broad representation of glioblastoma molecular subtypes. MGMT promoter methylation correlated with increased survival in response to temozolomide. WES of 24 matched patient tumors showed preservation of most genetic driver alterations, including EGFR amplification. However, in four patient–PDX pairs, driver alterations were gained or lost on engraftment, consistent with clonal selection. Conclusions: Our PDX panel captures the molecular heterogeneity of glioblastoma and recapitulates many salient genetic and phenotypic features. All models and genomic data are openly available to investigators.

Background: Pancreatic ductal adenocarcinoma (PDA) is a lethal cancer with complex genomes and dense fibrotic stroma. This study was designed to identify clinically relevant somatic aberrations in pancreatic cancer genomes of patients with primary and metastatic disease enrolled and treated in two clinical trials. Methods: Tumour nuclei were flow sorted prior to whole genome copy number variant (CNV) analysis. Targeted or whole exome sequencing was performed on most samples. We profiled biopsies from 68 patients enrolled in two Stand Up to Cancer (SU2C)-sponsored clinical trials. These included 38 resected chemoradiation naïve tumours (SU2C 20206-003) and metastases from 30 patients who progressed on prior therapies (SU2C 20206-001). Patient outcomes including progression-free survival (PFS) and overall survival (OS) were observed. Results: We defined: (a) CDKN2A homozygous deletions that included the adjacent MTAP gene, only its’ 3′ region, or excluded MTAP; (b) SMAD4 homozygous deletions that included ME2; (c) a pancreas-specific MYC super-enhancer region; (d) DNA repair-deficient genomes; and (e) copy number aberrations present in PDA patients with long-term (⩾ 40 months) and short-term (⩽ 12 months) survival after surgical resection. Conclusions: We provide a clinically relevant framework for genomic drivers of PDA and for advancing novel treatments.

Next-generation sequencing is enabling molecularly guided therapy for many cancer types, yet failure rates remain relatively high in pancreatic cancer (PC). The aim of this study is to investigate the feasibility of genomic profiling using endoscopic ultrasound (EUS) biopsy samples to facilitate personalised therapy for PC. Ninty-five patients underwent additional research biopsies at the time of diagnostic EUS. Diagnostic formalin-fixed (FFPE) and fresh frozen EUS samples underwent DNA extraction, quantification and targeted gene sequencing. Whole genome (WGS) and RNA sequencing was performed as proof of concept. Only 2 patients (2%) with a diagnosis of PC had insufficient material for targeted sequencing in both FFPE and frozen specimens. Targeted panel sequencing (n=54) revealed mutations in PC genes (KRAS, GNAS, TP53, CDKN2A, SMAD4) in patients with histological evidence of PC, including potentially actionable mutations (BRCA1, BRCA2, ATM, BRAF). WGS (n=5) of EUS samples revealed mutational signatures that are potential biomarkers of therapeutic responsiveness. RNA sequencing (n=35) segregated patients into clinically relevant molecular subtypes based on transcriptome. Integrated multi-omic analysis of PC using standard EUS guided biopsies offers clinical utility to guide personalized therapy and study the molecular pathology in all patients with PC.

Purpose: We sought to identify the genomic abnormalities in squamous cell carcinomas (SCC) arising in ovarian mature cystic teratoma (MCT), a rare gynecological malignancy of poor prognosis. Experimental design: We performed copy number, mutational state, and zygosity analysis of 151 genes in SCC arising in MCT (n = 25) using next-generation sequencing. The presence of high-/intermediate-risk HPV genotypes was assessed by quantitative PCR. Genomic events were correlated with clinical features and outcome. Results: MCT had a low mutation burden with a mean of only one mutation per case. Zygosity analyses of MCT indicated four separate patterns, suggesting that MCT can arise from errors at various stages of oogenesis. A total of 244 abnormalities were identified in 79 genes in MCT-associated SCC, and the overall mutational burden was high (mean 10.2 mutations per megabase). No SCC was positive for HPV. The most frequently altered genes in SCC were TP53 (20/25 cases, 80%), PIK3CA (13/25 cases, 52%), and CDKN2A (11/25 cases, 44%). Mutation in TP53 was associated with improved overall survival. In 8 of 20 cases with TP53 mutations, two or more variants were identified, which were bi-allelic. Conclusions: Ovarian SCC arising in MCT has a high mutational burden, with TP53 mutation the most common abnormality. The presence of TP53 mutation is a good prognostic factor. SCC arising in MCT share similar mutation profiles to other SCC. Given their rarity, they should be included in basket studies that recruit patients with SCC of other organs. Clin Cancer Res; 23(24); 7633–40. ©2017 AACR.

The aggressiveness of pancreatic cancer is associated with the acquisition of mesenchymal characteristics by a subset of pancreatic cancer cells. The factors driving the development of this subset are not well understood. In this study, we tested the hypothesis that acquisition of a mesenchymal phenotype occurs selectively in tumor cells that harbor specific enabling genetic alterations. We obtained whole-genome comparative genomic hybridization (CGH) measurements on pancreatic cancer cell lines that have either an epithelial-like (17 cell lines) or a mesenchymal-like (9 cell lines) phenotype in vitro. The total amounts of amplifications and deletions were equivalent between the epithelial and mesenchymal groups, but 20 genes showed a major difference between the groups in prevalence of alterations. All 20 alterations (18 deletions and 2 amplifications) were more prevalent in the mesenchymal group, confirming the advanced nature of this cellular subtype. CDKN2A was altered in more than 50% of both groups, but co-deletions in neighboring genes, and concomitant loss of gene expression, were more prevalent in the mesenchymal group, suggesting that the size of the loss around CDKN2A affects cell phenotype. Whole-genome CGH on 11 primary cancer tissues revealed that the 20 genes were altered at a higher prevalence (up to 55% of the cases for certain genes) than randomly selected sets of 20 genes, with the same direction of alteration as in the cell lines. These findings support the concept that specific genetic alterations enable phenotype plasticity and provide promising candidate genes for further research.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-β, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Liposarcoma is the most common soft tissue sarcoma, but little is known about the genomic basis of this disease. Given the low cell content of this tumor type, we utilized flow cytometry to isolate the diploid normal and aneuploid tumor populations from a well-differentiated liposarcoma prior to array comparative genomic hybridization and whole genome sequencing. This work revealed massive highly focal amplifications throughout the aneuploid tumor genome including MDM2, a gene that has previously been found to be amplified in well-differentiated liposarcoma. Structural analysis revealed massive rearrangement of chromosome 12 and 11 gene fusions, some of which may be part of double minute chromosomes commonly present in well-differentiated liposarcoma. We identified a hotspot of genomic instability localized to a region of chromosome 12 that includes a highly conserved, putative L1 retrotransposon element, LOC100507498 which resides within a gene cluster (NAV3, SYT1, PAWR) where 6 of the 11 fusion events occurred. Interestingly, a potential gene fusion was also identified in amplified DDR2, which is a potential therapeutic target of kinase inhibitors such as dastinib, that are not routinely used in the treatment of patients with liposarcoma. Furthermore, 7 somatic, damaging single nucleotide variants have also been identified, including D125N in the PTPRQ protein. In conclusion, this work is the first to report the entire genome of a well-differentiated liposarcoma with novel chromosomal rearrangements associated with amplification of therapeutically targetable genes such as MDM2 and DDR2.

Glioblastoma (GBM) is known to be a heterogeneous disease; however, the genetic composition of the cells within a given tumour is only poorly explored. In the advent of personalised medicine the understanding of intra-tumoural heterogeneity at the cellular and the genetic level is mandatory to improve treatment and clinical outcome. By combining ploidy-based flow sorting with array-comparative genomic hybridization we show that primary GBMs present as either mono- or polygenomic tumours (64 versus 36 %, respectively). Monogenomic tumours were limited to a pseudodiploid tumour clone admixed with normal stromal cells, whereas polygenomic tumours contained multiple tumour clones, yet always including a pseudodiploid population. Interestingly, pseudodiploid and aneuploid fractions carried the same aberrations as defined by identical chromosomal breakpoints, suggesting that evolution towards aneuploidy is a late event in GBM development. Interestingly, while clonal heterogeneity could be recapitulated in spheroid-based xenografts, we find that genetically distinct clones displayed different tumourigenic potential. Moreover, we show that putative cancer stem cell markers including CD133, CD15, A2B5 and CD44 were present on genetically distinct tumour cell populations. These data reveal the clonal heterogeneity of GBMs at the level of DNA content, tumourigenic potential and stem cell marker expression, which is likely to impact glioma progression and treatment response. The combined knowledge of intra-tumour heterogeneity at the genetic, cellular and functional level is crucial to assess treatment responses and to design personalized treatment strategies for primary GBM. Electronic supplementary material The online version of this article (doi:10.1007/s00401-013-1196-4) contains supplementary material, which is available to authorized users.

Formalin fixed paraffin embedded (FFPE) tissues are a vast resource of annotated clinical samples. As such, they represent highly desirable and informative materials for the application of high definition genomics for improved patient management and to advance the development of personalized therapeutics. However, a limitation of FFPE tissues is the variable quality of DNA extracted for analyses. Furthermore, admixtures of non-tumor and polyclonal neoplastic cell populations limit the number of biopsies that can be studied and make it difficult to define cancer genomes in patient samples. To exploit these valuable tissues we applied flow cytometry-based methods to isolate pure populations of tumor cell nuclei from FFPE tissues and developed a methodology compatible with oligonucleotide array CGH and whole exome sequencing analyses. These were used to profile a variety of tumors (breast, brain, bladder, ovarian and pancreas) including the genomes and exomes of matching fresh frozen and FFPE pancreatic adenocarcinoma samples.

Cancers frequently arise as a result of an acquired genomic instability and the subsequent clonal evolution of neoplastic cells with variable patterns of genetic aberrations. Thus, the presence and behaviors of distinct clonal populations in each patient's tumor may underlie multiple clinical phenotypes in cancers. We applied DNA content-based flow sorting to identify and isolate the nuclei of clonal populations from tumor biopsies, which was coupled with array CGH and targeted resequencing. The results produced high-definition genomic profiles of clonal populations from 40 pancreatic adenocarcinomas and a set of prostate adenocarcinomas, including serial biopsies from a patient who progressed to androgen-independent metastatic disease. The genomes of clonal populations were found to have patient-specific aberrations of clinical relevance. Furthermore, we identified genomic aberrations specific to therapeutically responsive and resistant clones arising during the evolution of androgen-independent metastatic prostate adenocarcinoma. We also distinguished divergent clonal populations within single biopsies and mapped aberrations in multiple aneuploid populations arising in primary and metastatic pancreatic adenocarcinoma. We propose that our high-definition analyses of the genomes of distinct clonal populations of cancer cells in patients in vivo can help guide diagnoses and tailor approaches to personalized treatment.

It is widely accepted that many cancers arise as a result of an acquired genomic instability and the subsequent evolution of tumor cells with variable patterns of selected and background aberrations. The presence and behaviors of distinct neoplastic cell populations within a patient's tumor may underlie multiple clinical phenotypes in cancers. A goal of many current cancer genome studies is the identification of recurring selected driver events that can be advanced for the development of personalized therapies. Unfortunately, in the majority of rare tumors, this type of analysis can be particularly challenging. Large series of specimens for analysis are simply not available, allowing recurring patterns to remain hidden. In this paper, we highlight the use of DNA content-based flow sorting to identify and isolate DNA-diploid and DNA-aneuploid populations from tumor biopsies as a strategy to comprehensively study the genomic composition and behaviors of individual cancers in a series of rare solid tumors: intrahepatic cholangiocarcinoma, anal carcinoma, adrenal leiomyosarcoma, and pancreatic neuroendocrine tumors. We propose that the identification of highly selected genomic events in distinct tumor populations within each tumor can identify candidate driver events that can facilitate the development of novel, personalized treatment strategies for patients with cancer.