Acute Myeloid Leukemia |
1 Institute of Hematology, University of Foggia, Foggia, Italy
2 Istituto Applicazioni del Calcolo "M. Picone", Consiglio Nazionale delle Ricerche (CNR), Rome, Italy
3 Dept. Surg. Med. Spec. and Public Health, University of Perugia, Italy
4 Department of Internal Medicine III, University of Ulm, Ulm, Germany
5 Institute of Hematology, University of Tor Vergata, Rome, Italy
6 Laboratory for Molecular Diagnostics, University Hospital Carl Gustav Carus, Dresden, Germany
7 MLL–Munich Leukemia Laboratory, Munich, Germany
8 Department of Hematology, Erasmus University Medical Center, Rotterdam, The Netherlands
9 Institute of Hematology, University of Perugia, Perugia, Italy
10 Department of Hematology, Hemostasis and Oncology, Hannover Medical School, Hannover, Germany
11 Institute of Hematology, University "La Sapienza", Rome, Italy
12 Computational Biology Center, Memorial Sloan Kettering Cancer Center, New York, NY, USA
Correspondence: Arcangelo Liso, Institute of Hematology, University of Foggia, Foggia, Italy. Brunangelo Falini, Institute of Hematology, University of Perugia, Italy. E-mail:a.liso{at}medicina.unifg.it or E-mail:faliniem{at}unipg.it
|
|
|---|
Key words: acute myeloid leukemia, nucleophosmin, mutation.
|
|
|---|
To investigate the age-specific incidence of AML with mutated NPM1, we adapted the one-mutation model that was originally designed to describe chronic myeloid leukemia age distribution.20 The model fits the NPMc+ AML age-specific incidence curve assuming plausible parameter values, supporting the hypothesis that a single genetic event, the NPM1 mutation, is sufficient to cause leukemia. The role of NPM1 mutations in AML development is discussed in the light of these findings.
|
|
|---|
A total of 1,444 AML patients (age range: 20–59; median 47) carrying a mutated NPM1 gene were included in this study (n=476 from AMLSG; n=354 from GIMEMA; n=251 from DSIL; n=223 from MLL; and n=140 from The Netherlands). Exclusion criteria were: i) patients under 20 years of age because few cases were available, due to the low frequency of NPM1 mutations in this age group22; and ii) patients over 59 years of age who are often treated in local hospitals. Consequently, those patients referred to major institutions for diagnosis and treatment may not be representative of the population of AML patients in this age group.
Information on FLT3 status was available in 1,386/1,444 AML patients with mutated NPM1 (96%). FLT3-ITD was detected in 553/1,386 cases (40%). For analysis, the 1,444 NPM1-mutated AML patients were stratified in 5-year age classes.
For this study, we assume that the mutational event needed to develop NPMc+ AML occurs independently of local exposure to environmental leukemogenic factors and that the age specific rates of NPM1-mutated AML patients 20–59 years in age reflect those of the general population in the three European countries included in the study.
Modeling age specific incidence
The mathematical model by Michor et al.20 was adapted to analyze the AML incidence data. Our model is based on the following considerations: (i) we consider a population of N hematopoietic stem cells. Initially, all cells are wild type and proliferate according to a stochastic process known as the Moran model:23 every
days, a cell is chosen at random proportional to fitness to divide; its offspring replaces another randomly chosen cell. The population size is strictly constant; (ii) a wild-type cell gives rise to a mutated cell at rate u per cell division. A mutated cell has a relative growth rate (fitness) of r. If r=1, the mutation is neutral as compared to wild type cells; if r<1, the mutant is disadvantageous, and if r>1, the mutant has a proliferation advantage over the wild type cell. We assume that an NPM1 mutation confers a fitness advantage to the cell, r>1; (iii) Our model adheres to standard Moran process until a surviving mutant cell appears; thereafter, clonal growth is initiated that continues until the mutated cell population reaches population size
. Unlike the model designed by Michor et al.,20 which assumes a constant population size of N cells, our model allows the mutant clone to expand until a maximal size,
. This change is intended to account for the marked expansion of the initial cell compartment which is peculiar to AML; (iv) the AML detection rate is proportional to the number of mutated cells present; if there are Nm mutated cells, the rate of diagnosis is q Nm. From assumptions (i) and (ii) it follows that the waiting time for the first successful (=surviving) mutation has a negative exponential distribution, b=Nu (1-1/r). Let a be the time since the occurrence of the first surviving mutation. Then assumption (3) states that the number of mutated cells, Nm, grows according to
![]() |
where c=(r-1)/
and Nm (0)=1/(1-1/r). To account for the significant expansion of the mutated clone, we assume
>>N. Finally, if (q) is the proportionality constant between the rate of detection and the number of mutated cells (assumption iv), then the probability of diagnosis20 at time t is given by
![]() | (1) |
We compared the predictions of equation (1) with the direct computer simulation of the stochastic process. The simulation is performed by first determining the time at which the first surviving mutated cell arises in a population of N wild type cells; this time follows a negative exponential distribution with mean 1/b. Once such a cell has emerged, the branching process of clonal expansion is simulated by choosing a cell for reproduction or for death at random at each time step. The probability that the number of wild type cells, N, increases by one is given by
![]() | (2a) |
where
= (1+d)N + (r+d)Nm. Here d denotes the death rate of both wild type and mutated cells. The probability that the number of mutated cells, Nm, increases by one is given by
![]() | (2b) |
The probabilities that the numbers of wild type and mutated cells decrease by one are respectively given by
![]() | (2c) |
![]() | (2d) |
A patient is diagnosed at rate qNm and is entered into the incidence data base of his age class. Online Supplementary Figure 1 shows the fit of equation (1) and system (2). Under particular circumstances, i.e. when the waiting time for the first successful mutation is long and clonal expansion occurs fast and reaches large cell numbers, the incidence data can be a kinked curve. A more detailed mathematical investigation of such situations is forthcoming (Michor F. et al., in preparation) but will not be discussed here since the experimentally determined incidence data is a straight line on a doubly logarithmic plot.
Finally, we compared equation (1) with the experimental data, which allowed us to quantify AML-specific parameters.
Statistical analysis
The
2 test (
<0.05) was used to assess independence of the age distribution of cases by center of diagnosis. The likelihood ratio test, comparing a Poisson regression model including age, country, and age x country interaction terms with the nested model without the interaction term was performed to evaluate dependence of age specific NPM1-mutated AML rates on the country.
|
|
|---|
Chi-square testing of the age distribution of cases on center of diagnosis was not significant (p=0.48) indicating that, although absolute incidence levels vary because they reflect the percentage of the general population that is covered by participating centers, number of cases by age class does not differ among study centers (Table 1). The likelihood ratio test comparing the Poisson model which includes a country x age class interaction with the simpler model without the interaction term (Online Supplementary Table 1) was non-significant (p=0.85). Together, these findings provide evidence that AML data from the three countries are comparable. In particular, AML incidence curves analyzed via linear regression all showed a slope of about 4 on a log-log scale (Figure 1).
|
View this table: [in a new window] [Download PPT slide] |
Table 1. Cumulative incidence of the NPM1-mutated acute myeloid leukemia in Germany (DE), the Netherlands (NL) and Italy (IT).*
|
![]() View larger version (12K): [in a new window] [Download PPT slide] |
Figure 1. Incidence data for the NPMc+ acute myeloid leukemia in Germany (DE), The Netherlands (NL) and Italy (IT). Linear regression shows a slope of about 4 on a doubly logarithmic scale.
|
2 and p values for the three countries were 0.02027 (p=0.9899), 0.00862 (p=0.9956), and 0.15275 (p=0.9264). The fitting procedure provided estimates of the parameters in each country (Table 2) generating numbers that are biologically plausible (see below).
![]() View larger version (10K): [in a new window] [Download PPT slide] |
Figure 2. The incidence data can be fit to the one-mutation model assuming plausible parameter values for Germany (A), The Netherlands (B) and Italy (C). 2 and corresponding p-values reported as labels in the corresponding figures show the quality of the fit.
|
|
View this table: [in a new window] [Download PPT slide] |
Table 2. Parameters of the one-mutation model for all NPM1 mutated acute myeloid leukemias.
|
The maximum number of mutated cells generated by the model was about
=1013. This number is consistent with the high tumor burden observed in leukemia patients, if one assumes that, under physiological conditions, the amount of human nucleated marrow cells per kg body weight has been calculated to be approximately 2.1x1010 (1.5x1012 in a subject of 70 kg).25 The relative fitness of mutated cells spanned the range 1.38–1.61. The mean cell generation time (i.e. the time needed for a cell to divide), was between 2.67 and three days, which concurs with early experimental findings26 and with clinical data.3 In the NPM1-mutated AML case, the rate of cancer detection per mutated cell was found to be in the range of 7.77x10–5–1.58x10–4 days. This implies that the total rate of detection (qN) is in the range of 0.26–1.78, which is higher than previous estimates in chronic myeloid leukemia.20 Leukemic clones are initiated by single NPM1 mutations occurring at rates ranging from 2.43x10–9 to 4.86x10–9 days per cell division. Taken together, these estimates imply that for a single individual the waiting time for the appearance of a surviving mutation is on average 1/(Nu(1–1/r)), which is about 5532, 9940 and 8779 days for Germany, The Netherlands and Italy respectively.
FLT3 gene status does not influence the age specific incidence of acute myeloid leukemia with mutated NPM1
Internal tandem duplication (ITD) at the FLT3 gene locus has been implicated as a cooperating genetic alteration in various AML subtypes.27,28 Since FLT3-ITD frequently associates with NPM1 mutations1 and appears to abrogate the favorable prognostic effect of NPM1 mutations in AML1,15,29 we determined whether the age incidence of NPM1-mutated AMLs with FLT3-ITD differs from cases with wild-type FLT3. No significant difference emerged in the slopes of FLT3-ITD-positive and -negative AML with mutated NPM1 (Figure 3). The quality of fit with the model-generated data was adequate and similar to the quality of fit for all AMLs with NPM1 mutations (Figure 4). The one-mutation model parameters for fitting FLT3-ITD positive and FLT3-ITD negative AML with mutated NPM1 are reported in Table 3. The slopes of the three groups (NPM1 mutated, NPM1 mutated/FLT3-ITD, NPM1-mutated/FLT3 wild-type) are not significantly different according to the Mann-Witney U test (p>0.05) (Online Supplementary Table 2).
![]() View larger version (13K): [in a new window] [Download PPT slide] |
Figure 3. Incidence data for the NPMc+ acute myeloid leukemia bearing FLT3-ITD (A) or FLT3 wild-type for Germany (DE), The Netherlands (NL) and Italy (IT) (B).
|
![]() View larger version (10K): [in a new window] [Download PPT slide] |
Figure 4. Fitting the incidence data with the one-mutation model (A) leads to plausible parameter values in all cases (NPMc+ acute myeloid leukemia with FLT3-ITD or FLT3 wild-type) for all three countries (Germany, The Netherlands (B) and Italy (C)). 2 and corresponding p-values reported as labels in the corresponding figures show the quality of the fit.
|
|
View this table: [in a new window] [Download PPT slide] |
Table 3. Parameters of the one-mutation model for the NPM1-mutated/FLT3-ITD and NPM1-mutated/FLT3 wild-type (wt) subsets.
|
|
|
|---|
Our findings add to the body of evidence that NPM1 mutation is a founder genetic lesion in NPMc+ AML: i) cytoplasmic mutated nucleophosmin is specific for AML1,30,31 and clinically shows close association with AML of de novo origin1,32–34; ii) all NPM1 mutations generate changes at the C-terminus of nucleophosmin protein which appear to maximise nuclear export of NPM leukemic mutants,3,35–37 pointing to cytoplasmic dislocation of the mutants as the central event for leukemogenesis; iii) NPM1 mutations are mutually exclusive with other recurrent genetic abnormalities,1,38 with the exception of rare cases in which both NPM1 and CEPBA (or FLT3-ITD) mutations are found;15 iv) they are stable during the course of the disease39,40 as the same type of NPM1 mutation is consistently detected at relapse in medullary and extramedullary sites;40 and v) quantitative real-time PCR shows that NPM1 mutations disappear at complete remission.41,42
The major finding in the present study is that the one-mutation mathematical model can explain the age specific incidence in NPMc+ AML. This hypothesis is in contrast to current concepts in AML development which, like other human cancers, is believed to be a consequence of more than one oncogenic hit.43 Indeed, several animal models of AML clearly point to leukemogenesis as a multi-step process.43 Moreover, in vitro findings that the NPM1 leukemic mutant specifically cooperates with the E1A adenovirus to transform primary MEFs in soft agar44 suggest that NPM1 mutations need to act in close concert with other oncogenic hits. In MEF cells, this mutual cooperation involves the NPM1 mutant inhibiting the E1A-elicited p19(Arf) induction and E1A overcoming NPM1 mutant-induced cellular senescence.44 Furthermore, an activating mutation of the FLT3 gene (FLT3-ITD) leading to an internal tandem duplication of the juxtamembrane portion of FLT3, a receptor which plays an important role in controlling proliferation and/or survival of hematopoietic progenitors, has been implicated as a cooperating genetic alteration in various AML subtypes.27,28 Since FLT3-ITD has been detected in about 40% of AML with mutated NPM1,1 it has been suggested that it may play an important role also in this leukemia subtype.
The findings of this paper suggest that the role of FLT3-ITD as a cooperative mutation in the pathogenesis of NPMc+ AML should be interpreted with caution. In fact, no difference can be detected between the slopes of the age specific incidence of FLT3-ITD-positive and -negative NPMc+ AML, supporting the view that NPMc+ AML is a homogeneous group irrespective of the FLT3 mutational status. This is consistent with the observation that the unique gene expression profile of AML with mutated NPM1, i.e. upregulation of HOX genes and downregulation of CD34,7,8 does not appear to be significantly influenced by the FLT3 gene status. This is also in keeping with the clinical observation that FLT3-ITD can appear or disappear in NPM1-mutated AML patients during the course of the disease.39 Moreover, in oncogenic cooperation tests, the NPM1 leukemic mutant and FLT3-ITD did not cooperate to transform mouse embryonic fibroblasts (MEFs).44 Hypothetically, FLT3-ITD may not be necessary for the development of AML but rather provide a selective advantage for leukemic cells that already harbor the NPM1 mutation. Unfortunately, there is as yet no experimental mouse model to prove or disprove this hypothesis. However, this interpretation would at least fit with the clinical observation that FLT3-ITD appears to abrogate the favorable prognostic impact of NPM1 mutations,29 suggesting that it may play a role at later stages of NPMc+ AML, leading to a more aggressive AML phenotype.
Thus, how can we reconcile the results of our one-mutation mathematical model with current evidence that favor the hypothesis that AML is the result of more than one oncogenic hit43 ? One possible explanation is that NPMc+ AML arises from the concerted action of an NPM1 mutation and another leukemogenic event occurring at the same time. Since the NPM1 mutant has intrinsic oncogenic properties44 and in knock-out mice NPM haploinsufficiency results in a MDS-like syndrome45 and in overt leukemia,46 an attractive hypothesis would be that these alterations act together to cause NPMc+ AML.3,47 Indeed, NPM1 mutations are associated with haploinsufficiency of wild-type NPM in leukemic cells, since mutations are always monoallelic3 and lead to dislocation of functionally active wild-type NPM from the nucleoli to the cytoplasm through formation of heterodimers with the NPM1 leukemic mutant.4
However, other scenarios cannot be excluded with certainty only on the basis of the mathematical model. NPM1 and yet undiscovered mutation(s) may act synergistically such that their actions cannot be discerned when investigating incidence data. Moreover, even though NPM1 mutations may be sufficient to cause leukemia, secondary mutations (e.g. FLT3-ITD) could increase the fitness of leukemic cells and/or result in the development of more aggressive AML stages. Finally, it is still possible that cancer incidence data cannot be used to identify the number of genetic changes necessary to cause cancer. Therefore, further experimental studies are warranted to clarify the oncogenic role of NPM1 mutations and other putative cooperating genetic lesions in NPMc+ AML.
This paper contains Supplementary Material. AL, FC and AC contributed equally to this work.
The online version of this article contains a supplemental appendix.
AL and BF had the original idea, coordinated the whole project and wrote the paper; FC and AC adapted the one-mutation mathematical model to the study of AML with mutated NPM1 and helped write the manuscript. FS performed the statistical analyses on incident cases before fitting the one-mutation model; RFS collected molecular and clinical data from patients of the AMLSG study and helped write the manuscript; SA was involved in designing the GIMEMA study and collecting clinical data from patients; CT performed molecular analyses of patients from DSIL and helped write the manuscript; SS performed mutational analysis in patients from the Munich Leukemia Laboratory (MLL) and helped write the manuscript; PJMV carried out molecular studies on AML patients from The Netherlands and reviewed the manuscript; KD collected molecular and clinical data from patients of AMLSG study and helped write the manuscript; MFM recruited patients in the GIMEMA study and reviewed the manuscript; MS designed and coordinated the clinical study (DSIL); JK collected molecular and clinical data from patients of the AMLSG study; AG collected clinical data from patients of the AMLSG study and coordinated the clinical study (AMLSG); MPM and NB performed immunohistochemical studies on the GIMEMA patients; BL recruited patients from The Netherlands and reviewed the manuscript; TH coordinated the study of patients from the Munich Leukemia Laboratory (MLL) and helped write the manuscript; GE designed and coordinated the clinical study (DSIL); FM designed and coordinated the clinical study (GIMEMA); HD designed and coordinated the clinical study (AMLSG); FM carried out computational simulation studies and helped write the manuscript.
The authors reported no potential conflicts of interest.
Received for publication April 9, 2008. Accepted for publication May 7, 2008.
|
|
|---|
This article has been cited by other articles:
![]() |
C. Haferlach, C. Mecucci, S. Schnittger, A. Kohlmann, M. Mancini, A. Cuneo, N. Testoni, G. Rege-Cambrin, A. Santucci, M. Vignetti, et al. AML with mutated NPM1 carrying a normal or aberrant karyotype show overlapping biologic, pathologic, immunophenotypic, and prognostic features Blood, October 1, 2009; 114(14): 3024 - 3032. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Facchetti, S. A. Pileri, C. Agostinelli, M. P. Martelli, M. Paulli, A. Venditti, M. F Martelli, and B. Falini Cytoplasmic nucleophosmin is not detected in blastic plasmacytoid dendritic cell neoplasm Haematologica, February 1, 2009; 94(2): 285 - 288. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Meloni, M. Mancini, V. Gianfelici, M. P. Martelli, R. Foa, and B. Falini Late relapse of acute myeloid leukemia with mutated NPM1 after eight years: evidence of NPM1 mutation stability Haematologica, February 1, 2009; 94(2): 298 - 300. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||