Blood Journal
Leading the way in experimental and clinical research in hematology

An FLT3 gene-expression signature predicts clinical outcome in normal karyotype AML

  1. Lars Bullinger1,
  2. Konstanze Döhner1,
  3. Raphael Kranz1,
  4. Christoph Stirner1,
  5. Stefan Fröhling1,
  6. Claudia Scholl1,
  7. Young H. Kim2,
  8. Richard F. Schlenk1,
  9. Robert Tibshirani3,
  10. Hartmut Döhner1,*, and
  11. Jonathan R. Pollack2,*
  1. 1Department of Internal Medicine III, University of Ulm, Ulm, Germany; and
  2. Departments of 2Pathology and
  3. 3Statistics and Health Research & Policy, Stanford University, CA

Abstract

Acute myeloid leukemia with normal karyotype (NK-AML) represents a cytogenetic grouping with intermediate prognosis but substantial molecular and clinical heterogeneity. Within this subgroup, presence of FLT3 (FMS-like tyrosine kinase 3) internal tandem duplication (ITD) mutation predicts less favorable outcome. The goal of our study was to discover gene-expression patterns correlated with FLT3-ITD mutation and to evaluate the utility of a FLT3 signature for prognostication. DNA microarrays were used to profile gene expression in a training set of 65 NK-AML cases, and supervised analysis, using the Prediction Analysis of Microarrays method, was applied to build a gene expression–based predictor of FLT3-ITD mutation status. The optimal predictor, composed of 20 genes, was then evaluated by classifying expression profiles from an independent test set of 72 NK-AML cases. The predictor exhibited modest performance (73% sensitivity; 85% specificity) in classifying FLT3-ITD status. Remarkably, however, the signature outperformed FLT3-ITD mutation status in predicting clinical outcome. The signature may better define clinically relevant FLT3 signaling and/or alternative changes that phenocopy FLT3-ITD, whereas the signature genes provide a starting point to dissect these pathways. Our findings support the potential clinical utility of a gene expression–based measure of FLT3 pathway activation in AML.

Introduction

In acute myeloid leukemia (AML), cytogenetics is used to stratify cases for appropriate risk-adapted therapy.1 The largest subset of patients, the 40% to 49% of AML cases with normal karyotype (NK-AML), comprises a heterogeneous group with intermediate prognosis.2 In the past few years, mutations in genes, such as FLT3, MLL, CEBPA, and NPM1, have been identified in NK-AML, and the presence of such mutations carries prognostic information.2,3 Despite these advances, there is still no clear consensus on how aggressively to treat patients with NK-AML.

Among these molecular alterations, FLT3 mutation has not only been linked to unfavorable outcome but is also a “drugable” target for molecularly directed therapy.4 FLT3 encodes a receptor tyrosine kinase that is preferentially expressed on hematopoietic progenitor cells and mediates stem cell differentiation and proliferation.5 In NK-AML, 28% to 33% of cases harbor an in-frame internal tandem duplication (ITD) mutation within the juxtamembrane (JM) domain of FLT3.6 ITD mutations disrupt the autoinhibitory conformation of the receptor7,8 and promote constitutive activation of FLT3 and downstream effectors, such as RAS and STAT5 (signal transducer and activator of transcription 5).9 The presence of FLT3-ITD mutation in NK-AML predicts shorter remission duration and overall survival.6 Within the group of FLT3-ITD mutations, the size of the ITD10 (but see Ponziani et al11), which varies from a few nucleotides to several hundred, and the ratio of FLT-ITD mutant to FLT3–wild-type (WT) alleles12,13 have each been positively correlated with less favorable outcome.

Approximately 5% to 15% of NK-AML cases carry an activating missense mutation of FLT3, mainly at Asp835 within the activation loop of the tyrosine kinase domain (TKD).14 The most commonly used diagnostic test (EcoRV digest) detects mutations of Asp835 and Ile836, although additional TKD mutations have been described (eg, residues 841 and 842).15,16 In comparison to FLT3-ITD, FLT3-TKD mutations are not strong activators of STAT517 and have not been definitively linked to poor prognosis.2,6,18,19

DNA microarray-based expression profiling also powerfully captures the molecular heterogeneity of cancers and has been applied to build classifiers and clinical-outcome predictors in AML.2023 Regarding FLT3, prior studies have defined gene-expression patterns associated with FLT3 activation in cultured cells24,25 and in clinical specimens,26,27 but none has assessed the clinical relevance of identified signatures. Here, we profile a large set of clinically annotated NK-AML specimens, apply supervised analysis to define gene-expression patterns characterizing FLT3-ITD mutation, and evaluate the utility of a FLT3 signature in NK-AML classification and prognostication.

