Mismatch repair-deficient (MMRd) colorectal cancers (CRCs) have high mutation burdens, which make these tumours immunogenic and many respond to immune checkpoint inhibitors. The MMRd hypermutator phenotype may also promote intratumour heterogeneity (ITH) and cancer evolution. We applied multiregion sequencing and CD8 and programmed death ligand 1 (PD-L1) immunostaining to systematically investigate ITH and how genetic and immune landscapes coevolve. All cases had high truncal mutation burdens. Despite pervasive ITH, driver aberrations showed a clear hierarchy. Those in WNT/β-catenin, mitogen-activated protein kinase, and TGF-β receptor family genes were almost always truncal. Immune evasion (IE) drivers, such as inactivation of genes involved in antigen presentation or IFN-γ signalling, were predominantly subclonal and showed parallel evolution. These IE drivers have been implicated in immune checkpoint inhibitor resistance or sensitivity. Clonality assessments are therefore important for the development of predictive immunotherapy biomarkers in MMRd CRCs. Phylogenetic analysis identified three distinct patterns of IE driver evolution: pan-tumour evolution, subclonal evolution, and evolutionary stasis. These, but neither mutation burdens nor heterogeneity metrics, significantly correlated with T-cell densities, which were used as a surrogate marker of tumour immunogenicity. Furthermore, this revealed that genetic and T-cell infiltrates coevolve in MMRd CRCs. Low T-cell densities in the subgroup without any known IE drivers may indicate an, as yet unknown, IE mechanism. PD-L1 was expressed in the tumour microenvironment in most samples and correlated with T-cell densities. However, PD-L1 expression in cancer cells was independent of T-cell densities but strongly associated with loss of the intestinal homeobox transcription factor CDX2. This explains infrequent PD-L1 expression by cancer cells and may contribute to a higher recurrence risk of MMRd CRCs with impaired CDX2 expression. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

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

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

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

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

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)

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 543K, Supplementary Table S2 shRNA constructs with Olaparib Drug Effect Z scores <-1.96

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.

Despite showing clinical activity in BRAF-mutant melanoma, the MEK inhibitor (MEKi) trametinib has failed to show clinical benefit in KRAS-mutant colorectal cancer. To identify mechanisms of resistance to MEKi, we employed a pharmacogenomic analysis of MEKi-sensitive versus MEKi-resistant colorectal cancer cell lines. Strikingly, interferon- and inflammatory-related gene sets were enriched in cell lines exhibiting intrinsic and acquired resistance to MEK inhibition. The bromodomain inhibitor JQ1 suppressed interferon-stimulated gene (ISG) expression and in combination with MEK inhibitors displayed synergistic effects and induced apoptosis in MEKi-resistant colorectal cancer cell lines. ISG expression was confirmed in patient-derived organoid models, which displayed resistance to trametinib and were resensitized by JQ1 co-treatment. In in vivo models of colorectal cancer, combination treatment significantly suppressed tumor growth. Our findings provide a novel explanation for the limited response to MEK inhibitors in KRAS-mutant colorectal cancer, known for its inflammatory nature. Moreover, the high expression of ISGs was associated with significantly reduced survival of colorectal cancer patients. Excitingly, we have identified novel therapeutic opportunities to overcome intrinsic and acquired resistance to MEK inhibition in colorectal cancer.

DNA somatic copy number aberrations (SCNAs) are key drivers in oesophagogastric adenocarcinoma (OGA). Whether minimally invasive SCNA analysis of circulating tumour (ct)DNA can predict treatment outcomes and reveal how SCNAs evolve during chemotherapy is unknown. We investigated this by low-coverage whole genome sequencing (lcWGS) of ctDNA from 30 patients with advanced OGA prior to first-line chemotherapy and on progression. SCNA profiles were detectable pretreatment in 23/30 (76.7%) patients. The presence of liver metastases, primary tumour in situ, or of oesophageal or junctional tumour location predicted for a high ctDNA fraction. A low ctDNA concentration associated with significantly longer overall survival. Neither chromosomal instability metrics nor ploidy correlated with chemotherapy outcome. Chromosome 2q and 8p gains before treatment were associated with chemotherapy responses. lcWGS identified all amplifications found by prior targeted tumour tissue sequencing in cases with detectable ctDNA as well as finding additional changes. SCNA profiles changed during chemotherapy, indicating that cancer cell populations evolved during treatment; however, no recurrent SCNA changes were acquired at progression. Tracking the evolution of OGA cancer cell populations in ctDNA is feasible during chemotherapy. The observation of genetic evolution warrants investigation in larger series and with higher resolution techniques to reveal potential genetic predictors of response and drivers of chemotherapy resistance. The presence of liver metastasis is a potential biomarker for the selection of patients with high ctDNA content for such studies.

Genome-wide association studies (GWAS) have identified approximately 100 breast cancer risk loci. Translating these findings into a greater understanding of the mechanisms that influence disease risk requires identification of the genes or non-coding RNAs that mediate these associations. Here, we use Capture Hi-C (CHi-C) to annotate 63 loci; we identify 110 putative target genes at 33 loci. To assess the support for these target genes in other data sources we test for associations between levels of expression and SNP genotype (eQTLs), disease-specific survival (DSS), and compare them with somatically mutated cancer genes. 22 putative target genes are eQTLs, 32 are associated with DSS and 14 are somatically mutated in breast, or other, cancers. Identifying the target genes at GWAS risk loci will lead to a greater understanding of the mechanisms that influence breast cancer risk and prognosis.

Intra-tumoral heterogeneity (ITH) is a major challenge for the detection of relevant genomic aberrations in primary tumor (TU) biopsies. Hence, there is an unmet need for alternative analytic procedures to obtain a more complete picture of the genomic landscape and to minimize sampling error in tumor patients. Analyzing cell-free tumor DNA (ctDNA) from more easily accessible patient biomaterial could surmount these challenges.

Clare Turnbull and colleagues report discovery of 19 new susceptibility loci for testicular germ cell tumor (TGCT) and provide evidence for a network of physical interactions between TGCT risk variants and candidate causal genes. Their findings implicate widespread disruption of developmental transcriptional regulators in TGCT susceptibility, consistent with failed primordial germ cell differentiation as an initiating step in oncogenesis. 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-analysis 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 oncogenesis1. 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.

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. Copyright © 2015 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

P<P Published online August 13, 2014 in advance of the print journal. Manuscript Accepted manuscript is likely to differ from the final, published version. Peer-reviewed and accepted for publication but not copyedited or typeset; accepted License Commons Creative . http://creativecommons.org/licenses/by-nc/4.0/ described at a Creative Commons License (Attribution-NonCommercial 4.0 International), as ). After six months, it is available under http://genome.cshlp.org/site/misc/terms.xhtml first six months after the full-issue publication date (see This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the

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.

AbstractSmall-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.

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.