Deciphering TP53 mutant Cancer Evolution with Single-Cell Multi-Omics

TP53 is the most commonly mutated gene in human cancer, typically occurring in association with complex cytogenetics and dismal outcomes. Understanding the genetic and non-genetic determinants of TP53-mutation driven clonal evolution and subsequent transformation is a crucial step towards the design of rational therapeutic strategies. Here, we carry out allelic resolution single-cell multi-omic analysis of haematopoietic stem/progenitor cells (HSPC) from patients with a myeloproliferative neoplasm who transform to TP53-mutant secondary acute myeloid leukaemia (AML), a tractable model of TP53-mutant cancer evolution. All patients showed dominant TP53 ‘multi-hit’ HSPC clones at transformation, with a leukaemia stem cell transcriptional signature strongly predictive of adverse outcome in independent cohorts, across both TP53-mutant and wild-type AML. Through analysis of serial samples and antecedent TP53-heterozygous clones, we demonstrate a hitherto unrecognised effect of chronic inflammation, which supressed TP53 wild-type HSPC whilst enhancing the fitness advantage of TP53 mutant cells. Our findings will facilitate the development of risk-stratification, early detection and treatment strategies for TP53-mutant leukaemia, and are of broader relevance to other cancer types.

[1]  Elisabeth F. Heuston,et al.  Transitions in lineage specification and gene regulatory networks in hematopoietic stem/progenitor cells over human development , 2021, Cell reports.

[2]  Allison P. Heath,et al.  The NCI Genomic Data Commons , 2021, Nature Genetics.

[3]  J. Marioni,et al.  Coordinated changes in gene expression kinetics underlie both mouse and human erythroid maturation , 2020, Genome Biology.

[4]  Jared S. Fowles,et al.  A Humanized Animal Model Predicts Clonal Evolution and Therapeutic Vulnerabilities in Myeloproliferative Neoplasms , 2020, bioRxiv.

[5]  Adam J. Mead,et al.  TARGET-Seq: A Protocol for High-Sensitivity Single-Cell Mutational Analysis and Parallel RNA Sequencing , 2020, STAR protocols.

[6]  V. Ugo,et al.  Leukemic evolution of polycythemia vera and essential thrombocythemia: genomic profiles predict time to transformation. , 2020, Blood advances.

[7]  A. Zebisch,et al.  Acute Myeloid Leukemia and Myelodysplastic Syndromes with TP53 Aberrations — A Distinct Stem Cell Disorder , 2020, Clinical Cancer Research.

[8]  R. Levine,et al.  Leukemia secondary to myeloproliferative neoplasms , 2020 .

[9]  Elisabeth F. Heuston,et al.  Single-Cell Analyses Reveal Megakaryocyte-Biased Hematopoiesis in Myelofibrosis and Identify Mutant Clone-Specific Targets , 2020, Molecular cell.

[10]  M. Wilkins,et al.  RNA Splicing Alterations Induce a Cellular Stress Response Associated with Poor Prognosis in Acute Myeloid Leukemia , 2020, Clinical Cancer Research.

[11]  P. A. Futreal,et al.  Clonal evolution of acute myeloid leukemia revealed by high-throughput single-cell genomics , 2020, Nature Communications.

[12]  M. Bashton,et al.  Concordance of copy number abnormality detection using SNP arrays and Multiplex Ligation-dependent Probe Amplification (MLPA) in acute lymphoblastic leukaemia , 2020, Scientific Reports.

[13]  Benjamin L. Ebert,et al.  Implications of TP53 allelic state for genome stability, clinical presentation and outcomes in myelodysplastic syndromes , 2019, Nature Medicine.

[14]  Howard Y. Chang,et al.  Single-cell multiomic analysis identifies regulatory programs in mixed-phenotype acute leukemia , 2019, Nature Biotechnology.

[15]  Fabian J Theis,et al.  Generalizing RNA velocity to transient cell states through dynamical modeling , 2019, Nature Biotechnology.

[16]  Jörg Menche,et al.  Mutational Landscape of the Transcriptome Offers Putative Targets for Immunotherapy of Myeloproliferative Neoplasms. , 2019, Blood.

[17]  S. Mccarroll,et al.  Monogenic and polygenic inheritance become instruments for clonal selection , 2019, bioRxiv.

[18]  R. Satija,et al.  Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression , 2019, Genome Biology.

[19]  Alice Giustacchini,et al.  Unravelling Intratumoral Heterogeneity through High-Sensitivity Single-Cell Mutational Analysis and Parallel RNA Sequencing , 2019, Molecular cell.

[20]  J. Tripodi,et al.  Advanced forms of MPNs are accompanied by chromosomal abnormalities that lead to dysregulation of TP53. , 2018, Blood advances.

[21]  I. Amit,et al.  Single cell dissection of plasma cell heterogeneity in symptomatic and asymptomatic myeloma , 2018, Nature Medicine.

[22]  Erik Sundström,et al.  RNA velocity of single cells , 2018, Nature.