Methods

Specimens

NK-AML samples were provided by the German-Austrian AML Study Group (AMLSG) and comprised Ficoll gradient purified mononuclear cells (> 80% blasts) from diagnostic (ie, pretreatment) peripheral blood or bone marrow. Written informed consent was obtained from all patients, and the expression-profiling study was approved by the institutional review board of each participating center. Patients received standard-of-care intensified treatment regimens (protocol AML HD98A), which included 2 courses of induction therapy with idarubicin, cytarabine, and etoposide, one consolidation cycle of high-dose cytarabine and mitoxantrone (HAM), followed by random assignment to a late consolidation cycle of HAM versus autologous hematopoietic stem cell transplantation in case no HLA identical family donor was available for allogeneic hematopoietic stem cell transplantation. The training (n = 65; 15 of which were previously profiled20) and test (n = 72) sets represented sequentially collected and profiled NK-AML sample sets. Conventional cytogenetic banding, fluorescence in situ hybridization, and analysis of FLT3 mutations were performed at the central reference laboratory of the AMLSG. FLT3-TKD mutations were identified by polymerase chain reaction (PCR) amplification of exon 20 followed by EcoRV digest. Determination of FLT3-ITD length and ITD/WT ratio were performed by GeneScan analysis according to published protocols.13 MLL-PTD, CEBPA, and NPM1 mutations were scored as previously described.2830 Clinicopathologic characteristics of specimens are summarized in Table 1.

View this table:
Table 1

Clinicopathologic characteristics of sample sets

Gene-expression profiling

RNA was isolated from stored, frozen mononuclear AML-cell pellets using Trizol reagent and RNA quality assessed by gel electrophoresis. We hybridized Cy5-labeled total RNA from AML samples, along with Cy3-labeled universal reference RNA (pooled from 11 cell lines), on cDNA microarrays (manufactured by the Stanford Functional Genomics Facility) containing 39 632 nonredundant cDNA clones, representing 14 613 named genes, 4925 additional UniGene clusters, and 3068 ESTs not mapping to UniGene clusters. Hybridized arrays were imaged using an Axon GenePix 4000B scanner, fluorescence ratios extracted using GenePix software, and data entered into the Stanford Microarray Database for subsequent analysis. The complete microarray dataset is accessible at SMD31 and at the Gene Expression Omnibus (GEO32; accession GSE8043).

Data analysis

Fluorescence ratios were normalized by mean-centering genes for each array and then mean-centering each gene across all arrays within different print runs, to minimize potential print run-specific bias. For all subsequent analyses, we included only the 11 233 cDNAs on the microarray whose expression was both well measured and highly variable among samples. Well-measured genes were defined as those with signal intensity/background greater than 2 for either the Cy5-labeled AML sample or the Cy3-labeled reference sample, in at least 75% of AML samples hybridized. We defined variably expressed genes as those whose expression was higher or lower by a factor of at least 4 from the average expression of all AML samples in at least 2 AML samples. The FLT3-ITD gene expression-based classifier was built using the Prediction Analysis of Microarrays (PAM) method, based on nearest shrunken centroids.33 An optimal predictor set of 25 cDNAs (representing 20 different genes) was determined by minimizing the cross-validation error rate in the training set. The predictor was then evaluated by classifying expression profiles into FLT3-ITD present or absent, from an independent test set of 72 NK-AML cases. Prognostic value of the FLT3 signature was assessed by comparing the clinical outcomes (remission duration, disease-free and overall survival) of the predicted FLT3 groups using Kaplan-Meier survival analysis. Multivariate proportional-hazards analysis was done using the coxph function in the R statistical package.

Quantitative RT-PCR

To validate differential expression of selected FLT3-ITD classifier genes, we performed quantitative reverse transcription–PCR (Q-RT-PCR) analysis using ABI Assays-on-Demand primer pairs and reagents, run on an ABI Prism 7700 instrument. Expression values are reported normalized to GAPDH levels.

Results

To discover gene-expression patterns associated with FLT3-ITD mutation, we used cDNA microarrays to profile expression across approximately 23 000 different genes in a training set of 65 NK-AML cases (Table 1) representing Ficoll-gradient purified blasts from bone marrow or peripheral blood. We then applied a supervised analysis method, PAM, to define a gene-expression classifier of FLT3-ITD status (hereafter, “FLT3 signature”). The optimal classifier comprised 25 cDNA, representing 20 different genes (Figure 1A).

