Search

Tumor behavior is intricately dependent on the oncogenic properties of cancer cells and their multi-cellular interactions. To understand these dependencies within the wider microenvironment, we studied over 270,000 single-cell transcriptomes and 100 microdissected whole exomes from 12 patients with kidney tumors, prior to validation using spatial transcriptomics. Tissues were sampled from multiple regions of the tumor core, the tumor-normal interface, normal surrounding tissues, and peripheral blood. We find that the tissue-type location of CD8+ T cell clonotypes largely defines their exhaustion state with intra-tumoral spatial heterogeneity that is not well explained by somatic heterogeneity. De novo mutation calling from single-cell RNA-sequencing data allows us to broadly infer the clonality of stromal cells and lineage-trace myeloid cell development. We report six conserved meta-programs that distinguish tumor cell function, and find an epithelial-mesenchymal transition meta-program highly enriched at the tumor-normal interface that co-localizes with IL1B-expressing macrophages, offering a potential therapeutic target.

Funder: E.P. is a Josie Robertson Investigator and is supported by the European Hematology Association, American Society of Hematology, Gabrielle’s Angels Foundation, V Foundation and The Geoffrey Beene Foundation and is a Damon Runyon Rachleff Innovator fellow. Work in the BJPH lab is funded by Cancer Research UK (C18680/A25508), the European Research Council (647685), MRC (MR-R009708-1), the Kay Kendall Leukaemia Fund (KKL1243), the Wellcome Trust (205254/Z/16/Z) and the Cancer Research UK Cambridge Major Centre (C49940/A25117). This research was supported by the NIHR Cambridge Biomedical Research Centre (BRC-1215-20014), and was funded in part, by the Wellcome Trust who supported the Wellcome - MRC Cambridge Stem Cell Institute (203151/Z/16/Z). The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care. L.B., H.D. and B.J.P.H. are supported by the HARMONY Alliance (IMI Project No. 116026; https://www.harmony-alliance.eu/). The UK-NCRI AML working group trials were supported with research grants from the Medical Research Council (MRC), Cancer Research UK (CRUK), Blood Cancer UK and Cardiff University. We would like to thank all patients and investigators for their participation in the trials and the study., Clinical recommendations for Acute Myeloid Leukemia (AML) classification and risk-stratification remain heavily reliant on cytogenetic findings at diagnosis, which are present in

The diversity in somatic variation across human tissues outshining human cancers has only recently been started to be appreciated. In particular, healthy individuals carry large numbers of cells and clonal expansions with mutations in driver genes commonly observed in cancer. We previously demonstrated that the majority of epithelial cells in the proximal airway accumulate DNA damage proportionally to an individual's smoking habits. However, the analysis of whole genomes indicated the abundance of a population of cells with a near-normal mutation burden, mirroring cellular populations found in never-smoking individuals. In contrast to the proximal airway, the genomics of cells in the distal airway of healthy individuals with diverse smoking histories, and thus the foundation of lung adenocarcinoma (LUAD), has not been investigated yet. To assess the influence of smoking on the distal airway and to compare the genomics of proximal and distal airway in the lung, we established an experimental procedure to create single-cell derived alveolar type II (AT2) organoids from primary human tissue of 9 patients. In particular, we focus on AT2 cells, given their hypothesised role as a potential cell of origin for LUAD. Whole genome sequencing (WGS) of >400 organoids was utilised to assess the landscape of somatic variation, including single nucleotide variants (SNVs), insertions and deletions (indels), copy number and structural variation (CNVs, SVs). Focusing on SNVs, we found that the mutation burden of never-smokers is significantly lower compared to ex- and current-smokers. Emphasising the effect of tobacco smoking on the distal airway, this difference was driven by the presence of mutations related to SBS4 and SBS92, mutational signatures associated with tobacco smoking. Interestingly, and in contrast to the proximal airway, we did not find a difference in the mutation burden of ex- and current-smokers, suggesting that smoking damage to alveolar cells is more severe and long-term. Lastly, to characterise the driver mutation landscape in depth, we used targeted single-molecule sequencing of 250 known cancer genes to profile 15 AT2 single-cell suspensions as well as 22 frozen parenchyma biopsies. To date, our preliminary results indicate positive selection on lineage defining mutations including SFTPB and SFTPC, coding for surfactant as well as TP53. As we begin to explore these data further, we expect this work to provide unprecedented insights into the contribution of somatic and driver mutations to the earliest stages of lung cancer development. Citation Format: Moritz Jakob Przybilla, Amany Ammar, Ryan Chuen Khaw, Andrew R. Lawson, Pantelis Nicola, Kate Davies, Zoe Frazer, Kate H. Gowers, Timothy M. Butler, Sarah E. Clarke, Inigo Martincorena, Sam M. Janes, Peter J. Campbell. Somatic mutations in single-cell derived alveolar organoids [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 218.

Tumour behaviour is dependent on the oncogenic properties of cancer cells and their multi-cellular interactions. These dependencies were examined through 270,000 single cell transcriptomes and 100 micro-dissected whole exomes obtained from 12 patients with kidney tumours. Tissue was sampled from multiple regions of tumour core, tumour-normal interface, normal surrounding tissues, and peripheral blood. We found the principal spatial location of CD8+ T cell clonotypes largely defined exhaustion state, with clonotypic heterogeneity not explained by somatic intra-tumoural heterogeneity. De novo mutation calling from single cell RNA sequencing data allows us to lineage-trace and infer clonality of cells. We discovered six meta-programmes that distinguish tumour cell function. An epithelial-mesenchymal transition meta-programme, enriched at the tumour-normal interface appears modulated through macrophage expressed IL1B, potentially forming a therapeutic target.

Most cancers are characterized by the somatic acquisition of genomic rearrangements during tumour evolution that eventually drive the oncogenesis. Here, using multiplatform sequencing technologies, we identify and characterize a remarkable mutational mechanism in human hepatocellular carcinoma caused by Hepatitis B virus, by which DNA molecules from the virus are inserted into the tumour genome causing dramatic changes in its configuration, including non-homologous chromosomal fusions, dicentric chromosomes and megabase-size telomeric deletions. This aberrant mutational mechanism, present in at least 8% of all HCC tumours, can provide the driver rearrangements that a cancer clone requires to survive and grow, including loss of relevant tumour suppressor genes. Most of these events are clonal and occur early during liver cancer evolution. Real-time timing estimation reveals some HBV-mediated rearrangements occur as early as two decades before cancer diagnosis. Overall, these data underscore the importance of characterising liver cancer genomes for patterns of HBV integration., Hepatitis B virus (HBV) infection and DNA integration is a frequent cause of hepatocellular carcinoma (HCC), but the consequences of this process are not fully understood. Here the authors use whole-genome and long-read sequencing data from HCC patient samples to study the timing and alterations induced by HBV insertions.

Fanconi anemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink (ICL) repair resulting in chromosome breakage1–3. The FA repair pathway comprises at least 22 FANC proteins including BRCA1 and BRCA24–6, and protects against carcinogenic endogenous and exogenous aldehydes7–10. Individuals with FA are hundreds to thousands-fold more likely to develop head and neck (HNSCC), esophageal and anogenital squamous cell carcinomas (SCCs) with a median onset age of 31 years11. The aggressive nature of these tumors and poor patient tolerance of platinum and radiation-based therapy have been associated with short survival in FA11–16. Molecular studies of SCCs from individuals with FA (FA SCCs) have been limited, and it is unclear how they relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or human papillomavirus (HPV) infection17. Here, by sequencing FA SCCs, we demonstrate that the primary genomic signature of FA-deficiency is the presence of a high number of structural variants (SVs). SVs are enriched for small deletions, unbalanced translocations, and fold-back inversions that arise in the context of TP53 loss. The SV breakpoints preferentially localize to early replicating regions, common fragile sites, tandem repeats, and SINE elements. SVs are often connected forming complex rearrangements. Resultant genomic instability underlies elevated copy number alteration (CNA) rates of key HNSCC-associated genes, including PIK3CA, MYC, CSMD1, PTPRD, YAP1, MXD4, and EGFR. In contrast to sporadic HNSCC, we find no evidence of HPV infection in FA HNSCC, although positive cases were identified in gynecologic tumors. A murine allograft model of FA pathway-deficient SCC was enriched in SVs, exhibited dramatic tumor growth advantage, more rapid epithelial-to-mesenchymal transition (EMT), and enhanced autonomous inflammatory signaling when compared to an FA pathway-proficient model. In light of the protective role of the FA pathway against SV formation uncovered here, and recent findings of FA pathway insufficiency in the setting of increased formaldehyde load resulting in hematopoietic stem cell failure and carcinogenesis18–20, we propose that high copy-number instability in sporadic HNSCC may result from functional overload of the FA pathway by endogenous and exogenous DNA crosslinking agents. Our work lays the foundation for improved FA patient treatment and demonstrates that FA SCC is a powerful model to study tumorigenesis resulting from DNA crosslinking damage.

The rate of nonalcoholic fatty liver disease (NAFLD)-associated hepatocellular carcinoma (HCC) is increasing worldwide, but the steps in precancerous hepatocytes which lead to HCC driver mutations are not well understood. Here we provide evidence that metabolically-driven histone hyperacetylation in steatotic hepatocytes can increase DNA damage to initiate carcinogenesis. Genome-wide histone acetylation is increased in steatotic livers of rodents fed high fructose or high fat diet. In vitro, steatosis relaxes chromatin and increases DNA damage marker γH2AX, which is reversed by inhibiting acetyl-CoA production. Steatosis-associated acetylation and γH2AX are enriched at gene clusters in telomere-proximal regions which contain HCC tumor suppressors in hepatocytes and human fatty livers. Regions of metabolically-driven epigenetic change also have increased levels of DNA mutation in non-cancerous tissue from NAFLD patients. Finally, genome-scale network modelling indicates that redox balance is a key contributor to this mechanism. Thus abnormal histone hyperacetylation is a potential initiating event in HCC carcinogenesis.

Most cancers are characterized by the somatic acquisition of genomic rearrangements during tumour evolution that eventually drive the oncogenesis. There are different mutational mechanisms causing structural variation, some of which are specific to particular cancer types. Here, using multiplatform sequencing technologies, we identify and characterize a remarkable mutational mechanism in human hepatocellular carcinoma caused by Hepatitis B virus, by which DNA molecules from the virus are inserted into the tumour genome causing dramatic changes in its configuration, including non-homologous chromosomal fusions and megabase-size telomeric deletions. This aberrant mutational process, present in at least 8% of all HCC tumours, is active early during liver cancer evolution and can provide the driver rearrangements that a cancer clone requires to survive and grow.

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n = 15) are characterized by later initiation in the patient’s life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable., The factors that are associated with myeloma precursor condition progression are not well understood. Here the authors find that monoclonal gammopathies of undetermined significance and smoldering myelomas that did not progress to multiple myelomas have a distinct genomic profile and emerge later in the patient’s life.

