Systematic computation with functional gene-sets among leukemic and hematopoietic stem cells reveals a favorable prognostic signature for acute myeloid leukemia

BackgroundGenes that regulate stem cell function are suspected to exert adverse effects on prognosis in malignancy. However, diverse cancer stem cell signatures are difficult for physicians to interpret and apply clinically. To connect the transcriptome and stem cell biology, with potential clinical applications, we propose a novel computational “gene-to-function, snapshot-to-dynamics, and biology-to-clinic” framework to uncover core functional gene-sets signatures. This framework incorporates three function-centric gene-set analysis strategies: a meta-analysis of both microarray and RNA-seq data, novel dynamic network mechanism (DNM) identification, and a personalized prognostic indicator analysis. This work uses complex disease acute myeloid leukemia (AML) as a research platform.ResultsWe introduced an adjustable “soft threshold” to a functional gene-set algorithm and found that two different analysis methods identified distinct gene-set signatures from the same samples. We identified a 30-gene cluster that characterizes leukemic stem cell (LSC)-depleted cells and a 25-gene cluster that characterizes LSC-enriched cells in parallel; both mark favorable-prognosis in AML. Genes within each signature significantly share common biological processes and/or molecular functions (empirical p = 6e-5 and 0.03 respectively). The 25-gene signature reflects the abnormal development of stem cells in AML, such as AURKA over-expression. We subsequently determined that the clinical relevance of both signatures is independent of known clinical risk classifications in 214 patients with cytogenetically normal AML. We successfully validated the prognosis of both signatures in two independent cohorts of 91 and 242 patients respectively (log-rank p < 0.0015 and 0.05; empirical p < 0.015 and 0.08).ConclusionThe proposed algorithms and computational framework will harness systems biology research because they efficiently translate gene-sets (rather than single genes) into biological discoveries about AML and other complex diseases.

[1]  J. Boultwood,et al.  The 5q-syndrome. , 1994, Blood.

[2]  R. Foà,et al.  Interleukin-2 may induce prolonged remissions in advanced acute myelogenous leukemia. , 1994, Blood.

[3]  M. Caligiuri,et al.  A cell initiating human acute myeloid leukaemia after transplantation into SCID mice , 1994, Nature.

[4]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[5]  A. Leong,et al.  Manual of Diagnostic Antibodies for Immunohistology , 1998 .

[6]  Dekang Lin,et al.  An Information-Theoretic Definition of Similarity , 1998, ICML.

[7]  E. Jones.  Manual of Diagnostic Antibodies for Immunohistology , 1999 .

[8]  N. Zhao,et al.  Transcript map and comparative analysis of the 1.5-Mb commonly deleted segment of human 5q31 in malignant myeloid diseases with a del(5q). , 2001, Genomics.

[9]  J. Wainscoat,et al.  The 5 q-Syndrome , 2002 .

[10]  Jeanne Kowalski,et al.  Microarray and Serial Analysis of Gene Expression Analyses Identify Known and Novel Transcripts Overexpressed in Hematopoietic Stem Cells , 2004, Cancer Research.

[11]  D. Koller,et al.  A module map showing conditional activity of expression modules in cancer , 2004, Nature Genetics.

[12]  Andrew I Su,et al.  Uncovering regulatory pathways that affect hematopoietic stem cell function using 'genetical genomics' , 2005, Nature Genetics.

[13]  G. Ehninger,et al.  Activation of the RAS Pathway Is Predictive for a Chemosensitive Phenotype of Acute Myelogenous Leukemia Blasts , 2005, Clinical Cancer Research.

[14]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[15]  A. Levine,et al.  The p53 pathway: positive and negative feedback loops , 2005, Oncogene.

[16]  Derick R. Peterson,et al.  The sesquiterpene lactone parthenolide induces apoptosis of human acute myelogenous leukemia stem and progenitor cells. , 2005, Blood.

[17]  Ulrich Pfeffer,et al.  The Transforming Growth Factor-β Family Members Bone Morphogenetic Protein-2 and Macrophage Inhibitory Cytokine-1 as Mediators of the Antiangiogenic Activity of N-(4-Hydroxyphenyl)Retinamide , 2005, Clinical Cancer Research.

