Genetic and immune landscape evolution in MMR-deficient colorectal cancer.

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.

The Mutational Concordance of Fixed Formalin Paraffin Embedded and Fresh Frozen Gastro-Oesophageal Tumours Using Whole Exome Sequencing.

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.

Correction: Suppression of interferon gene expression overcomes resistance to MEK inhibition in KRAS-mutant colorectal cancer.

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

Detecting and Tracking Circulating Tumour DNA Copy Number Profiles during First Line Chemotherapy in Oesophagogastric Adenocarcinoma.

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.

Suppression of interferon gene expression overcomes resistance to MEK inhibition in KRAS-mutant colorectal cancer.

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.

Capture Hi-C identifies putative target genes at 33 breast cancer risk loci.

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.

Identification of 19 new risk loci and potential regulatory mechanisms influencing susceptibility to testicular germ cell tumor.

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

Whole-exome DNA sequence analysis of Brca2- and Trp53-deficient mouse mammary gland tumours.

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.

Unbiased analysis of potential targets of breast cancer susceptibility loci by Capture Hi-C.

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.

Genome-wide profiling of genetic synthetic lethality identifies CDK12 as a novel determinant of PARP1/2 inhibitor sensitivity.

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.

The genomic landscape of oesophagogastric junctional adenocarcinoma.

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.

Efficacy of chemotherapy in BRCA1/2 mutation carrier ovarian cancer in the setting of PARP inhibitor resistance: a multi-institutional study.

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.

A genetic screen using the PiggyBac transposon in haploid cells identifies Parp1 as a mediator of olaparib toxicity.

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.

Secondary mutations in BRCA2 associated with clinical resistance to a PARP inhibitor.

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.

Genome-wide analysis of p63 binding sites identifies AP-2 factors as co-regulators of epidermal differentiation.

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.

A modified method for whole exome resequencing from minimal amounts of starting DNA.

Next generation DNA sequencing (NGS) technologies have revolutionized the pace at which whole genome and exome sequences can be generated. However, despite these advances, many of the methods for targeted resequencing, such as the generation of high-depth exome sequences, are somewhat limited by the relatively large amounts of starting DNA that are normally required. In the case of tumour analysis this is particularly pertinent as many tumour biopsies often return submicrogram quantities of DNA, especially when tumours are microdissected prior to analysis. Here, we present a method for exome capture and resequencing using as little as 50 ng of starting DNA. The sequencing libraries generated by this minimal starting amount (MSA-Cap) method generate datasets that are comparable to standard amount (SA) whole exome libraries that use three micrograms of starting DNA. This method, which can be performed in most laboratories using commonly available reagents, has the potential to enhance large scale profiling efforts such as the resequencing of tumour exomes.

Genome-wide functional screen identifies a compendium of genes affecting sensitivity to tamoxifen.

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.

High-throughput RNA interference screening using pooled shRNA libraries and next generation sequencing.

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.

Comprehensive genomic analysis of a BRCA2 deficient human pancreatic cancer.

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.