Introduction: Clonal evolution drives cancer development due to the emergence and/or selection of proliferatively advantageous subclones. Its understanding may facilitate the design of anticipation-based management strategies. Richter transformation (RT) is a paradigmatic tumor evolution in which chronic lymphocytic leukemia (CLL), an indolent neoplasia of mature B-cells, transforms into a high-grade lymphoma, usually diffuse large B-cell lymphoma (DLBCL), conferring a dismal prognosis. The evolutionary trajectories of RT and its driving (epi)genomic mechanisms remain largely unknown. Aims: To reconstruct the evolutionary history of RT and to reveal the molecular processes underlying this transformation. Methods: We characterized the whole genome (WGS), epigenome (DNA methylation, H3K27ac, ATAC-seq), and transcriptome (RNA-seq), combined with single-cell DNA and RNA sequencing analyses, of 19 CLL patients developing RT before (n=3) or after treatment with chemoimmunotherapy (n=6) and targeted therapies (BCR or BCL2 inhibitors, n=10). We analyzed 54 longitudinal samples covering up to 19 years of disease course. Results: Our WGS analyses uncovered that RT is characterized by a remarkable structural complexity. We also identified a novel treatment-independent RT-specific mutational process, which we named SBS-RT. The genetic driver landscape of RT is a compendium of alterations in genes involved in cell cycle, MYC, and NF-κB pathways, frequently targeted in single catastrophic events including chromothripsis and chromoplexy. The WGS-based phylogenic reconstruction and single-cell DNA/RNA-seq analyses identified a very early diversification of CLL leading to emergence of RT-cells carrying specific genetic drivers and transcriptomic profiles of RT already at CLL diagnosis. These small subclones were dormant for 6-19 years until rapid expansion associated with the clinical transformation. While the DNA methylome kept track of the cell of origin and proliferative history of RT cells, their chromatin configuration and transcriptional program converged into the overexpression of cell cycle regulators, Toll-like receptors, MYC, MTORC1, and OXPHOS related transcripts, as well as downregulation of BCR pathway. This phenotypic shift was related to de novo activation of key transcription factors. In vitro experiments confirmed that RT cells have a 4-fold higher oxygen consumption at routine respiration and electron transfer system capacity compared to CLL. The resistance of RT to BCR inhibition is consistent with its high OXPHOS and low BCR signaling, which mimics de novo DLBCL-OXPHOS insensitive to BCR inhibition. This OXPHOShigh-BCRlow transcriptional axis of RT can be exploited therapeutically. Conclusions: These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover therapeutic targets for the, once expanded, lethal Richter transformation. Citation Format: Ferran Nadeu, Romina Royo, Ramon Massoni-Badosa, Beatriz Garcia-Torre, Martí Duran-Ferrer, Kevin J. Dawson, Marta Kulis, Ander Diaz-Navarro, Neus Villamor, Juan L. Melero, Vicente Chapaprieta, Ana Dueso-Barroso, Julio Delgado, Riccardo Moia, Sara Ruiz-Gil, Domenica Marchese, Núria Verdaguer-Dot, Mónica Romo, Maria Rozman, Gerard Frigola, Alfredo Rivas-Delgado, Tycho Baumann, Miguel Alcoceba, Marcos González, Fina Climent, Pau Abrisqueta, Josep Castellví, Francesc Bosch, Marta Aymerich, Anna Enjuanes, Sílvia Ruiz-Gaspà, Armando López-Guillermo, Pedro Jares, Sílvia Beà, Dolors Colomer, Núria López-Bigas, Josep LlGelpí, David Torrents, Peter J. Campbell, Ivo Gut, Pablo M. Garcia-Roves, Davide Rossi, Gianluca Gaidano, Xose S. Puente, Holger Heyn, Francesco Maura, José I. Martín-Subero, Elías Campo. Early seeding of Richter transformation in chronic lymphocytic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3795.

Fanconi anemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink (ICL) repair resulting in chromosome breakage. The FA repair pathway comprises at least 22 FANC proteins including BRCA1 and BRCA2 and protects against carcinogenic endogenous and exogenous aldehydes. Individuals with FA are hundreds to thousands-fold more likely to develop head and neck (HNSCC), esophageal and anogenital squamous cell carcinomas (SCCs) with a median onset age of 31 years. The aggressive nature of these tumors and poor patient tolerance of platinum and radiation-based therapy have been associated with short survival in FA. Molecular studies of SCCs from individuals with FA (FA SCCs) have been limited, and it is unclear how they relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or human papillomavirus (HPV) infection. Here, by sequencing FA SCCs, we demonstrate that the primary genomic signature of FA-deficiency is the presence of a high number of structural variants (SVs). SVs are enriched for small deletions, unbalanced translocations, and fold-back inversions that arise in the context of TP53 loss. The SV breakpoints preferentially localize to early replicating regions, common fragile sites, tandem repeats, and SINE elements. SVs are often connected forming complex rearrangements. Resultant genomic instability underlies elevated copy number alteration (CNA) rates of key HNSCC-associated genes, including PIK3CA, MYC, CSMD1, PTPRD, YAP1, MXD4, and EGFR. In contrast to sporadic HNSCC, we find no evidence of HPV infection in FA HNSCC, although positive cases were identified in gynecologic tumors. A murine allograft model of FA pathway-deficient SCC was enriched in SVs, exhibited dramatic tumor growth advantage, more rapid epithelial-to-mesenchymal transition, and enhanced autonomous inflammatory signaling when compared to an FA pathway-proficient model. In light of the protective role of the FA pathway against SV formation uncovered here, and recent findings of FA pathway insufficiency in the setting of increased formaldehyde load resulting in hematopoietic stem cell failure and carcinogenesis, we propose that high copy-number instability in sporadic HNSCC may result from functional overload of the FA pathway by endogenous and exogenous DNA crosslinking agents. Our work lays the foundation for improved FA patient treatment and demonstrates that FA SCC is a powerful model to study tumorigenesis resulting from DNA crosslinking damage. Citation Format: Andrew L. Webster, Mathijs A. Sanders, Krupa Patel, Ralf Dietrich, Raymond J. Noonan, Francis P. Lach, Ryan R. White, Audrey M. Goldfarb, Kevin Hadi, Matthew M. Edwards, Frank X. Donovan, Moonjung Jung, Sunandini Sridhar, Olivier Fedrigo, Huasong Tian, Joel Rosiene, Thomas Heineman, Jennifer Kennedy, Lorenzo Bean, Rasim O. Rosti, Rebecca Tryon, Ashlyn-Maree Gonzalez, Allana Rosenberg, Ji-Dung Luo, Thomas Carrol, Eunike Velleuer, Jeff C. Rastatter, Susanne I. Wells, Jordi Surrallés, Grover Bagby, Margaret L. MacMillan, John E. Wagner, Maria Cancio, Farid Boulad, Theresa Scognamiglio, Roger Vaughan, Amnon Koren, Marcin Imielinski, Settara Chandrasekharappa, Arleen D. Auerbach, Bhuvanesh Singh, David Kutler, Peter J. Campbell, Agata Smogorzewska. Fanconi anemia pathway deficiency drives copy number variation in squamous cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6196.