[18]  O. Yli-Harja,et al.  Global Gene Expression Profile of Human Cord Blood–Derived CD133+ Cells , 2006, Stem cells.

[19]  M. Nagasawa,et al.  Phospho-regulation of human protein kinase Aurora-A: analysis using anti-phospho-Thr288 monoclonal antibodies , 2006, Oncogene.

[20]  S. Knuutila,et al.  DNA copy number amplification profiling of human neoplasms , 2006, Oncogene.

[21]  G. Wang,et al.  Identification of genes targeted by the androgen and PKA signaling pathways in prostate cancer cells , 2006, Oncogene.

[22]  S. Takada,et al.  Epigenetic silencing of AXIN2 in colorectal carcinoma with microsatellite instability , 2006, Oncogene.

[23]  Rainer Spang,et al.  Similarities of Ordered Gene Lists , 2006, J. Bioinform. Comput. Biol..

[24]  Holger Fröhlich,et al.  GOSim – an R-package for computation of information theoretic GO similarities between terms and gene products , 2007, BMC Bioinformatics.

[25]  Rainer Spang,et al.  OrderedList - a bioconductor package for detecting similarity in ordered gene lists , 2006, Bioinform..

[26]  C. Goolsby,et al.  NUP98-HOXA9 induces long-term proliferation and blocks differentiation of primary human CD34+ hematopoietic cells. , 2006, Cancer research.

[27]  K. Mills,et al.  Transcriptional dysregulation mediated by RUNX1-RUNX1T1 in normal human progenitor cells and in acute myeloid leukaemia , 2007, Leukemia.

[28]  Satoshi Tanaka,et al.  Chemotherapy-resistant human AML stem cells home to and engraft within the bone-marrow endosteal region , 2007, Nature Biotechnology.

[29]  Kai Wang,et al.  Comparative analysis of microarray normalization procedures: effects on reverse engineering gene networks , 2007, ISMB/ECCB.

[30]  T. Holyoake,et al.  Transcriptional Analysis of Quiescent and Proliferating CD34+ Human Hemopoietic Cells from Normal and Chronic Myeloid Leukemia Sources , 2007, Stem cells.

[31]  Vincent C. Manganiello,et al.  Foxp3-dependent programme of regulatory T-cell differentiation , 2007, Nature.

[32]  J. Gribben,et al.  Anti-CD38 antibody-mediated clearance of human repopulating cells masks the heterogeneity of leukemia-initiating cells. , 2008, Blood.

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

[34]  B. Williams,et al.  Mapping and quantifying mammalian transcriptomes by RNA-Seq , 2008, Nature Methods.

[35]  Bas J. Wouters,et al.  Brief Report Results and Discussion , 2022 .

[36]  Nicolas Mermod,et al.  Nuclear Factor I-C Links Platelet-Derived Growth Factor and Transforming Growth Factor β1 Signaling to Skin Wound Healing Progression , 2009, Molecular and Cellular Biology.

[37]  L. Hood,et al.  Dysregulated gene expression networks in human acute myelogenous leukemia stem cells , 2009, Proceedings of the National Academy of Sciences.

[38]  M. Andreeff,et al.  Cyclin‐dependent kinase 1 inhibitor RO‐3306 enhances p53‐mediated Bax activation and mitochondrial apoptosis in AML , 2009, Cancer science.

[39]  David Haussler,et al.  Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM , 2010, Bioinform..

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

[41]  R. Arceci.  FLT3-ITD up-regulates MCL-1 to promote survival of stem cells in acute myeloid leukemia via FLT3-ITD–specific STAT5 activation , 2010 .

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

[43]  Keji Zhao,et al.  Regulation of nucleosome landscape and transcription factor targeting at tissue-specific enhancers by BRG1. , 2011, Genome research.

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

[45]  Yong-mei Zhu,et al.  Gene mutation patterns and their prognostic impact in a cohort of 1185 patients with acute myeloid leukemia. , 2011, Blood.

[46]  Craig T Jordan,et al.  Leukemia stemness signatures step toward the clinic. , 2011, Cell stem cell.

