Search

Supplementary Table 1 - PDF file 57K, Supplementary Table 1. Details of platinum sensitivity and best response to olaparib and post-PARPi chemotherapy of the six patients whose tumor samples were analyzed by massively parallel sequencing. (NE: not evaluable)

Purpose: Preclinical data suggest that exposure to PARP inhibitors (PARPi) may compromise benefit to subsequent chemotherapy, particularly platinum-based regimens, in patients with BRCA1/2 mutation carrier ovarian cancer (PBMCOC), possibly through the acquisition of secondary BRCA1/2 mutations. The efficacy of chemotherapy in the PARPi-resistant setting was therefore investigated.Experimental Design: We conducted a retrospective review of PBMCOC who received chemotherapy following disease progression on olaparib, administered at ≥200 mg twice daily for one month or more. Tumor samples were obtained in the post-olaparib setting where feasible and analyzed by massively parallel sequencing.Results: Data were collected from 89 patients who received a median of 3 (range 1–11) lines of pre-olaparib chemotherapy. The overall objective response rate (ORR) to post-olaparib chemotherapy was 36% (24 of 67 patients) by Response Evaluation Criteria in Solid Tumors (RECIST) and 45% (35 of 78) by RECIST and/or Gynecologic Cancer InterGroup (GCIG) CA125 criteria with median progression-free survival (PFS) and overall survival (OS) of 17 weeks [95% confidence interval (CI), 13–21] and 34 weeks (95% CI, 26–42), respectively. For patients receiving platinum-based chemotherapy, ORRs were 40% (19 of 48) and 49% (26/53), respectively, with a median PFS of 22 weeks (95% CI, 15–29) and OS of 45 weeks (95% CI, 15–75). An increased platinum-to-platinum interval was associated with an increased OS and likelihood of response following post-olaparib platinum. No evidence of secondary BRCA1/2 mutation was detected in tumor samples of six PARPi-resistant patients [estimated frequency of such mutations adjusted for sample size: 0.125 (95%-CI: 0–0.375)].Conclusions: Heavily pretreated PBMCOC who are PARPi-resistant retain the potential to respond to subsequent chemotherapy, including platinum-based agents. These data support the further development of PARPi in PBMCOC. Clin Cancer Res; 19(19); 5485–93. ©2013 AACR.

PDF file 543K, Supplementary Table S2 shRNA constructs with Olaparib Drug Effect Z scores

PDF file 32K, Supplementary Table S1 Antibody reagent list used in study

PDF file 58K, Supplementary Table S4. DNA repair involvement of PARP inhibitor sensitisation genes found in the screen. PMID numbers refer to pubmed entries

Small-molecule inhibitors of PARP1/2, such as olaparib, have been proposed to serve as a synthetic lethal therapy for cancers that harbor BRCA1 or BRCA2 mutations. Indeed, in clinical trials, PARP1/2 inhibitors elicit sustained antitumor responses in patients with germline BRCA gene mutations. In hypothesizing that additional genetic determinants might direct use of these drugs, we conducted a genome-wide synthetic lethal screen for candidate olaparib sensitivity genes. In support of this hypothesis, the set of identified genes included known determinants of olaparib sensitivity, such as BRCA1, RAD51, and Fanconi's anemia susceptibility genes. In addition, the set included genes implicated in established networks of DNA repair, DNA cohesion, and chromatin remodeling, none of which were known previously to confer sensitivity to PARP1/2 inhibition. Notably, integration of the list of candidate sensitivity genes with data from tumor DNA sequencing studies identified CDK12 deficiency as a clinically relevant biomarker of PARP1/2 inhibitor sensitivity. In models of high-grade serous ovarian cancer (HGS-OVCa), CDK12 attenuation was sufficient to confer sensitivity to PARP1/2 inhibition, suppression of DNA repair via homologous recombination, and reduced expression of BRCA1. As one of only nine genes known to be significantly mutated in HGS-OVCa, CDK12 has properties that should confirm interest in its use as a biomarker, particularly in ongoing clinical trials of PARP1/2 inhibitors and other agents that trigger replication fork arrest. Cancer Res; 74(1); 287–97. ©2013 AACR.

PDF file 470K, Supplementary Table S3 shRNA constructs with Olaparib Drug Effect Z scores >1.96