Introduction: Airway basal cells isolated from histologically normal epithelium in current and former smokers demonstrate a significantly higher mutation burden than those in never-smokers. However, a proportion of basal cells with lower (“near normal”) mutation burden are observed, which show preferential expansion in response to smoking cessation(1). We believe this expansion represents a “cancer protective” process, reducing the risk of future malignancy. However, mechanisms determining the dynamic balance of high versus low mutant basal cells throughout the bronchial tree are unclear. We have developed protocols that enable systematic scanning of airway epithelium in resected lung specimens, through which we aim to elucidate the mutation burden and spatial clonal architecture across contiguous regions of normal epithelium through the bronchial tree. Study design: 5-10 current or former smokers undergoing curative lobectomy for peripheral lung cancer with histologically normal airways will be consented to the study. To date, two participants have donated their samples. Fresh lobectomies were collected at the point of resection and kept on wet ice until pathological examination. All samples undergo examination within 1 hour of resection. The lobe is orientated, and the tumour and segmental bronchi identified. A sterile metal cannula is inserted through the airway lumen and sterile pathology scissors are used to carefully isolate the airways from the surrounding parenchyma. Upon macro-dissection, airway branches are immediately placed into transport media containing antibiotics. Samples are gently washed using transport media to remove residual blood, then sectioned into 1 cm3 segments, which are orientated, embedded and frozen in room-temperature Optimal Cutting Temperature compound. An 8 μM section is cut and stained with haematoxylin and eosin from each block immediately after a “full-face” tissue confirmation is observed.These will be used as guiding slides for laser capture microdissection of the regions of epithelium prior to downstream sequencing pathway. Conclusion: We have developed a pathological approach in which we can define somatic mutations in normal epithelium along the bronchial tree. We will use these data to determine basal cell clonal relationships and interpret these in a spatial context, creating a more detailed map of the smoking-induced changes that occur in the airway epithelium and how these may contribute to early carcinogenesis. References1.Yoshida K, Gowers KHC, Lee-Six H, Chandrasekharan DP, Coorens T, Maughan EF, et al. Tobacco smoking and somatic mutations in human bronchial epithelium. Nature. 2020;578(7794):266-72. Citation Format: Kate E. Davies, Helen Hall, Zoe Frazer, Kate Gowers, David Moore, Fleur Monk, Moritz J. Przybilla, Peter J. Campbell, Sam M. Janes. Contiguous airway sampling for basal cell isolation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2202.

It is well established that tobacco smoking causes lung cancer and this is thought to be the result of a progressive accumulation of genomic mutations induced by tobacco smoke exposure; however, until recently a detailed understanding of the effects of smoke exposure on the genomes of cells in the normal bronchial epithelium was lacking. To address this, we defined the genomic architecture of normal airway basal cells in children and non-smoking, smoking and ex-smoking adults. We discovered two populations of basal cells in individuals with a smoking history: those with high mutational burden, driver gene mutations and tobacco-associated mutational signatures, and those with low mutational burden, similar to basal cells from never smokers and with an absence of mutational signatures associated with tobacco exposure. These cells expand preferentially in ex-smokers, accounting for 20-40% of basal cells compared with 2-5% in current smokers. We hypothesise that these cells are cancer-protective and their expansion reduces the risk of lung cancer development after smoking cessation. Understanding the functional differences between high- and low-mutant airway basal cells and the mechanisms by which some cells resist the accumulation of smoking-induced mutations will be crucial to understanding the earliest stages of lung carcinogenesis. To address this, we performed RNAseq on a subset of whole-genome sequenced airway basal cells. Preliminary data show that expression of markers of stemness is significantly different between high- and low-mutant basal cells. In addition, initial analysis shows that pathways such as carcinogen metabolism and MHC class I antigen presentation are higher in high-mutant basal cells. Further analysis in more patients is ongoing. To complement these data, we expanded whole-genome sequenced basal cell clones in culture and performed a range of assays to assess their progenitor and differentiation capacity. We assessed proliferation, colony-forming efficiency, longevity and differentiation and found no apparent differences between low- and high-mutant basal cells. Initial analysis suggests that carcinogen metabolism and MHC class I antigen presentation may be key pathways in establishing heterogeneity in mutational burden and clonal dynamics in the normal airway epithelium. Additional assays and analysis are ongoing and will be the focus of future research. Citation Format: Kate H. Gowers, Sarah E. Clarke, Ayu Hutami Syaraf, Kenichi Yoshida, Moritz J. Przybilla, Hugh Selway, Adam Pennycuick, Peter J. Campbell, Sam M. Janes. Defining the mechanisms that lead to mutational heterogeneity in the normal respiratory epithelium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3158.

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate, for the first time, the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n=15) are characterized by later initiation in the patient’s life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes(1,2). The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution(1,3,4). Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such asBCL9, MCL1,ARNT(also known asHIF1B),TERTandMYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompassesBCL9, MCL1andARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassingMYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassingCCND1) in HER2(+)breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.Chromosomal instability enables the continuous selection of somatic copy number alterations, which are established as ordered events that often occur in parallel, throughout tumour evolution and metastasis.

Clonal dynamics and mutation burden in healthy human prostate epithelium are relevant to prostate cancer. We sequenced whole genomes from 409 microdissections of normal prostate epithelium across 8 donors, using phylogenetic reconstruction with spatial mapping in a 59-year-old man's prostate to reconstruct tissue dynamics across the lifespan. Somatic mutations accumulate steadily at ∼16 mutations/year/clone, with higher rates in peripheral than peri-urethral regions. The 24-30 independent glandular subunits are established as rudimentary ductal structures during fetal development by 5-10 embryonic cells each. Puberty induces formation of further side and terminal branches by local stem cells disseminated throughout the rudimentary ducts during development. During adult tissue maintenance, clonal expansions have limited geographic scope and minimal migration. Driver mutations are rare in aging prostate epithelium, but the one driver we did observe generated a sizable intraepithelial clonal expansion. Leveraging unbiased clock-like mutations, we define prostate stem cell dynamics through fetal development, puberty, and aging. ispartof: CELL STEM CELL vol:28 issue:7 pages:1262-+ ispartof: location:United States status: published

The rare benign giant cell tumour of bone (GCTB) is defined by an almost unique mutation in the H3.3 family of histone genes H3-3A or H3-3B; however, the same mutation is occasionally found in primary malignant bone tumours which share many features with the benign variant. Moreover, lung metastases can occur despite the absence of malignant histological features in either the primary or metastatic lesions. Herein we investigated the genetic events of 17 GCTBs including benign and malignant variants and the methylation profiles of 122 bone tumour samples including GCTBs. Benign GCTBs possessed few somatic alterations and no other known drivers besides the H3.3 mutation, whereas all malignant tumours harboured at least one additional driver mutation and exhibited genomic features resembling osteosarcomas, including high mutational burden, additional driver event(s), and a high degree of aneuploidy. The H3.3 mutation was found to predate the development of aneuploidy. In contrast to osteosarcomas, malignant H3.3-mutated tumours were enriched for a variety of alterations involving TERT, other than amplification, suggesting telomere dysfunction in the transformation of benign to malignant GCTB. DNA sequencing of the benign metastasising GCTB revealed no additional driver alterations; polyclonal seeding in the lung was identified, implying that the metastatic lesions represent an embolic event. Unsupervised clustering of DNA methylation profiles revealed that malignant H3.3-mutated tumours are distinct from their benign counterpart, and other bone tumours. Differential methylation analysis identified CCND1, encoding cyclin D1, as a plausible cancer driver gene in these tumours because hypermethylation of the CCND1 promoter was specific for GCTBs. We report here the genomic and methylation patterns underlying the rare clinical phenomena of benign metastasising and malignant transformation of GCTB and show how the combination of genomic and epigenomic findings could potentially distinguish benign from malignant GCTBs, thereby predicting aggressive behaviour in challenging diagnostic cases. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland. ispartof: JOURNAL OF PATHOLOGY vol:252 issue:4 pages:433-440 ispartof: location:England status: published

