Prediction of molecular subtypes in acute myeloid leukemia based on gene expression profiling (original) (raw)

Abstract

We examined the gene expression profiles of two independent cohorts of patients with acute myeloid leukemia [n=247 and n=214 (younger than or equal to 60 years)] to study the applicability of gene expression profiling as a single assay in prediction of acute myeloid leukemia-specific molecular subtypes. The favorable cytogenetic acute myeloid leukemia subtypes, i.e., acute myeloid leukemia with t(8;21), t(15;17) or inv(16), were predicted with maximum accuracy (positive and negative predictive value: 100%). Mutations in NPM1 and CEBPA were predicted less accurately (positive predictive value: 66% and 100%, and negative predictive value: 99% and 97% respectively). Various other characteristic molecular acute myeloid leukemia subtypes, i.e., mutant FLT3 and RAS, abnormalities involving 11q23, −5/5q-, −7/7q-, abnormalities involving 3q (abn3q) and t(9;22), could not be correctly predicted using gene expression profiling. In conclusion, gene expression profiling allows accurate prediction of certain acute myeloid leukemia subtypes, e.g. those characterized by expression of chimeric transcription factors. However, detection of mutations affecting signaling molecules and numerical abnormalities still requires alternative molecular methods.

Introduction

Acute myeloid leukemia (AML) is not a singledisease but a group of neoplasms with various genetic abnormalities and variable responses to treatment. The pre-treatment karyotype is still essential in therapy decision-making in AML.1–3 In recent years, a number of novel molecular markers have been associated with AML prognostics.2,3 Several attempts have been made to investigate whether genome-wide GEP could be valuable for prediction of certain subtypes of AML.4–12 Although there was concordance in predictive signatures in the various studies, none of those studies validated the derived signatures to predict the recurrent molecular markers using independent representative AML cohorts. The question, therefore, remains whether GEP could substitute current diagnostic techniques and could be applied as a reliable single test to simultaneously detect known cytogenetic and molecular abnormalities. The aim of this study was to validate GEP as preferred single assay to predict prognostically relevant AML subtypes using two large independent cohorts of young adults with AML.

Design and Methods

Bone marrow aspirates or peripheral blood samples of two independent representative cohorts of de novo AML patients (lower or equal than 60 years), consisting of 247 and 214 patients, were collected (Table 1). The first cohort represents a subset of 285 patients previously studied,8 while the second cohort has not yet been described.

Blast cell purification and RNA isolation were carried out as previously described.8 All samples were analyzed using Affymetrix Human Genome U133Plus2.0 GeneChips (Affymetrix, Santa Clara, CA, USA). Labeling, hybridization, scanning and data normalization were performed as previously described.8 The variation between the scaling/normalization factors of the GeneChips in both cohorts was less than 3-fold [cohort1: 0.53(±0.15); cohort2: 0.73(±0.20)]. Also, the percentage of genes present [cohort1: 39.1(±3.1); cohort2: 40.6(±3.7)], GAPDH 3′/5′ ratio [cohort1: 1.07(±0,13); cohort2: 1.08(±0.16)] and actin 3′/5′ ratio [cohort1: 1.26(±0.21); cohort2: 1.33(±0.29)] were indicative for high overall quality and consistency between both AML sample populations. Mutational analyses to detect recurrent mutations in AML were performed as previously described.13–16 All supervised class prediction analyses were performed with Prediction Analysis for Microarrays (PAM) software version 1.28 in R version 2.1.0.17

Clinical, cytogenetic and molecular information as well as the gene expression profiles of all primary AML cases is available at the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo , accession number GSE6891).

Results and Discussion

In this study of 461 clinically and molecularly well-characterized cases of AML (Table 1), we were able to comprehensively validate the application of GEP to predict therapeutically relevant molecular subtypes in AML.

We applied PAM to investigate whether karyotypic and mutational abnormalities with prognostic or therapeutic value in AML were accurately predictable based on GEP. PAM allows the selection of the minimal number of genes required for optimal prediction, which may be beneficial in a diagnostic setting. The AML cohort1 (n=247) was used as training set to derive predictive signatures that were subsequently validated on AML cohort2 (n=214). The deduced expression signatures are available in the Online Supplementary Tables S1–18.