PDF file 3292K, Supplementary Figure S1. Genome wide shRNA screen information. Supplementary Figure S2. Validation of ATAD5 as a genetic determinant of olaparib response. Supplementary Figure S3. A working model of PARP1/2 inhibitor-induced DNA repair. Supplementary Figure S4. Inhibition of CDK12 sensitises serous ovarian cancer cells to olaparib and cisplatin. Supplementary Figure S5. CDK12 silencing and the impact on expression of DNA repair proteins. Supplementary Figure S6. Targeting of CCNK sensitises serous ovarian cancer cells to olaparib and cisplatin. Supplementary Figure S7. Animal body weights from in vivo study

PDF file 54K, Supplementary Table S5. Tumor-specific CDK12 mutations described in publically available sequencing data

PDF file 120K, Supplementary figure legends and supplementary materials and methods

Mismatch repair deficient colorectal cancers have high mutation loads and many respond to immune checkpoint-inhibitors. We investigated how genetic and immune landscapes co-evolve in these tumors. All cases had high truncal mutation loads. Driver aberrations showed a clear hierarchy despite pervasive intratumor heterogeneity: Those in WNT/βCatenin, mitogen-activated protein kinase and TGFβ receptor family genes were almost always truncal. Immune evasion drivers were predominantly subclonal and showed parallel evolution. Pan-tumor evolution, subclonal evolution, and evolutionary stasis of genetic immune evasion drivers defined three MMRd CRC subtypes with distinct T-cell infiltrates. These immune evasion drivers have been implicated in checkpoint-inhibitor resistance. Clonality and subtype assessments are hence critical for predictive immunotherapy biomarker development. Cancer cell PD-L1 expression was conditional on loss of the intestinal homeobox transcription factor CDX2. This explains infrequent PD-L1 expression by cancer cells and likely contributes to the high recurrence risk of MMRd CRCs with impaired CDX2 expression.

1. Background: The application of massively parallel sequencing has led to the identification of aberrant druggable pathways and somatic mutations within therapeutically relevant genes in gastro-oesophageal cancer. Given the widespread use of formalin-fixed paraffin-embedded (FFPE) samples in the study of this disease, it would be beneficial, especially for the purposes of biomarker evaluation, to assess the concordance between comprehensive exome-wide sequencing data from archival FFPE samples originating from a prospective clinical study and those derived from fresh-frozen material. 2. Methods: We analysed whole-exome sequencing data to define the mutational concordance of 16 matched fresh-frozen and FFPE gastro-oesophageal tumours (N = 32) from a prospective clinical study. We assessed DNA integrity prior to sequencing and then identified coding mutations in genes that have previously been implicated in other cancers. In addition, we calculated the mutant-allele heterogeneity (MATH) for these samples. 3. Results: Although there was increased degradation of DNA in FFPE samples compared with frozen samples, sequencing data from only two FFPE samples failed to reach an adequate mapping quality threshold. Using a filtering threshold of mutant read counts of at least ten and a minimum of 5% variant allele frequency (VAF) we found that there was a high median mutational concordance of 97% (range 80.1–98.68%) between fresh-frozen and FFPE gastro-oesophageal tumour-derived exomes. However, the majority of FFPE tumours had higher mutant-allele heterogeneity (MATH) scores when compared with corresponding frozen tumours (p < 0.001), suggesting that FFPE-based exome sequencing is likely to over-represent tumour heterogeneity in FFPE samples compared to fresh-frozen samples. Furthermore, we identified coding mutations in 120 cancer-related genes, including those associated with chromatin remodelling and Wnt/β-catenin and Receptor Tyrosine Kinase signalling. 4. Conclusions: These data suggest that comprehensive genomic data can be generated from exome sequencing of selected DNA samples extracted from archival FFPE gastro-oesophageal tumour tissues within the context of prospective clinical trials.