Myeloproliferative neoplasms (MPNs) are characterized by deregulation of mature blood cell production and increased risk of myelofibrosis (MF) and leukemic transformation. Numerous driver mutations have been identified but substantial disease heterogeneity remains unexplained, implying the involvement of additional as yet unidentified factors. The inflammatory microenvironment has recently attracted attention as a crucial factor in MPN biology, in particular whether inflammatory cytokines and chemokines contribute to disease establishment or progression. Here we present a large-scale study of serum cytokine profiles in more than 400 MPN patients and identify an essential thrombocythemia (ET)-specific inflammatory cytokine signature consisting of Eotaxin, GRO-α, and EGF. Levels of 2 of these markers (GRO-α and EGF) in ET patients were associated with disease transformation in initial sample collection (GRO-α) or longitudinal sampling (EGF). In ET patients with extensive genomic profiling data (n = 183) cytokine levels added significant prognostic value for predicting transformation from ET to MF. Furthermore, CD56+CD14+ pro-inflammatory monocytes were identified as a novel source of increased GRO-α levels. These data implicate the immune cell microenvironment as a significant player in ET disease evolution and illustrate the utility of cytokines as potential biomarkers for reaching beyond genomic classification for disease stratification and monitoring.

Summary Tobacco smoking causes lung cancer1–3, driven by the 60+ carcinogens in cigarette smoke that directly damage and mutate DNA4,5. The profound effects of tobacco on the lung cancer genome have been well documented6–10, but we lack equivalent data for normal bronchial cells. We sequenced whole genomes of 632 colonies derived from single bronchial epithelial cells across 16 subjects. Tobacco smoking was the major influence on mutation burden, adding 1000-10,000+ mutations/cell, massively increasing both within-subject and between-subject variance, and generating several distinct signatures of substitutions and indels. A population of cells in subjects with smoking history had mutation burdens equivalent to that expected for never-smokers: these cells had less damage from tobacco-specific mutational processes, were four-fold more frequent in ex-smokers than current smokers, and had significantly longer telomeres than their more mutated counterparts. Driver mutations increased in frequency with age, affecting 4-14% of cells in middle-aged never-smokers. In current smokers, ≥25% of cells carried driver mutations and 0-6% cells had 2 or even 3 drivers. Thus, tobacco smoking increases mutation burden, cell-to-cell heterogeneity and driver mutations, but quitting promotes replenishment of bronchial epithelium from mitotically quiescent cells that have avoided tobacco mutagenesis.

Soft tissue tumors of infancy encompass an overlapping spectrum of diseases that pose unique diagnostic and clinical challenges. We studied genomes and transcriptomes of cryptogenic congenital mesoblastic nephroma (CMN), and extended our findings to five anatomically or histologically related soft tissue tumors: infantile fibrosarcoma (IFS), nephroblastomatosis, Wilms tumor, malignant rhabdoid tumor, and clear cell sarcoma of the kidney. A key finding is recurrent mutation of EGFR in CMN by internal tandem duplication of the kinase domain, thus delineating CMN from other childhood renal tumors. Furthermore, we identify BRAF intragenic rearrangements in CMN and IFS. Collectively these findings reveal novel diagnostic markers and therapeutic strategies and highlight a prominent role of isolated intragenic rearrangements as drivers of infant tumors.

This abstract was withdrawn by the authors.

Allogeneic haematopoietic stem cell transplantation (HSCT) is a potentially curative treatment for over 40,000 patients/year in Europe and the US alone. However, substantial treatment-related mortality and morbidity, as well as risks of disease relapse give a survival rate of about 50% and leave considerable room for improvement. Despite being an established treatment for over 50 years, fundamental questions remain regarding its biology. For example, how many of the frequently >100 million transplanted CD34+ cells are true haematopoietic stem cells (HSCs), determined by long-term engraftment and contribution to multi-lineage hematopoiesis (long-term engrafting HSCs [LTE-HSCs])? What are the mutational consequences for transplanted HSCs given their proliferation and potential mutagenic insults in the post-transplant period? Most recently, the discovery of clonal haematopoiesis (CH) has raised interest in the interaction between this and HSCT. Do such clones further expand during HSCT? This may potentially lead to the devastating complication of donor-cell leukemia or other CH-related risks, e.g. cardiovascular disease. Recently, some studies have addressed this using targeted sequencing panels for myeloid cancer genes. However, many clonal expansions in normal blood are not driven by mutations in such genes, with evidence suggesting that the set of potential 'driver' genes is much larger than currently recognized (Poon et al, bioRXiv 2020). Advances in HSC tracking methodologies - using naturally-occurring somatic mutations as clonal markers (Lee-Six et al, Nature 2018) - provide a powerful tool to simultaneously address these questions. Whole-genome sequencing of hundreds of single-cell derived haematopoietic stem and progenitor cell (HSPC) colonies from a single individual is used to compile a complete set of somatic mutations in each colony founder cell, and the pattern of shared mutations amongst cells used to infer their phylogeny or 'family tree'. The constant rate of mutation acquisition during post-development life allows estimation of the timing of mutation acquisition. Using phylodynamic approaches borrowed from pathogen biology, patterns of branching points can be used to infer important parameters such as the size of population 'bottlenecks' (in this context the number of LTE-HSCs), and the growth dynamics of expanded clones. We selected 7 donor/ recipient (D/ R) pairs who had undergone HSCT 9-31 years previously. For each individual (D and R), peripheral blood CD34+ HSPC-derived colonies were grown on methylcellulose medium. Whole-genome sequencing (WGS) was performed on 100-300 colonies per individual - a total of 2,278 genomes. Mutations were called using established pipelines, then filtered to remove artefacts, germline variants, and in vitro mutations leaving only somatically-acquired mutations. Phylogenies for each D/R pair were inferred, using a maximum parsimony algorithm. Mutational signatures were extracted using a hierarchical dirichlet process. D/R phylogenies were compared using metrics of phylogenetic diversity. Clonal fractions of expanded clones in D/R were compared. Approximate Bayesian computation was used to estimate numbers of LTE-HSCs. Our results reveal that HSCT engraftment is remarkably polyclonal, with thousands of transplanted HSCs (in most cases >5,000) actively contributing to haematopoiesis decades after transplant. HSCs suffer little consequence in terms of their somatic mutation burden. Recipient haematopoiesis showed decreased clonal diversity compared to their donors with a mean 20% decrease of the Shannon's Diversity Index. This may partly result from increased selective pressures during HSCT. Intriguingly, several DNMT3A-driven expansions seen in donors had lower clonal fractions in recipients. Conversely, clones with >1 driver mutation (e.g. DNMT3A/CHEK2) showed larger expansions in recipients compared to donors, despite originating in the donor. DNMT3A mutations frequently originated in early development - in one case occurring in utero. We demonstrate the power of applying a novel clonal tracking approach to HSCT, for the first time giving a detailed picture of the clonal dynamics of engraftment. Overall, our findings are reassuring from a safety perspective, but the different clonal composition in recipients merits further investigation to better understand the factors involved. Disclosures Manz: University of Zurich: Patents & Royalties: CD117xCD3 TEA; CDR-Life Inc: Consultancy, Current holder of stock options in a privately-held company. Campbell: Mu Genomics: Current holder of individual stocks in a privately-held company, Membership on an entity's Board of Directors or advisory committees.

BACKGROUND The authors describe a family with a high penetrance of plasma cell dyscrasias, suggesting inheritance of an autosomal dominant risk allele. METHODS The authors performed whole-exome sequencing and reported on a combined approach aimed at the identification of causative variants and risk loci, using the wealth of data provided by this approach. RESULTS The authors identified gene mutations and single-nucleotide polymorphisms of potential significance, and pinpointed a known risk locus for myeloma as a potential area of transmissible risk in the family. CONCLUSIONS To the authors' knowledge, the current study is the first to provide a whole-exome sequencing approach to such cases, and a framework analysis that could be applied to further understanding of the inherited risk of developing plasma cell dyscrasias. Cancer 2017. © 2017 American Cancer Society.

