Targeting DNA Damage Response and Replication Stress in Pancreatic Cancer.

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.

Genomic analyses identify molecular subtypes of pancreatic cancer.

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.

Clinical study of genomic drivers in pancreatic ductal adenocarcinoma.

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.

Glioblastomas are composed of genetically divergent clones with distinct tumourigenic potential and variable stem cell-associated phenotypes.

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.

Advancing a clinically relevant perspective of the clonal nature of cancer.

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.

Genomic analysis and selected molecular pathways in rare cancers.

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.

HNF4A and GATA6 Loss Reveals Therapeutically Actionable Subtypes in Pancreatic Cancer.

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.

Genomic and Phenotypic Characterization of a Broad Panel of Patient-Derived Xenografts Reflects the Diversity of Glioblastoma.

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.

Feasibility and clinical utility of endoscopic ultrasound guided biopsy of pancreatic cancer for next-generation molecular profiling.

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.

The Driver Mutational Landscape of Ovarian Squamous Cell Carcinomas Arising in Mature Cystic Teratoma.

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.

Mesenchymal-like pancreatic cancer cells harbor specific genomic alterations more frequently than their epithelial-like counterparts.

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.

STAG2 is a clinically relevant tumor suppressor in pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDA) is a highly lethal cancer characterized by complex aberrant genomes. A fundamental goal of current studies is to identify those somatic events arising in the variable landscape of PDA genomes that can be exploited for improved clinical outcomes. METHODS: We used DNA content flow sorting to identify and purify tumor nuclei of PDA samples from 50 patients. The genome of each sorted sample was profiled by oligonucleotide comparative genomic hybridization and targeted resequencing of STAG2. Transposon insertions within STAG2 in a KRAS (G12D)-driven genetically engineered mouse model of PDA were screened by RT-PCR. We then used a tissue microarray to survey STAG2 protein expression levels in 344 human PDA tumor samples and adjacent tissues. Univariate Kaplan Meier analysis and multivariate Cox Regression analysis were used to assess the association of STAG2 expression relative to overall survival and response to adjuvant therapy. Finally, RNAi-based assays with PDA cell lines were used to assess the potential therapeutic consequence of STAG2 expression in response to 18 therapeutic agents. RESULTS: STAG2 is targeted by somatic aberrations in a subset (4%) of human PDAs. Transposon-mediated disruption of STAG2 in a KRAS (G12D) genetically engineered mouse model promotes the development of PDA and its progression to metastatic disease. There was a statistically significant loss of STAG2 protein expression in human tumor tissue (Wilcoxon-Rank test) with complete absence of STAG2 staining observed in 15 (4.3%) patients. In univariate Kaplan Meier analysis nearly complete STAG2 positive staining (>95% of nuclei positive) was associated with a median survival benefit of 6.41 months (P = 0.031). The survival benefit of adjuvant chemotherapy was only seen in patients with a STAG2 staining of less than 95% (median survival benefit 7.65 months; P = 0.028). Multivariate Cox Regression analysis showed that STAG2 is an independent prognostic factor for survival in pancreatic cancer patients. Finally, we show that RNAi-mediated knockdown of STAG2 selectively sensitizes human PDA cell lines to platinum-based therapy. CONCLUSIONS: Based on these iterative findings we propose that STAG2 is a clinically significant tumor suppressor in PDA.

Whole genome analyses of a well-differentiated liposarcoma reveals novel SYT1 and DDR2 rearrangements.

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.

Elucidating potentially significant genomic regions involved in the initiation and progression of undifferentiated pleomorphic sarcoma.

Sarcomas are cancers that arise in soft tissues or bone and make up a small percentage of malignancies. In an effort to identify potential genetic targets for therapy, this study explores the genomic landscape of a metastatic undifferentiated pleomorphic sarcoma (UPS) with spindle cell morphology. Thick sections (50 µm) of formalin-fixed, paraffin-embedded tissue from a primary, recurrent, and metastatic tumor were collected and processed from a single patient for DNA content-based flow-sorting and analyses. Nuclei of diploid and aneuploid populations were sorted from the malignant tissues and their genomes interrogated with array comparative genomic hybridization. The third sample was highly degraded and did not contain any intact ploidy peaks in our flow assays. A 2.5N aneuploid population was identified in the primary and recurrent sample. We detected a series of shared and unique genomic aberrations in the sorted aneuploid populations. The patterns of aberrations suggest that two similar but independent clonal populations arose during the clinical history of this rare tumor. None of these aberrations were detected in the matching sorted diploid samples. The targeted regions of interest might play a role in UPS and may lead to clinical significance with further investigation.

Clonal evolution and therapeutic resistance in solid tumors.

Tumors frequently arise as a result of an acquired genomic instability and the subsequent evolution of neoplastic populations with variable genomes. A barrier to the study of the somatic genetics of human solid tumors in vivo is the presence of admixtures of non-neoplastic cells with normal genomes in patient samples. These can obscure the presence of somatic aberrations including mutations, homozygous deletions, and breakpoints in biopsies of interest. Furthermore, clinical samples frequently contain multiple neoplastic populations that cannot be distinguished by morphology. Consequently, it is difficult to determine whether mutations detected in a sample of interest are concurrent in a single clonal population or if they occur in distinct cell populations in the same sample. The advent of targeted therapies increases the selection for preexisting populations. However the asymmetric distribution of therapeutic targets in clonal populations provides a mechanism for the rapid evolution of resistant disease. Thus, there is a need to not only isolate tumor from normal cells, but to also enrich distinct populations of clonal neoplastic cells in order to apply genome technologies to identify clinically relevant genomic aberrations that drive disease in patients in vivo. To address this we have applied single and multiparameter DNA content based flow assays to the study of solid tumors. Our work has identified examples of clonal resistance to effective therapies. This includes androgen withdrawal in advanced prostate cancer. In addition we demonstrate examples of co-existing clonal populations with highly aberrant genomes and ploidies in a wide variety of solid tumors. We propose that clonal analysis of tumors, based on flow cytometry and high resolution genome analyses of purified neoplastic populations, provides a unique approach to the study of therapeutic responses and the evolution of resistance.

Deep clonal profiling of formalin fixed paraffin embedded clinical samples.

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.