Genome-wide association studies (GWAS) have transformed understanding of susceptibility to testicular germ cell tumors (TGCTs), but much of the heritability remains unexplained. Here we report a new GWAS, a meta- A nalysis with previous GWAS and a replication series, totaling 7,319 TGCT cases and 23,082 controls. We identify 19 new TGCT risk loci, roughly doubling the number of known TGCT risk loci to 44. By performing in situ Hi-C in TGCT cells, we provide evidence for a network of physical interactions among all 44 TGCT risk SNPs and candidate causal genes. Our findings implicate widespread disruption of developmental transcriptional regulators as a basis of TGCT susceptibility, consistent with failed primordial germ cell differentiation as an initiating step in oncogenesis. Defective microtubule assembly and dysregulation of KIT-MAPK signaling also feature as recurrently disrupted pathways. Our findings support a polygenic model of risk and provide insight into the biological basis of TGCT.

Small-molecule inhibitors of PARP1/2, such as olaparib, have been proposed to serve as a synthetic lethal therapy for cancers that harbor BRCA1 or BRCA2 mutations. Indeed, in clinical trials, PARP1/2 inhibitors elicit sustained antitumor responses in patients with germline BRCA gene mutations. In hypothesizing that additional genetic determinants might direct use of these drugs, we conducted a genome-wide synthetic lethal screen for candidate olaparib sensitivity genes. In support of this hypothesis, the set of identified genes included known determinants of olaparib sensitivity, such as BRCA1, RAD51, and Fanconi's anemia susceptibility genes. In addition, the set included genes implicated in established networks of DNA repair, DNA cohesion, and chromatin remodeling, none of which were known previously to confer sensitivity to PARP1/2 inhibition. Notably, integration of the list of candidate sensitivity genes with data from tumor DNA sequencing studies identified CDK12 deficiency as a clinically relevant biomarker of PARP1/2 inhibitor sensitivity. In models of high-grade serous ovarian cancer (HGS-OVCa), CDK12 attenuation was sufficient to confer sensitivity to PARP1/2 inhibition, suppression of DNA repair via homologous recombination, and reduced expression of BRCA1. As one of only nine genes known to be significantly mutated in HGS-OVCa, CDK12 has properties that should confirm interest in its use as a biomarker, particularly in ongoing clinical trials of PARP1/2 inhibitors and other agents that trigger replication fork arrest. Cancer Res; 74(1); 287–97. ©2013 AACR.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

The incidence of oesophagogastric junctional (OGJ) adenocarcinoma is rising rapidly in western countries, in contrast to the declining frequency of distal gastric carcinoma. Treatment options for adenocarcinomas involving the oesophagogastric junction are limited and the overall prognosis is extremely poor. To determine the genomic landscape of OGJ adenocarcinoma, exomes of eight tumours and matched germline DNA were subjected to massively parallel DNA sequencing. Microsatellite instability was observed in three tumours which coincided with an elevated number of somatic mutations. In total, 117 genes were identified that had predicted coding alterations in more than one tumour. Potentially actionable coding mutations were identified in 67 of these genes, including those in CR2, HGF , FGFR4, and ESRRB. Twenty-nine genes harbouring somatic coding mutations and copy number changes in the MSS OGJ dataset are also known to be altered with similar predicted functional consequence in other tumour types. Compared with the published mutational profile of gastric cancers, 49% (57/117) of recurrently mutated genes were unique to OGJ tumours. TP53, SYNE1, and ARID1A were amongst the most frequently mutated genes in a larger OGJ cohort. Our study provides an insight into the mutational landscape of OGJ adenocarcinomas and confirms that this is a highly mutated and heterogeneous disease. Furthermore, we have uncovered somatic mutations in therapeutically relevant genes which may represent candidate drug targets.