Figure 1

Expression-based classification of FLT3-ITD status in NK-AML training and test sets. Heatmap representations of gene-expression levels of the PAM-derived 20-gene FLT3-ITD predictor, shown for the NK-AML (A) training and (B) independent test sets. Mean-centered log2 gene-expression ratios are depicted by ratio color scale, indicated. Gray denotes poorly measured data. Genes are ordered by fractional contribution to the PAM predictor, and samples by PAM probabilities of FLT3-ITD classification, also displayed graphically above. Molecular alterations are coded by shaded boxes: FLT3-ITD (■), not FLT3-ITD (Embedded Image; FLT3-ITD length > 48 nt (median; ■); ≤ 48 nt (Embedded Image); not FLT3-ITD (Embedded Image); FLT3-ITD/WT ratio > 0.587 (median; ■), FLT3-ITD/WT ratio ≤ 0.587 (Embedded Image); FLT3-TKD (■), not FLT3-TKD (Embedded Image); not available (Embedded Image). Misclassified samples are indicated by gray arrowhead (misclassified as not FLT3-ITD) or black arrowhead (misclassified as FLT3-ITD).

As technical validation of the gene-expression differences, we carried out Q-RT-PCR of 4 representative signature genes (CFH, SOCS2, STON2, and MYO1B) on 20 samples exhibiting a range of expression levels for these genes. For each of the 4 genes, expression levels quantified by microarray and by Q-RT-PCR exhibited good correlation (r = 0.55-0.82; Figure 2).

Figure 2

Validation of expression levels of selected FLT3-ITD classifier genes by Q-RT-PCR. Heatmap representation of gene-expression levels of CFH, SOCS2, STON2, and MYO1B, quantified by microarray or Q-RT-PCR, measured for 20 NK-AML samples (10 each from training and test sets). Samples are ordered by PAM probabilities of FLT3-ITD classification. Q-RT-PCR expression levels shown are normalized to GAPDH levels, and both microarray and Q-RT-PCR expression ratios are log2-transformed and mean centered (ratio color scale indicated). Correlations (r) between microarray and Q-RT-PCR measured expression levels are indicated.

To evaluate classification performance of the FLT3 signature, we profiled gene expression in an additional 72 cases of NK-AML (Table 1). In this independent test set, the FLT3 signature classified cases with FLT3-ITD mutations with an overall accuracy of 81%, with 73% sensitivity and 85% specificity (Figure 1B). None of the cases misclassified as FLT3-ITD carried FLT3-TKD mutations. For the subset of FLT3-ITD cases misclassified as not FLT3-ITD, there was no significant difference in the length of the ITD (P = 1, Fisher exact test; combined training and test set, comparing proportion > vs ≤ median) or the FLT3-ITD/WT ratio (P = .52). Expression levels for FLT3, measured by microarray (Figure 1), were also not significantly different in cases correctly or incorrectly classified as FLT3-ITD (P = .58; Mann-Whitney U test), or as not FLT3-ITD (P = .16).

Although classification performance was modest, we nevertheless assessed the potential prognostic value of the FLT3 signature by Kaplan-Meier analysis. The FLT3 signature was a significant predictor of early relapse (P = .003, log-rank test; hazard ratio = 3.6; 1.5-8.8; Figure 3B) as well as unfavorable disease-free (P < .001, hazard ratio = 3.3; 1.6-6.6; Figure 3D) and overall survival (P = .001, hazard ratio = 3.9; 1.5-10.6; Figure 3F), and, remarkably, for all 3 outcome measures outperformed FLT3-ITD mutation status (Figure 3A,C,E). In multivariate analysis along with other prognostic factors, the FLT3 signature was the strongest prognosticator (Table 2).

Figure 3

FLT3-ITD classifier predicts clinical outcome. Kaplan-Meier survival analysis of NK-AML independent test set for (A) FLT3-ITD mutation and remission duration, (B) PAM-predicted FLT3-ITD mutation and remission duration, (C) FLT3-ITD mutation and disease-free survival, (D) PAM-predicted FLT3-ITD mutation and disease-free survival, (E) FLT3-ITD mutation and overall survival, and (F) PAM-predicted FLT3-ITD mutation and overall survival. For the PAM prediction groups, samples not evaluated for FLT3-ITD status were excluded. P values (log rank test) are indicated.

View this table:
Table 2

Multivariate proportional hazards analysis

Of note, none of the 6 cases misclassified as FLT3-ITD carried MLL-partial tandem duplication (PTD) (an alteration associated with poor prognosis), and cases misclassified (compared with correctly classified) as not FLT3-ITD were not enriched for CEBPA or NPM1 mutations (associated with good prognosis; P = 1, Fisher exact test). For misclassified cases, there were also no significant differences in patient age and serum lactate dehydrogenase levels (P values .08-.59, Mann-Whitney U test), nor in the proportion of patients treated by allogeneic stem cell transplant (P = 1, Fisher exact test). The only significant association found was between cases misclassified as FLT3-ITD and lower (ie, nonelevated) white blood cell (WBC) counts (P = .04, Mann-Whitney U test). Interestingly, the FLT3 signature prognostically stratified both the FLT3-ITD-positive and FLT3-ITD-negative subgroups, although only the latter reached statistical significance (Figure 4).

Figure 4

FLT3 signature prognostically stratifies both FLT3-ITD+ and FLT3-ITD subgroups. Kaplan-Meier overall survival analysis of NK-AML independent test set for PAM-predicted FLT3-ITD mutation in (A) FLT3-ITD+ subgroup and (B) FLT3-ITD subgroup. P values (log rank test) are indicated.

To assess the robustness of our findings, we performed 50 random training test splits of the combined 137 NK-AML samples. For each random split, an optimal PAM FLT3-ITD classifier was built using the training set and evaluated on the test set. On average, the FLT3 signature (median P = .001, hazard ratio = 7.4) outperformed FLT3-ITD mutation status (median P = .003, hazard ratio = 2.7) in predicting overall survival. The FLT3 signature also exhibited relative stability, with 8 of the 20 original signature genes occurring within the FLT3 classifier of at least 50% of the 50 random training test splits and 14 of the original 20 genes found in at least 40% of the splits.

Discussion

By expression profiling and PAM analysis, we have defined a 20-gene FLT3-ITD classifier signature. Notably, despite modest FLT3-ITD classification accuracy, the 20-gene classifier outperformed FLT3-ITD mutation status in predicting disease relapse as well as disease-free and overall survival for patients with NK-AML, a finding not explainable by the presence of other known prognostic mutations (MLL-PTD, CEBP1, or NPM1), by other prognostic factors (age, lactate dehydrogenase, WBC) or by a spurious association with treatment. (Note, the one significant association found between cases misclassified as FLT3-ITD and lower WBC count (P = .04) would lose significance after correcting for multiple hypothesis testing; regardless, because lower WBC counts are associated with good prognosis any confounding effect would be in the wrong direction.) Rather, we speculate that the superior prognostic performance of the FLT3 signature reflects one or both of 2 possibilities. First, FLT3-ITD cases classified falsely by the signature as not having ITD might be associated with quantitatively or qualitatively less FLT3 pathway activation. Because systematic differences in the length of ITD, the ITD/WT ratio, and FLT3 mRNA levels were not observed, other alterations in FLT3 or downstream pathway genes would presumably account for the observed decreased FLT3 signature.

A second possibility is that cases classified falsely as FLT3-ITD might harbor alternative genetic or epigenetic alterations that phenocopy FLT3-ITD and have similar downstream gene expression read-outs and clinical impact. Because FLT3-TKD mutations were not observed among these cases and there was no significant difference in FLT3 mRNA levels, other alterations in FLT3 or other genes would have to account for the FLT3 pathway activation signature. Other alterations in FLT3 might include JM domain point mutations3436 as well as a recently identified mutation in the N-terminal lobe of the TKD (K663Q),37 which have similar signaling characteristics, at least in vitro. Alterations in other genes might include recently described CBL-inactivating mutations, resulting in up-regulated wild-type FLT3 receptor signaling in AML.38,39 Because the FLT3 signature prognostically stratified both the FLT3-ITD-positive and -negative subgroups, probably both the above classes of possibility are operative. Regardless, our data suggest that a gene expression-based read-out might be a more accurate measure of clinically meaningful FLT3 pathway activation.

The genes included in the FLT3 signature may provide insight into FLT3-ITD–associated leukemogenesis. Some of the genes have known links to AML. For example, SOCS2 (Suppressor of cytokine signaling 2) was previously identified among genes up-regulated in 32Dc23 murine leukemia cells stably transfected with FLT3-ITD.24 SOCS2 is a known target gene induced by STAT5 and functions as a negative regulator of cytokine signaling.40 HOX genes, such as FLT3 signature genes HOXB2, HOXB5, and PBX, are frequently overexpressed in AML and have been linked to leukemogenesis.41 Other signature genes have no prior link to leukemia but have plausible leukemogenic roles based on known functions. For example, DPPA4 was recently identified in a shRNA screen for genes required for self-renewal in embryonic stem cells42 and may have a similar function in leukemic cells. LAPTM4b (lysosome-associated protein transmembrane 4 beta), up-regulated in several solid tumor types,43 has been shown to transform NIH3T3 cells.44 Other signature genes have no obvious link to FLT3 pathway activation but may suggest new areas of investigation. For example, CFH (complement factor H), one of the top FLT3-ITD discriminator genes, encodes a serum glycoprotein that regulates the alternative complement pathway but has not previously been linked to cancer pathogenesis. CFHR1 is a transcribed pseudogene of CFH.45

Other prognostic signatures have been reported for NK-AML. By unsupervised analysis, Valk et al described clinically relevant AML subgroups with distinct gene-expression patterns, including 3 subgroups comprising mostly NK-AML cases.21 In our own prior studies, we used a semisupervised approach46 to train a clinical-outcome predictor in AML (including cases with cytogenetic aberration).20 The resultant 133-gene signature predicted outcome independent of FLT3-ITD and was also prognostic within the NK-AML subgroup.20 The 133-gene signature was recently validated in another NK-AML cohort but was found not independent of FLT3-ITD.47 We have applied a similar semisupervised approach (training on outcomes rather than FLT3-ITD) in the current NK-AML dataset, finding prognostic signatures (L.B., J.R.P., unpublished observation, December 2006) with substantial overlap with the FLT3 signature, suggesting that FLT3-ITD is a major component of prognostically relevant gene expression in NK-AML. Interestingly, in breast cancer, where multiple prognostic signatures have also been described, a recent study revealed concordance among individual patient outcome predictions, suggesting the different signatures are tracking a shared underlying biology.48 Future studies should clarify which of the proposed NK-AML outcome signatures might ultimately translate to clinical utility, using either microarray-based testing or other platforms.

Our finding that a gene expression–based FLT3-ITD mutation classifier is a better prognosticator than mutation status itself is reminiscent of a recent study of TP53 mutation. Miller et al49 defined a gene-expression signature of TP53 mutation in breast tumors. Although TP53 mutation is known to be prognostic, the TP53 signature was a stronger prognosticator and predictor of therapeutic response. Our study provides further support that gene-expression signatures might in general be more effective reporters of pathway status, by integrating downstream transcriptional responses of pathway activation.

Our study indicates a potential utility of the FLT3 signature for improved prognostication and risk-assessed treatment stratification for patients with NK-AML. Validation of prognostic performance in additional patient cohorts is warranted. In addition, and although speculative, the FLT3 signature might provide utility in guiding therapy with selective FLT3 inhibitors, several of which are now being evaluated in clinical trial.4 Specifically, the FLT3 signature might inform the selection of patients likely (or unlikely) to respond, where patients with the activated signature might benefit from treatment in the absence of FLT3 mutation and vice versa. Measuring the FLT3 signature might also be useful in monitoring the therapeutic response in patients treated with FLT3 inhibitors. Our findings provide a starting point to further dissect FLT3-ITD activation pathways and support the potential clinical utility of a gene expression–based measure of FLT3 pathway activation in patients with NK-AML.

Authorship

Contribution: L.B. and J.R.P. designed and performed research, analyzed and interpreted data, and wrote the paper; K.D. designed the research, contributed vital new reagents, and analyzed and interpreted data; R.K., C. Stirner, C. Scholl, and Y.H.K. performed research and analyzed and interpreted data; S.F. performed research, analyzed and interpreted data, and wrote the paper; R.F.S. contributed vital new reagents, analyzed and interpreted data, and wrote the paper; R.T. and H.D. designed the research, analyzed and interpreted data, and wrote the paper.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Jonathan R. Pollack, Department of Pathology, Stanford University School of Medicine, 269 Campus Drive, CCSR-3245A, Stanford, CA 94305-5176; e-mail: pollack1{at}stanford.edu.

Acknowledgments

The authors thank the Stanford Functional Genomics Facility for microarray manufacture, the Stanford Microarray Database for database support, the members of the German-Austrian AMLSG for their continuous support of the treatment protocols, and Keyan Salari for helpful discussion.

This work was supported by grants from the Leukemia & Lymphoma Society (6151-06; J.R.P., H.D.) and Else-Kröner-Fresenius-Stiftung (L.B., K.D.).

Footnotes

  • *H.D. and J.R.P. share senior authorship.

  • The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

  • Submitted September 25, 2007.
  • Accepted February 25, 2008.

References

View Abstract