Search

PURPOSE Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with a 10% annual risk of progression. Various prognostic models exist for risk stratification; however, those are based on solely clinical metrics. The discovery of genomic alterations that underlie disease progression to MM could improve current risk models. METHODS We used next-generation sequencing to study 214 patients with SMM. We performed whole-exome sequencing on 166 tumors, including 5 with serial samples, and deep targeted sequencing on 48 tumors. RESULTS We observed that most of the genetic alterations necessary for progression have already been acquired by the diagnosis of SMM. Particularly, we found that alterations of the mitogen-activated protein kinase pathway ( KRAS and NRAS single nucleotide variants [SNVs]), the DNA repair pathway (deletion 17p, TP53, and ATM SNVs), and MYC (translocations or copy number variations) were all independent risk factors of progression after accounting for clinical risk staging. We validated these findings in an external SMM cohort by showing that patients who have any of these three features have a higher risk of progressing to MM. Moreover, APOBEC associated mutations were enriched in patients who progressed and were associated with a shorter time to progression in our cohort. CONCLUSION SMM is a genetically mature entity whereby most driver genetic alterations have already occurred, which suggests the existence of a right-skewed model of genetic evolution from monoclonal gammopathy of undetermined significance to MM. We identified and externally validated genomic predictors of progression that could distinguish patients at high risk of progression to MM and, thus, improve on the precision of current clinical models.