Purpose: Preclinical data suggest that exposure to PARP inhibitors (PARPi) may compromise benefit to subsequent chemotherapy, particularly platinum-based regimens, in patients with BRCA1/2 mutation carrier ovarian cancer (PBMCOC), possibly through the acquisition of secondary BRCA1/2 mutations. The efficacy of chemotherapy in the PARPi-resistant setting was therefore investigated. Experimental Design: We conducted a retrospective review of PBMCOC who received chemotherapy following disease progression on olaparib, administered at ≥200 mg twice daily for one month or more. Tumor samples were obtained in the post-olaparib setting where feasible and analyzed by massively parallel sequencing. Results: Data were collected from 89 patients who received a median of 3 (range 1–11) lines of pre-olaparib chemotherapy. The overall objective response rate (ORR) to post-olaparib chemotherapy was 36% (24 of 67 patients) by Response Evaluation Criteria in Solid Tumors (RECIST) and 45% (35 of 78) by RECIST and/or Gynecologic Cancer InterGroup (GCIG) CA125 criteria with median progression-free survival (PFS) and overall survival (OS) of 17 weeks [95% confidence interval (CI), 13–21] and 34 weeks (95% CI, 26–42), respectively. For patients receiving platinum-based chemotherapy, ORRs were 40% (19 of 48) and 49% (26/53), respectively, with a median PFS of 22 weeks (95% CI, 15–29) and OS of 45 weeks (95% CI, 15–75). An increased platinum-to-platinum interval was associated with an increased OS and likelihood of response following post-olaparib platinum. No evidence of secondary BRCA1/2 mutation was detected in tumor samples of six PARPi-resistant patients [estimated frequency of such mutations adjusted for sample size: 0.125 (95%-CI: 0–0.375)]. Conclusions: Heavily pretreated PBMCOC who are PARPi-resistant retain the potential to respond to subsequent chemotherapy, including platinum-based agents. These data support the further development of PARPi in PBMCOC. Clin Cancer Res; 19(19); 5485–93. ©2013 AACR.

PARP inhibitors (PARPi) for the treatment of BRCA1 or BRCA2 deficient tumours are currently the focus of seminal clinical trials exploiting the concept of synthetic lethality. Although clinical resistance to PARPi has been described, the mechanism underlying this has not been elucidated. Here, we investigate tumour material from patients who had developed resistance to the PARPi olaparib, subsequent to showing an initial clinical response. Massively parallel DNA sequencing of treatment-naive and post-olaparib treatment biopsies identified tumour-specific BRCA2 secondary mutations in olaparib-resistant metastases. These secondary mutations restored full-length BRCA2 protein, and most likely cause olaparib resistance by re-establishing BRCA2 function in the tumour cells.

The p63 transcription factor (TP63) is critical in development, growth and differentiation of stratifying epithelia. This is highlighted by the severity of congenital abnormalities caused by TP63 mutations in humans, the dramatic phenotypes in knockout mice and de-regulation of TP63 expression in neoplasia altering the tumour suppressive roles of the TP53 family. In order to define the normal role played by TP63 and provide the basis for better understanding how this network is perturbed in disease, we used chromatin immunoprecipitation combined with massively parallel sequencing (ChIP-seq) to identify >7500 high-confidence TP63-binding regions across the entire genome, in primary human neonatal foreskin keratinocytes (HFKs). Using integrative strategies, we demonstrate that only a subset of these sites are bound by TP53 in response to DNA damage. We identify a role for TP63 in transcriptional regulation of multiple genes genetically linked to cleft palate and identify AP-2alpha (TFAP2A) as a co-regulator of a subset of these genes. We further demonstrate that AP-2gamma (TFAP2C) can bind a subset of these regions and that acute depletion of either TFAP2A or TFAP2C alone is sufficient to reduce terminal differentiation of organotypic epidermal skin equivalents, indicating overlapping physiological functions with TP63.

Germline mutations in the tumour suppressor BRCA2 predispose to breast, ovarian and a number of other human cancers. Brca2-deficient mouse models are used for preclinical studies but the pattern of genomic alterations in these tumours has not yet been described in detail. We have performed whole-exome DNA sequencing analysis of mouse mammary tumours from Blg-Cre Brca2(f/f) Trp53(f/f) animals, a model of BRCA2-deficient human cancer. We also used the sequencing data to estimate DNA copy number alterations in these tumours and identified a recurrent copy number gain in Met, which has been found amplified in other mouse mammary cancer models. Through a comparative genomic analysis, we identified several mouse Blg-Cre Brca2(f/f) Trp53(f/f) mammary tumour somatic mutations in genes that are also mutated in human cancer, but few of these genes have been found frequently mutated in human breast cancer. A more detailed analysis of these somatic mutations revealed a set of genes that are mutated in human BRCA2 mutant breast and ovarian tumours and that are also mutated in mouse Brca2-null, Trp53-null mammary tumours. Finally, a DNA deletion surrounded by microhomology signature found in human BRCA1/2-deficient cancers was not common in the genome of these mouse tumours. Although a useful model, there are some differences in the genomic landscape of tumours arising in Blg-Cre Brca2(f/f) Trp53(f/f) mice compared to human BRCA-mutated breast cancers. Therefore, this needs to be taken into account in the use of this model.