Background Under normal physiological conditions, colonic crypts accrue ~40 somatic mutations for every year of life. That somatic mutations contribute to the development of cancer is well established, but their patterns, burden and functional consequences in diseases other than cancer have not been extensively studied and our understanding of the effects of chronic inflammation on the mutational profile and clonal structure of the colon is limited. Here, we investigated how the recurrent cycles of inflammation, ulceration and regeneration seen in IBD impact the mutational and clonal structure of intestinal epithelia. Methods We isolated and whole-genome sequenced ~400 individual colonic crypts from 46 IBD patients and compared these to ~400 crypts from 41 non-IBD controls. We compared the mutation burden, mutational signature exposure, clonal structure and cancer driver mutation landscape in crypts from actively and previously inflamed regions with crypts dissected from controls. Results We estimated the base substitution rate of affected colonic epithelial cells to be doubled after IBD onset. This change was primarily driven by acceleration of mutational processes ubiquitously observed in normal colon (Figure 1), and we did not detect an IBD-specific mutational process. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimetre-scale clonal expansions, even in the absence of mutations in KRAS, TP53 and APC (Figure 2). We discovered that mutations in ARID1A, PIGR and ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, were under positive selection in colonic crypts from IBD patients (Figure 3). With the exception of ARID1A, these genes and pathways have not been previously associated with cancer risk. A previously published mouse model of ZC3H12A suggests that LoF mutations in this gene may facilitate healing of affected mucosa without promoting tumorigenesis. This could make the encoded protein an attractive drug target. The observed enrichment of mutations in PIGR and IL17 and TLR pathways suggests that somatic mutations may initiate, maintain or perpetuate IBD pathogenesis through disruption of microbe-epithelial homeostasis. Conclusion Our results provide new insights into the consequences of chronic intestinal inflammation on the mutational profile and clonal structure of colonic epithelia. We identify the mutagens driving the increase in mutation burden and mutations which are under positive selection in the context of inflammation. Our results suggest that studying somatic mutations in the colon can reveal putative drug targets and pathogenic mechanisms for IBD.

A childhood tumor—from the beginning Many adult cancers arise from clonal expansions of mutant cells in normal tissue. These premalignant expansions are defined by somatic mutations shared by the cancers. Whether pediatric cancers originate in a similar way is unknown. Coorens et al. studied Wilms tumor, a childhood kidney cancer. Phylogenetic analyses revealed large clones of mutant cells in histologically and functionally normal kidney tissue long before tumor development. Thus, like adult tumors, Wilms tumor appears to arise from a premalignant tissue bed. Science , this issue p. 1247

Summary Inflammatory bowel disease (IBD) is a chronic inflammatory disease associated with increased risk of gastrointestinal cancers. We whole-genome sequenced 446 colonic crypts from 46 IBD patients and compared these to 412 crypts from 41 non-IBD controls from our previous publication on the mutation landscape of the normal colon. The average mutation rate of affected colonic epithelial cells is 2.4-fold that of healthy colon, and this increase is mostly driven by acceleration of mutational processes ubiquitously observed in normal colon. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimeter-scale clonal expansions. We discovered non-synonymous mutations in ARID1A, FBXW7, PIGR, ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, under positive selection in IBD. These results suggest distinct selection mechanisms in the colitis-affected colon and that somatic mutations potentially play a causal role in IBD pathogenesis., Graphical Abstract, Highlights • IBD-affected colons accrue substitutions and indels 2.4 and 7 times faster than normal • 17 signatures of mutational processes, including treatment • Millimeter-scale clonal expansions late in molecular time • Distinct mechanisms of positive selection of mutations in immune-related genes, Whole-genome sequencing of inflammatory bowel disease patient samples allows insight into mutational burdens and processes associated with disease, including putative driver mutations positively selected in the diseased colon.

Before squamous cell lung cancer develops, pre-cancerous lesions can be found in the airways. From longitudinal monitoring, we know that only half of such lesions become cancer, whereas a third spontaneously regress. While recent studies have described the presence of an active immune response in high-grade lesions, the mechanisms underpinning clinical regression of pre-cancerous lesions remain unknown. Here, we show that host immune surveillance is strongly implicated in lesion regression. Using bronchoscopic biopsies from human subjects, we find that regressive carcinoma in-situ lesions harbour more infiltrating immune cells than those that progress to cancer. Moreover, molecular profiling of these lesions identifies potential immune escape mechanisms specifically in those that progress to cancer: antigen presentation is impaired by genomic and epigenetic changes, TGF-beta signalling is overactive, and the immunomodulator TNFSF9 is downregulated. Changes appear intrinsic to the CIS lesions as the adjacent stroma of progressive and regressive lesions are transcriptomically similar. This study identifies mechanisms by which pre-cancerous lesions evade immune detection during the earliest stages of carcinogenesis and forms a basis for new therapeutic strategies that treat or prevent early stage lung cancer.

The colorectal adenoma-carcinoma sequence has provided a paradigmatic framework for understanding the successive somatic genetic changes and consequent clonal expansions that lead to cancer1. However, our understanding of the earliest phases of colorectal neoplastic changes-which may occur in morphologically normal tissue-is comparatively limited, as for most cancer types. Here we use whole-genome sequencing to analyse hundreds of normal crypts from 42 individuals. Signatures of multiple mutational processes were revealed; some of these were ubiquitous and continuous, whereas others were only found in some individuals, in some crypts or during certain periods of life. Probable driver mutations were present in around 1% of normal colorectal crypts in middle-aged individuals, indicating that adenomas and carcinomas are rare outcomes of a pervasive process of neoplastic change across morphologically normal colorectal epithelium. Colorectal cancers exhibit substantially increased mutational burdens relative to normal cells. Sequencing normal colorectal cells provides quantitative insights into the genomic and clonal evolution of cancer.

INTRODUCTION: Cancer pathogenesis is usually characterized by a long evolutionary process where genomic driver events accumulate over time, conferring advantage to distinct subclones, allowing their expansion and progression. METHODS: To investigate the multiple myeloma (MM) evolutionary history, we characterized the mutational processes' landscape and activity over time utilizing a large cohort of 89 whole genomes and 973 exomes. To improve the accuracy of mutational signatures analysis, we analyzed both the 3' and 5' nucleotide context of each mutation and we developed the novel fitting algorithm mmSig, which fits the entire mutational catalogue of each patient with the mutational signatures involved in MM pathogenesis. The contribution of each mutational signature was then corrected based on the cosine similarity between the original 96-mutational profile and the reconstructed profile generated without that signature. To reconstruct the genetic evolutionary history of each patient's cancer, we integrated two approaches. First dividing all mutations into clonal (early) or subclonal (late), then subdivided the clonal mutations into duplicated mutations (present on two alleles and therefore acquired before the duplication) or non-duplicated mutations (detected on a single allele), reflecting either pre-gain and post-gain mutations on the minor allele, or post-gain mutations acquired on one of the duplicated alleles. RESULTS Eight mutational signatures were identified, seven of which showed significant similarity with the most recent mutational signature catalogue (i.e SBS1, SBS2, SBS5, SBS8, SBS9, SBS13 and SBS18). The new mutational signature (named SBS-MM1) was observed only among relapsed patients exposed to alkylating agents (i.e melphalan). The etiology of this specific signature was further confirmed by analyzing recent whole genomes public data from human-induced pluripotent stem cells exposed to melphalan (Kucab et al, Cell 2019). Reconstructing the chronological activity of each mutational signature, we identified four different routes to acquire the full mutational spectrum in MM based on the differential temporal activity of AID (SBS9) and APOBEC (SBS2 and SBS13). Our data indicate that AID activity is not limited to the first contact with the GC, but persists in the majority of patients, behaving similarly to a B-memory cells, capable of re-entering the germinal center upon antigen stimulation to undergo clonal expansion several times before MM diagnosis. Next, we confirmed the clock-like nature (i.e constant mutation rate) of SBS5 in MM and other post-germinal center disorders such as chronic lymphocytic leukemia and B-cell lymphomas. Based on the SBS5 mutation rates and the corrected ratio between duplicated and non-duplicated mutations within large chromosomal gains, we could time the acquisition of the first copy number gain during the life history of each MM patient. Intriguingly, the first MM chromosomal duplication was acquired on average 38 years (ranges 11-64) before sample collection. In 23/27 (85%) cases the first multi gain event occurred before 30 years of age, and in 13/27 (48%) before 20 years reflecting a long and slow process potentially influenced and accelerated by extrinsic and intrinsic factors. DISCUSSION Our analysis provides a glimpse into the early stages of myelomagenesis, where acquisition of the first key drivers precedes cancer diagnosis by decades. Defining the time window when transformation occurs opens up for new avenues of research: to identify causal mechanisms of disease initiation and evolution, to better define the optimal time to start therapy, and ultimately develop early prevention strategies. Disclosures Bolli: CELGENE: Honoraria; JANSSEN: Honoraria; GILEAD: Other: Travel expenses. Corradini:Janssen: Honoraria, Other: Travel Costs; Jazz Pharmaceutics: Honoraria; KiowaKirin: Honoraria; Servier: Honoraria; Takeda: Honoraria, Other: Travel Costs; Kite: Honoraria; Novartis: Honoraria, Other: Travel Costs; Gilead: Honoraria, Other: Travel Costs; Roche: Honoraria; Sanofi: Honoraria; BMS: Other: Travel Costs. Anderson:Janssen: Consultancy, Speakers Bureau; Takeda: Consultancy, Speakers Bureau; Celgene: Consultancy, Speakers Bureau; Bristol-Myers Squibb: Other: Scientific Founder; Oncopep: Other: Scientific Founder; Amgen: Consultancy, Speakers Bureau; Sanofi-Aventis: Other: Advisory Board. Moreau:Celgene: Consultancy, Honoraria; Janssen: Consultancy, Honoraria; Amgen: Consultancy, Honoraria; Takeda: Consultancy, Honoraria; AbbVie: Consultancy, Honoraria. Papaemmanuil:Celgene: Research Funding. Avet-Loiseau:takeda: Consultancy, Other: travel fees, lecture fees, Research Funding; celgene: Consultancy, Other: travel fees, lecture fees, Research Funding. Munshi:Adaptive: Consultancy; Amgen: Consultancy; Celgene: Consultancy; Janssen: Consultancy; Takeda: Consultancy; Oncopep: Consultancy; Abbvie: Consultancy. Landgren:Karyopharm: Membership on an entity's Board of Directors or advisory committees; Theradex: Other: IDMC; Adaptive: Honoraria, Membership on an entity's Board of Directors or advisory committees; Abbvie: Membership on an entity's Board of Directors or advisory committees; Sanofi: Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Takeda: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Merck: Other: IDMC; Amgen: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding.