The cytogenetic status of all AML patients with favorable risk, i.e. those with t(8;21), t(15;17) or inv(16) abnormalities, was predicted with 100 percent accuracy (Table 2). In fact, among these predicted AML cases, there were cases with favorable cytogenetics that had previously been missed by routine cytogenetics (4 out of 37 inv(16) and 4 out of 25 t(15;17)). The presence of the translocation-related fusion transcripts in these specific cases was confirmed by real-time quantitative PCR. Thus, GEP is a reliable alternative to discriminate these three AML subtypes,2,3 which represent approximately 20% of all cases.2,3 Prediction of t(15;17) and inv(16) required only few genes, as seen previously.8 For the t(8;21) cases, 76 probe sets were needed to correctly classify all samples. However, as few as two probe sets, including one associated with the RUNX1T1 (ETO) gene, were sufficient to accurately classify all but one t(8;21) cases, which is also consistent with earlier studies8 (Online Supplementary Figure S3).

AML cases with mutations in the transcription factor CCAAT/enhancer binding protein α (CEBPA), which are associated with a relatively favorable treatment outcome, were predicted with positive and negative predictive values of 100% and 97% respectively. Six out of 15 CEBPA mutant cases were missed in the validation set (sensitivity 60%; Table 2). Of note, the misclassified cases all carried a single heterozygous CEBPA mutation, whereas samples with biallelic mutations (either homo- or heterozygous) were all correctly recognized (data not shown). In the training cohort, all but two (14/16) samples carried biallelic mutations14,18 and in cross-validation in the training cohort the two heterozygous mutants were the only misclassified samples as well.

Previous work has shown that mutations in nucleophosmin (NPM1) are strongly associated with a discriminative _HOX_- and TALE gene-specific signature.16,19 In this study, AML cases carrying a NPM1 mutation were indeed recognized with high accuracy based on such a signature (Table 2 and Online Supplementary Table S5). However, a relatively high number of AML cases without NPM1 mutations was incorrectly predicted positive (32 out of 151), suggesting the presence of genetic alterations resulting in a similar upregulation of the HOX- and TALE genes in those cases. Among these false positives were several AMLs carrying 11q23 abnormalities, which is in line with the role of the mixed lineage leukemia (MLL) protein as an important regulator of HOX gene expression.16 Of note, all t(6;9) AML cases in the training and validation cohort (n=6) were predicted to also carry an NPM1 mutation, raising the possibility that the DEK-CAN fusion protein also induces _HOX_-related gene expression. Interestingly, prediction of t(6;9) translocation was partly feasible using a unique signature (Table2 and Online Supplementary Table S14), although these results are based on a relatively low number of cases.

Table 1.Clinical and molecular data.

Table 2.Class prediction using prediction analysis for microarrays.

NPM1 mutations are associated with relatively favorable survival parameters in patients with a normal karyotype and standard risk AML.16,20–22 The favorable risk is particularly associated with AMLs lacking internal tandem duplications (ITD) in the fms-related tyrosine kinase (FLT3) gene.16,20–22 Analyses of AML subsets defined by combined presence or absence of NPM1 and FLT3 ITD abnormalities demonstrated that only patients carrying both mutations could be moderately predicted, whereas the remaining subtypes could not be discriminated (Table 2). Restriction of these analyses to normal karyotype cases only did not result in a significant improvement in prediction accuracy (Online Supplementary Table S19). Of note, prediction of NPM1 mutation in preselected normal karyotype samples led to a slightly increased positive predictive value (83 vs. 66%), which may be consistent with the lack of interfering 11q23 positive samples. The remaining cytogenetic and molecular subgroups we studied were not associated with strong predictive signatures. Whereas the positive predictive value for FLT3 ITD aberrations was relatively high (77%), the high number of false predictions eliminates GEP, with the currently available analyses tools, as a reliable test to determine the FLT3 ITD status. Restriction to the normal karyotype group did not lead to a marked improvement (Online Supplementary Table S19). Likewise, the low positive predictive values for FLT3 tyrosine kinase domain (TKD) or RAS mutations, abnormalities involving 11q23, −5/5q-, −7/7q- and abn3q, and the translocation t(9;22), disqualify GEP as single detection method for these abnormalities. Similarly, 3q aberrations were not readily predictable. Nevertheless, the most discriminative gene for abn3q abnormalities was the oncogenic transcription factor ecotropic viral integration site1 (EVI1) (Online Supplementary Table S15), which is frequently involved in 3q26 abnormalities. Of note, in these predictions we included the cases carrying a cryptic abn3q recently identified by gene expression analyses and fluorescence in situ hybridization.23