Multiple myeloma (MM) is a plasma-cell neoplasm that is treated with high-dose chemotherapy, autologous stem cell transplant (ASCT) and long-term immunomodulatory drug (IMiD) maintenance. The presence of somatic mutations in the peripheral blood is termed clonal hematopoiesis of indeterminate potential (CHIP) and is associated with adverse outcomes. Targeted sequencing of the stem cell product from 629 MM patients treated by ASCT at the Dana-Farber Cancer Institute (2003–2011) detects CHIP in 136/629 patients (21.6%). The most commonly mutated genes are DNMT3A, TET2, TP53, ASXL1 and PPM1D. Twenty-one from fifty-six patients (3.3%) receiving first-line IMiD maintenance develop a therapy-related myeloid neoplasm (TMN). However, regardless of CHIP status, the use of IMiD maintenance associates with improved PFS and OS. In those not receiving IMiD maintenance, CHIP is associated with decreased overall survival (OS) (HR:1.34, p = 0.02) and progression free survival (PFS) (HR:1.45, p, Multiple myeloma (MM) is treated with induction chemotherapy, autologous stem cell transplant (ASCT) and long-term immunomodulatory drug (IMiD) maintenance. Here, the authors show that the presence of clonal haematopoiesis of indeterminate potential (CHIP) at time of ASCT is associated with adverse outcomes in MM patients.

Clonal hematopoiesis results from somatic mutations in cancer driver genes in hematopoietic stem cells. We sought to identify novel drivers of clonal expansion using an unbiased analysis of sequencing data from 84,683 persons and identified common mutations in the 5-methylcytosine reader, ZBTB33, as well as in YLPM1, SRCAP, and ZNF318. We also identified these mutations at low frequency in myelodysplastic syndrome patients. Zbtb33 edited mouse hematopoietic stem and progenitor cells exhibited a competitive advantage in vivo and increased genome-wide intron retention. ZBTB33 mutations potentially link DNA methylation and RNA splicing, the two most commonly mutated pathways in clonal hematopoiesis and MDS.

Cancer immunotherapy with checkpoint blockade (CPB) leads to improved outcomes in melanoma and other tumor types, but a majority of patients do not respond. High tumor mutation burden (TMB) and high levels of tumor-infiltrating T cells have been associated with response to immunotherapy, but integrative models to predict clinical benefit using DNA or RNA alone have not been comprehensively explored. We sequenced DNA and RNA from melanoma patients receiving CPB, and aggregated previously published data, yielding whole exome sequencing data for 189 patients and bulk RNA sequencing data for 178 patients. Using these datasets, we derived genomic and transcriptomic factors that predict overall survival (OS) and response to immunotherapy. Using whole-exome DNA data alone, we calculated T cell burden (TCB) and B cell burden (BCB) based on rearranged TCR/Ig DNA sequences and found that patients whose melanomas have high TMB together with either high TCB or high BCB survived longer and had higher response rates as compared to patients with either low TMB or TCB/BCB. Next, using bulk RNA-Seq data, differential expression analysis identified 83 genes associated with high or low OS. By combining pairs of immune-expressed genes with tumor-expressed genes, we identified three gene pairs associated with response and survival (Bonferroni PPMAP4K1 (Mitogen- Activated Protein Kinase Kinase Kinase Kinase 1) combined with the transcription factor TBX3 (T-Box Transcription Factor 3) which is overexpressed in poorly differentiated melanomas. We conclude that RNA-based (MAP4K1&TBX3) or DNA-based (TCB&TMB) models combining immune and tumor measures improve predictions of outcome after checkpoint blockade in melanoma.

Purpose: Leiomyosarcoma and liposarcoma are common subtypes of soft tissue sarcoma (STS). Patients with metastatic leiomyosarcoma or dedifferentiated liposarcoma (DDLPS) typically have worse outcomes compared with localized leiomyosarcoma or well-differentiated liposarcoma (WDLPS). A better understanding of genetic changes between primary/metastatic leiomyosarcoma and between WDLPS/DDLPS may provide insight into their genetic evolution. Experimental Design: We interrogated whole-exome sequencing (WES) from “trios” of normal tissue, primary tumor, and metastatic tumor from individual patients with leiomyosarcoma (n = 9), and trios of normal tissue, well-differentiated tumor, and dedifferentiated tumor from individual patients with liposarcoma (n = 19). Specifically, we performed mutational, copy number, and tumor evolution analyses on these cohorts and compared patterns among leiomyosarcoma and liposarcoma trios. Results: Leiomyosarcoma cases harbored shared drivers through a typical parent/child relationship where the metastatic tumor was derived from the primary tumor. In contrast, while all liposarcoma cases shared the characteristic focal chromosome 12 amplicon, most paired liposarcoma cases did not share additional mutations, suggesting a divergent evolutionary pattern from a common precursor. No highly recurrent genomic alterations from WES were identified that could be implicated as driving the progression of disease in either sarcoma subtype. Conclusions: From a genomic perspective, leiomyosarcoma metastases contain genetic alterations that are also found in primary tumors. WDLPS and DDLPS, however, appear to divergently evolve from a common precursor harboring 12q amplification, rather than as a transformation to a higher-grade tumor. Further efforts to identify specific drivers of these distinct evolutionary patterns may inform future translational and clinical research in STS.

How the genomic features of a patient’s cancer relate to individual disease kinetics remains poorly understood. Here we used the indolent growth dynamics of chronic lymphocytic leukaemia (CLL) to analyse the growth rates and corresponding genomic patterns of leukaemia cells from 107 patients with CLL, spanning decades-long disease courses. We found that CLL commonly demonstrates not only exponential expansion but also logistic growth, which is sigmoidal and reaches a certain steady-state level. Each growth pattern was associated with marked differences in genetic composition, the pace of disease progression and the extent of clonal evolution. In a subset of patients, whose serial samples underwent next-generation sequencing, we found that dynamic changes in the disease course of CLL were shaped by the genetic events that were already present in the early slow-growing stages. Finally, by analysing the growth rates of subclones compared with their parental clones, we quantified the growth advantage conferred by putative CLL drivers in vivo. Analysis of growth dynamics in a dataset from 107 patients with chronic lymphocytic leukaemia (CLL) reveals both exponential and logistic patterns of growth, which are associated with differences in genetic attributes and clinical outcomes.

Genomics of radiation-induced damage The potential adverse effects of exposures to radioactivity from nuclear accidents can include acute consequences such as radiation sickness, as well as long-term sequelae such as increased risk of cancer. There have been a few studies examining transgenerational risks of radiation exposure but the results have been inconclusive. Morton et al. analyzed papillary thyroid tumors, normal thyroid tissue, and blood from hundreds of survivors of the Chernobyl nuclear accident and compared them against those of unexposed patients. The findings offer insight into the process of radiation-induced carcinogenesis and characteristic patterns of DNA damage associated with environmental radiation exposure. In a separate study, Yeager et al. analyzed the genomes of 130 children and parents from families in which one or both parents had experienced gonadal radiation exposure related to the Chernobyl accident and the children were conceived between 1987 and 2002. Reassuringly, the authors did not find an increase in new germline mutations in this population. Science , this issue p. eabg2538 , p. 725

Diffuse large B-cell lymphoma (DLBCL), the most common lymphoid malignancy in adults, is a clinically and genetically heterogeneous disease that is further classified into transcriptionally defined activated B-cell (ABC) and germinal center B-cell (GCB) subtypes. We carried out a a comprehensive genetic analysis of 304 primary DLBCLs that identifies low-frequency alterations, captured recurrent mutations, somatic copy number alterations (SCNAs) and structural variants (SVs), and defined coordinate signatures in patients with available outcome data. We integrated these genetic drivers using consensus clustering and identified five robust DLBCL subsets, including a previously unrecognized group of low-risk ABC-DLBCLs of extrafollicular/marginal zone origin; two distinct subsets of GCB-DLBCLs with different outcomes and targetable alterations; and an ABC/GCB-independent group with biallelic inactivation of TP53, CDKN2A loss and associated genomic instability. The genetic features of the newly characterized subsets, their mutational signatures and the temporal ordering of identified alterations provide new insights into DLBCL pathogenesis. The coordinate genetic signatures also predict outcome independent of the clinical International Prognostic Index and suggest new combination treatment strategies. More broadly, our results provide a roadmap for an actionable DLBCL classification.

Introduction: Multiple Myeloma (MM) is an incurable plasma cell malignancy commonly preceded by the asymptomatic stage smoldering multiple myeloma (SMM). MM is characterized with significant genomic heterogeneity of chromosomal gains and losses (CNVs), translocations, and point mutations (SNVs); alterations that are also observed in SMM patients. However, current SMM risk models rely solely on clinical markers and do not accurately capture progression risk. While incorporating some genomic biomarkers improves prediction, using all MM genomic features to comprehensively stratify patients may increase risk stratification precision in SMM. Methods: We obtained a total of 214 patient samples at SMM diagnosis. We performed whole-exome sequencing on 166 tumors; of these, RNA sequencing was performed on 100. Targeted capture was done on 48 additional tumors. Upon binarization of DNA features, we performed consensus non-negative matrix factorization to identify distinct molecular clusters. We then trained a random forest classifier on translocations, SNVs, and CNVs. The predicted clinical outcomes for the molecular subtypes were further validated in an independent SMM cohort of 74 patients. Results: We identified six genomic subtypes, four with hyperdiploidy (>48 chromosomes, HMC, HKR, HNT, HNF) and two with IgH translocations (FMD, CND) (Table 1). In multivariate analysis accounting for IMWG (20-2-20) clinical risk stages, high-risk (HMC, FMD, HKR) and intermediate-risk (HNT, HNF) genetic subtypes were independent predictors of progression (Hazards ratio [HR]: 3.8 and 5.5, P = 0.016 and 0.001, respectively). The low-risk, CND subtype harboring translocation (11;14) was enriched for the previously defined CD-2 MM signature defined by the B cell markers CD20 and CD79A (FDR = 0.003 ), showed upregulation of CCND1, E2F1, and E2F7 (FDR = 0.01, 0.0004, 0.08), and was enriched for G2M checkpoint, heme metabolism, and monocyte cell signature (FDR = 0.003, 0.003, 0.003, respectively). The FMD subtype with IgH translocations (4;14) and (14;16) was enriched for P53, mTORC1, unfolded protein signaling pathways and plasmacytoid dendritic cell signatures (FDR = 0.01, 0.005, 0.008, respectively). The HKR tumors were enriched for inflammatory cytokine signaling, MYC target genes, T regulatory cell signature, and the MM proliferative (PR) signatures (FDR = 0.02, 0.03, 0.007, 0.02, respectively). The APOBEC mutational signature was enriched in HMC and FMD tumors (P = 0.005), while there was no statistical difference across subtypes in the AID signature. The median follow-up for the primary cohort is 7.1 years. Median TTP for patients in HMC, FMD, and HKR was 3.8, 2.6, and 2.2 years, respectively; TTP for HNT and HNF was 4.3 and 5.2, respectively, while it was 11 years in CND patients (P = 0.007). Moreover, by analyzing the changes in MM clinical biomarkers over time, we found that patients from high-risk subgroups had higher odds of developing evolving hemoglobin and monoclonal protein levels over time (P = 0.01 and 0.002, respectively); Moreover, the absolute increase in M-protein was significantly higher in patients from the high-risk genetic subtypes at one, two, and five years from diagnosis (P = 0.001, 0.03, and 0,01, respectively). Applying the classifier to the external cohort replicated our findings where intermediate and high-risk genetic subgroups conferred increased risk of progression to MM in multivariate analysis after accounting for IMWG staging (HR: 5.5 and 9.8, P = 0.04 and 0.005, respectively). Interestingly, within the intermediate-risk clinical group in the primary cohort, patients in the high-risk genetic subgroups had increased risk of progression (HR: 5.2, 95% CI 1.5 - 17.3, P = 0.007). In the validation cohort, these patients also had an increased risk of progression to MM (HR: 6.7, 95% CI 1.2 - 38.3, P = 0.03), indicating that molecular classification improves the clinical risk-stratification models. Conclusion: We identified and validated in an independent dataset six SMM molecular subgroups with distinct DNA alterations, transcriptional profiles, dysregulated pathways, and risks of progression to active MM. Our results underscore the importance of molecular classification in addition to clinical evaluation in better identifying high-risk SMM patients. Moreover, these subgroups may be used to identify tumor vulnerabilities and target them with precision medicine efforts. Figure 1 Figure 1. Disclosures Bustoros: Janssen, Bristol Myers Squibb: Honoraria, Speakers Bureau; Takeda: Consultancy, Honoraria. Casneuf: Janssen: Current Employment. Kastritis: Amgen: Consultancy, Honoraria, Research Funding; Takeda: Honoraria; Pfizer: Consultancy, Honoraria, Research Funding; Genesis Pharma: Honoraria; Janssen: Consultancy, Honoraria, Research Funding. Walker: Bristol Myers Squibb: Research Funding; Sanofi: Speakers Bureau. Davies: Takeda: Consultancy, Honoraria; Amgen: Consultancy, Honoraria; Abbvie: Consultancy, Honoraria; BMS: Consultancy, Honoraria; Roche: Consultancy, Honoraria; Janssen: Consultancy, Honoraria. Dimopoulos: Amgen: Honoraria; BMS: Honoraria; Takeda: Honoraria; Beigene: Honoraria; Janssen: Honoraria. Bergsagel: Genetech: Consultancy, Honoraria; Oncopeptides: Consultancy, Honoraria; Janssen: Consultancy, Honoraria; Pfizer: Consultancy, Honoraria; Novartis: Consultancy, Honoraria, Patents & Royalties: human CRBN mouse; GSK: Consultancy, Honoraria; Celgene: Consultancy, Honoraria. Yong: BMS: Research Funding; Autolus: Research Funding; Takeda: Honoraria; Janssen: Honoraria, Research Funding; Sanofi: Honoraria, Research Funding; GSK: Honoraria; Amgen: Honoraria. Morgan: BMS: Membership on an entity's Board of Directors or advisory committees; Jansen: Membership on an entity's Board of Directors or advisory committees; Karyopharm: Membership on an entity's Board of Directors or advisory committees; Oncopeptides: Membership on an entity's Board of Directors or advisory committees; GSK: Membership on an entity's Board of Directors or advisory committees. Getz: IBM, Pharmacyclics: Research Funding; Scorpion Therapeutics: Consultancy, Current holder of individual stocks in a privately-held company, Membership on an entity's Board of Directors or advisory committees. Ghobrial: AbbVie, Adaptive, Aptitude Health, BMS, Cellectar, Curio Science, Genetch, Janssen, Janssen Central American and Caribbean, Karyopharm, Medscape, Oncopeptides, Sanofi, Takeda, The Binding Site, GNS, GSK: Consultancy.

Background: Multiple Myeloma (MM) is an incurable plasma cell malignancy with significant genomic heterogeneity. It is usually preceded by the asymptomatic stage known as smoldering multiple myeloma (SMM). SMM patients have a 10% annual risk of progression to MM. Genomic alterations that are observed in SMM patients include chromosomal gains and losses, translocations, and point mutations. However, current SMM risk models rely solely on clinical markers that do not accurately capture the progression risk. While incorporating some genomic biomarkers improves prediction, using all MM genomic features to comprehensively stratify patients may increase the precision of risk models. Methods: We obtained a total of 214 patients' samples at SMM diagnosis in the US and Europe. We performed whole exome sequencing on 166 tumors; of these, RNA sequencing was performed on 100. Targeted capture with a MM gene panel was done on an additional 48 tumors. We identified subgroups using binarized DNA features and performing consensus binary non-negative matrix factorization. Results: We identified six clusters (C1-C6) with the following features: four with a hyperdiploidy (HD) (>48 chromosomes) and two with IgH translocations. These subgroups have unique transcriptomic profiles overlapping with known MM signatures and biological pathways. One of the clusters harboring translocation (11;14), which we call C4-CCND1, was enriched with the previously defined CD-2 MM signature that uniquely expresses B cell markers CD20 and CD79A; shows upregulation of CCND1 and E2F7; and is enriched with pathways like DNA replication, heme metabolism, and NFkB signaling. The C3-MS_MF cluster with the IgH translocations (4;14) and (14;16) shows downregulation of ribosomal genes, TRAF2, and DUSP2. The MYC oncogene was highly expressed in the four HD clusters: C1-HD_NRAS, C2-HD_MAFB, C5-HD_KRAS, and C6-HD_1q (BH-P = 0.037). The clusters also showed different outcomes in terms of time to progression (TTP) to active MM (P = 0.005). Median TTP for patients in C2-HD_MAFB, C3-MS_MF, and C5-HD_KRAS was 3.7, 2.6, and 2.2 years, respectively; TTP for C1-HD_NRAS, C4-CCND1, and C6-HD_1q was 4.3, 11, and not reached, respectively. In multivariate analysis, C2-HD_MAFB, C3-MS_MF, and C5-HD_KRAS were independent predictors of progression after accounting for the clinical risk stage. Moreover, the odds of having evolving hemoglobin and monoclonal protein levels in these three clusters were 3.5 and 12.3 times higher than the other clusters, respectively (P = 0.01 and 0.002). Conclusion: We identified six distinct SMM molecular groups with corresponding transcription profiles and dysregulated pathways. These groups have different progression risks to active MM, with three groups being independent predictors of progression. Our results underscore the importance of molecular classification in MM to better understand and target various tumor vulnerabilities. Citation Format: Shankara Anand, Mark Bustoros, François Aguet, Romanos Sklavenitis-Pistofidis, Robert Redd, Benny Zhitomirsky, Andrew J. Dunford, Yu-Tzu Tai, Selina J. Chavda, Cody Boehner, Carl J. Neuse, Tineke Casneuf, Lorenzo Trippa, Chip Stewart, Kwee Yong, Irene Ghobrial, Gad Getz. Genomic profiling of smoldering multiple myeloma classifies distinct molecular groups [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2240.

The advent of PD-1/PD-L1 agents has transformed the therapeutic landscape of many advanced cancers, including non-small cell lung cancer (NSCLC). However, our understanding of the genomic biomarkers underlying effective treatment response remain limited. Here we provide updated results from our ongoing effort, the Stand Up To Cancer Lung (SU2C-Lung)/Mark Foundation EXTOLConsortium, a multi-institution collaboration to expand our understanding of the molecular determinants of immunotherapy response in NSCLC. Comprising a set of nearly 400 patients, this cohort enables evaluation of both genomic and transcriptomic factors associated with checkpoint blockade response. In addition to validating previously known associations including TMB and neoantigen burden, we examined transcriptional predictors associated with response. We generated a list of differentially expressed genes with respect to best overall response (BOR), and performed dimensionality reduction using semi-supervised Bayesian Non-Negative Factorization (ssBNMF). We identified 3 distinct clusters with strong sample membership. Characterization of these subtypes revealed varying levels of immune infiltrate, histologic composition, and response rates to checkpoint blockade. Of these three subtypes, two were associated with low response rates to PD-1/PD-L1 blockade, suggesting the existence of distinct avenues toward resistance. To further characterize these transcriptional subtypes, we used ssBNMF marker genes to classify publicly available NSCLC samples from The Cancer Genome Atlas (TCGA), along with gene expression from a smaller cohort of large-cell neuroendocrine (LCNE) samples. Histologic composition showed good concordance with our SU2C samples, and redemonstrated observations within our smaller SU2C cohort of distinct immuno-suppressive and immuno-depleted milieus associated with resistance. Citation Format: Monica B. Arniella, Arvind Ravi, Justin Gainor, Chip Stewart, Sam Freeman, Mark Awad, Patrick Forde, Valsamo Anagnostou, Brian Henick, Jonathan W. Riess, Don Gibbons, Nathan Pennell, Vamsidhar Velcheti, Ignaty Leshchiner, Jaegil Kim, Subba Digumarthy, Mari Mino-Kenudson, John Heymach, Nir Hacohen, Naiyer Rizvi, Roy Herbst, Victor E. Velculescu, Julie Brahmer, Kurt Schalper, Pasi Jänne, Jedd Wolchok, Alice Shaw, Gad Getz, Matthew D. Hellman. An SU2C-Mark Foundation Lung collaborative update: integrative genomics identifies distinct transcriptional states associated with checkpoint blockade resistance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB197.

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid cancer (PTC) incidence in surrounding regions, particularly for 131I-exposed children. To investigate the contribution of environmental radiation to PTC characteristics and improve understanding of radiation-induced carcinogenesis, we analyzed genomic, transcriptomic, and epigenomic characteristics of 440 pathologically-confirmed fresh-frozen PTCs from Ukraine (359 with estimated childhood or in utero 131I exposure and 81 from unexposed children born after March 1987) and matched normal tissue (non-tumor thyroid tissue and/or blood). Mean age at PTC was 28.0 years (range: 10.0-45.6). Among 131I-exposed individuals, mean radiation dose was 250 mGy (range: 11.0-8,800). In multivariable models adjusted for age at PTC and sex, we observed radiation dose-dependent enrichment of fusion drivers (P=6.6 × 10−8), nearly all occurring in the MAPK pathway, as well as increases in small deletions (P=8.0 × 10−9) and simple/balanced structural variants (P=1.2 × 10−14). Further analyses demonstrated even stronger associations for those small deletions and simple/balanced structural variants that were clonal and bore hallmarks of non-homologous end-joining repair (deletions: P=4.9 × 10−31; simple/balanced structural variants: P=5.5 × 10−19). In contrast, radiation dose was not associated with subclonal small deletions (P=0.82) or subclonal simple/balanced structural variants (P=0.91). Additionally, radiation dose was not associated with TINS (locally templated insertions), which are characteristic of alt-end-joining repair (P=0.69). The effects of radiation on genomic alterations with more pronounced for those younger at exposure. Analyses generally were consistent with a linear radiation dose-response for all molecular characteristics except clonal small deletions. Analyses of transcriptomic and epigenomic features demonstrated strong associations with the PTC driver gene but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth following environmental radiation exposure. Citation Format: Lindsay M. Morton, Danielle Karyadi, Chip Stewart, Tetiana Bogdanova, Eric Dawson, Mia Steinberg, Jieqiong Dai, Stephen Hartley, Sara Schonfeld, Joshua Sampson, Yosi Maruvka, Vidushi Kapoor, Dale Ramsden, Juan Carvajal-Garcia, Chuck Perou, Joel Parker, Marko Krznaric, Meredith Yeager, Joseph Boland, Amy Hutchinson, Belynda Hicks, Casey Dagnall, Julie Gastier-Foster, Jay Bowen, Olivia Lee, Mitchell Machiela, Elizabeth Cahoon, Alina Brenner, Kiyohiko Mabuchi, Vladimir Drozdovitch, Sergii Masiuk, Mykola Chepurny, Liudmyla Yu Zurnadzhy, Maureen Hatch, Amy Berrington de Gonzalez, Gerry Thomas, Mykola Tronko, Gad Getz, Stephen Chanock. Molecular characterization of papillary thyroid cancer in relation to ionizing radiation dose following the Chernobyl accident [abstract]. In: Proceedings of the AACR Virtual Special Conference on Radiation Science and Medicine; 2021 Mar 2-3. Philadelphia (PA): AACR; Clin Cancer Res 2021;27(8_Suppl):Abstract nr PO-055.

Classical Hodgkin lymphoma (cHL) is composed of rare malignant Hodgkin Reed-Sternberg (HRS) cells within an extensive, but ineffective, inflammatory/immune cell infiltrate. HRS cells exhibit near-universal somatic copy gains of chromosome 9p/9p24.1, which increase expression of the programmed cell death protein 1 (PD-1) ligands. To define genetic mechanisms of response and resistance to PD-1 blockade and identify complementary treatment targets, we performed whole-exome sequencing of flow cytometry–sorted HRS cells from 23 excisional biopsies of newly diagnosed cHLs, including 8 Epstein-Barr virus–positive (EBV(+)) tumors. We identified significantly mutated cancer candidate genes (CCGs) as well as somatic copy number alterations and structural variations and characterized their contribution to disease-defining immune evasion mechanisms and nuclear factor κB (NF-κB), JAK/STAT, and PI3K signaling pathways. EBV(–) cHLs had a higher prevalence of genetic alterations in the NF-κB and major histocompatibility complex class I antigen presentation pathways. In this young cHL cohort (median age, 26 years), we identified a predominant mutational signature of spontaneous deamination of cytosine- phosphate-guanines (“Aging”), in addition to apolipoprotein B mRNA editing catalytic polypeptide-like, activation-induced cytidine deaminase, and microsatellite instability (MSI)–associated hypermutation. In particular, the mutational burden in EBV(–) cHLs was among the highest reported, similar to that of carcinogen-induced tumors. Together, the overall high mutational burden, MSI-associated hypermutation, and newly identified genetic alterations represent additional potential bases for the efficacy of PD-1 blockade in cHL. Of note, recurrent cHL alterations, including B2M, TNFAIP3, STAT6, GNA13, and XPO1 mutations and 2p/2p15, 6p21.32, 6q23.3, and 9p/9p24.1 copy number alterations, were also identified in >20% of primary mediastinal B-cell lymphomas, highlighting shared pathogenetic mechanisms in these diseases.

Driver mutations alter cells from normal to cancer through several evolutionary epochs: premalignancy, early malignancy, subclonal diversification, metastasis and resistance to therapy. Later stages of disease can be explored through analyzing multiple samples collected longitudinally, on or between successive treatments, and finally at time of autopsy. It is also possible to study earlier stages of cancer development through probabilistic reconstruction of developmental trajectories based on mutational information preserved in the genome. Here we present a suite of tools, called Phylogic N-Dimensional with Timing (PhylogicNDT), that statistically model phylogenetic and evolutionary trajectories based on mutation and copy-number data representing samples taken at single or multiple time points. PhylogicNDT can be used to infer: (i) the order of clonal driver events (including in pre-cancerous stages); (ii) subclonal populations of cells and their phylogenetic relationships; and (iii) cell population dynamics. We demonstrate the use of PhylogicNDT by applying it to whole-exome and whole-genome data of 498 lung adenocarcinoma samples (434 previously available and 64 of newly generated data). We identify significantly different progression trajectories across subtypes of lung adenocarcinoma (EGFR mutant, KRAS mutant, fusion-driven and EGFR/KRAS wild type cancers). In addition, we study the progression of fusion-driven lung cancer in 21 patients by analyzing samples from multiple timepoints during treatment with 1st and next generation tyrosine kinase inhibitors. We characterize their subclonal diversification, dynamics, selection, and changes in mutational signatures and neoantigen load. This methodology will enable a systematic study of tumour initiation, progression and resistance across cancer types and therapies.

Summary Large-scale genomic characterization of tumors from prospective cohort studies may yield new insights into cancer pathogenesis. We performed whole-exome sequencing of 619 incident colorectal cancers (CRCs) and integrated the results with tumor immunity, pathology, and survival data. We identified recurrently mutated genes in CRC, such as BCL9L, RBM10, CTCF, and KLF5, that were not previously appreciated in this disease. Furthermore, we investigated the genomic correlates of immune-cell infiltration and found that higher neoantigen load was positively associated with overall lymphocytic infiltration, tumor-infiltrating lymphocytes (TILs), memory T cells, and CRC-specific survival. The association with TILs was evident even within microsatellite-stable tumors. We also found positive selection of mutations in HLA genes and other components of the antigen-processing machinery in TIL-rich tumors. These results may inform immunotherapeutic approaches in CRC. More generally, this study demonstrates a framework for future integrative molecular epidemiology research in colorectal and other malignancies., Graphical Abstract, Highlights • Whole-exome sequencing of 619 colorectal cancers with clinicopathologic annotations • Discovery of significantly mutated genes in colorectal cancer • Neoantigen load correlation with infiltrating lymphocytes and memory T cells • Positive selection for HLA mutations in immune-cell-infiltrated tumors, Through whole-exome sequencing of annotated colorectal tumors, Giannakis et al. identify additional colorectal cancer driver genes and correlate high neoantigen load with increased lymphocytic infiltration and improved survival. They also find positive selection for HLA mutations in immune-cell-infiltrated tumors. These results may inform immunotherapeutic approaches in colorectal cancer.

Primary central nervous system lymphomas (PCNSLs) and primary testicular lymphomas (PTLs) are extranodal large B-cell lymphomas (LBCLs) with inferior responses to current empiric treatment regimens. To identify targetable genetic features of PCNSL and PTL, we characterized their recurrent somatic mutations, chromosomal rearrangements, copy number alterations (CNAs), and associated driver genes, and compared these comprehensive genetic signatures to those of diffuse LBCL and primary mediastinal large B-cell lymphoma (PMBL). These studies identify unique combinations of genetic alterations in discrete LBCL subtypes and subtype-selective bases for targeted therapy. PCNSLs and PTLs frequently exhibit genomic instability, and near-uniform, often biallelic, CDKN2A loss with rare TP53 mutations. PCNSLs and PTLs also use multiple genetic mechanisms to target key genes and pathways and exhibit near-uniform oncogenic Toll-like receptor signaling as a result of MYD88 mutation and/or NFKBIZ amplification, frequent concurrent B-cell receptor pathway activation, and deregulation of BCL6. Of great interest, PCNSLs and PTLs also have frequent 9p24.1/PD-L1/PD-L2 CNAs and additional translocations of these loci, structural bases of immune evasion that are shared with PMBL.

Juvenile myelomonocytic leukemia (JMML) is a myeloproliferative neoplasm (MPN) of childhood with a poor prognosis. Mutations in NF1, NRAS, KRAS, PTPN11 and CBL occur in 85% of patients, yet there are currently no risk stratification algorithms capable of predicting which patients will be refractory to conventional treatment and therefore be candidates for experimental therapies. In addition, there have been few other molecular pathways identified aside from the Ras/MAPK pathway to serve as the basis for such novel therapeutic strategies. We therefore sought to genomically characterize serial samples from patients at diagnosis through relapse and transformation to acute myeloid leukemia in order to expand our knowledge of the mutational spectrum in JMML. We identified recurrent mutations in genes involved in signal transduction, gene splicing, the polycomb repressive complex 2 (PRC2) and transcription. Importantly, the number of somatic alterations present at diagnosis appears to be the major determinant of outcome.

PD-1/PD-L1 checkpoint blockade has been a landmark advance for many patients suffering from advanced non-small cell lung cancer (NSCLC). However, detailed biomarkers of response beyond tumor mutational burden (TMB) are still poorly understood. As part of the effort to elucidate these additional signatures, we describe our progress on the Stand Up to Cancer Lung (SU2C-Lung) cohort. Initial characterization of exomes recapitulates mutational and copy number profiles seen in The Cancer Genome Atlas (TCGA) project. To better define expression subtypes using RNA sequencing, we performed non-negative matrix factorization (NMF) across an aggregated set of publicly available NSCLC expression data (including adenocarcinoma, squamous cell carcinoma, and large cell neuroendocrine histologies), and demonstrate good concordance in the SU2C-Lung cohort between this expression-based classifier and clinically annotated histology. To gain further insight into how immune cell infiltrates vary across our cohort, we additionally tested two common deconvolution algorithms, EPIC and CIBERSORT. While these two methods agree for some prominent cell types, such as B cells and CD4 T cells, discrepancies in minor infiltrating components such as NK cells may suggest a limit to the inference of rare subpopulations from bulk sequencing data.Finally, we describe a novel approach for determining single-gene predictors of response. Using the method, which is based on comparison of top single transcriptional features identified from random bootstraps of the full cohort as compared to a set of background shuffles, we are able to show that we remain powered for discovery of RNA response biomarkers despite the typically burdensome toll of multiple hypothesis correction at genome wide scale. Acknowledgment: Supported by Stand Up To Cancer-American Cancer Society Lung Cancer Dream Team Translational Research Grant SU2C-AACR-DT17-15. Citation Format: Monica Arniella, Arvind Ravi, Chip Stewart, Sam Freeman, Mark Awad, Patrick Forde, Valsamo Anagnostou, Brian Henick, Jonathan Riess, Don Gibbons, Nathan Pennell, Vamsidhar Velcheti, Ignaty Leshchiner, Jaegil Kim, Subba Digumarthy, Mari Mino-Kenudson, John Heymach, Nir Hacohen, Naiyer Rizvi, Roy Herbst, Victor Velculescu, Julie Brahmer, Kurt Schalper, Pasi Jänne, Jedd Wolchok, Alice Shaw, Justin Gainor, Matthew Hellmann, Gad Getz. Integrative genomic analysis of checkpoint blockade in lung cancer: A multi-institution SU2C collaborative [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5902.

Cancer immunotherapy with checkpoint blockade has improved survival and outcomes in melanoma, but still a majority of patients do not respond. Both high tumor mutation burden (TMB) and high T cell infiltration have been associated with response, but integrative models based on DNA or RNA assays have not been comprehensively explored and validated. Focusing on melanomas from patients receiving checkpoint blockade, we generated new and aggregated existing datasets of whole exome sequencing (WES) (n = 189 total) and bulk RNA sequencing (n = 154 total) to derive genomic and transcriptomic factors that predict survival and response to immunotherapy in melanoma. We quantified T and B cell infiltrates using rearranged T cell receptor (TCR) and immunoglobulin (Ig) sequences, respectively, from DNA or RNA sequencing. High levels of rearranged TCR reads or rearranged Ig reads in RNA-seq were associated with survival (P = 0.0046, P = 0.015) and response (P = 0.0034, P = 0.047). We created RNA-based metrics of T and B cell burden (TCBRNA or BCBRNA) by normalizing the number of rearranged TCR reads by the total number of mapped reads. When we analyzed WES data in patients for whom DNA and RNA were extracted from the same region, we found that the TCBDNA correlated with TCBRNA (rho = 0.73) and BCBDNA with BCBRNA (rho = 0.41), demonstrating that the level of lymphocyte infiltration can be estimated using rearranged TCR or Ig reads from tumor WES alone. We found that TCBDNA and BCBDNA both associated with survival (P = 0.0023 and 0.0089). In a combined model, patients with high TMB and high TCB DNA survived longer (P = 2.4e-4, HR = 2.68) and had a higher response rate (Fisher P = 0.028). This combined model was superior to models with TMB or TCBDNA alone. Similarly, patients with high TMB and high BCBDNA had longer survival and higher response rates (log-rank P = 0.0029, HR = 2.64, Fisher P = 0.015). We reanalyzed stage III/IV melanomas from TCGA and found that the TMB high, TCBDNA high subgroup had increased survival (P = 0.007). Next, clustering of tumor transcriptomes identified 5 tumor subtypes based on melanocyte differentiation, immune infiltration and keratin levels. These melanoma subtypes were associated with survival outcomes after immunotherapy (P = 0.019). We found that TBX3, a tumor-expressed transcription factor enriched in poorly differentiated melanomas, was over-expressed among non-responders within the immune-infiltrated subtype and among all patients (P = 3.9e-4, P = 8.7e-5). Patients whose tumors had high immune infiltrate and low expression of TBX3 had longer survival (P = 1.6e-5, HR = 3.39), however this subgroup did not have longer survival in an independent cohort (n = 73, P = 0.10, HR = 2.63). In conclusion, we demonstrate both RNA-based (immune infiltrate and tumor subtype) and DNA-based metrics (TMB/TCB or TMB/BCB) can be used as pre-treatment predictors of survival after checkpoint blockade in melanoma. Citation Format: Samuel S. Freeman, Moshe Sade-Feldman, Jaegil Kim, Chip Stewart, Arvind Ravi, Monica Arniella, Keren Yizhak, Ignaty Leshchiner, Liudmila Elagina, Oliver Spiro, Dimitri Livitz, Daniel Rosebrock, François Aguet, Jian Carrot-Zhang, Anna Gonye, Gavin Ha, Ziao Lin, Jonathan H. Chen, Dennie T. Frederick, Michal Barzily-Rokni, Marc R. Hammond, Hans Vitzthum, Shauna M. Blackmon, Yunxin J. Jiao, Donald P. Lawrence, Lyn M. Duncan, Anat Stemmer-Rachamimov, Jennifer A. Wargo, Keith T. Flaherty, Genevieve M. Boland, Ryan J. Sullivan, Matthew Meyerson, Gad Getz, Nir Hacohen. Combined signals from tumor and immune cells predict outcomes of checkpoint inhibition in melanoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6670.

Chronic lymphocytic leukemia (CLL) accounts for one third of all diagnosed leukemias. It exhibits tremendous clinical variability, ranging from indolent disease kinetics to a fast progressing aggressive disease. Two large-scale studies by our groups, each analyzing more than 500 CLL subjects, have shown that a large part of this clinical variability can be explained by underlying molecular features of the cancerous B-cells, including genetic and epigenetic alterations. Despite this progress, patients do not yet fully benefit from these discoveries in their treatment plans, and prognosis is not tailored to them based on the molecular features of their cancer. To achieve precision medicine, we compile a large dataset of multi-omic and clinical data from over 1100 CLL patients. Our extensive cohort confers power to detect novel driver genes, doubling the number of known drivers (114 total) and identifying drivers specific to major CLL subtypes, such as IGHV-mutated and IGHV-unmutated. The analysis identifies recurrent copy number alterations, order-of-event timing and mutational processes acting in CLL. We integrate these with expression data from 700 CLL RNA-seqs and epigenetic data to define molecular subtypes that present differential clinical outcomes. Altogether, this work delineates the molecular landscape of CLL and sets the stage for improved classification, prognostics and treatment of CLL patients. Citation Format: Binyamin A. Knisbacher, Ziao Lin, Chip Stewart, Cynthia K. Hahn, Kristen E. Stevenson, Donna S. Neuberg, José I. Martín-Subero, Xose S. Puente, Elias Campo, Catherine J. Wu, Gad Getz. Discovery of clinically distinct CLL subgroups by integrative mapping of large-scale genetic, epigenetic, expression and clinical data [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 169.

Classical Hodgkin lymphoma (cHL) and primary mediastinal large B-cell lymphoma (PMBL) are aggressive tumors with distinct cells of origin and pathomorphological features. However, these lymphomas share certain transcriptional signatures and aberrant signaling pathways. CHLs and PMBLs both exhibit constitutive activation of NF-κB and JAK/STAT signaling and genetic bases of PD-1 mediated immune evasion including frequent 9p24.1/PD-L1/PD-L2 copy gains. In both lymphomas, PD-1 blockade is a FDA-approved therapy for relapsed/refractory disease. To characterize genetic bases of response to PD-1 blockade and identify complementary treatment targets in cHL and PMBL, we defined the comprehensive genetic signatures of both diseases. First, we obtained flow cytometry-sorted Hodgkin Reed Sternberg (HRS) cells from 23 biopsies of newly diagnosed cHLs and intact tumor biopsy specimens from 37 newly diagnosed PMBLs. The isolated HRS cells and paired normal DNAs and PMBL biopsy specimens were subjected to whole exome sequencing using an optimized workflow for low input samples and an expanded bait set to capture structural variants (SVs), including translocations. We used newly developed and established analytical pipelines to analyze tumor samples without paired normals (PMBLs) and identify significantly mutated genes (candidate cancer genes [CCGs], MutSig2CV, CLUMPS), SCNAs (GISTIC2.0) and SVs(4 algorithms) in both cHL and PMBL. In cHL, we identified 15 CCGs, 13 recurrent SCNAs, SVs in ETV6 and CIITA, complementary alterations of JAK/STAT, NF-κB and PI3K signaling pathway components and a median number of 11 genetic drivers per tumor. Previously unappreciated aspects of the cHL genetic signature included the increased incidence of driver mutational events in cHLs with ARID1A alterations (p=0.012). Analyses of co-occurring genetic events in EBV+ and EBV- cHLs confirmed that EBV- cHLs were significantly more likely to exhibit alterations of specific NF-κB signaling intermediaries (such as TNFAIP3 mutation and/or focal copy loss, p=0.006) and perturbations of MHC class I antigen presentation pathway components (inactivating B2M mutations, HLA-B mutations or focal copy loss of 6p21.32/HLA-B, p=0.008). The latter findings provide genetic bases for the reported differences in cell surface expression of MHC class I in EBV+ and EBV- cHLs. In PMBL, we defined 15 CCGs and more selective perturbations of specific epigenetic modifiers (ZNF217 and EZH2), transcription factors (PAX5 and IRF2BP2) and TP53, in comparison with cHL. The majority of these alterations were clonal supporting their role as early drivers. We identified 18 SCNAs and additional SVs in CIITA and PD-1 ligands, recurrent alterations of JAK/STAT and NF-κB signaling pathway components and a median of 9 genetic drivers per PMBL. Antigen presentation pathways in PMBL were perturbed by multiple recurrent alterations, including B2M mutations, focal copy losses of B2M and the MHCI/II loci, SVs of CTIIA and EZH2 mutations. There was a significant correlation between genetic perturbations of MHC class I pathway components and absence of MHC class I expression in PMBL, as previously described in cHL. Recurrent cHL alterations including B2M, TNFAIP3, STAT6, GNA13 and XPO1 CCGs and 2p/2p15/2p16.1, 6p21.32, 6q23.2 and 9p/9p24.1 SCNAs were also identified in >20% of PMBLs, highlighting shared pathogenetic mechanisms in these diseases. These tumors of predominantly young adults (median age: cHL 26 yrs; PMBL 34 yrs) both had a high rate of spontaneous deamination of CpGs, a clock-like mutational signature that is typically associated with aging. CHLs and PMBLs both exhibited previously uncharacterized molecular features that may increase sensitivity to PD-1 blockade, including high mutational burdens, in comparison with other lymphoid and solid tumors. In particular, the mutational burden in EBV- cHLs was among the highest reported, similar to that in carcinogen-induced cancers (melanoma and NSCLC). Additionally, both cHLs and PMBLs had an increased incidence of microsatellite instability and APOBEC mutational signatures, features associated with a more favorable response to PD-1 blockade. Taken together, these data define genetic similarities and differences in cHL and PMBL and establish a framework to comprehensively assess molecular bases of response to PD-1 blockade and develop rational combination therapies in these diseases. Disclosures Armand: Merck: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding, Speakers Bureau; Bristol-Myers Squibb: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Otsuka: Research Funding; Sigma Tau: Research Funding; Adaptive: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Affimed: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Roche: Research Funding; Pfizer: Consultancy; ADC Therapeutics: Consultancy; Infinity: Consultancy; Genentech: Research Funding; Tensha: Research Funding. Rodig:Merck: Research Funding; Affirmed: Research Funding; Kite, a Gilead Company: Research Funding; Bristol Myers Squib: Consultancy, Honoraria, Other: Travel Expenses, Speakers Bureau. Fromm:Merck, Inc.: Research Funding. Getz:Pharmacyclics: Research Funding; IBM: Research Funding; MuTect, ABSOLTUE, MutSig and POLYSOLVER: Patents & Royalties: MuTect, ABSOLTUE, MutSig and POLYSOLVER. Shipp:AstraZeneca: Honoraria, Membership on an entity's Board of Directors or advisory committees; Gilead Sciences: Honoraria, Membership on an entity's Board of Directors or advisory committees; Takeda Pharmaceuticals: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bayer: Research Funding; Merck & Co.: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; BMS: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding.

Diffuse large B-cell lymphoma (DLBCL) is a clinically and molecularly heterogeneous disease with recognized transcriptional subtypes associated with normal cells of origin, activated B-cell (ABC) and germinal center B-cell (GCB) tumors. Emerging data suggested that additional heterogeneity existed, prompting us to comprehensively characterize genomic signatures of 304 newly diagnosed DLBCLs from patients treated with state-of-the-art therapy. We integrated recurrent mutations, somatic copy number alterations (SCNAs) and structural variants (SVs) and identified 5 genetically distinct DLBCL clusters (C1- C5 DLBCLs; Chapuy, Stewart, Dunford, et al. Nat Med 2018). Specifically, we identified two genetically distinct ABC subtypes, including favorable-risk C1 DLBCLs with features of extrafollicular origin and alterations also seen in transformed marginal zone lymphomas (NOTCH2 and NF-κB pathway member mutations and BCL6 SVs). Unfavorable-risk C5 ABC DLBCLs harbored frequent 18q/BCL2 copy gain and co-occurring CD79B and MYD88L265P mutations. We also identified two genetically distinct GCB subtypes, including unfavorable-risk C3 DLBCLs with frequent BCL2 SVs, mutations in chromatin-modifying enzymes (CREBBP, MLL2, EZH2) and BCR/PI3K signaling pathway members (including inactivating PTEN mutations and copy loss). Favorable-risk C4 GCB DLBCLs had frequent mutations in core and linker histones and signaling intermediates (SGK1, BRAF and STAT3). Additionally, we identified an ABC/GCB-independent subtype, C2 DLBCLs, characterized by frequent bi-allelic TP53 inactivation, 9p21.23/CDKN2A copy loss and associated genomic instability reflected in recurrent SCNAs, increased genome doublings and a distinct outcome following induction therapy. A next step in utilizing the characterized genetic substructure was to confirm it in an independent series and develop a molecular classifier that allows prospective identification of C1-C5 DLBCLs. To this end, we accessed whole exome sequencing, copy number and SV data from a recent cohort of newly diagnosed DLBCLs (39 tumor-normal pairs, 462 tumor-only samples; Schmitz et al. NEJM 2018). All samples were re-analyzed using our mutational and SCNA pipelines and our newly generated tumor-only algorithm (Chapuy, Stewart, Dunford, et al. Nat Med 2018) to avoid batch effects and harmonize the datasets. SVs were used as reported. Purity and ploidy were inferred using ABSOLUTE and samples with missing data or low purity were removed. For the combined cohort (579 samples), we assessed our previously characterized 158 genetic drivers (Chapuy, Stewart, Dunford, et al. Nat Med 2018) and confirmed equal distribution of their marginal frequencies (R=0.88, p=1.5e-51), excluding batch effects. Next, we applied non-negative matrix factorization (NNF) consensus clustering to the combined dataset (158 genetic drivers vs. 579 tumors) and confirmed the C1-C5 DLBCL genetic clusters. Notably, tumors from both series contributed at comparable frequencies to the respective C1-C5 DLBCLs. We also noted an enrichment of the alternative genetic labels from Schmitz et al. in 3 of our C1-C5 DLBCL subtypes (B2N in C1 DLBCLs, p Next, we developed a molecular classifier that prospectively identified C1-C5 DLBCLs using a minimum number of easy-to-measure features. The NMF-defined classes of the combined cohort were used as gold-standard training and validation datasets. We tested different models for classification and selected an artificial neural network approach which provides accurate classification of individual samples and well-calibrated confidence metrics. To minimize potential overtraining, we developed a reduced input feature set of the 22 most discriminating features, constructed confidence metrics for each sample and trained an ensemble of Feed-Forward Neural Networks via 10-fold cross validation. With this approach, our classifier had 84% accuracy for the total set and 94% accuracy for the high-confidence samples (70% of all samples). The newly developed parsimonious classifier will allow prospective identification of the independently confirmed C1-C5 DLBCL subtypes in newly diagnosed patients, a necessity for clinical application. Disclosures Getz: Pharmacyclics: Research Funding; IBM: Research Funding; MuTect, ABSOLTUE, MutSig and POLYSOLVER: Patents & Royalties: MuTect, ABSOLTUE, MutSig and POLYSOLVER. Shipp:BMS: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Merck & Co.: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; AstraZeneca: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bayer: Research Funding; Gilead Sciences: Honoraria, Membership on an entity's Board of Directors or advisory committees; Takeda Pharmaceuticals: Honoraria, Membership on an entity's Board of Directors or advisory committees.

Summary Hurthle cell carcinoma of the thyroid (HCC) is a form of thyroid cancer recalcitrant to radioiodine therapy that exhibits an accumulation of mitochondria. We performed whole-exome sequencing on a cohort of primary, recurrent, and metastatic tumors, and identified recurrent mutations in DAXX, TP53, NRAS, NF1, CDKN1A, ARHGAP35, and the TERT promoter. Parallel analysis of mtDNA revealed recurrent homoplasmic mutations in subunits of complex I of the electron transport chain. Analysis of DNA copy-number alterations uncovered widespread loss of chromosomes culminating in near-haploid chromosomal content in a large fraction of HCC, which was maintained during metastatic spread. This work uncovers a distinct molecular origin of HCC compared with other thyroid malignancies.

BACKGROUND: Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas. METHODS: We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes. RESULTS: Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma. CONCLUSIONS: The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.)

Cancer cells can acquire profound alterations to the structure of their genomes, including rearrangements that fuse distant DNA breakpoints. We analyze the distribution of somatic rearrangements across the cancer genome, using whole-genome sequencing data from 2,693 tumor-normal pairs. We observe substantial variation in the density of rearrangement breakpoints, with enrichment in open chromatin and sites with high densities of repetitive elements. After accounting for these patterns, we identify significantly recurrent breakpoints (SRBs) at 52 loci, including novel SRBs near BRD4 and AKR1C3. Taking into account both loci fused by a rearrangement, we observe different signatures resembling either single breaks followed by strand invasion or two separate breaks that become joined. Accounting for these signatures, we identify 90 pairs of loci that are significantly recurrently juxtaposed (SRJs). SRJs are primarily tumor-type specific and tend to involve genes with tissue-specific expression. SRJs were frequently associated with disruption of topology-associated domains, juxtaposition of enhancer elements, and increased expression of neighboring genes. Lastly, we find that the power to detect SRJs decreases for short rearrangements, and that reliable detection of all driver SRJs will require whole-genome sequencing data from an order of magnitude more cancer samples than currently available.

Genomic analysis of tumours has led to the identification of hundreds of cancer genes on the basis of the presence of mutations in protein-coding regions. By contrast, much less is known about cancer-causing mutations in non-coding regions. Here we perform deep sequencing in 360 primary breast cancers and develop computational methods to identify significantly mutated promoters. Clear signals are found in the promoters of three genes. FOXA1, a known driver of hormone-receptor positive breast cancer, harbours a mutational hotspot in its promoter leading to overexpression through increased E2F binding. RMRP and NEAT1, two non-coding RNA genes, carry mutations that affect protein binding to their promoters and alter expression levels. Our study shows that promoter regions harbour recurrent mutations in cancer with functional consequences and that the mutations occur at similar frequencies as in coding regions. Power analyses indicate that more such regions remain to be discovered through deep sequencing of adequately sized cohorts of patients.

Small cell lung carcinoma (SCLC) is a highly lethal, smoking-associated cancer with few known targetable genetic alterations. Using genome sequencing, we characterized the somatic evolution of a genetically engineered mouse model (GEMM) of SCLC initiated by loss of Trp53 and Rb1. We identified alterations in DNA copy number and complex genomic rearrangements and demonstrated a low somatic point mutation frequency in the absence of tobacco mutagens. Alterations targeting the tumor suppressor Pten occurred in the majority of murine SCLC studied, and engineered Pten deletion accelerated murine SCLC and abrogated loss of Chr19 in Trp53; Rb1; Pten compound mutant tumors. Finally, we found evidence for polyclonal and sequential metastatic spread of murine SCLC by comparative sequencing of families of related primary tumors and metastases. We propose a temporal model of SCLC tumorigenesis with implications for human SCLC therapeutics and the nature of cancer-genome evolution in GEMMs, National Cancer Institute (U.S.) (K08CA160658)

Craniopharyngiomas are epithelial tumors that typically arise in the suprasellar region of the brain. Patients experience substantial clinical sequelae from both extension of the tumors and therapeutic interventions that damage the optic chiasm, the pituitary stalk and the hypothalamic area. Using whole-exome sequencing, we identified mutations in CTNNB1 (β-catenin) in nearly all adamantinomatous craniopharyngiomas examined (11/12, 92%) and recurrent mutations in BRAF (resulting in p.Val600Glu) in all papillary craniopharyngiomas (3/3, 100%). Targeted genotyping revealed BRAF p.Val600Glu in 95% of papillary craniopharyngiomas (36 of 39 tumors) and mutation of CTNNB1 in 96% of adamantinomatous craniopharyngiomas (51 of 53 tumors). The CTNNB1 and BRAF mutations were clonal in each tumor subtype, and we detected no other recurrent mutations or genomic aberrations in either subtype. Adamantinomatous and papillary craniopharyngiomas harbor mutations that are mutually exclusive and clonal. These findings have important implications for the diagnosis and treatment of these neoplasms.

Purpose Lung squamous cell carcinoma (SCC) is the second most prevalent type of lung cancer. Currently, no targeted therapeutics are approved for treatment of this cancer, largely because of a lack of systematic understanding of the molecular pathogenesis of the disease. To identify therapeutic targets and perform comparative analyses of lung SCC, we probed somatic genome alterations of lung SCC by using samples from Korean patients. Patients and Methods We performed whole-exome sequencing of DNA from 104 lung SCC samples from Korean patients and matched normal DNA. In addition, copy-number analysis and transcriptome analysis were conducted for a subset of these samples. Clinical association with cancer-specific somatic alterations was investigated. Results This cancer cohort is characterized by a high mutational burden with an average of 261 somatic exonic mutations per tumor and a mutational spectrum showing a signature of exposure to cigarette smoke. Seven genes demonstrated statistical enrichment for mutation: TP53, RB1, PTEN, NFE2L2, KEAP1, MLL2, and PIK3CA). Comparative analysis between Korean and North American lung SCC samples demonstrated a similar spectrum of alterations in these two populations in contrast to the differences seen in lung adenocarcinoma. We also uncovered recurrent occurrence of therapeutically actionable FGFR3-TACC3 fusion in lung SCC. Conclusion These findings provide new steps toward the identification of genomic target candidates for precision medicine in lung SCC, a disease with significant unmet medical needs.

In the version of this article originally published, an asterisk was omitted from Fig. 1a. The asterisk has been added to the figure. Additionally, a “NOTCH2” label was erroneously included in Fig. 4a. The label has been removed. The errors have been corrected in the PDF and HTML versions of this article.

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

Primary mediastinal large B-cell lymphomas (PMBL) typically occur in young women who present with localized, large mediastinal masses. These tumors share certain clinical, pathomorphological and transcriptional features with classical Hodgkin lymphoma (cHL). To date, PMBL genetic analyses focused on limited sets of genes and recurrent somatic copy number alterations (SCNAs). Previously, we identified frequent 9p24.1/PD-L1/PD-L2 copy gains and increased expression of the PD-1 ligands as a genetically-defined immune escape mechanism in PMBL. The demonstrated efficacy of PD-1 blockade in relapsed/refractory PMBL led to recent FDA approval and underscored the importance of characterizing targetable genetic vulnerabilities in this disease. For these reasons, we obtained diagnostic biopsy specimens from 37 patients with PMBL (median age 34; female 70%) and performed whole exome sequencing (WES) with an expanded bait set to capture structural variants (SVs). Somatic alterations (mutations, SCNAs and SVs) were determined using established analytical pipelines including our algorithm for evaluating tumors without paired normal samples. Genes more frequently mutated than by chance, Candidate Cancer Genes (CCGs), were identified with MutSig2CV and recurrent SCNAs were defined with GISTIC2.0. SVs were characterized with a recently described 4-algorithm pipeline (Nature Medicine, 2018;24(5):679-690). First, we identified 15 CCGs (q-value We next analyzed the PMBL mutational signatures and identified 3 cases as hypermutators with MSI signatures, including 2 with MLH1 frameshift mutations and 1 with a nonsense PMS2 mutation. Despite the young age of the PMBL patient cohort, the majority of remaining mutations were caused by spontaneous deamination at CpGs, a genetic signature usually associated with aging. The next most prevalent mutational signatures were APOBEC and, infrequently, AID. We observed a higher median mutational density in PMBL (7.56 mutations/MB), compared to diffuse large B-cell lymphoma (DLBCL) and most solid cancers, providing a potential basis for increased neoantigen production and responsiveness to PD-1 blockade. Next, we identified 18 recurrent SCNAs, including 10 copy gains (2 focal and 8 arm level) and 8 copy losses (7 focal and 1 arm level). Copy gains of 9p24.1/PD-L1/PD-L2 were detected in 70% of cases. SVs were defined at base-pair resolution and included infrequent (2/37) tandem duplications of both PD-1 ligands and inactivating CTIIA SVs (deletions and inversions) in 10% (4/37) of cases. Although PMBL had a higher mutational density than DLBCL, the PMBL alterations involved a smaller number of median genetic drivers (9 [PMBL] vs 17 [DLBCL], respectively). Combined analyses of recurrent CCGs, SCNAs and SVs revealed that certain candidate driver genes were perturbed by multiple mechanisms. Examples include: TNFAIP3 (59% overall, 41% mutations, 24% copy loss, 6% biallelic); and B2M (51% overall, 30% mutations, 27% copy loss, 6% biallelic). Concurrent analyses of the 3 types of genetic alterations also revealed multiple bases of perturbing specific signaling pathways. In this PMBL series, 73% (27/33) of tumors exhibited one or more alterations of JAK/STAT pathway components: IL4R mutations (32%), JAK2 (9p24.1 focal copy gain [70%]) and STAT6 mutations (43%). Additionally, 59% of PMBLs had alterations of antigen presentation pathway components including B2M copy loss or mutations, copy loss of 6q21.33 (which includes the HLA class I/II loci) and SVs of CTIIA. These findings provide a genetic framework for analyzing the precise mechanism of action of PD-1 blockade in PMBL. Taken together, these findings underscore the importance of a comprehensive genomic analysis in PMBL and define additional candidate treatment targets and pathogenetic mechanisms in this disease. ____ BC, CS and AD contributed equally. GG and MAS contributed equally. Disclosures Rodig: Merck: Research Funding; KITE: Research Funding; Affimed: Research Funding; Bristol Myers Squibb: Research Funding. Shipp:Merck: Research Funding; AstraZeneca: Honoraria; Bristol-Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Bayer: Research Funding.

Classical Hodgkin lymphoma (cHL) is composed of rare malignant Hodgkin Reed Sternberg (HRS) cells within an extensive, but ineffective, inflammatory/immune cell infiltrate. Emerging data suggests that cHLs use multiple genetic mechanisms to evade immune recognition. We previously found that HRS cells exhibit near-universal somatic copy number alterations (SCNAs) involving chromosome 9p24.1/PD-1-L1/PD-L2 and rare chromosomal rearrangements of PD-L1 or PD-L2. The 9p24.1 amplicon also includes JAK2, which increases JAK2 copy numbers, augments JAK2/STAT signaling and further induces PD-1 ligand expression. However, HRS cells also have inactivating mutations of B2M and decreased or absent MHC class I expression. In cHL, clinical responses to PD-1 blockade are unrelated to HRS cell expression of MHC class I but closely associated with HRS cell expression of MHC class II, highlighting the potential role of CD4+ T-cell effectors (J Clin Oncol 2018;36:942-50). To define genetic bases of response and resistance to PD-1 blockade and identify complementary treatment targets, we performed whole exome sequencing (WES) of HRS cells. We first used a previously described multi-color flow cytometric sorting protocol (Methods 2012; 57:368-75) to obtain highly purified CD30+ HRS cells and normal B cells from the excisional biopsies of 25 newly diagnosed cHLs. The isolated HRS cells and paired normal B cells were then subjected to WES using an optimized workflow for low input samples and an expanded bait set to capture structural variants (SVs). We used established analytical pipelines to identify significantly mutated genes (candidate cancer genes [CCGs], MutSig2CV), SCNAs (GISTIC2.0) and SVs (4 algorithms). With improved methodology and purity (median of 80%) of the isolated HRS cells, we defined 15 significantly mutated CCGs, 21 recurrent SCNAs, including 6 CN gains (4 focal and 2 arm level) and 15 CN losses (14 focal and 1 arm level), and low frequency SVs. We identified 2 cHLs as hypermutators with MSI signatures due to splice site mutations in MSH2 or missense mutations in POLE. Excluding the 2 hypermutators, the analyzed cHLs had a median mutational density of 6.4 mutations/Mb, that falls within the top quartile of reported cancer mutational frequencies (Nature 2013 499:214). We also identified a previously unappreciated high incidence of ARID1A mutations (24%) in cHL. This is noteworthy because ARID1A deficiency increases mutational load and augments the efficacy of PD-1 blockade in murine models (Nature Med 2018;24:556). Together, the observed MSI signatures, relatively high mutational burden and newly identified ARID1A mutations in cHL represent additional potential genetic bases for the efficacy of PD-1 blockade. Notably, these cHLs also exhibited recurrent 9p24.1 copy gain (80%) and multiple genetic bases of enhanced JAK/STAT signaling including JAK2 copy gain (80%), STAT6 mutations (32%) involving known hotspots (D419 and N421) in the DNA-binding domain and frequent inactivating SOCS1 mutations (68%). We also identified multiple genetic bases for immune evasion, including B2M inactivating mutations (36%), HLA-B mutations (16%) and 6p21.32/HLA-B copy loss (28%), copy loss of the larger 6p21.32 region and inactivating CIITA SVs (8%). Additional signaling pathways were perturbed by multiple genetic mechanisms in these cHLs. For example, NF-κB pathway alterations included: TNFAIP3 mutations (24%) and 6q23.2/TNFAIP3 copy loss (56%), 12% biallelic; NFKBIE mutations (24%) and 6q21.32/NFKBIE copy loss (12%); and NFKBIA mutations (16%). The gene encoding the nuclear export protein, XPO1, was perturbed by E571K mutations (24%) and frequent 2p15/XPO1 copy gain (72%). Additionally, GNA13, an activator of RHOA and modifier of PI3K signaling, was mutated in 24% of cases. Of interest, cHL recurrent alterations including B2M, TNFAIP3, STAT6, and GNA13 mutations and 6q23.2 and 9p24.1 SCNAs were also identified in > 20% of examined primary mediastinal B-cell lymphomas, highlighting shared pathogenetic mechanisms in these diseases. In summary, comprehensive genomic analyses of purified HRS cells reveal new genetic bases of immune evasion, potential mechanisms of response and resistance to PD-1 blockade and additional targetable alterations. KW, BC, CS, AD and DW contributed equally. JF, GG and MS contributed equally. Disclosures Rodig: Affimed: Research Funding; KITE: Research Funding; Merck: Research Funding; Bristol Myers Squibb: Research Funding. Shipp:Merck: Research Funding; Bayer: Research Funding; Bristol-Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; AstraZeneca: Honoraria.

Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.

Genetic changes underlying clear cell renal cell carcinoma (ccRCC) include alterations in genes controlling cellular oxygen sensing (for example, VHL) and the maintenance of chromatin states (for example, PBRM1). We surveyed more than 400 tumours using different genomic platforms and identified 19 significantly mutated genes. The PI(3)K/AKT pathway was recurrently mutated, suggesting this pathway as a potential therapeutic target. Widespread DNA hypomethylation was associated with mutation of the H3K36 methyltransferase SETD2, and integrative analysis suggested that mutations involving the SWI/SNF chromatin remodelling complex (PBRM1, ARID1A, SMARCA4) could have far-reaching effects on other pathways. Aggressive cancers demonstrated evidence of a metabolic shift, involving downregulation of genes involved in the TCA cycle, decreased AMPK and PTEN protein levels, upregulation of the pentose phosphate pathway and the glutamine transporter genes, increased acetyl-CoA carboxylase protein, and altered promoter methylation of miR-21 (also known as MIR21) and GRB10. Remodelling cellular metabolism thus constitutes a recurrent pattern in ccRCC that correlates with tumour stage and severity and offers new views on the opportunities for disease treatment.

The incidence of esophageal adenocarcinoma (EAC) has risen 600% over the last 30 years. With a 5-year survival rate of ~15%, the identification of new therapeutic targets for EAC is greatly important. We analyze the mutation spectra from whole-exome sequencing of 149 EAC tumor-normal pairs, 15 of which have also been subjected to whole-genome sequencing. We identify a mutational signature defined by a high prevalence of A>C transversions at AA dinucleotides. Statistical analysis of exome data identified 26 significantly mutated genes. Of these genes, five (TP53, CDKN2A, SMAD4, ARID1A and PIK3CA) have previously been implicated in EAC. The new significantly mutated genes include chromatin-modifying factors and candidate contributors SPG20, TLR4, ELMO1 and DOCK2. Functional analyses of EAC-derived mutations in ELMO1 identifies increased cellular invasion. Therefore, we suggest the potential activation of the RAC1 pathway as a contributor to EAC tumorigenesis., National Human Genome Research Institute (U.S.) (Large Scale Sequencing Program Grant U54 HG003067)

SummaryClonal evolution is a key feature of cancer progression and relapse. We studied intratumoral heterogeneity in 149 chronic lymphocytic leukemia (CLL) cases by integrating whole-exome sequence and copy number to measure the fraction of cancer cells harboring each somatic mutation. We identified driver mutations as predominantly clonal (e.g., MYD88, trisomy 12, and del(13q)) or subclonal (e.g., SF3B1 and TP53), corresponding to earlier and later events in CLL evolution. We sampled leukemia cells from 18 patients at two time points. Ten of twelve CLL cases treated with chemotherapy (but only one of six without treatment) underwent clonal evolution, predominantly involving subclones with driver mutations (e.g., SF3B1 and TP53) that expanded over time. Furthermore, presence of a subclonal driver mutation was an independent risk factor for rapid disease progression. Our study thus uncovers patterns of clonal evolution in CLL, providing insights into its stepwise transformation, and links the presence of subclones with adverse clinical outcomes.

Neuroblastoma is a malignancy of the developing sympathetic nervous system that often presents with widespread metastatic disease, resulting in survival rates of less than 50%. To determine the spectrum of somatic mutation in high-risk neuroblastoma, we studied 240 affected individuals (cases) using a combination of whole-exome, genome and transcriptome sequencing as part of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative. Here we report a low median exonic mutation frequency of 0.60 per Mb (0.48 nonsilent) and notably few recurrently mutated genes in these tumors. Genes with significant somatic mutation frequencies included ALK (9.2% of cases), PTPN11 (2.9%), ATRX (2.5%, and an additional 7.1% had focal deletions), MYCN (1.7%, causing a recurrent p.Pro44Leu alteration) and NRAS (0.83%). Rare, potentially pathogenic germline variants were significantly enriched in ALK, CHEK2, PINK1 and BARD1. The relative paucity of recurrent somatic mutations in neuroblastoma challenges current therapeutic strategies that rely on frequently altered oncogenic drivers., National Human Genome Research Institute (U.S.) (Grant U54HG003067), National Cancer Institute (U.S.) (Contract HHSN261200800001E)

We performed genomic, epigenomic, transcriptomic, and proteomic characterizations of uterine carcinosarcomas (UCSs). Cohort samples had extensive copy-number alterations and highly recurrent somatic mutations. Frequent mutations were found in TP53, PTEN, PIK3CA, PPP2R1A, FBXW7, and KRAS, similar to endometrioid and serous uterine carcinomas. Transcriptome sequencing identified a strong epithelial-to-mesenchymal transition (EMT) gene signature in a subset of cases that was attributable to epigenetic alterations at microRNA promoters. The range of EMT scores in UCS was the largest among all tumor types studied via The Cancer Genome Atlas. UCSs shared proteomic features with gynecologic carcinomas and sarcomas with intermediate EMT features. Multiple somatic mutations and copy-number alterations in genes that are therapeutic targets were identified.

Resistance to the Bruton’s tyrosine kinase (BTK) inhibitor ibrutinib has been attributed solely to mutations in BTK and related pathway molecules. Using whole-exome and deep-targeted sequencing, we dissect evolution of ibrutinib resistance in serial samples from five chronic lymphocytic leukaemia patients. In two patients, we detect BTK-C481S mutation or multiple PLCG2 mutations. The other three patients exhibit an expansion of clones harbouring del(8p) with additional driver mutations (EP300, MLL2 and EIF2A), with one patient developing trans-differentiation into CD19-negative histiocytic sarcoma. Using droplet-microfluidic technology and growth kinetic analyses, we demonstrate the presence of ibrutinib-resistant subclones and estimate subclone size before treatment initiation. Haploinsufficiency of TRAIL-R, a consequence of del(8p), results in TRAIL insensitivity, which may contribute to ibrutinib resistance. These findings demonstrate that the ibrutinib therapy favours selection and expansion of rare subclones already present before ibrutinib treatment, and provide insight into the heterogeneity of genetic changes associated with ibrutinib resistance., University of Texas M.D. Anderson Cancer Center (Support Grant CA016672), National Science Foundation (U.S.) (DMR-1310266), Harvard University. Materials Research Science and Engineering Center (DMR-1420570), National Natural Science Foundation (China) (81372496)

Recurrence rates after breast-conserving therapy may depend on genomic characteristics of cancer-adjacent, benign-appearing tissue. Studies have not evaluated recurrence in association with multiple genomic characteristics of cancer-adjacent breast tissue. To estimate the prevalence of DNA defects and RNA expression subtypes in cancer-adjacent, benign-appearing breast tissue at least 2 cm from the tumor margin, cancer-adjacent, pathologically well-characterized, benign-appearing breast tissue specimens from The Cancer Genome Atlas project were analyzed for DNA sequence, copy-number variation, DNA methylation, messenger RNA (mRNA) sequence, and mRNA/microRNA expression. Additional samples were also analyzed by at least one of these genomic data types and associations between genomic characteristics of normal tissue and overall survival were assessed. Approximately 40% of cancer-adjacent, benign-appearing tissues harbored genomic defects in DNA copy number, sequence, methylation, or in RNA sequence, although these defects did not significantly predict 10-year overall survival. Two mRNA/microRNA expression phenotypes were observed, including an active mRNA subtype that was identified in 40% of samples. Controlling for tumor characteristics and the presence of genomic defects, this active subtype was associated with significantly worse 10-year survival among estrogen receptor (ER)-positive cases. This multi-platform analysis of breast cancer-adjacent samples produced genomic findings consistent with current surgical margin guidelines, and provides evidence that extratumoral RNA expression patterns in cancer-adjacent tissue predict overall survival among patients with ER-positive disease. RNA expression patterns in the breast tissue adjacent to where tumors are surgically removed are linked to patient survival outcomes. A research team led by Melissa Troester from the University of North Carolina at Chapel Hill, USA, looked for genomic alterations in the healthy, benign-looking breast tissue found more than 2 centimeters from the margin of where cancer had been cut out and treated with radiation. Approximately 40% of the 142 tissue samples analyzed harbored defects in DNA copy number, sequence, methylation, or RNA sequence. However, the only characteristic significantly associated with 10-year survival rates in patients was a particular expression profile of messenger RNA known as “active” and this connection held up only in women with a subtype of breast cancer called ER-positive. In this population, RNA profiling could help predict survival outcomes.

There is substantial heterogeneity among primary prostate cancers, evident in the spectrum of molecular abnormalities and its variable clinical course. As part of The Cancer Genome Atlas (TCGA), we present a comprehensive molecular analysis of 333 primary prostate carcinomas. Our results revealed a molecular taxonomy in which 74% of these tumors fell into one of seven subtypes defined by specific gene fusions (ERG, ETV1/4, and FLI1) or mutations (SPOP, FOXA1, and IDH1). Epigenetic profiles showed substantial heterogeneity, including an IDH1 mutant subset with a methylator phenotype. Androgen receptor (AR) activity varied widely and in a subtype-specific manner, with SPOP and FOXA1 mutant tumors having the highest levels of AR-induced transcripts. 25% of the prostate cancers had a presumed actionable lesion in the PI3K or MAPK signaling pathways, and DNA repair genes were inactivated in 19%. Our analysis reveals molecular heterogeneity among primary prostate cancers, as well as potentially actionable molecular defects., National Institutes of Health (U.S.). (grant 5U24CA143799), National Institutes of Health (U.S.). (grant 5U24CA143835), National Institutes of Health (U.S.). (grant 5U24CA143840), National Institutes of Health (U.S.). (grant 5U24CA143843), National Institutes of Health (U.S.). (grant 5U24CA143845), National Institutes of Health (U.S.). (grant 5U24CA143848), National Institutes of Health (U.S.). (grant 5U24CA143858), National Institutes of Health (U.S.). (grant 5U24CA143866), National Institutes of Health (U.S.). (grant 5U24CA143867), National Institutes of Health (U.S.). (grant 5U24CA143882), National Institutes of Health (U.S.). (grant 5U24CA143883), National Institutes of Health (U.S.). (grant 5U24CA144025), National Institutes of Health (U.S.). (grant U54HG003067), National Institutes of Health (U.S.). (grant U54HG003079), National Institutes of Health (U.S.). (grant U54HG003273), National Institutes of Health (U.S.). (grant P30CA16672)

In the version of this article originally published, some text above the “Tri–nucleotide sequence motifs” label in Fig. 2a appeared incorrectly. The text was garbled and should have appeared as nucleotide codes. Additionally, the labels on the bars in Fig. 2c were not italicized in the original publication. These are gene symbols, and they should have been italicized. The colored labels above the graphs in Fig. 4b were also erroneously not italicized. These labels represent gene names and loci, and they should have been italicized.

Barrett's esophagus is thought to progress to esophageal adenocarcinoma (EAC) through a stepwise progression with loss of CDKN2A followed by TP53 inactivation and aneuploidy. Here we present whole-exome sequencing from 25 pairs of EAC and Barrett's esophagus and from 5 patients whose Barrett's esophagus and tumor were extensively sampled. Our analysis showed that oncogene amplification typically occurred as a late event and that TP53 mutations often occurred early in Barrett's esophagus progression, including in non-dysplastic epithelium. Reanalysis of additional EAC exome data showed that the majority (62.5%) of EACs emerged following genome doubling and that tumors with genomic doubling had different patterns of genomic alterations, with more frequent oncogenic amplification and less frequent inactivation of tumor suppressors, including CDKN2A. These data suggest that many EACs emerge not through the gradual accumulation of tumor-suppressor alterations but rather through a more direct path whereby a TP53-mutant cell undergoes genome doubling, followed by the acquisition of oncogenic amplifications., National Human Genome Research Institute (U.S.) Large-Scale Sequencing Program (Grant U54 HG0003067)

Which genetic alterations drive tumorigenesis and how they evolve over the course of disease and therapy are central questions in cancer biology. Here we identify 44 recurrently mutated genes and 11 recurrent somatic copy number variations through whole-exome sequencing of 538 chronic lymphocytic leukaemia (CLL) and matched germline DNA samples, 278 of which were collected in a prospective clinical trial. These include previously unrecognized putative cancer drivers (RPS15, IKZF3), and collectively identify RNA processing and export, MYC activity, and MAPK signalling as central pathways involved in CLL. Clonality analysis of this large data set further enabled reconstruction of temporal relationships between driver events. Direct comparison between matched pre-treatment and relapse samples from 59 patients demonstrated highly frequent clonal evolution. Thus, large sequencing data sets of clinically informative samples enable the discovery of novel genes associated with cancer, the network of relationships between the driver events, and their impact on disease relapse and clinical outcome., National Human Genome Research Institute (U.S.) (Grant U54HG003067)

Lymphoma is the most common hematological malignancy in developed countries. Outcome is strongly determined by molecular subtype, reflecting a need for new and improved treatment options. Dogs spontaneously develop lymphoma, and the predisposition of certain breeds indicates genetic risk factors. Using the dog breed structure, we selected three lymphoma predisposed breeds developing primarily T-cell (boxer), primarily B-cell (cocker spaniel), and with equal distribution of B- and T-cell lymphoma (golden retriever), respectively. We investigated the somatic mutations in B- and T-cell lymphomas from these breeds by exome sequencing of tumor and normal pairs. Strong similarities were evident between B-cell lymphomas from golden retrievers and cocker spaniels, with recurrent mutations in TRAF3-MAP3K14 (28% of all cases), FBXW7 (25%), and POT1 (17%). The FBXW7 mutations recurrently occur in a specific codon; the corresponding codon is recurrently mutated in human cancer. In contrast, T-cell lymphomas from the predisposed breeds, boxers and golden retrievers, show little overlap in their mutation pattern, sharing only one of their 15 most recurrently mutated genes. Boxers, which develop aggressive T-cell lymphomas, are typically mutated in the PTEN-mTOR pathway. T-cell lymphomas in golden retrievers are often less aggressive, and their tumors typically showed mutations in genes involved in cellular metabolism. We identify genes with known involvement in human lymphoma and leukemia, genes implicated in other human cancers, as well as novel genes that could allow new therapeutic options.

SummaryIn Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the “living biobank” agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design.PaperClip

Summary We performed integrated genomic, transcriptomic, and proteomic profiling of 150 pancreatic ductal adenocarcinoma (PDAC) specimens, including samples with characteristic low neoplastic cellularity. Deep whole-exome sequencing revealed recurrent somatic mutations in KRAS , TP5 3, CDKN2A , SMAD4 , RNF43 , ARID1A , TGFβR2 , GNAS , RREB1 , and PBRM1 . KRAS wild-type tumors harbored alterations in other oncogenic drivers, including GNAS , BRAF , CTNNB1 , and additional RAS pathway genes. A subset of tumors harbored multiple KRAS mutations, with some showing evidence of biallelic mutations. Protein profiling identified a favorable prognosis subset with low epithelial-mesenchymal transition and high MTOR pathway scores. Associations of non-coding RNAs with tumor-specific mRNA subtypes were also identified. Our integrated multi-platform analysis reveals a complex molecular landscape of PDAC and provides a roadmap for precision medicine.