Assessing smoke damage in cancer genomes We have known for over 60 years that smoking tobacco is one of the most avoidable risk factors for cancer. Yet the detailed mechanisms by which tobacco smoke damages the genome and creates the mutations that ultimately cause cancer are still not fully understood. Alexandrov et al. examined mutational signatures and DNA methylation changes in over 5000 genome sequences from 17 different cancer types linked to smoking (see the Perspective by Pfeifer). They found a complex pattern of mutational signatures. Only cancers originating in tissues directly exposed to smoke showed a signature characteristic of the known tobacco carcinogen benzo[ a ]pyrene. One mysterious signature was shared by all smoking-associated cancers but is of unknown origin. Smoking had only a modest effect on DNA methylation. Science , this issue p. 618 ; see also p. 549

Background: At present, conventional clinical and histopathological evaluations are not sufficient to distinguish biologically indolent cancers from those that will exhibit aggressive behavior. We hypothesize that global transcriptomic activity of tumor cells reflects the end cumulative result of somatic, germline, and epigenetic alterations, as well as additional transcriptional regulatory events. Therefore, it may be more directly associated with clinical outcomes. However, the total number of mRNA molecules is not directly measurable, either in bulk or single-cell RNA sequencing data. To this end, we develop a novel metric: the transcriptional activity score (TAS), to measure the relative global tumor-cell specific transcriptional activity in heterogeneous tumor samples. Materials and Methods: We propose TAS as the ratio of average total transcript proportion over the count proportion of tumor cells versus surrounding non-tumor cells. The transcript proportions are estimated using RNAseq deconvolution method DeMixT and the count proportions are estimated using DNAseq deconvolution methods such as ASCAT and ABSOLUTE. Using matching bulk RNA and DNA sequencing data from the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC), we calculated TAS for a total of 5,031 patient samples across 15 cancer types. For validation, we obtained TAS for two genomic studies: 1) from patients with early-onset prostate cancer (n=99) as part of the ICGC, and 2) from patients with localized prostate cancer as part of the Canadian Prostate Cancer Genome Network (CPC-GENE, n=144). Results: We found that higher TAS corresponds to a more aggressive state of cancer, as characterized by MYC dysregulation, genome instability, known marker genes, and molecular subtypes. By examining the association between TAS and survival outcomes across cancer types, we also found that TAS refines the prognostic ability of pathologic stage, identifying aggressive early-stage tumors associated with poor survival as well as late-stage tumors with favorable outcomes. In prostate cancer, TAS is linearly associated with progression-free probabilities, useful to rank patients within the median risk group (Gleason score = 7). This added prediction power is consistent in TCGA and two independent validation data (ICGC and CPC-GENE). Conclusion: We have developed a new summary metric using matched DNA and RNA sequencing data from tumor samples, to compute, in vivo and using deconvolution, the relative global gene expression level of tumor cells. The TAS metric evaluates global transcriptional activity, an intrinsic behavior of cells that is well-known, but now for the first time is shown through TAS to be associated with prognosis. TAS may serve as a tractable phenotype to help elucidate the biology that underlies metastasis, prognosis and response to treatment in cancer patients. Citation Format: Shaolong Cao, Jennifer R. Wang, Jonas Demeulemeester, Jingxiao Chen, Kaixian Yu, Peng Yang, Bora Lim, Alfonso Urbanucci, Peter Campbell, Hongtu Zhu, Peter Van loo, Wenyi Wang. Global tumor transcriptional activity reveals aggressiveness across multiple cancers [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 829.

Background Lung cancer is the leading cause of cancer death, of which 80-90% are attributed to tobacco smoking. Our understanding of how tobacco exposure affects mutational burden, mutational signatures, driver mutations and clonal dynamics in normal lung tissue is limited. Also, genetic differences between normal cells and cancer cells in lung has not been fully elucidated. Methods To access the landscape of somatic mutations in normal bronchial epithelia, we sequenced whole genomes of 632 single-cell-derived colonies of lung basal cells from 16 patients, including 3 children, 4 never-smokers, 6 ex-smokers and 3 current smokers. Five patients had squamous cell carcinomas or carcinoma in situ, which we also sequenced to compare the genetic alteration between normal and cancer cells. Results A positive correlation between the number of base substitutions and age was observed, and 22 somatic mutations accumulated per year (95% confidence interval:20-25; P=10−8). Previous or current smoking significantly increased mutational burden: 2330 substitutions in ex-smokers and 5300 in current smokers. In addition, tobacco smoking is massively increasing both between-subject and within-subject variance of mutation burden. A population of cells in subjects with smoking history had mutation burdens equivalent to that expected for never-smokers: these cells had less damage from tobacco-specific mutational processes, and were four-fold more frequent in ex-smokers than current smokers. Signature analysis revealed that the same mutational signatures seen in lung cancers operate both in patients with and without smoking history: endogeneous signatures, including COSMIC signatures 1 and 5, and APOBEC-related signatures. Three mutational signatures were largely restricted to current or ex-smokers, including known smoking-related COSMIC signature 4 characterized by C>A transversions, signature 16 characterized by T>C mutations with extremely strong transcription strand bias and a new signature characterized by T>A and T>C mutations. dN/dS method identified driver genes in normal bronchial epithelium, including NOTCH1, TP53 and FAT1, which were overlapped with those seen in squamous cell lung cancers and other normal squamous tissues such as esophagus and skin. Driver mutations increased in frequency with age, affecting 4-14% of cells in middle-aged never-smokers. In current smokers, ≥25% of cells carried driver mutations and 0-6% cells had 2 or even 3 drivers. Layering driver mutations onto phylogenetic trees revealed that driver mutations occurred in early life. Compared to the normal bronchial epithelial cells, lung cancers and precancerous lesions were characterized by extensive copy number changes or structural variants and distinct set of driver mutations. Conclusions Tobacco smoking increases mutation burden, cell-to-cell heterogeneity and driver mutations. Our data of genetic lesions in normal bronchial cells provides insights into genetic alterations that drive carcinogenesis in lung. Citation Format: Kenichi Yoshida, Kate HC Gowers, Henry Lee-Six, Deepak P. Chandrasekharan, Tim Coorens, Elizabeth F. Maughan, Kathryn Beal, Andrew Menzies, Fraser R. Millar, Elizabeth Anderson, Sarah E. Clarke, Adam Pennycuick, Ricky M. Thackeray, Colin R. Butler, Nobuyuki Kakiuchi, Tomonori Hirano, Robert E. Hynds, Michael R. Stratton, Inigo Martincorena, Sam M. Janes, Peter J. Campbell. Tobacco exposure and somatic mutations in normal bronchial epithelia [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 2342.

The congenital bone marrow failure syndrome Shwachman-Diamond Syndrome (SDS) is characterized by defective ribosome biogenesis, with a significant risk of developing a severe myelodysplastic syndrome or leukemia. The most common causal germ line mutations are in the SBDS gene. It remains unclear how the germline ribosomal dysfunction/stress affects individual hematopoietic stem cells and contributes to the aberrant population dynamics that drive specific clonal expansions and eventual leukemic transformation. To address this, we adopted a clonal tracking approach using whole genome sequencing of single hematopoietic stem and progenitor cell (HSPC)-derived colonies followed by phylogenetic tree construction to assess clonal relatedness in 9 SDS patients. The mean mutation burden in SDS HSPCs increased linearly with age and the majority of mutations were C >T variants, consistent with mutational signatures in ageing. Single nucleotide variants in annotated genes across patients were found in known drivers TP53, EIF6 and other ribosomal proteins: all of which were mutually exclusive in different HSPC clones. Distinct mutations in the same genes were found within the same patients, demonstrating that individual clones must acquire independent genetic rescue mechanisms. Whole genome sequencing permitted us to uncover a high degree of intra-patient heterogeneity, charting multiple independent paths that individual cells can take to combat defective ribosome deficiency caused by SBDS loss. This study sets the stage for functional analysis of potential leukemia driver mutations and longitudinal tracking of clonal dynamics throughout disease progression and re-enforces the value of using whole genome sequencing to track individual stem cell clones in disease.

We analyzed whole genomes of unique paired samples from smoldering multiple myeloma (SMM) patients progressing to multiple myeloma (MM). We report that the genomic landscape, including mutational profile and structural rearrangements at the smoldering stage is very similar to MM. Paired sample analysis shows two different patterns of progression: a “static progression model”, where the subclonal architecture is retained as the disease progressed to MM suggesting that progression solely reflects the time needed to accumulate a sufficient disease burden; and a “spontaneous evolution model”, where a change in the subclonal composition is observed. We also observe that activation-induced cytidine deaminase plays a major role in shaping the mutational landscape of early subclinical phases, while progression is driven by APOBEC cytidine deaminases. These results provide a unique insight into myelomagenesis with potential implications for the definition of smoldering disease and timing of treatment initiation., Smoldering MM (SMM) is a premalignant stage of multiple myeloma (MM). Here the authors perform whole genome sequencing of unique paired samples of SMM progressing to MM, and show that the genomic landscape at the SMM stage is very similar to MM, but trajectories of evolution can vary from patient to patient.

The colorectal adenoma-carcinoma sequence has provided a paradigmatic framework for understanding the successive somatic genetic changes and consequent clonal expansions leading to cancer. As for most cancer types, however, understanding of the earliest phases of colorectal neoplastic change, which may occur in morphologically normal tissue, is comparatively limited because of the difficulty of detecting somatic mutations in normal cells. Each colorectal crypt is a small clone of cells derived from a single recently-existing stem cell. Here, we whole genome sequenced hundreds of normal crypts from 42 individuals. Signatures of multiple mutational processes were revealed, some ubiquitous and continuous, others only found in some individuals, in some crypts or during some phases of the cell lineage from zygote to adult cell. Likely driver mutations were present in ∼1% of normal colorectal crypts in middle-aged individuals, indicating that adenomas and carcinomas are rare outcomes of a pervasive process of neoplastic change across morphologically normal colorectal epithelium.

The molecular alterations that occur in cells before cancer is manifest are largely uncharted. Lung carcinoma in situ (CIS) lesions are the pre-invasive precursor to squamous cell carcinoma. Although microscopically identical, their future is in equipoise, with half progressing to invasive cancer and half regressing or remaining static. The cellular basis of this clinical observation is unknown. Here, we profile the genomic, transcriptomic, and epigenomic landscape of CIS in a unique patient cohort with longitudinally monitored pre-invasive disease. Predictive modeling identifies which lesions will progress with remarkable accuracy. We identify progression-specific methylation changes on a background of widespread heterogeneity, alongside a strong chromosomal instability signature. We observed mutations and copy number changes characteristic of cancer and chart their emergence, offering a window into early carcinogenesis. We anticipate that this new understanding of cancer precursor biology will improve early detection, reduce overtreatment, and foster preventative therapies targeting early clonal events in lung cancer. Postprint

The transcription factor FOS has long been implicated in the pathogenesis of bone tumours, following the discovery that the viral homologue, v-fos, caused osteosarcoma in laboratory mice. However, mutations of FOS have not been found in human bone-forming tumours. Here, we report recurrent rearrangement of FOS and its paralogue, FOSB, in the most common benign tumours of bone, osteoblastoma and osteoid osteoma. Combining whole-genome DNA and RNA sequences, we find rearrangement of FOS in five tumours and of FOSB in one tumour. Extending our findings into a cohort of 55 cases, using FISH and immunohistochemistry, provide evidence of ubiquitous mutation of FOS or FOSB in osteoblastoma and osteoid osteoma. Overall, our findings reveal a human bone tumour defined by mutations of FOS and FOSB., FOS has been linked to bone tumour pathogenesis, and viral homologue v-fos causes osteosarcoma in mice. Here, the authors report rearrangement of FOS and its paralogue FOSB in osteoblastoma and osteoid osteoma, revealing human bone tumours that are defined by mutations of FOS and FOSB.

Every cancer originates from a single cell. During expansion of the neoplastic cell population, individual cells acquire genetic and phenotypic differences from each other. Here, to investigate the nature and extent of intra-tumour diversification, we characterized organoids derived from multiple single cells from three colorectal cancers as well as from adjacent normal intestinal crypts. Colorectal cancer cells showed extensive mutational diversification and carried several times more somatic mutations than normal colorectal cells. Most mutations were acquired during the final dominant clonal expansion of the cancer and resulted from mutational processes that are absent from normal colorectal cells. Intra-tumour diversification of DNA methylation and transcriptome states also occurred; these alterations were cell-autonomous, stable, and followed the phylogenetic tree of each cancer. There were marked differences in responses to anticancer drugs between even closely related cells of the same tumour. The results indicate that colorectal cancer cells experience substantial increases in somatic mutation rate compared to normal colorectal cells, and that genetic diversification of each cancer is accompanied by pervasive, stable and inherited differences in the biological states of individual cancer cells.

Summary Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5′ UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor’s most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention., Graphical Abstract, Highlights • Novel hotspot of driver mutations in 5′-UTR repressor of TERT, expanding telomeres • Most common cause of 3p loss is a chromothripsis event, generating concurrent 5q gain • t(3;5) event occurs in childhood or adolescence, decades before tumor diagnosed • Initial clonal expansion after 3p loss starts from only a few hundred cells, Combination of whole-genome sequencing analysis and a multi-region sampling approach provides insights into the nature and timing of key oncogenic events in clear cell renal cell carcinoma, depicts the evolutionary trajectories of tumors in patients and highlights the opportunity for early intervention.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with > 10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance. We thank Aida Murra, Naheed Shaikh, Justine Korteweg, Jeremy Tai, Eleanor Carlyle, Leonora Conneely, KimEdmonds, Karla Lingard, Karen O'Meara, Helen Breeze, Sarah Sarker, Lesley Cooper, Linda Shephard, Susie Slater, and Catherine Rogers for study support. We thank the patients and their families. S.T. and H.X. are funded by Cancer Research UK (CRUK) (C50947/A18176). S.T., T.C., J.L., and M.G. are funded by the National Institute for Health Research (NIHR) Biomedical Research Centre (BRC) at the Royal Marsden Hospital and Institute of Cancer Research (A109). J.I.L. is funded by the Ministerio de Economia y Competitividad (MINECO, SAF2016-79847-R). M.S., A.S., J. Lynch, R.F., L.A., and L.S. are funded by the Royal Marsden Cancer Charity. K.L. is funded by UK Medical Research Council (MR/P014712/1). T.B.K.W. is funded by the European Union Seventh Framework Programme (FP7-People-2013-ITN). M.J.-H. is funded by the NIHR. N.M.L. is a Winton Group Leader in recognition of the Winton Charitable Foundation's support towards the establishment of The Francis Crick Institute. N.M.L. is additionally funded by a Wellcome Trust Joint Investigator Award (103760/Z/14/Z) and the MRC eMedLab Medical Bioinformatics Infrastructure Award (MR/L016311/1). M.K. is funded by the Danish Cancer Society grant (R90-A6213). P.C. is funded by the Wellcome Trust (WT088340MA). N.M. receives funding from CRUK, Rosetrees, and the NIHR BRC at University College London Hospitals. C.S. is a Royal Society Napier Research Professor. C.S. is funded by Cancer Research UK (TRACERx and CRUK Cancer Immunotherapy Catalyst Network), the CRUK Lung Cancer Centre of Excellence, Stand Up 2 Cancer (SU2C), the Rosetrees and Stoneygate Trusts, NovoNordisk Foundation (THESEUS), Marie Curie Network PloidyNet, the NIHR BRC at University College London Hospitals, and the CRUK University College London Experimental Cancer Medicine Centre. C.S., O.P., D.R., I.C., and Z.S. are funded by NovoNordisk Foundation (16584). O.P., D.R., and I.C. are funded by the National Research, Development and Innovation Office of Hungary (NVKP_16-1-20160004). This work was supported by CRUK (Clinical Scientist Fellowship to S.T., C50947/A18176), The Francis Crick Institute, which receives its core funding from Cancer Reseach UK (FC010110), the UK Medical Research Council (FC010110), the Wellcome Trust (FC010110), and NIHR BRC at the Royal Marsden Hospital and Institute of Cancer Research (A109). High-performance computing was supported by eMedLab. The results published here are in whole or part based upon data generated by the TCGA Research Network: http://cancergenome.nih.gov/.

Cancers emerge from an on-going Darwinian evolutionary process, often leading to multiple competing subclones within a single primary tumour1-4. This evolutionary process culminates in the formation of metastases, which is the cause of 90% of cancer-related deaths5. However, despite its clinical importance, little is known about the principles governing the dissemination of cancer cells to distant organs. Although the hypothesis that each metastasis originates from a single tumour cell is generally supported6-8, recent studies using mouse models of cancer demonstrated the existence of polyclonal seeding from and inter-clonal cooperation between multiple subclones9,10. In this study, we sought definitive evidence for the existence of polyclonal seeding in human malignancy and to establish the clonal relationship among different metastases in the context of androgen-deprived metastatic prostate cancer. Using whole genome sequencing, we characterised multiple metastases arising from prostate tumours in ten patients. Integrated analyses of subclonal architecture revealed the patterns of metastatic spread in unprecedented detail. Metastasis-to-metastasis spread was found to be common, either through de novo monoclonal seeding of daughter metastases or, in five cases, through the transfer of multiple tumour clones between metastatic sites. Lesions affecting tumour suppressor genes usually occur as single events, whereas mutations in genes involved in androgen receptor signalling commonly involve multiple, convergent events in different metastases. Our results elucidate in detail the complex patterns of metastatic spread and further our understanding of the development of resistance to androgen deprivation therapy in prostate cancer.

International audience; In multiple myeloma, next generation sequencing (NGS) has expanded our knowledge of genomic lesions, and highlighted a dynamic and heterogeneous composition of the tumor. Here, we used NGS to characterize the genomic landscape of 418 multiple myeloma cases at diagnosis and correlate this with prognosis and classification. Translocations and copy number changes (CNAs) had a preponderant contribution over gene mutations in defining the genotype and prognosis of each case. Known and novel independent prognostic markers were identified in our cohort of proteasome inhibitor and IMiD-treated patients with long follow-up, including events with context-specific prognostic value, such as deletions of the PRDM1 gene. Taking advantage of the comprehensive genomic annotation of each case, we used innovative statistical approaches to identify potential novel myeloma subgroups. We observed clusters of patients stratified based on the overall number of mutations and number/type of CNAs, with distinct effects on survival, suggesting that extended genotype of multiple myeloma at diagnosis may lead to improved disease classification and prognostication.

Chordoma is a malignant, often incurable bone tumour showing notochordal differentiation. Here, we defined the somatic driver landscape of 104 cases of sporadic chordoma. We reveal somatic duplications of the notochordal transcription factor brachyury (T) in up to 27% of cases. These variants recapitulate the rearrangement architecture of the pathogenic germline duplications of T that underlie familial chordoma. In addition, we find potentially clinically actionable PI3K signalling mutations in 16% of cases. Intriguingly, one of the most frequently altered genes, mutated exclusively by inactivating mutation, was LYST (10%), which may represent a novel cancer gene in chordoma., Chordoma is a rare often incurable malignant bone tumour. Here, the authors investigate driver mutations of sporadic chordoma in 104 cases, revealing duplications in notochordal transcription factor brachyury (T), PI3K signalling mutations, and mutations in LYST, a potential novel cancer gene in chordoma.

With a few exceptions, cancers typically carry more than one driver mutation, sometimes five, ten, or more, and these driver mutations do not necessarily assort randomly. In this issue of Cancer Cell, Mina et al. systematically characterize patterns of co-mutation and mutual exclusivity in 6,456 cancers across 23 tumor types.

Summary Patterns of genomic evolution between primary and metastatic breast cancer have not been studied in large numbers, despite patients with metastatic breast cancer having dismal survival. We sequenced whole genomes or a panel of 365 genes on 299 samples from 170 patients with locally relapsed or metastatic breast cancer. Several lines of analysis indicate that clones seeding metastasis or relapse disseminate late from primary tumors, but continue to acquire mutations, mostly accessing the same mutational processes active in the primary tumor. Most distant metastases acquired driver mutations not seen in the primary tumor, drawing from a wider repertoire of cancer genes than early drivers. These include a number of clinically actionable alterations and mutations inactivating SWI-SNF and JAK2-STAT3 pathways., Graphical Abstract, Highlights • Metastases mostly disseminate late from primary breast tumors, keeping most drivers • Drivers at relapse sample from a wider range of cancer genes than in primary tumors • Mutations in SWI-SNF complex and inactivated JAK-STAT signaling enriched at relapse • Mutational processes similar in primary and relapse; radiotherapy can damage genome, By sequencing primary, locally relapsed, and metastatic breast cancers, Yates et al. show that clones seeding metastasis or relapse disseminate late from primary tumors but continue to acquire mutations, including clinically actionable alterations and mutations inactivating the SWI/SNF and JAK2-STAT3 pathways.

International audience; Approximately 1–5% of breast cancers are attributed to inherited mutations in BRCA1 or BRCA2 and are selectively sensitive to poly(ADP-ribose) polymerase (PARP) inhibitors. In other cancer types, germline and/or somatic mutations in BRCA1 and/or BRCA2 (BRCA1/BRCA2) also confer selective sensitivity to PARP inhibitors. Thus, assays to detect BRCA1/BRCA2-deficient tumors have been sought. Recently, somatic substitution, insertion/deletion and rearrangement patterns, or 'mutational signatures', were associated with BRCA1/BRCA2 dysfunction. Herein we used a lasso logistic regression model to identify six distinguishing mutational signatures predictive of BRCA1/BRCA2 deficiency. A weighted model called HRDetect was developed to accurately detect BRCA1/BRCA2-deficient samples. HRDetect identifies BRCA1/BRCA2-deficient tumors with 98.7% sensitivity (area under the curve (AUC) = 0.98). Application of this model in a cohort of 560 individuals with breast cancer, of whom 22 were known to carry a germline BRCA1 or BRCA2 mutation, allowed us to identify an additional 22 tumors with somatic loss of BRCA1 or BRCA2 and 47 tumors with functional BRCA1/BRCA2 deficiency where no mutation was detected. We validated HRDetect on independent cohorts of breast, ovarian and pancreatic cancers and demonstrated its efficacy in alternative sequencing strategies. Integrating all of the classes of mutational signatures thus reveals a larger proportion of individuals with breast cancer harboring BRCA1/BRCA2 deficiency (up to 22%) than hitherto appreciated (~1–5%) who could have selective therapeutic sensitivity to PARP inhibition.

Our knowledge of the genetic basis of myelodysplastic syndromes (MDS) and myelodysplastic/myeloproliferative neoplasms (MDS/MPN) has considerably improved. To define genotype/phenotype relationships of clinical relevance, we studied 308 patients with MDS, MDS/MPN, or acute myeloid leukemia evolving from MDS. Unsupervised statistical analysis, including the World Health Organization classification criteria and somatic mutations, showed that MDS associated with SF3B1-mutation (51 of 245 patients, 20.8%) is a distinct nosologic entity irrespective of current morphologic classification criteria. Conversely, MDS with ring sideroblasts with nonmutated SF3B1 segregated in different clusters with other MDS subtypes. Mutations of genes involved in DNA methylation, splicing factors other than SF3B1, and genes of the RAS pathway and cohesin complex were independently associated with multilineage dysplasia and identified a distinct subset (51 of 245 patients, 20.8%). No recurrent mutation pattern correlated with unilineage dysplasia without ring sideroblasts. Irrespective of driver somatic mutations, a threshold of 5% bone marrow blasts retained a significant discriminant value for identifying cases with clonal evolution. Comutation of TET2 and SRSF2 was highly predictive of a myeloid neoplasm characterized by myelodysplasia and monocytosis, including but not limited to, chronic myelomonocytic leukemia. These results serve as a proof of concept that a molecular classification of myeloid neoplasms is feasible.