Classifiers were also deduced using a number of other approaches, i.e. compound covariate predictor, linear discriminant analysis, 1-nearest neighbor and 3-nearest neighbors, nearest centroid and support vector machines (probe set selection at 0.001 significance level). These alternative analyses were carried out in BRB-ArrayTools, version 3.7.0 β2 release, developed by Dr. Richard Simon and Amy Peng Lam. Overall, this comparative analysis yielded highly similar results, i.e. the favorable cytogenetic subclasses were predictable with (close to) 100% accuracy, whereas other subtypes showed a similar prediction pattern as depicted in Table 2 (data not shown). One exception was NPM1 mutation status, for which prediction accuracy was better using an approach based on support vector machines (positive predictive value 91% with a negative predictive value of 99%). Several general causes for the inability to predict specific recurrent abnormalities could apply: (i) if different recurrent genetic aberrations affect similar pathways, their GEP signatures may overlap; (ii) mutations affecting signaling pathways may not result in strong discriminative mRNA expression signatures; (iii) the expression of differentiation-related genes may affect accurate prediction; (iv) secondary mutations, or biallelic versus monoallelic mutations as in the case of CEBPA, may prohibit reliable prediction. More specifically, (v) the various partners of the MLL gene may affect reliable prediction of 11q23 abnormalities, and (vi) the numerical changes in (part of) the chromosomes 5 and 7 may only result in minor changes in gene expression that are insufficient for GEP prediction. Of note, still almost all discriminative genes with decreased expression in the deduced signature for 7(q) abnormalities were located on chromosome 7, including FASTK, GSTK1, LSM8 and_ZNF746_ (Online Supplementary Table S17).

Altogether, we conclude that AML cases with favorable cytogenetics are predictable with high accuracy with the currently available genome-wide gene expression technology and analyses tools. All other prognostically and therapeutically known abnormalities in AML still require additional molecular methods for detection.

Footnotes