[47]  E. Vellenga,et al.  Gene expression profiling in the leukemic stem cell-enriched CD34+ fraction identifies target genes that predict prognosis in normal karyotype AML , 2011, Leukemia.

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

[49]  L. Bullinger,et al.  High BRE expression predicts favorable outcome in adult acute myeloid leukemia, in particular among MLL-AF9-positive patients. , 2011, Blood.

[50]  David Venet,et al.  Most Random Gene Expression Signatures Are Significantly Associated with Breast Cancer Outcome , 2011, PLoS Comput. Biol..

[51]  Daniel R. Salomon,et al.  Strategies for aggregating gene expression data: The collapseRows R function , 2011, BMC Bioinformatics.

[52]  Helga Thorvaldsdóttir,et al.  Molecular signatures database (MSigDB) 3.0 , 2011, Bioinform..

[53]  K. Hansen,et al.  Removing technical variability in RNA-seq data using conditional quantile normalization , 2012, Biostatistics.

[54]  Yong Huang,et al.  Towards Mechanism Classifiers: Expression-anchored Gene Ontology Signature Predicts Clinical Outcome in Lung Adenocarcinoma Patients , 2012, AMIA.

[55]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[56]  Yong Huang,et al.  Single Sample Expression-Anchored Mechanisms Predict Survival in Head and Neck Cancer , 2012, PLoS Comput. Biol..

[57]  Q. Cui,et al.  Corrigendum: Identification of high-quality cancer prognostic markers and metastasis network modules , 2012, Nature Communications.

[58]  T. Furey,et al.  Integrating genetic and gene expression evidence into genome-wide association analysis of gene sets. , 2011, Genome research.

[59]  Atul J. Butte,et al.  Ten Years of Pathway Analysis: Current Approaches and Outstanding Challenges , 2012, PLoS Comput. Biol..

[60]  B. Huntly,et al.  Recent advances in acute myeloid leukemia stem cell biology , 2012, Haematologica.

[61]  Kazuyuki Aihara,et al.  Detecting early-warning signals for sudden deterioration of complex diseases by dynamical network biomarkers , 2012, Scientific Reports.

[62]  O. Abdel-Wahab,et al.  The role of mutations in epigenetic regulators in myeloid malignancies , 2012, Nature Reviews Cancer.

[63]  Michal Sheffer,et al.  Pathway-based personalized analysis of cancer , 2013, Proceedings of the National Academy of Sciences.

[64]  Xinan Yang,et al.  Bridging Cancer Biology with the Clinic: Relative Expression of a GRHL2-Mediated Gene-Set Pair Predicts Breast Cancer Metastasis , 2013, PloS one.

[65]  T. Ikezoe,et al.  CD34+/CD38− acute myelogenous leukemia cells aberrantly express Aurora kinase A , 2011, International journal of cancer.

[66]  Ulrich Mansmann,et al.  Identification of a 24-gene prognostic signature that improves the European LeukemiaNet risk classification of acute myeloid leukemia: an international collaborative study. , 2013, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[67]  M. Minden,et al.  Preferential eradication of acute myelogenous leukemia stem cells by fenretinide , 2013, Proceedings of the National Academy of Sciences.

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

[69]  Roberto Romero,et al.  A Comparison of Gene Set Analysis Methods in Terms of Sensitivity, Prioritization and Specificity , 2013, PloS one.

[70]  Rui Liu,et al.  Detecting tissue-specific early warning signals for complex diseases based on dynamical network biomarkers: study of type 2 diabetes by cross-tissue analysis , 2014, Briefings Bioinform..

[71]  T. Somervaille,et al.  The variety of leukemic stem cells in myeloid malignancy , 2014, Oncogene.

[72]  D. Sasaki,et al.  Expression of myeloperoxidase in acute myeloid leukemia blasts mirrors the distinct DNA methylation pattern involving the downregulation of DNA methyltransferase DNMT3B , 2014, Leukemia.

[73]  Xinan Yang,et al.  Computational prognostic indicators for breast cancer , 2014, Cancer management and research.

[74]  K. Yanagihara,et al.  Title Expression of myeloperoxidase in acute myeloid leukemia blasts mirrors the distinct DNA methylation pattern involving the downregulation of DNA methyltransferase DNMT 3 B , 2017 .