[23]  K. Ballman,et al.  Somatic mutations precede acute myeloid leukemia years before diagnosis , 2018, Nature Medicine.

[24]  Yakir A Reshef,et al.  Insights about clonal hematopoiesis from 8,342 mosaic chromosomal alterations , 2018, Nature.

[25]  A. Tefferi,et al.  Targeted next-generation sequencing in blast phase myeloproliferative neoplasms. , 2018, Blood advances.

[26]  Alice Giustacchini,et al.  Ezh2 and Runx1 Mutations Collaborate to Initiate Lympho-Myeloid Leukemia in Early Thymic Progenitors. , 2018, Cancer cell.

[27]  Scott W. Lowe,et al.  Putting p53 in Context , 2017, Cell.

[28]  M. Ladanyi,et al.  Therapy-Related Clonal Hematopoiesis in Patients with Non-hematologic Cancers Is Common and Associated with Adverse Clinical Outcomes. , 2017, Cell stem cell.

[29]  Hannah A. Pliner,et al.  Reversed graph embedding resolves complex single-cell trajectories , 2017, Nature Methods.

[30]  Š. Pospíšilová,et al.  Low-burden TP53 mutations in chronic phase of myeloproliferative neoplasms: association with age, hydroxyurea administration, disease type and JAK2 mutational status , 2017, Leukemia.

[31]  R. Sandberg,et al.  Single-cell transcriptomics uncovers distinct molecular signatures of stem cells in chronic myeloid leukemia , 2017, Nature Medicine.

[32]  A. Mead,et al.  Myeloproliferative neoplasm stem cells. , 2017, Blood.

[33]  M. Kurokawa,et al.  Loss of p53 induces leukemic transformation in a murine model of Jak2 V617F-driven polycythemia vera , 2017, Oncogene.

[34]  D. Birnbaum,et al.  Genomic analysis of myeloproliferative neoplasms in chronic and acute phases , 2017, Haematologica.

[35]  T. Druley,et al.  Clonal haematopoiesis harbouring AML-associated mutations is ubiquitous in healthy adults , 2016, Nature Communications.

[36]  Nicola D. Roberts,et al.  Genomic Classification and Prognosis in Acute Myeloid Leukemia. , 2016, The New England journal of medicine.

[37]  N. Beerenwinkel,et al.  Tree inference for single-cell data , 2016, Genome Biology.

[38]  W. Vainchenker,et al.  P53 activation inhibits all types of hematopoietic progenitors and all stages of megakaryopoiesis , 2016, Oncotarget.

[39]  Jinghui Zhang,et al.  Reply to Artifacts in the data of Hu et al. , 2015, Nature Genetics.

[40]  M. Nørgaard,et al.  Epidemiology and Clinical Significance of Secondary and Therapy-Related Acute Myeloid Leukemia: A National Population-Based Cohort Study. , 2015, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[41]  A. Regev,et al.  Spatial reconstruction of single-cell gene expression data , 2015 .

[42]  D. Neuberg,et al.  Acute myeloid leukemia ontogeny is defined by distinct somatic mutations. , 2015, Blood.

[43]  P. Vyas,et al.  Development and Evaluation of the Clinical Utility of a Next Generation Sequencing (NGS) Tool for Myeloid Disorders , 2014 .

[44]  S. Armstrong,et al.  Genomic and functional analysis of leukemic transformation of myeloproliferative neoplasms , 2014, Proceedings of the National Academy of Sciences.

[45]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[46]  Robert Gentleman,et al.  Software for Computing and Annotating Genomic Ranges , 2013, PLoS Comput. Biol..

[47]  Benjamin E. Gross,et al.  The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. , 2012, Cancer discovery.

[48]  G. Getz,et al.  GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers , 2011, Genome Biology.

[49]  P. Vyas,et al.  Coexistence of LMPP-like and GMP-like leukemia stem cells in acute myeloid leukemia. , 2011, Cancer cell.

[50]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[51]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[52]  P. Campbell,et al.  Mutation of JAK2 in the myeloproliferative disorders: timing, clonality studies, cytogenetic associations, and role in leukemic transformation. , 2006, Blood.

[53]  R. Mesa,et al.  Leukemic transformation in myelofibrosis with myeloid metaplasia: a single-institution experience with 91 cases. , 2005, Blood.

[54]  T. Jacks,et al.  Mutant p53 Gain of Function in Two Mouse Models of Li-Fraumeni Syndrome , 2004, Cell.

[55]  M. Wigler,et al.  Circular binary segmentation for the analysis of array-based DNA copy number data. , 2004, Biostatistics.

[56]  Mark Coles,et al.  Transgenic mice with hematopoietic and lymphoid specific expression of Cre , 2003, European journal of immunology.

[57]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[58]  Ira M. Hall,et al.  BEDTools: a flexible suite of utilities for comparing genomic features , 2010, Bioinform..

[59]  A. Olshen,et al.  A Faster Circular Binary Segmentation Algorithm for the Analysis of Array Cgh Data , 2022 .