References

  1. Mrozek K, Heerema NA, Bloomfield CD. Cytogenetics in acute leukemia. Blood Rev. 2004; 18:115-36. Google Scholar
  2. Mrozek K, Bloomfield CD. Chromosome aberrations, gene mutations and expression changes, and prognosis in adult acute myeloid leukemia. Hematology Am Soc Hematol Educ Program. 2006;169-77. Google Scholar
  3. Estey E, Dohner H. Acute myeloid leukaemia. Lancet. 2006; 368:1894-907. Google Scholar
  4. Virtaneva K, Wright FA, Tanner SM, Yuan B, Lemon WJ, Caligiuri MA. Expression profiling reveals fundamental biological differences in acute myeloid leukemia with isolated trisomy 8 and normal cytogenetics. Proc Natl Acad Sci USA. 2001; 98:1124-9. Google Scholar
  5. Schoch C, Kohlmann A, Schnittger S, Brors B, Dugas M, Mergenthaler S. Acute myeloid leukemias with reciprocal rearrangements can be distinguished by specific gene expression profiles. Proc Natl Acad Sci USA. 2002; 99:10008-13. Google Scholar
  6. Debernardi S, Lillington DM, Chaplin T, Tomlinson S, Amess J, Rohatiner A. Genome-wide analysis of acute myeloid leukemia with normal karyotype reveals a unique pattern of homeobox gene expression distinct from those with translocation-mediated fusion events. Genes Chromosomes Cancer. 2003; 37:149-58. Google Scholar
  7. Kohlmann A, Schoch C, Schnittger S, Dugas M, Hiddemann W, Kern W. Molecular characterization of acute leukemias by use of microarray technology. Genes Chromosomes Cancer. 2003; 37:396-405. Google Scholar
  8. Valk PJ, Verhaak RG, Beijen MA, Erpelinck CA, Barjesteh van Waalwijk van Doorn-Khosrovani S, Boer JM. Prognostically useful gene-expression profiles in acute myeloid leukemia. N Engl J Med. 2004; 350:1617-28. Google Scholar
  9. Bullinger L, Döhner K, Bair E, Fröhling S, Schlenk RF, Tibshirani R. Use of gene-expression profiling to identify prognostic subclasses in adult acute myeloid leukemia. N Engl J Med. 2004; 350:1605-16. Google Scholar
  10. Haferlach T, Kohlmann A, Schnittger S, Dugas M, Hiddemann W, Kern W. Global approach to the diagnosis of leukemia using gene expression profiling. Blood. 2005; 106:1189-98. Google Scholar
  11. Vey N, Mozziconacci MJ, Groulet-Martinec A, Debono S, Finetti P, Carbuccia N. Identification of new classes among acute myelogenous leukaemias with normal karyotype using gene expression profiling. Oncogene. 2004; 23:9381-91. Google Scholar
  12. Bullinger L, Rücker FG, Kurz S, Du J, Scholl C, Sander S. Gene-expression profiling identifies distinct subclasses of core binding factor acute myeloid leukemia. Blood. 2007; 110:1291-300. Google Scholar
  13. Valk PJM, Bowen DT, Frew ME, Goodeve AC, Löwenberg B, Reilly JT. Second hit mutations in the RTK/RAS signalling pathway in acute myeloid leukaemia and inv(16). Haematologica. 2004; 89:106. Google Scholar
  14. Barjesteh van Waalwijk van Doorn-Khosrovani S, Erpelinck C, Meijer J, van Oosterhoud S, van Putten WL, Valk PJ. Biallelic mutations in the CEBPA gene and low CEBPA expression levels as prognostic markers in intermediate-risk AML. Hematol J. 2003; 4:31-40. Google Scholar
  15. Care RS, Valk PJ, Goodeve AC, Abu-Duhier FM, Geertsma-Kleinekoort WM, Wilson GA. Incidence and prognosis of c-KIT and FLT3 mutations in core binding factor (CBF) acute myeloid leukaemias. Br J Haematol. 2003; 121:775-7. Google Scholar
  16. Verhaak RG, Goudswaard CS, van Putten W, Bijl MA, Sanders MA, Hugens W. Mutations in nucleophosmin NPM1 in acute myeloid leukemia (AML): association with other gene abnormalities and previously established gene expression signatures and their favorable prognostic significance. Blood. 2005; 106:3747-54. Google Scholar
  17. Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci USA. 2002; 99:6567-72. Google Scholar
  18. Wouters BJ, Jordà MA, Keeshan K, Louwers I, Erpelinck-Verschueren CA, Tielemans D. Distinct gene expression profiles of acute myeloid/T-lymphoid leukemia with silenced CEBPA and mutations in NOTCH1. Blood. 2007; 110:3706-14. Google Scholar
  19. Alcalay M, Tiacci E, Bergomas R, Bigerna B, Venturini E, Minardi SP. Acute myeloid leukemia bearing cytoplasmic nucleophosmin (NPMc+AML) shows a distinct gene expression profile characterized by up-regulation of genes involved in stem cell maintenance. Blood. 2005; 106:899-902. Google Scholar
  20. Schnittger S, Schoch C, Kern W, Mecucci C, Tschulik C, Martelli MF. Nucleophosmin gene mutations are predictors of favorable prognosis in acute myelogenous leukemia with a normal karyotype. Blood. 2005; 106:3733-9. Google Scholar
  21. Falini B, Mecucci C, Tiacci E, Alcalay M, Rosati R, Pasqualucci L. Cytoplasmic nucleophosmin in acute myelogenous leukemia with a normal karyotype. GIMEMA Acute Leukemia Working Party. N Engl J Med. 2005; 352:254-66. Google Scholar
  22. Döhner K, Schlenk RF, Habdank M, Scholl C, Rücker FG, Corbacioglu A. Mutant nucleophosmin (NPM1) predicts favorable prognosis in younger adults with acute myeloid leukemia and normal cytogenetics: interaction with other gene mutations. Blood. 2005; 106:3740-6. Google Scholar
  23. Lugthart S, Drunen EV, Norden YV, van Hoven A, Erpelinck CA, Valk PJ. High EVI1 levels predict adverse outcome in acute myeloid leukemia: prevalence of EVI1 overexpression and chromosome 3q26 abnormalities underestimated. Blood. 2008; 111:4329-37. Google Scholar