A widespread approach to modern cancer therapy is to identify a single oncogenic driver gene and target its mutant-protein product (for example, EGFR-inhibitor treatment in EGFR-mutant lung cancers). ...However, genetically driven resistance to targeted therapy limits patient survival. Through genomic analysis of 1,122 EGFR-mutant lung cancer cell-free DNA samples and whole-exome analysis of seven longitudinally collected tumor samples from a patient with EGFR-mutant lung cancer, we identified critical co-occurring oncogenic events present in most advanced-stage EGFR-mutant lung cancers. We defined new pathways limiting EGFR-inhibitor response, including WNT/β-catenin alterations and cell-cycle-gene (CDK4 and CDK6) mutations. Tumor genomic complexity increases with EGFR-inhibitor treatment, and co-occurring alterations in CTNNB1 and PIK3CA exhibit nonredundant functions that cooperatively promote tumor metastasis or limit EGFR-inhibitor response. This study calls for revisiting the prevailing single-gene driver-oncogene view and links clinical outcomes to co-occurring genetic alterations in patients with advanced-stage EGFR-mutant lung cancer.
Colorectal cancer (CRC) shows variable response to immune checkpoint blockade, which can only partially be explained by high tumor mutational burden (TMB). We conducted an integrated study of the ...cancer tissue and associated tumor microenvironment (TME) from patients treated with pembrolizumab (KEYNOTE 177 clinical trial) or nivolumab to dissect the cellular and molecular determinants of response to anti- programmed cell death 1 (PD1) immunotherapy.
We selected multiple regions per tumor showing variable T-cell infiltration for a total of 738 regions from 29 patients, divided into discovery and validation cohorts. We performed multiregional whole-exome and RNA sequencing of the tumor cells and integrated these with T-cell receptor sequencing, high-dimensional imaging mass cytometry, detection of programmed death-ligand 1 (PDL1) interaction in situ, multiplexed immunofluorescence, and computational spatial analysis of the TME.
In hypermutated CRCs, response to anti-PD1 immunotherapy was not associated with TMB but with high clonality of immunogenic mutations, clonally expanded T cells, low activation of Wnt signaling, deregulation of the interferon gamma pathway, and active immune escape mechanisms. Responsive hypermutated CRCs were also rich in cytotoxic and proliferating PD1+CD8 T cells interacting with PDL1+ antigen-presenting macrophages.
Our study clarified the limits of TMB as a predictor of response of CRC to anti-PD1 immunotherapy. It identified a population of antigen-presenting macrophages interacting with CD8 T cells that consistently segregate with response. We therefore concluded that anti-PD1 agents release the PD1-PDL1 interaction between CD8 T cells and macrophages to promote cytotoxic antitumor activity.
Display omitted
Colorectal cancers responsive to anti-programmed cell death 1 immunotherapy show clonal immunogenic mutations, low Wnt activation, beta-2-microglobulin deregulation, and high infiltration of antigen presenting macrophages interacting with programmed cell death 1-positive cluster of differentiation 8 T cells.
As tumors grow, they acquire mutations, some of which create neoantigens that influence the response of patients to immune checkpoint inhibitors. We explored the impact of neoantigen intratumor ...heterogeneity (ITH) on antitumor immunity. Through integrated analysis of ITH and neoantigen burden, we demonstrate a relationship between clonal neoantigen burden and overall survival in primary lung adenocarcinomas. CD8⁺ tumor-infiltrating lymphocytes reactive to clonal neoantigens were identified in early-stage non–small cell lung cancer and expressed high levels of PD-1. Sensitivity to PD-1 and CTLA-4 blockade in patients with advanced NSCLC and melanoma was enhanced in tumors enriched for clonal neoantigens. T cells recognizing clonal neoantigens were detectable in patients with durable clinical benefit. Cytotoxic chemotherapy–induced subclonal neoantigens, contributing to an increased mutational load, were enriched in certain poor responders. These data suggest that neoantigen heterogeneity may influence immune surveillance and support therapeutic developments targeting clonal neoantigens.
DNA methylation is an epigenetic mark that has a crucial role in many biological processes. To understand the functional consequences of DNA methylation on phenotypic plasticity, a genome-wide ...analysis should be embraced. This in turn requires a technique that balances accuracy, genome coverage, resolution and cost, yet is low in DNA input in order to minimize the drain on precious samples. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) fulfils these criteria, combining MeDIP with massively parallel DNA sequencing. Here we report an improved protocol using 100-fold less genomic DNA than that commonly used. We show comparable results for specificity (>97%) and enrichment (>100-fold) over a wide range of DNA concentrations (5,000-50 ng) and demonstrate the utility of the protocol for the generation of methylomes from rare bone marrow cells using 160-300 ng of starting DNA. The protocol described here, i.e., DNA extraction to generation of MeDIP-seq library, can be completed within 3-5 d.
Esophageal adenocarcinomas are associated with a dismal prognosis. Deciphering the evolutionary history of this disease may shed light on therapeutically tractable targets and reveal dynamic ...mutational processes during the disease course and following neoadjuvant chemotherapy (NAC). We exome sequenced 40 tumor regions from 8 patients with operable esophageal adenocarcinomas, before and after platinum-containing NAC. This revealed the evolutionary genomic landscape of esophageal adenocarcinomas with the presence of heterogeneous driver mutations, parallel evolution, early genome-doubling events, and an association between high intratumor heterogeneity and poor response to NAC. Multiregion sequencing demonstrated a significant reduction in thymine to guanine mutations within a CpTpT context when comparing early and late mutational processes and the presence of a platinum signature with enrichment of cytosine to adenine mutations within a CpC context following NAC. Esophageal adenocarcinomas are characterized by early chromosomal instability leading to amplifications containing targetable oncogenes persisting through chemotherapy, providing a rationale for future therapeutic approaches.
This work illustrates dynamic mutational processes occurring during esophageal adenocarcinoma evolution and following selective pressures of platinum exposure, emphasizing the iatrogenic impact of therapy on cancer evolution. Identification of amplifications encoding targetable oncogenes maintained through NAC suggests the presence of stable vulnerabilities, unimpeded by cytotoxics, suitable for therapeutic intervention.
Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined ...genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10(-4), permutation p = 1.0×10(-3)). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10(-7)). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within disease-associated loci, thus providing a novel route to aid unravelling common complex diseases.
Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes
. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation ...necessary for tumour evolution
. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2
breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.
The integration of genomic and epigenomic data is an increasingly popular approach for studying the complex mechanisms driving cancer development. We have developed a method for evaluating both ...methylation and copy number from high-density DNA methylation arrays. Comparing copy number data from Infinium HumanMethylation450 BeadChips and SNP arrays, we demonstrate that Infinium arrays detect copy number alterations with the sensitivity of SNP platforms. These results show that high-density methylation arrays provide a robust and economic platform for detecting copy number and methylation changes in a single experiment. Our method is available in the ChAMP Bioconductor package: http://www.bioconductor.org/packages/2.13/bioc/html/ChAMP.html.
Malignant Pleural Mesothelioma (MPM) is typically diagnosed 20-50 years after exposure to asbestos and evolves along an unknown evolutionary trajectory. To elucidate this path, we conducted ...multi-regional exome sequencing of 90 tumour samples from 22 MPMs acquired at surgery. Here we show that exomic intratumour heterogeneity varies widely across the cohort. Phylogenetic tree topology ranges from linear to highly branched, reflecting a steep gradient of genomic instability. Using transfer learning, we detect repeated evolution, resolving 5 clusters that are prognostic, with temporally ordered clonal drivers. BAP1/-3p21 and FBXW7/-chr4 events are always early clonal. In contrast, NF2/-22q events, leading to Hippo pathway inactivation are predominantly late clonal, positively selected, and when subclonal, exhibit parallel evolution indicating an evolutionary constraint. Very late somatic alteration of NF2/22q occurred in one patient 12 years after surgery. Clonal architecture and evolutionary clusters dictate MPM inflammation and immune evasion. These results reveal potentially drugable evolutionary bottlenecking in MPM, and an impact of clonal architecture on shaping the immune landscape, with potential to dictate the clinical response to immune checkpoint inhibition.
Isocitrate dehydrogenase (IDH) genes 1 and 2 are frequently mutated in acute myeloid leukaemia (AML), low-grade glioma, cholangiocarcinoma (CC) and chondrosarcoma (CS). For AML, low-grade glioma and ...CC, mutant IDH status is associated with a DNA hypermethylation phenotype, implicating altered epigenome dynamics in the aetiology of these cancers. Here we show that the IDH variants in CS are also associated with a hypermethylation phenotype and display increased production of the oncometabolite 2-hydroxyglutarate, supporting the role of mutant IDH-produced 2-hydroxyglutarate as an inhibitor of TET-mediated DNA demethylation. Meta-analysis of the acute myeloid leukaemia, low-grade glioma, cholangiocarcinoma and CS methylation data identifies cancer-specific effectors within the retinoic acid receptor activation pathway among the hypermethylated targets. By analysing sequence motifs surrounding hypermethylated sites across the four cancer types, and using chromatin immunoprecipitation and western blotting, we identify the transcription factor EBF1 (early B-cell factor 1) as an interaction partner for TET2, suggesting a sequence-specific mechanism for regulating DNA methylation.