Genome-wide association studies have identified more than 70 common variants that are associated with breast cancer risk. Most of these variants map to non-protein-coding regions and several map to gene deserts, regions of several hundred kilobases lacking protein-coding genes. We hypothesized that gene deserts harbor long-range regulatory elements that can physically interact with target genes to influence their expression. To test this, we developed Capture Hi-C (CHi-C), which, by incorporating a sequence capture step into a Hi-C protocol, allows high-resolution analysis of targeted regions of the genome. We used CHi-C to investigate long-range interactions at three breast cancer gene deserts mapping to 2q35, 8q24.21, and 9q31.2. We identified interaction peaks between putative regulatory elements (“bait fragments”) within the captured regions and “targets” that included both protein-coding genes and long noncoding (lnc) RNAs over distances of 6.6 kb to 2.6 Mb. Target protein-coding genes were IGFBP5, KLF4, NSMCE2, and MYC; and target lncRNAs included DIRC3, PVT1, and CCDC26. For one gene desert, we were able to define two SNPs (rs12613955 and rs4442975) that were highly correlated with the published risk variant and that mapped within the bait end of an interaction peak. In vivo ChIP-qPCR data show that one of these, rs4442975, affects the binding of FOXA1 and implicate this SNP as a putative functional variant.

The incidence of oesophagogastric junctional (OGJ) adenocarcinoma is rising rapidly in western countries, in contrast to the declining frequency of distal gastric carcinoma. Treatment options for adenocarcinomas involving the oesophagogastric junction are limited and the overall prognosis is extremely poor. To determine the genomic landscape of OGJ adenocarcinoma, exomes of eight tumours and matched germline DNA were subjected to massively parallel DNA sequencing. Microsatellite instability was observed in three tumours which coincided with an elevated number of somatic mutations. In total, 117 genes were identified that had predicted coding alterations in more than one tumour. Potentially actionable coding mutations were identified in 67 of these genes, including those in CR2, HGF , FGFR4, and ESRRB. Twenty-nine genes harbouring somatic coding mutations and copy number changes in the MSS OGJ dataset are also known to be altered with similar predicted functional consequence in other tumour types. Compared with the published mutational profile of gastric cancers, 49% (57/117) of recurrently mutated genes were unique to OGJ tumours. TP53, SYNE1, and ARID1A were amongst the most frequently mutated genes in a larger OGJ cohort. Our study provides an insight into the mutational landscape of OGJ adenocarcinomas and confirms that this is a highly mutated and heterogeneous disease. Furthermore, we have uncovered somatic mutations in therapeutically relevant genes which may represent candidate drug targets.

Background: Prostate cancer (PrCa) is one of the most common diseases to affect men worldwide and among the leading causes of cancer-related death. The purpose of this study was to use second-generation sequencing technology to assess the frequency of deleterious mutations in 22 tumour suppressor genes in familial PrCa and estimate the relative risk of PrCa if these genes are mutated. Methods: Germline DNA samples from 191 men with 3 or more cases of PrCa in their family were sequenced for 22 tumour suppressor genes using Agilent target enrichment and Illumina technology. Analysis for genetic variation was carried out by using a pipeline consisting of BWA, Genome Analysis Toolkit (GATK) and ANNOVAR. Clinical features were correlated with mutation status using standard statistical tests. Modified segregation analysis was used to determine the relative risk of PrCa conferred by the putative loss-of-function (LoF) mutations identified. Results: We discovered 14 putative LoF mutations in 191 samples (7.3%) and these mutations were more frequently associated with nodal involvement, metastasis or T4 tumour stage (P=0.00164). Segregation analysis of probands with European ancestry estimated that LoF mutations in any of the studied genes confer a relative risk of PrCa of 1.94 (95% CI: 1.56–2.42). Conclusions: These findings show that LoF mutations in DNA repair pathway genes predispose to familial PrCa and advanced disease and therefore warrants further investigation. The clinical utility of these findings will become increasingly important as targeted screening and therapies become more widespread.

