A 17-gene stemness score for rapid determination of risk in acute leukaemia

Refractoriness to induction chemotherapy and relapse after achievement of remission are the main obstacles to cure in acute myeloid leukaemia (AML). After standard induction chemotherapy, patients are assigned to different post-remission strategies on the basis of cytogenetic and molecular abnormalities that broadly define adverse, intermediate and favourable risk categories. However, some patients do not respond to induction therapy and another subset will eventually relapse despite the lack of adverse risk factors. There is an urgent need for better biomarkers to identify these high-risk patients before starting induction chemotherapy, to enable testing of alternative induction strategies in clinical trials. The high rate of relapse in AML has been attributed to the persistence of leukaemia stem cells (LSCs), which possess a number of stem cell properties, including quiescence, that are linked to therapy resistance. Here, to develop predictive and/or prognostic biomarkers related to stemness, we generated a list of genes that are differentially expressed between 138 LSC+ and 89 LSC− cell fractions from 78 AML patients validated by xenotransplantation. To extract the core transcriptional components of stemness relevant to clinical outcomes, we performed sparse regression analysis of LSC gene expression against survival in a large training cohort, generating a 17-gene LSC score (LSC17). The LSC17 score was highly prognostic in five independent cohorts comprising patients of diverse AML subtypes (n = 908) and contributed greatly to accurate prediction of initial therapy resistance. Patients with high LSC17 scores had poor outcomes with current treatments including allogeneic stem cell transplantation. The LSC17 score provides clinicians with a rapid and powerful tool to identify AML patients who do not benefit from standard therapy and who should be enrolled in trials evaluating novel upfront or post-remission strategies.

[1]  M. Carroll,et al.  Human acute myelogenous leukemia stem cells are rare and heterogeneous when assayed in NOD/SCID/IL2Rγc-deficient mice. , 2011, The Journal of clinical investigation.

[2]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[3]  G. Ehninger,et al.  Long-term prognosis of acute myeloid leukemia according to the new genetic risk classification of the European LeukemiaNet recommendations: evaluation of the proposed reporting system. , 2011, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[4]  Kai Fu,et al.  Determining cell-of-origin subtypes of diffuse large B-cell lymphoma using gene expression in formalin-fixed paraffin-embedded tissue. , 2014, Blood.

[5]  Mei Yu,et al.  PERT: A Method for Expression Deconvolution of Human Blood Samples from Varied Microenvironmental and Developmental Conditions , 2012, PLoS Comput. Biol..

[6]  Bob Löwenberg,et al.  Resistance Prediction in AML: Analysis of 4,601 Patients from MRC/NCRI, HOVON/SAKK, SWOG, and MD Anderson Cancer Center , 2014, Leukemia.

[7]  Benjamin J. Raphael,et al.  Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. , 2013, The New England journal of medicine.

[8]  Benjamin M. Bolstad,et al.  affy - analysis of Affymetrix GeneChip data at the probe level , 2004, Bioinform..

[9]  R. Hills,et al.  Refinement of cytogenetic classification in acute myeloid leukemia: determination of prognostic significance of rare recurring chromosomal abnormalities among 5876 younger adult patients treated in the United Kingdom Medical Research Council trials. , 2010, Blood.

[10]  Bob Löwenberg,et al.  Review Articles (434 articles) , 2008 .

[11]  Andrew P. Feinberg,et al.  An LSC epigenetic signature is largely mutation independent and implicates the HOXA cluster in AML pathogenesis , 2015, Nature Communications.

[12]  Christopher A. Miller,et al.  Association Between Mutation Clearance After Induction Therapy and Outcomes in Acute Myeloid Leukemia. , 2015, JAMA.

[13]  Y. Kanda,et al.  Investigation of the freely available easy-to-use software ‘EZR' for medical statistics , 2012, Bone Marrow Transplantation.

[14]  N. Friedman,et al.  Densely Interconnected Transcriptional Circuits Control Cell States in Human Hematopoiesis , 2011, Cell.

[15]  Angelo J. Canty,et al.  Stem cell gene expression programs influence clinical outcome in human leukemia , 2011, Nature Medicine.

[16]  Trevor Hastie,et al.  Regularization Paths for Cox's Proportional Hazards Model via Coordinate Descent. , 2011, Journal of statistical software.

[17]  G. Sauvageau,et al.  RNA-Seq Reveals Spliceosome and Proteasome Genes as Most Consistent Transcripts in Human Cancer Cells , 2013, PloS one.

[18]  Ulrich Mansmann,et al.  An 86-probe-set gene-expression signature predicts survival in cytogenetically normal acute myeloid leukemia. , 2008, Blood.

