- Uri Rozovski1,2,3,
- Srdan Verstovsek1,
- Taghi Manshouri1,
- Vilma Dembitz1,4,
- Ksenija Bozinovic1,
- Kate Newberry1,
- Ying Zhang1,
- Joseph E. Bove IV1,
- Sherry Pierce1,
- Hagop Kantarjian1 and
- Zeev Estrov1⇑
- 1Department of Leukemia, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
- 2Davidoff Medical Center, Beilinson Hospital, Petah Tikva, Israel
- 3Tel Aviv University, Sackler School of Medicine, Israel
- 4University of Zagreb School of Medicine, Croatian Institute for Brain Research, Croatia
In most patients with primary myelofibrosis, one of three mutually exclusive somatic mutations is detected. In approximately 60% of patients, the Janus kinase 2 gene is mutated, in 20%, the calreticulin gene is mutated, and in 5%, the myeloproliferative leukemia virus gene is mutated. Although patients with mutated calreticulin or myeloproliferative leukemia genes have a favorable outcome, and those with none of these mutations have an unfavorable outcome, prognostication based on mutation status is challenging due to the heterogeneous survival of patients with mutated Janus kinase 2. To develop a prognostic model based on mutation status, we screened primary myelofibrosis patients seen at the MD Anderson Cancer Center, Houston, USA, between 2000 and 2013 for the presence of Janus kinase 2, calreticulin, and myeloproliferative leukemia mutations. Of 344 primary myelofibrosis patients, Janus kinase 2V617F was detected in 226 (66%), calreticulin mutation in 43 (12%), and myeloproliferative leukemia mutation in 16 (5%); 59 patients (17%) were triple-negatives. A 50% cut-off dichotomized Janus kinase 2-mutated patients into those with high Janus kinase 2V617F allele burden and favorable survival and those with low Janus kinase 2V617F allele burden and unfavorable survival. Patients with a favorable mutation status (high Janus kinase 2V617F allele burden/myeloproliferative leukemia/calreticulin mutation) and aged 65 years or under had a median survival of 126 months. Patients with one risk factor (low Janus kinase 2V617F allele burden/triple-negative or age >65 years) had an intermediate survival duration, and patients aged over 65 years with an adverse mutation status (low Janus kinase 2V617F allele burden or triple-negative) had a median survival of only 35 months. Our simple and easily applied age- and mutation status-based scoring system accurately predicted the survival of patients with primary myelofibrosis.
Primary myelofibrosis (PMF) is a myeloproliferative neoplasm characterized by bone marrow fibrosis and extramedullary hematopoiesis, resulting in variable degrees of splenomegaly, leukocytosis, anemia, thrombocytopenia, and impaired quality of life.1 The median survival of patients with PMF is five years from diagnosis,2,3 but the clinical course is variable. Some patients succumb to the disease within one year, whereas others survive for more than ten years.2–4
Several prognostic scoring systems have been developed for PMF that are based on clinical characteristics and blood counts.2,5 The international prognostic scoring system (IPSS) stratifies patients into 4 risk groups (low, intermediate-1, intermediate-2, and high) based on age (>65 years), the presence of constitutional symptoms, hemoglobin less than 10 g/dL, white blood cell (WBC) count of more than 25×109/L, and circulating blast cells of 1% or more at time of diagnosis.5 Based on the IPSS, a dynamic IPSS (DIPSS) was developed, which accounts for acquisition of risk factors over time.6 A refinement of the DIPPS that incorporates adverse karyotype, transfusion dependency, and thrombocytopenia has been suggested by the Mayo Clinic.7
In most patients with PMF, one of three mutually exclusive hematopoietic cell somatic mutations is commonly identified.8–10 In approximately 60% of PMF patients, an activating substitution mutation at position 617 of the pseudo kinase domain of Janus kinase 2 (JAK2V617F) is detected.11–13 In 20%–25% of patients, frameshift mutations caused either by deletions or insertions in the last exon of the calreticulin (CALR) gene are detected. CALR encodes a Ca++ binding protein that is primarily localized to the endoplasmic reticulum (ER). When CALR is mutated, the ER C-terminal ER retention signal (KDEL) is lost and, as a result, the protein is no longer localized to the ER.8,9 In 5% of patients, an activating mutation in the myeloproliferative leukemia virus (MPL) (thrombopoietin receptor) gene is found.14
Patients with a mutated JAK2 present a more aggressive disease than patients with mutated CALR. However, the overall survival (OS) of patients with high JAK2V617F allele burden is better than that of patients with a low JAK2V617F allele burden.15,16
In this study, we developed an easy-to-use scoring system that integrates age and mutation status, and accurately predicts the survival of patients with newly diagnosed PMF.
Included in our study were patients with PMF who were referred to MD Anderson Cancer Center, Houston, USA, between June 2000 and July 2013; the diagnosis was established in accordance with World Health Organization (WHO) criteria.17 Demographic and clinical information at the time of presentation was obtained from patients’ medical records by using a retrospective chart review protocol that was approved by the MD Anderson Institutional Review Board. The IPSS5 and DIPSS scores were assigned to each patient, as previously described.6
After obtaining patients’ informed consent, we analyzed residual blood and/or bone marrow cells which had been obtained from the patients for diagnostic purposes and stored, in accordance with a research protocol that was approved by the MD Anderson Institutional Review Board. Before freezing, all samples were fractionated with use of the Ficoll Hypaque 1077 (Sigma-Aldrich, St. Louis, MO, USA). Low-density cells were recovered from the Ficoll interface and collected by centrifugation. Genomic DNA was extracted by using Puregene DNA purification reagents (Gentra, Minneapolis, MN, USA).
To detect JAK2V617F mutation and measure the JAK2V617F allele burden, we extracted 50 ng of total genomic DNA and performed quantitative allele-specific suppressive polymerase chain reaction (PCR) with the use of the 7900HT FAST platform sequence detection system (Applied Biosystems, Foster City, CA, USA), as previously described.18
Detection of frame shift mutations in exon 9 of CALR was performed as previously described8 by using the following primer pairs: Forward: 5′ -FAM-GGCAAGGCCCTGAGGTGT; reverse: GGCCTCAGTCCAGCCCTG. This reaction captures the two frameshift type 1 (52 bp deletion) and type 2 (5-bp deletions) mutations. To detect mutations in exon 10 of MPL, we amplified genomic DNA by using the following primer set: MPL13474-F; GTGACCGCTCTGCATCTAGTG, MPL13726-R; GTGGGCGTGTTAGAG TGT. The resulting 250-bp PCR product was purified with use of a Qiagen PCR purification kit (Qiagen, Valencia, CA, USA) and was subjected to Sanger sequencing by using the above primers on a 3300 Genetic Analyzer (Applied Biosystems). The DNA sequencing fragments were analyzed with the use of Lasergene 11 (DNASTAR, Madison, WI, USA).
Patients’ characteristics were summarized by using frequencies (percentages) for categorical variables and median and range for continuous variables. To compare patients on the basis of categorical variables, we used the χ2 test. To compare medians, we used the Mann-Whitney test. To determine the optimal survival cut-off point that dichotomized patients according to their JAK2V617F allele burden, we used the X-Tile statistical software (http://www.tissuearray.org/rimmlab). The cut-off point used corresponds to the maximum χ2 value of the Mantel-Cox test for OS between groups above and below the cut-off point threshold.19 The probability of OS was estimated by the Kaplan-Meier method. The log-rank test was used to compare patients’ survival. Univariable and multivariable Cox proportional hazard regression models were fit to assess the association between mutation status and OS. The Wald test was used to assess the significance of covariates in Cox models. To compare competing models, we used the log-likelihood ratio. The replicability of the prognostic scoring system was tested by bootstrap resampling. One thousand samples, the same size as the original series, were built through random extraction with reposition. To predict the risk of transformation based on mutation status, we applied a logistic regression model and used the Exp (β) to estimate the odds ratio and the 95% confidence interval (CI) around it. Statistical analyses were performed with SPSS software (version 21, SPSS Inc., Chicago, IL, USA) and Graph Pad Prism (version 6.0, San Diego, CA, USA).
JAK2V617F, CALR, and MPL mutation frequency
A total of 344 PMF patients, aged 26 to 86 years (median: 65 years; 64% males) were included. Patients’ characteristics are shown in Table 1. Of the 344 patients, 226 (66%) had a JAK2V617F mutation, 43 (12%) had a CALR mutation (40 patients had 50–52-bp deletions and 3 patients had 5–10 bp insertion), and 16 (5%) had an MPL mutation. Fifty-nine patients (17%) had none of these mutations and were designated ‘triple-negative’.
JAK2V617 allele burden and survival
When used as a continuous variable, JAK2V617F allele burden (ranging from 0% to 98%) had only marginal power to predict OS [Hazard Ratio (HR): 0.997, 95% Confidence Interval (CI): 0.990–1.00]. However, a 50% cut off of the JAK2V617F allele burden dichotomized patients into two groups with different survival outcomes. Patients with a JAK2V617F allele burden of 50% or over had a median OS of 80 months (95%CI: 51–109 months), whereas patients with a JAK2V617F allele burden of less than 50% had a median OS of 50 months (95%CI: 40–60 months) (P=0.01) (Figure 1A). Remarkably, patients with a high (≥50%) JAK2V617F allele burden had a larger spleen, a higher hemoglobin level, and a higher WBC count than did patients with a low (<50%) JAK2V617F allele burden (Table 2).
Mutation status and survival outcome
The longest OS was observed in patients with mutated MPL (median survival 221 months, 95%CI: 40–401 months), followed by patients with mutated CALR (median survival 131 months, 95%CI: 100–160 months), high JAK2V617F burden (median survival 80 months, 95%CI: 51–109 months), triple-negatives (median survival 56 months, 95%CI: 35–77 months), and low JAK2V617F burden (median survival 50 months, 95%CI: 38–62 months) (Figure 1B). The incorporation of high- and low-JAK2V617F mutation status divides PMF patients into two groups: patients with either high JAK2V617F allele burden, mutated CALR, or mutated MPL had a median survival of 104 months (95%CI: 86–122 months), whereas patients with low JAK2V617F allele burden or triple-negative mutation status had a median survival of 48 months (95%CI: 39–57 months) (Figure 1C).
Development of an age- and mutation status-based prognostic model for survival
When mutation status and DIPSS variables were included as covariates in a multivariable analysis, only unfavorable mutation status, older age, and a high percentage of peripheral blood blasts predicted a shorter survival (Table 3). Age dichotomized patients into two groups with different survival outcomes. Patients aged 65 years or under had a median OS of 97 months (95%CI: 67–127 months; n=175), whereas patients older than 65 had a median OS of 47 months (95% CI: 39–55 months; n=169) (P<0.0001) (Figure 1D). However, a model that included 2 parameter estimates (age and mutation status) was superior in fitting the data, as had the largest log-likelihood ratio, and divided the cohort into 4 groups of almost equal size. Bootstrapping resampling procedure confirmed the stability of the model. Patients with a favorable mutation status (high JAK2V617F allele burden, CALR, or MPL mutations) and aged 65 years or under had a median OS of 126 months (95%CI: 91–161 months; n=82). Patients with one risk factor, either age over 65 years (n=88) or adverse mutation status (n=87) had an intermediate OS of 72 months, whereas patients with two risk factors, e.g. age over 65 years and an adverse mutation status (low JAK2V617F allele burden or triple-negative; n=87) had a median OS of 35 months (95%CI: 31–113 months) (Table 4 and Figure 1E). In comparison, the DIPSS uses 5 risk factors to classify patients into one of 4 groups. In our cohort, most patients (n=228, 78%) were classified as either intermediate-1 or -2 (Table 5) and had similar survival outcome.
Mutation status and the risk of transformation
Thirty-two patients (9%) transformed to acute myeloid leukemia. Median time to transformation was 33 months (range 1–271 months). None of the mutations predicted transformation to acute myeloid leukemia.
The clinical outcome of patients with PMF is partially dictated by mutually exclusive driver mutations in the genes JAK2, CALR, or MPL.20 Here we show that patients’ mutation status can be integrated into a prognostic model.
Although the survival of patients with a JAK2V617F mutation is heterogeneous, dividing these patients into subgroups of high and low JAK2V617F allele burden enabled the development of a prognostic model that integrates genetic information. Patients with low JAK2V617F allele burden or a triple-negative mutation status had a shorter OS than patients in the other groups. Because we found that patients’ age is of prognostic significance, similar to other investigators’ findings, we integrated the patients’ mutation status and age, and analyzed 4 equally-sized cohorts. Patients with no risk factors (aged 65 years or under with a favorable mutation status) had the longest survival (median OS 126 months); patients with a single risk factor (over 65 years of age or an adverse mutation status) had an intermediate survival duration (median OS 72 months); and patients with two risk factors (over 65 years of age and an adverse mutation such as JAK2V617F allele burden or triple-negative mutation status) had the worst prognosis (median OS 35 months). Although the percentage of circulating blasts emerged as a prognostic indicator, it did not contribute to the overall variance and was not included in the final model. Given that our hospital is a tertiary care cancer center, our PMF patient cohort has a high proportion of high-risk patients; this PMF patient population is particularly suitable for this analysis. Notably, only 7% of our patients were low-risk patients according to the DIPSS, compared with 44% in the IPSS. Hence, while our prognostic scale divides our patient cohort into 4 groups of equal size, it is possible that lower-risk patients were under-represented. Nevertheless, because prognostic scales for patients with PMF are routinely used to identify high-risk patients who are suitable for allogeneic transplantation, this scale might prove to be very useful.
In our cohort, 13% of PMF patients had triple-negative mutation status. It is possible that these patients carry mutations in yet unidentified genes or that triple-negative status may represent a late event in clonal evolution that gives proliferation and/or survival advantage to a dominant neoplastic clone that is no longer dependent on the initial ‘driver’ mutagenic event.
The survival of patients with a high (≥50%) JAK2V617F allele burden was significantly better than that of patients with a low JAK2V617F allele burden. It has been reported that PMF patients with a homozygous JAK2 mutation have distinct clinical features such as splenomegaly.21 Here we show that the clinical features on presentation of patients with a high JAK2V617F allele burden were reminiscent of patients with polycythemia vera (PV); they had discernible splenomegaly, leukocytosis, and higher hemoglobin levels compared with the group with low JAK2V617F allele burden. Therefore, it is possible that this group consists, at least in part, of patients with post-PV myelofibrosis that evolved from an undiagnosed PV. Interestingly, in a study of 68 patients with post-PV myelofibrosis, all patients carried a high JAK2V617F allele burden, and 78% had an allele burden of more than 50%.6,22 CALR mutations have been divided into two types. In PMF the type 1/type 1-like mutations are the most common ones.23–25 In our patient cohort only 3 patients had type 2 mutations. Therefore, our study was not powered to determine the prognostic value of CALR mutation subtypes.
Since the initial publication of the IPSS prognostic score,5 several refinements have been proposed, most of which attempt to incorporate recurrent gene mutations that have been identified in patients with PMF.26 Some mutations, such as those in DNTM327 or TET2,28 have not been shown to correlate with survival outcome. Conversely, mutations in ASXL1, SRSF2, and EZH2 predicted short survival in a large cohort of patients, and only the ASXL1 mutation remained statistically significant when added to the IPSS prognostic score.29 A report by Tefferi et al.10 points to the CALR−/ASXL1+ profile as the most detrimental mutation profile in PMF.
The applicability of our prognostic scale depends on screening for mutations in CALR and MPL and quantification of the JAK2V617F allele burden. Recently, the WHO added CALR and MPL mutations to the PMF diagnostic criteria30 and, as a result, most diagnostic laboratories perform these tests. Moreover, most diagnostic laboratories assess the presence of JAK2 mutations by using quantitative PCR. Although the JAK2V617F allele burden is readily available, it is not routinely reported, although various assays yield similar quantification results.31
Here we present a prognostic model that is based on a relatively large cohort. The internal validation of this model was confirmed by bootstrap resampling. By using only 2 variables, we developed a simple, easily applied model with excellent discrimination power for survival outcome of patients with newly diagnosed PMF.
Although this prognostic model needs to be validated with a large, independent patient population, it has a larger log-likelihood ratio than that of the IPSS or DIPSS, suggesting that it has superior, clinically applicable value.
We thank Tamara Locke for editing our manuscript.
Check the online version for the most updated information on this article, online supplements, and information on authorship & disclosures: www.haematologica.org/content/102/1/79
This study was supported by the NIH/NCI under award number P30 CA016672.
- Received May 23, 2016.
- Accepted September 14, 2016.
- Copyright© Ferrata Storti Foundation