Genetic perturbation screens have the potential to dissect a wide range of cellular phenotypes. Such screens have historically been difficult in diploid mammalian cells. The recent derivation of haploid embryonic stem cells provides an opportunity to cause loss of function mutants with a random mutagen in a mammalian cell with a normal genetic background. We describe an approach to genetic screens that exploits the highly active piggyBac transposon in haploid mammalian cells. As an example of haploid transposon (HTP) screening, we apply this approach to identifying determinants of cancer drug toxicity and resistance. In a screen for 6-thioguanine resistance we recovered components of the DNA mismatch repair pathway, a known requirement for toxicity. In a further screen for resistance to the clinical poly(ADP-ribose) polymerase (PARP) inhibitor olaparib we recovered multiple Parp1 mutants. Our results show that olaparib toxicity to normal cells is mediated predominantly via Parp1, and suggest that the clinical side effects of olaparib may be on target. The transposon mutant libraries are stable and can be readily reused to screen other drugs. The screening protocol described has several advantages over other methods such as RNA interference: it is rapid and low cost, and mutations can be easily reverted to establish causality.

PARP inhibitors (PARPi) for the treatment of BRCA1 or BRCA2 deficient tumours are currently the focus of seminal clinical trials exploiting the concept of synthetic lethality. Although clinical resistance to PARPi has been described, the mechanism underlying this has not been elucidated. Here, we investigate tumour material from patients who had developed resistance to the PARPi olaparib, subsequent to showing an initial clinical response. Massively parallel DNA sequencing of treatment-naive and post-olaparib treatment biopsies identified tumour-specific BRCA2 secondary mutations in olaparib-resistant metastases. These secondary mutations restored full-length BRCA2 protein, and most likely cause olaparib resistance by re-establishing BRCA2 function in the tumour cells.

Therapies that target estrogen signaling have made a very considerable contribution to reducing mortality from breast cancer. However, resistance to tamoxifen remains a major clinical problem. Here we have used a genome-wide functional profiling approach to identify multiple genes that confer resistance or sensitivity to tamoxifen. Combining whole-genome shRNA screening with massively parallel sequencing, we have profiled the impact of more than 56,670 RNA interference reagents targeting 16,487 genes on the cellular response to tamoxifen. This screen, along with subsequent validation experiments, identifies a compendium of genes whose silencing causes tamoxifen resistance (including BAP1 , CLPP , GPRC5D , NAE1 , NF1 , NIPBL , NSD1 , RAD21 , RARG , SMC3 , and UBA3 ) and also a set of genes whose silencing causes sensitivity to this endocrine agent ( C10orf72 , C15orf55/NUT , EDF1 , ING5 , KRAS , NOC3L , PPP1R15B , RRAS2 , TMPRSS2 , and TPM4 ). Multiple individual genes, including NF1 , a regulator of RAS signaling, also correlate with clinical outcome after tamoxifen treatment.

Capan-1 is a well-characterised BRCA2-deficient human cell line isolated from a liver metastasis of a pancreatic adenocarcinoma. Here we report a genome-wide assessment of structural variations and high-depth exome characterization of single nucleotide variants and small insertion/deletions in Capan-1. To identify potential somatic and tumour-associated variations in the absence of a matched-normal cell line, we devised a novel method based on the analysis of HapMap samples. We demonstrate that Capan-1 has one of the most rearranged genomes sequenced to date. Furthermore, small insertions and deletions are detected more frequently in the context of short sequence repeats than in other genomes. We also identify a number of novel mutations that may represent genetic changes that have contributed to tumour progression. These data provide insight into the genomic effects of loss of BRCA2 function.

RNA interference (RNAi) screening is a state-of-the-art technology that enables the dissection of biological processes and disease-related phenotypes. The commercial availability of genome-wide, short hairpin RNA (shRNA) libraries has fueled interest in this area but the generation and analysis of these complex data remain a challenge. Here, we describe complete experimental protocols and novel open source computational methodologies, shALIGN and shRNAseq, that allow RNAi screens to be rapidly deconvoluted using next generation sequencing. Our computational pipeline offers efficient screen analysis and the flexibility and scalability to quickly incorporate future developments in shRNA library technology.