[19]  H. Dombret,et al.  Effect of gemtuzumab ozogamicin on survival of adult patients with de-novo acute myeloid leukaemia (ALFA-0701): a randomised, open-label, phase 3 study , 2012, The Lancet.

[20]  L. Bullinger,et al.  Gene expression profiling in AML with normal karyotype can predict mutations for molecular markers and allows novel insights into perturbed biological pathways , 2010, Leukemia.

[21]  R. Gray A Class of $K$-Sample Tests for Comparing the Cumulative Incidence of a Competing Risk , 1988 .

[22]  R Simon,et al.  A non-parametric graphical representation of the relationship between survival and the occurrence of an event: application to responder versus non-responder bias. , 1984, Statistics in medicine.

[23]  Yurii S. Aulchenko,et al.  PredictABEL: an R package for the assessment of risk prediction models , 2011, European Journal of Epidemiology.

[24]  Xavier Robin,et al.  pROC: an open-source package for R and S+ to analyze and compare ROC curves , 2011, BMC Bioinformatics.

[25]  Sean Ferree,et al.  Analytical validation of the PAM50-based Prosigna Breast Cancer Prognostic Gene Signature Assay and nCounter Analysis System using formalin-fixed paraffin-embedded breast tumor specimens , 2014, BMC Cancer.

[26]  Nathan Mantel,et al.  Evaluation of Response-Time Data Involving Transient States: An Illustration Using Heart-Transplant Data , 1974 .

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

[28]  Jian-Bing Fan,et al.  The transcriptional architecture of early human hematopoiesis identifies multilevel control of lymphoid commitment , 2013, Nature Immunology.

[29]  R. Schlenk,et al.  The European LeukemiaNet AML Working Party consensus statement on allogeneic HSCT for patients with AML in remission: an integrated-risk adapted approach , 2012, Nature Reviews Clinical Oncology.

[30]  E. Estey,et al.  Addition of gemtuzumab ozogamicin to induction chemotherapy in adult patients with acute myeloid leukaemia: a meta-analysis of individual patient data from randomised controlled trials. , 2014, The Lancet. Oncology.

[31]  R. Gill,et al.  Cox's regression model for counting processes: a large sample study : (preprint) , 1982 .

[32]  Jennifer L. Osborn,et al.  Direct multiplexed measurement of gene expression with color-coded probe pairs , 2008, Nature Biotechnology.

[33]  Gary D Bader,et al.  miR-126 Regulates Distinct Self-Renewal Outcomes in Normal and Malignant Hematopoietic Stem Cells , 2016, Cancer cell.

[34]  Robert Gray,et al.  A Proportional Hazards Model for the Subdistribution of a Competing Risk , 1999 .

[35]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[36]  Bas J. Wouters,et al.  Prediction of molecular subtypes in acute myeloid leukemia based on gene expression profiling , 2009, Haematologica.

[37]  O. Ohara,et al.  Identification of Therapeutic Targets for Quiescent, Chemotherapy-Resistant Human Leukemia Stem Cells , 2010, Science Translational Medicine.

[38]  C. Bloomfield,et al.  Revised recommendations of the International Working Group for Diagnosis, Standardization of Response Criteria, Treatment Outcomes, and Reporting Standards for Therapeutic Trials in Acute Myeloid Leukemia. , 2003, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[39]  Sean C. Bendall,et al.  Data-Driven Phenotypic Dissection of AML Reveals Progenitor-like Cells that Correlate with Prognosis , 2015, Cell.

[40]  F. Ferrara,et al.  Acute myeloid leukaemia in adults , 2013, The Lancet.

[41]  D. Grimwade,et al.  Molecular landscape of acute myeloid leukemia in younger adults and its clinical relevance. , 2016, Blood.

[42]  Ø. Bruserud,et al.  SIRT1 activation by a c-MYC oncogenic network promotes the maintenance and drug resistance of human FLT3-ITD acute myeloid leukemia stem cells. , 2014, Cell stem cell.

[43]  Ash A. Alizadeh,et al.  Association of a leukemic stem cell gene expression signature with clinical outcomes in acute myeloid leukemia. , 2010, JAMA.

[44]  Pan Du,et al.  lumi: a pipeline for processing Illumina microarray , 2008, Bioinform..

[45]  A. Krešo,et al.  Evolution of the cancer stem cell model. , 2014, Cell stem cell.

[46]  R. Myers,et al.  Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data , 2005, Nucleic acids research.

[47]  Mark A. Dawson,et al.  BET inhibitor resistance emerges from leukaemia stem cells , 2015, Nature.