Figures
Abstract
Background
myelodysplastic syndromes (MDS) are a heterogeneous group of hematopoietic clonal disorders. So, prognostic variables are important to separate patients with a similar biology and clinical outcome. We compared the importance of risk stratification in primary MDS of IPSS and WPSS with the just described revision of IPSS (IPSS-R), and examined if variables obtained by bone marrow immunophenotyping could add prognostic information to any of the scores.
Methods
In this prospective study of 101 cases of primary MDS we compared the relation of patients’ overall survival with WHO types, IPSS, IPSS-R, WPSS and phenotypic abnormalities of hematopoietic precursors. We examined aberrancies in myelomonocytic precursors and CD34+ cells. Patients were censored when receiving chemotherapy or BM transplantation. Survival analysis was made by Cox regressions and stability of the models was examined by bootstrap resampling.
Results
median age: 64 years (15–93). WHO types: 2 cases of 5q- syndrome, 7 of RA, 64 of RCDM and 28 of RAEB. In the univariate Cox analysis, increasing risk category of all scores, degree of anemia, higher percentage of BM blasts, higher number of CD34+ cells and their myeloid fractions besides increasing number of phenotypic abnormalities detected were significantly associated with a shorter survival. In the multivariate analysis comparing the three scores, IPSS-R was the only independent risk factor. Comparing WPSS with phenotypic variables (CD34+/CD13+ cells, CD34+/CD13− cells and “total alterations”) the score and “CD34+/CD13+ cells” remained in the model. When IPSS was tested together with these phenotypic variables, only “CD34+/CD13+ cells”, and “total alterations” remained in the model. Testing IPSS-R with the phenotypic variables studied, only the score and “CD34+/CD13+ cells” entered the model.
Citation: Reis-Alves SC, Traina F, Harada G, Campos PM, Saad STO, Metze K, et al. (2013) Immunophenotyping in Myelodysplastic Syndromes Can Add Prognostic Information to Well-Established and New Clinical Scores. PLoS ONE 8(12): e81048. https://doi.org/10.1371/journal.pone.0081048
Editor: Robert K. Hills, Cardiff University, United Kingdom
Received: March 27, 2013; Accepted: October 17, 2013; Published: December 6, 2013
Copyright: © 2013 Reis-Alves et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: Grants: FAPESP (Research Funding Agency of the State of São Paulo) proc. 04/08882-3, CNPq/INCTS (National Reseach Council/Funding for Research Excellence Centers) proc.2008-57895/1, FAEPEX (Research Fund of the University of Campinas) proc 1208/11 and MDS Foundation (Tito Bastianello Young Investigator grant 2009). Irene Lorand-Metze, Konradin Metze and Sara T.O. Saad have research grants from CNPq (302277/2009-9, proc 307270/2010-6, and 301483/2008-6 respectively). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Author Konradin Metze is an Academic Editor for PLOS ONE. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.
Introduction
Myelodysplastic Syndromes (MDS) constitute a wide spectrum of hematopoietic clonal disorders with a variable clinical course [1]–[5]. In this context, the study of features that are able to predict patients’ survival and progression to AML is very important [4]–[10], in order to assign patients correctly in clinical trials which test new treatment approaches.
The severity of PB cytopenias, the percentage of BM blasts and the kind of cytogenetic abnormalities found have long been recognized as independent prognostic factors in MDS, and have been included in the currently used prognostic scores, such as IPSS and WPSS [4]–[10].
The IPSS, [6] that is currently the most frequently used risk stratification score for MDS, is based on the number of cytopenias found in peripheral blood counts (PB), percentage of bone marrow (BM) blasts and kind of cytogenetic abnormalities. This score has recently been revised (IPSS-R) [9] as the importance of karyotype abnormalities was underscored in the former classification system. Cytogenetic findings were re-analyzed in a large international multicentric study and their importance for risk stratification was revised [7]. Besides, IPSS-R risk categories are based not only on these revised cytogenetic groups, but also on a more detailed categorization of the peripheral blood values and BM blast counts [9]. Therefore, a considerable proportion of cases with IPSS intermediate risk switched to a higher risk category.
Karyotype abnormalities are found in 30%–80% of cases with primary MDS. They are very heterogeneous, and some abnormalities are very rare, precluding the assessment of their real prognostic value [8]. On the other hand, several frequent point mutations such as TP53, EHZ2, ETV6, RUNX1 and ASXL1 have been considered independent prognostic factors when compared to age, sex and IPSS classification [11], [12]. TET2 mutations were frequent in cases with a normal karyotype and those of TP53 were associated with abnormalities of chromosome 17 or a complex karyotype, but mutations of EZH2 (localized in chromosome 7) were not associated with 7q deletion. So, the presence of multiple gene mutations in MDS may help to explain the clinical heterogeneity of these clonal disorders and could help to improve the prediction of the patients’ prognosis.
In 2005, a WHO classification-based Prognostic Scoring System (WPSS), considering WHO categories, transfusion requirement, and karyotype abnormalities (risk categories as in IPSS) was described. As transfusion dependency is difficult to standardize, this parameter was substituted by the hemoglobin level of the patient: <9 g/dL for males and <8 g/dL for females [10]. This score is dynamic, and can be applied at every point in the course of the disease.
In the last decade, multiparametric flow cytometric analysis of BM hematopoietic precursors has been extensively studied in MDS and is nowadays recognized as a useful diagnostic tool, especially in cases with a normal karyotype [13]–[22]. Immunophenotyping in MDS is based on the knowledge that antigen expression during maturation of normal hematopoiesis is tightly controlled. In MDS, deviations of the normal pattern, with over- or underexpression of antigens, as well as maturation asynchrony are indicative of a clonal abnormality [13]–[22]. Increased number and aberrant antigen expression of CD34+ cells, as well as the total number of phenotypic abnormalities in BM precursors, [13], [14], [17]–[22] have shown to be independent risk factors for survival.
Immunophenotyping in MDS is feasible in all patients, nowadays reasonably well standardized, and has been recommended as an ancillary diagnostic tool for the differential diagnosis of MDS with a normal karyotype and non-clonal reactive PB cytopenias [13]–[16], [18], [23].
So, many independent prognostic variables and risk stratification scores have been described in MDS. However, a head-to-head comparison of them has seldom been performed. Especially it is not clear at the moment, to which degree the IPSS-R increases the prognostic power when compared with the old IPSS classification. Therefore, the aim of our prospective study was to compare in patients with primary MDS the two well-established scores (WPSS and IPSS) with the newly created IPSS-R classification and to investigate whether variables obtained by BM flow cytometric analysis could add prognostic information to any of these scores.
Methods
Patients
We analyzed prospectively 87 cases with adult primary MDS diagnosed at our Institution between 2006 and 2012 and 14 patients who had already been diagnosed at other Institutions and were referred to us, but had not received any treatment. The project was approved by the Ethics Committee of the Faculty of Medical Sciences, University of Campinas (Proc 0652.0.146.000-08).
Diagnosis was based on clinical data, PB counts, BM cytology and histology, as well as cytogenetics, as recommended by the WHO criteria. Stability of PB counts for at least 3 months was required for all cases. Non-clonal disorders presenting PB cytopenias and a cellular marrow with hematopoietic cell atypias (HIV, hepatitis and other viral infections, autoimmune disorders, etc) were excluded. At the time of the patients’ evaluation none of them was being treated.
For cytogenetic examination, nucleated cells from aspirated bone marrow were separated by Ficoll-Hypaque centrifugation and then cultured for 48 hours in RPMI1640 supplemented with calf fetal serum without stimulation. Finally, G-banding was performed. Results were reported using the ISCN (International Standing Committee on Human Cytogenetic Nomenclature) criteria [24].
The cases were classified by WHO criteria and the risk was accessed according to IPSS, IPSS-R and WPSS using the hemoglobin level instead of “transfusion dependency” [10]. Only patients with a high risk disease according to WHO criteria (refractory anemia with excess of blasts - RAEB) were considered for specific treatment with cytotoxic therapy. The other patients received only supportive care.
All PB and BM samples for the study were obtained from each patient after written informed consent.
Flow Cytometric Analysis
Flow cytometric analysis was performed as previously described [15], [17], [20]. In brief, aspirated BM was collected in EDTA and diluted to obtain a concentration of 5–7×106 cells in 100 µl per test. Immunofluorescence staining was made using a standardized direct lyse-and-wash technique within 24 hours after sample collection following the European LeukemiaNet (ELN) recommendations [15], [23]. Antigenic expression of the myelomonocytic series and CD34+ cell subsets were analyzed using four color combinations of MoAbs: CD64/CD14/CD45/HLA-DR; CD16/CD11b/CD45/CD13; CD13/CD117/CD45/CD34; CD19/CD10/CD45/CD34 and CD7/CD56/CD45/CD34 (Table 1).
This platform is able to reveal up to 16 alterations in the granulocytic and monocytic cell lines as well as in CD34+ progenitors [14], [15], [17], [20], [23]:
Granulocytic precursors (4 abnormalities): decreased SSC, abnormal maturation pattern in the CD16/CD13 and CD16/CD11b combinations (deficiency or increase in antigen expression or asynchronous maturation) and expression of CD34 in >10% of mature granulocytes (Fig. 1).
The granulocytic population was selected on the SSC/CD45 plot (A). The four differentiation stages are distinguishable based on the expression of CD11b, CD16 and CD13 (B–D: orange: promyelocytes; green: myelocytes; purple: metamyelocytes and blue: mature neutrophils). Maturing granulocytes show an abnormal decreased granularity demonstrated by a low SSC (A). Increase of a uniform population with the phenotype CD11blow/CD13low/CD16− (promyelocytes) with a gap of maturation between them and myelocytes. Abnormally high expression of CD13 in myelocytes and metamyelocytes (C and D).
Monocytic cell line (4 abnormalities): increased number, increase in CD16+ monocytes, asynchronous/abnormal maturation pattern and expression of CD34 in >10% of the cells (Fig. 2).
Left: normal BM. Right: a case of RCMD (monocytes = 2.4%). Abnormal continuous acquisition of CD16 in the mature monocytes (B).
CD34+ cells (6 abnormalities): increase of CD34+ cells, decrease of B-cell precursors (CD34+/CD19+/CD10+), increase of myeloblasts (SSCint/CD13+/CD117+/CD34+) (Fig. 3), increase of immature non-lymphoid precursors (SSCint/CD34+/CD117− or SSCint/CD34+/CD13−), and abnormal expression of CD7 or CD56 in CD34+ cells.
(combination CD19/CD10/CD45/CD34 in A and B and CD13/CD117/CD45/CD34 in C and D). In a first step, (A) all CD34+ cells were selected in the SSC/CD34 plot. Two populations can be observed: SSClow (in brown) that are CD19+/CD10+ (B-cell precursors) on the CD19/CD10 plot (B), and SSCint/CD34+/CD19− (light blue). These cells are characterized as myeloid subsets of CD34+ cells (C and D) when analysed in the CD13/CD117/CD45/CD34 combination. Among SSCint/CD34+ cells, the CD34+/CD117+ subset (dark blue) and CD34+/CD117− subset (red) were defined by their expression of CD117 and CD13. Myeloblasts SSCint/CD34+/CD117+/CD13+ (dark blue) and immature cells SSCint/CD34+/CD117−/CD13− (red).
Immediately after staining, samples were acquired in a FACSCalibur flow cytometer (Becton Dickinson – BD Biosciences, San José, CA, USA) using the CellQuest software (BD Biosciences). Information of least 100,000 events was acquired for each sample. Data analysis was made using the Infinicyt software (Cytognos, Salamanca, Spain). We computed the total number of granulocytic and monocytic abnormalities and those of CD34+ cells. The sum of all these abnormalities was called “total alterations” Concerning CD34+ cells, we calculated the total percentage of positive cells among all nucleated cells, and the subsets: B-precursor cells (CD34+/CD19+/CD10+), SSCint/CD34+/CD13+ and SSCint/CD34+/CD13− cells [17], [20], [23].
Statistics
Differences among groups and relations between phenotypic features and other risk factors were analyzed by non-parametric tests (Mann-Whitney, Kruskall-Wallis and Spearman’s test). Overall survival of the patients according to WPSS, IPSS and IPSS-R was calculated by the Kaplan-Meier method followed by the log-rank test. Cluster analysis for cytometric data was performed according to the Ward algorithm [25]. We tested the WHO classification, WPSS, IPSS and IPSS-R as categorical variables and all other variables entered the models as continuous variables. Patients were censored if they underwent chemotherapy or bone marrow transplantation.
The prognostic relevance of PB counts, WHO classification, the three prognostic scores analyzed, revised cytogenetic categories, % BM blasts, total CD34+ cells, SSCint/CD34+/CD13+ and SSCint/CD34+/CD13− cells, as well as the “total alterations” was analyzed by univariate and multivariate Cox regressions. The limits for input were p = 0.05 and p = 0.1 for output. Backward conditional step-wise selection was applied to all calculations with the exception of the multivariate comparison of the three risk stratification scores, where the forward stepwise conditional mode was used.
Since our aim was to examine if variables derived from flow cytometry analysis which had impact on patients’ survival could add prognostic information to the three scores analyzed we created three multivariate models, testing always the variables with p<0.05 in the univariate Cox models together with one of the three risk classifications systems. Furthermore, we run the three risk stratification systems together in a multivariate proportional hazard model.
The internal stability of the models was tested by bootstrap resampling [26]–[39]. In brief, new data sets with the same size of the original one were created by random sampling with replacement. In each new bootstrap data set, a patient may be represented once, several times or not at all. Cox regressions with the same conditions as in the original data set were performed for each of these new bootstrap data sets.
The 95% confidence intervals of the R-squared values were calculated based on the values obtained by bootstrapping. In each bootstrap set we analyzed the consistency of the three classifications, i.e. whether with increasing risk group the hazard ratio increased. SPSS 8.0 (to calculate the R values), SPSS 15.0 and Winstat softwares were used for calculations.
Results
A total of 101 patients (64 males and 37 females) were included in this prospective study. The median age at diagnosis was 64 years (15–93). Among the WHO types, the majority of the cases had refractory cytopenias with multilineage dysplasias (RCMD) and had a good risk cytogenetics (Table 2). Concerning the IPSS, most patients were classified as low or intermediate-1 risk.
Comparing the cytogenetic risk classification used in IPSS and the revised one, only 17 cases changed category. The change in risk categories after reclassification of the patients according to IPSS-R is shown in Table 3. The patients with “IPSS low risk” were divided into “IPSS-R very low” and “low risk”, but those with “IPSS intermediate risk” were distributed into many categories in IPSS-R, while the majority of those with “IPSS high risk” were classified as “IPSS-R very high risk”
Concerning immunophenotypic abnormalities, all cases presented at least 2 alterations (Table 4). Cross-lineage aberrant expressions (CD7 and CD56) in myeloid CD34+ cells were found in 54% of the cases. Percentage of BM blasts was correlated with total CD34+ cells (r = 0.66; p<0.0001); “CD34+/CD13+”, and “CD34+/CD13− cells” (r = 0.51; p<0.0001 and r = 0.42; p<0.0001 respectively).
The total number of phenotypic abnormalities increased with the increasing risk in the classifications: WHO type (r = 0.27; p = 0.005), WPSS (r = 0.40; p<0.0001), IPSS (r = 0.36; p<0.0001) and IPSS-R (r = 0.49; p<0.0001). The total number of phenotypic abnormalities was negatively correlated with the hemoglobin levels (r = - 0.40; p<0.0001), but there were no significant correlations with PB neutrophil and platelet counts.
Survival Analysis
Median duration of the patients’ follow-up was 28 months (1–134). At the end of the observation period 29 patients had died. IPSS, IPSS-R and WPSS were able to stratify patients with a different survival (p<0.0001 for all three scores) (Fig. 4).
All three scores were able to separate patients with a significantly different survival (p<0.0001).
In the univariate Cox analysis (Table 5) all the classification systems and prognostic scores analyzed had a significant association with the overall survival of the patients.
Age was not related to survival but diminished hemoglobin values, revised cytogenetic risk categories and the percentage of BM blasts were risk factors. The same was true for the total number of immunophenotypic abnormalities found, number of abnormalities in monocytes and CD34+ cells as well as total number of CD34+ cells and their myeloid subsets (SSCint/CD34+/CD13+ and SSCint/CD34+/CD13−) (Table 6).
In the multivariate analysis, comparing the three prognostic scores studied, only IPSS-R remained as an independent risk factor of overall survival (p<0.0001; R2 = 0.07952). In the bootstrap resampling, IPSS-R was present in 88% of the new models, IPSS in appeared in 20% and WPSS in 18% of the new data sets. The variation of the R-squared values in the bootstrap resampling sets was: median 0.0462 (CI 95%: 0.0148–0.0960) for IPSS-R, median 0.0433 (CI 95%: 0.0106–0.1179) for IPSS and median 0.06993 (CI 95%: 0.0139–0.1156) for WPSS.
In the univariate analysis, the risk stratifications were consistent (B values continuously increasing with increasing risk) in 84% of the bootstrap sets analyzed with the IPSS or WPSS score and in 77% of the sets analyzed with the IPSS-R score.
When joining together WPSS and phenotypic variables (“CD34+/CD13+ cells”, “CD34+/CD13− cells” and “total alterations”), only WPSS and “CD34+/CD13+ cells” remained in the final multivariate Cox model (p = 0.02 and p = 0.004 respectively) (Table 7). In the bootstrap resampling stability test, WPSS remained in 88%, “CD34+/CD13+ cells” in 79%, “total alterations” in 44% while “CD34+/CD13− cells” remained only in 5% of the models.
In the multivariate analysis of IPSS and the three immunophenotypic variables considered (Table 7), only “CD34+/CD13+ cells”, and “total alterations” remained in the model (p = 0.0005 and p = 0.01 respectively), whereas the IPSS score was eliminated. In the stability test, IPSS appeared in 69% “CD34+/CD13+ cells” in 79%, “total alterations” in 55% and “CD34+/CD13− cells” in only 3% of the models.
Finally, in a multivariate analysis containing IPSS-R and flow cytometric variables, only IPSS-R and “CD34+/CD13+ cells” entered the model (p = 0.0012 and p = 0.01 respectively). In the bootstrap resampling, IPSS-R was found in 96%, “CD34+/CD13+ cells” in 76%, “total alterations” in 17% and “CD34+/CD13− cells” in only 5% of the models.
After normalization by logarithmic transformation, the variable “CD34+/CD13+ cells” was submitted to a cluster analysis with the Ward algorithm, which suggested 3 groups: group 1 (16 cases) median CD34+/CD13+ cells = 0.07 (0.02–0.11); group 2 (69 cases) median CD34+/CD13+ cells = 0.46 (0.14–2.36) and group 3 (15 cases) median CD34+/CD13+ cells = 5.07 (2.9–27.9). These groups had a significantly different survival in the log-rank test (p<0.0001) (Fig. 5).
Group 1 median: 0.07 (0.02–0.11), group 2 median: 0.46 (0.14–2.36) and group 3: median 5.07 (2.90–27.9) (p<0.0001).
Discussion
Myelodysplastic syndromes comprise a heterogeneous group of clonal hematopoietic disorders concerning patients’ clinical behavior and laboratory features [4]–[10]. In recent years, important advances in the understanding of its pathophysiology have been achieved and some disease-altering treatment strategies have been established [11], [12], [40]–[43]. In this context, new prognostic factors may enable us to improve the choice of the best therapeutic options.
In the present work, which is prospective and mono-institutional, from an Institution not participating in the International Working Group for Prognosis in MDS (IWG-PM) project [9] we compared the prognostic value of two well-established risk scores (IPSS and WPSS) and the newly described revised IPSS (IPSS-R). Furthermore, we examined if variables obtained by immunophenotyping could improve the risk stratification obtained by each of the clinical score systems.
IPSS is based on cytogenetic risk stratification, number of PB cytopenias and a rough categorization of BM blasts. These features have long been recognized as independent risk factors in MDS. WPSS incorporated the WHO classification of MDS that is also, per se, able to separate groups of patients with a different outcome [5], [10]. It also includes the concept that anemia and transfusion dependency are an important hallmark of disease severity. Anemia is the most common peripheral cytopenia in MDS and is associated with the most important morbidity such as transfusion dependency and iron overload [4], [10]. This was also observed in our patients, where the hemoglobin level was inversely related to survival, but neutrophil and platelet counts were not. This is also in keeping with the findings of the IWG-PM [9], and with another previous study of our group [5].
The recent IPSS-R classification system is not only based on a more refined cytogenetic risk classification, but also on a more detailed splitting of BM blast count and PB values. In the present work, we could validate this new score system in a patient cohort which had not participated in the IWG-PM study.
Concerning the internal validation of our model, different solutions could be chosen. An excellent method is data splitting, where a sub-sample will serve for model selection and the remaining cases will be used for the internal validation. But in our study, the relatively low number of cases in some subgroups precluded the use of this methodology. Alternative methods would be the jack-knife procedure [29]–[35], or the bootstrap resampling, which can be considered an elaboration of the jack-knife analysis. Bootstrap resampling is interesting, since it allows to investigate the consistency of the inclusion of each variable in the regression model [29]–[39] and to estimate “internal validity”, which means that the model closely mirrors the survival patterns of the original data [38]. Before a model is applied in routine clinical practice, it should however, be tested for external validity, which implies testing its generalizability with other groups of patients, preferentially in other Institutions.
In linear regressions, R2 helps to estimate the explanatory power of a model [33], [36], [37]. In a similar way, R-squared values can be calculated for survival regressions, where they asses the proportion of variability in the patients’ outcome explained by each variable. But here, the R-squared values decrease with increasing percentage of censored patients and earlier censoring. Moreover, the generalized R-squared does not represent exactly the proportion of variance of the dependent variable explained by the independent variable, but provides rather an estimate of association between these variables [38]. Therefore, the R-squared value may be small even in good fitting models, and is able to summarize considerable information.
When comparing the three classifications in the bootstrap sets in the univariate analyses, we found that the median R-squared values for IPSS and IPSS-R were rather similar, but the variation among the 100 new sets was larger for IPSS. In the multivariate analysis joining together the three scores, IPSS-R was present in most new models but the two other scores were present in only a small number of the bootstrap sets. Besides, when joining together each score with the same immunophenotypic variables in multivariate analyses, the IPSS score entered only in 69% of the bootstrap sets, whereas IPSS-R in 96%. For all these reasons the IPSS-R score may be regarded as a more consolidated risk classification system. This may be due to the fact that in IPSS-R the degree of the peripheral cytopenias and the percentage of bone marrow blasts are stratified more in detail.
A head-to-head comparison among the three scores, as done in our study has not been reported so far, to the best of our knowledge. So, a confirmation of our results in larger cohorts of patients should be performed.
Immunophenotyping is now a well-established ancillary technique for the diagnosis of MDS [13]–[18] and has been standardized by the European LeukemiaNet [15], [23]. This technology is able to demonstrate, in a simple and rapid way, the degree of functional damage of the abnormal clone, revealing maturational abnormalities in the myelomonocytic [14], [15]–[18], [22], [23] and erythroid lineages [14], [15] through abnormal antigen expression, asynchronous antigen co-expression as well as cross-lineage aberrant expressions. Moreover, increase in number of CD34+ progenitors, as well as their phenotypic abnormalities corresponding to leukemia-associated phenotypes (LAIPs), are also frequently found [15], [22], [23]. Most consensus reports recommend the analysis of the myelomonocytic lineages. However, the analysis of the erythroid precursors is not mandatory, because only few antibodies are available for this lineage. Besides, sample preparation may lyse a part of the erythroblasts, which can affect the measures obtained [13], [15], [18], [23].
Accordingly, several publications have emphasized the prognostic relevance of some of the features obtained, such as the total number of CD34+ cells and their subsets corresponding to early myeloid progenitors, aberrant expressions of CD34+ cells as well as the total number of phenotypic abnormalities detected [13]–[23]. An increase of the latter has been observed during disease progression [18], [20], [44], [45], [47] and has been related to a worse outcome after BM transplantation [47]. Some study groups reported a relation between the total number of phenotypic abnormalities and IPSS [14], [18], [43], [46] or WPSS [17], [20]–[22]. Aberrant cross-lineage antigen expressions in CD34+ cells, may be found in up to 38% of the cases [22], [45], can be detected in some patients with a low BM blast count and have been associated with disease progression, transfusion dependency and a shorter survival. Several flow scores have also been described, but there is no agreement about their use in clinical practice [18], [22], [43].
For these reasons, our main question was whether flow cytometric features could refine the prognostic value not only of IPSS-R, but also IPSS and WPSS. This comparison has not yet been done to the best of our knowledge. In our study, the number of phenotypic abnormalities found in monocytes and in CD34+ cells but not those of maturing granulocytes presented a significant relation to the patients’ survival in the univariate analysis. Consequently, also the total number of abnormalities detected was related to the overall survival of the patients, but had a lower relevance. The total number of CD34+ cells and their myeloid subsets were the most important features capable to predict the patients’ outcome. So we decided to examine in multivariate analyses if one of these variables could retain an independent prognostic value together with IPSS, WPSS or IPSS-R. We chose to test “CD34+/CD13+ cells”, “CD34+/CD13− cells” and “total alterations” as these features are easy to obtain and standardize and can be measured in all cases, including those lacking cytogenetic information or having a normal karyotype.
In the multivariate analysis comparing all scores examined with flow cytometric features, the number of SSCint/CD34+/CD13+ cells was the strongest independent feature that could add prognostic information to the clinical scores examined. These cells represent the major fraction of CD34+ cells and correspond to an objective measure of BM myeloid blasts [17], [20], [42], [48], [49]. They may be quantified in all patients, differently from LAIPs that are present in only about half of the cases, as was also observed in our patients [18], [23], [45].
Although all clinical scores consider the percentage of BM blasts in a categorized way or indirectly in the WHO types, our results demonstrate that enumeration of myeloid progenitors by FCM is able to increase risk stratification of the patients.
Different reasons may explain this phenomenon. First, we have to remind that every form of stratification, as done for the BM blast numbers in the score systems, will diminish considerably the test power [50]–[52]. Therefore, a part of the information “lost” by categorizing the blast numbers will be regained by adding the continuous flow cytometric variables in the Cox regression.
Moreover, blasts are counted by the human observer under the microscope [41]. This kind of measure is subject to considerable inter-observer variability and “expert opinion”. Besides, a relatively low number of blasts are counted in the specimens. So, the “manual” blast count cannot be regarded as a very precise measure (a maximum of total 500 cells are counted in routine work), especially when compared with flow cytometric data which are based on a large number of cells (usually at least 100000 events) [53]. Moreover, especially in MDS cases, the reproducibility of cell typing is low, due to the high degree of cell atypias which hampers considerably cell classification [41], [53]. Both kinds of samples can have varying degrees of hemodilution when compared to BM histological specimens [23], [54] resulting in different counts when assessed by both techniques. For all the reasons above one can consider both features as related but independent measures [23].
In order to illustrate the importance of the CD34+/CD13+ cell counts we performed a cluster analysis and stratified the data into three non-overlapping groups of patients with a different survival. These groups separated mainly patients with a low number of SSCint/CD34+/CD13+ cells and a low blast count in cytology, emphasizing again the importance of considering subtle abnormalities in the hematopoietic progenitors for the evaluation of the aggressiveness of the abnormal clone.
It was interesting to note that comparing IPSS with flow cytometric features, these last ones were more powerful in predicting patients’ outcome. This may be due to the fact that IPSS considers only the number of PB cytopenias and not the degree of BM failure. Besides, the categorization of BM blasts assigns more points to counts above 10%. However, in our cluster analysis we could separate three groups with a different survival mainly among cases with low numbers of SSCint/CD34+/CD13+ cells and therefore immunophenotypic variables were able to add important prognostic information to IPSS. This underlines the importance to refine stratification of patients with a low BM blast count in cytology, which was made in the revised IPSS-R, where more points are assigned also low BM blast counts.
Comparing WPSS with immunophenotypic variables, the clinical score remained in the model together with “SSCint/CD34+/CD13+ cells”, which were stable in the bootstrap resampling, although “total alterations” entered a large number of new data sets. However, in the comparison using IPSS-R, the best clinical score in our study, only “SSCint/CD34+/CD13+ cells” remained as a stable independent risk factor. However, all these results should be validated further in a larger cohort of patients.
We could confirm previous reports that maturation abnormalities in the granulocytic series and decrease in B-cell precursors are rather an adjuvant criterion for the differential diagnosis of MDS with a normal karyotype and non-clonal cytopenias, but have a limited impact on patients’ survival [13], [18], [20], [55], [56]. Interestingly, the number of abnormalities and especially the percentage of monocytes were significant risk factors, at least in the univariate analysis. It has been shown that, with the progression to AML, the number of mature monocytes decrease [14], [16], [20].
In conclusion, our data corroborate that FCM is able to add important prognostic information to all the prognostic scores studied. It is feasible in all patients, is especially useful in cases with a normal karyotype and gives more detailed information about the maturation abnormalities and the number of myeloblasts. Immunophenotyping is complementary, and should be used together with the clinical scores for more precise risk stratification of MDS patients.
Author Contributions
Conceived and designed the experiments: SCRA KM ILM. Performed the experiments: SCRA FT GH PMC ILM. Analyzed the data: SCRA KM ILM. Contributed reagents/materials/analysis tools: STOS FT SCRA. Wrote the paper: SCRA FT KM ILM.
References
- 1. Bennett JM, Catovsky D, Daniel MT, Flandrin G, Galton DAG, et al. (1982) Proposals for the classification of the myelodysplastic syndromes. Br J Haematol 51: 189–199.
- 2.
Swerdlow S, Camp E, Harris N, Jaffe ES, Pileri SA, et al. (2008) WHO classification of tumors of haematopoietic and lymphoid tissues. Lyon: IARC.
- 3. Yin CC, Medeiros LJ, Bueso-Ramos CE (2010) Recent advances in the diagnosis and classification of myeloid neoplasms-comments on the 2008 WHO classification. Int J Lab Hematol 32: 461–476.
- 4. Valent P, Horny HP, Bennett JM, Fonatsch C, Germing U, et al. (2007) Definitions and standards in the diagnosis and treatment of the myelodysplastic syndromes: Consensus statements and report from a working conference. Leuk Res 31: 727–739.
- 5. Lorand-Metze I, Pinheiro MP, Ribeiro E, de Paula EV, Metze K (2004) Factors influencing survival in myelodysplastic syndromes in a Brazilian population: Comparison of FAB and WHO classifications. Leuk Res 28: 587–594.
- 6. Greenberg P, Cox C, LeBeau MM, Fenaux P, Morel P, et al. (1997) International scoring system for evaluation in MDS. Blood 89: 2079–2088.
- 7. Schanz J, Tüchler H, Solé F, Mallo M, Luño E, et al. (2012) New comprehensive cytogenetic scoring system for primary myelodysplastic syndromes (MDS) and oligoblastic acute myeloid leukemia after MDS derived from an international database merge. J Clin Oncol 30: 820–829.
- 8. Pardanani A, Tefferi AC (2012) Cytogenetic risk stratification in myelodysplastic syndromes: are we there yet? J Clin Oncol 30: 2703–2704.
- 9. Greenberg P, Tuechler H, Schanz J, Sanz G, Garcia-Manero G, et al. (2012) Revised International Prognostic Scoring System for myelodysplastic syndromes. Blood 120: 2454–2465.
- 10. Malcovati L, Della Porta MG, Strupp C, Ambaglio I, Kuendgen A, et al. (2011) Impact of the degree of anemia on the outcome of patients with myelodysplastic syndrome and its integration into the WHO classification-based Prognostic Scoring System (WPSS). Haematologica 96: 1433–1440.
- 11. Thol F, Yun H, Sonntang AK, Damm F, Weissinger EM, et al. (2012) Prognostic significance of combined MN1, ERG, BAALC, and EV11 (MEBE) expression in patients with myelodysplastic syndromes. Ann Hematol 91: 1221–1233.
- 12. Bejar R, Stevenson K, Abdel-Wahab O, Galili N, Nilsson B, et al. (2012) Clinical effect of point mutations in myelodysplastic syndromes. N Engl J Med 364: 2496–2506.
- 13. Stetler-Stevenson M (2009) Flow cytometric immunophenotyping: emerging as an important diagnostic tool in the evaluation of cytopenic patients. Leuk Res 33: 1020–1021.
- 14. Lorand-Metze I, Ribeiro E, Lima CSP, Batista LS, Metze K (2007) Detection of hematopoietic maturation abnormalities by flow cytometry in myelodysplastic syndromes and its utility for the differential diagnosis with non-clonal disorders. Leuk Res 31: 147–155.
- 15. van de Loosdrecht AA, Alhan C, Béné MC, Della Porta MG, Dräger AM, et al. (2009) Standardization of flow cytometry in myelodysplastic syndromes: report from the first European LeukemiaNet working conference on flow cytometric in myelodysplastic syndromes. Haematologica 94: 1124–1134.
- 16. Matarraz S, Lopez A, Barrena S, Fernandez C, Jensen E, et al. (2008) The immunophenotype of different immature, myeloid and B-cell lineage-committed CD34+ hematopoietic cells allows discrimination between normal/reactive and myelodysplastic syndrome precursors. Leukemia 22: 1175–1183.
- 17. Reis SC, Traina F, Metze K, Saad ST, Lorand-Metze I (2009) Variation of bone marrow CD34+ cell subsets in myelodysplastic syndromes according WHO types. Neoplasma 56: 435–440.
- 18. Chu SC, Wang TF, Li CC, Kao RH, Li DK, et al. (2011) Flow cytometric scoring system as a diagnostic and prognostic tool in myelodysplastic syndromes. Leuk Res 35: 868–873.
- 19. Ribeiro E, Matarraz SS, Santiago M, Lima CSP, Metze K, et al. (2006) Maturation-associated immnophenotypic abnormalities in bone marrow B-lymphocytes in myelodysplastic syndromes. Leuk Res 30: 9–16.
- 20. Reis-Alves SC, Traina F, Saad ST, Metze K, Lorand-Metze I (2010) The impact of several phenotypic features at diagnosis on survival of patients with myelodysplastic syndromes. Neoplasma 57: 530–536.
- 21. Lorand-Metze I, Califani SMV, Ribeiro E, Lima CSP, Metze K (2008) The prognostic value of maturation-associated phenotypic abnormalities in myelodysplastic syndromes. Leuk Res 32: 211–213.
- 22. Ossenkoppele GJ, van de Loosdrecht AA, Schuurhuis GJ (2011) Review of the relevance of aberrant antigen expression by flow cytometry in myeloid neoplasms. Br J Haematol 153: 421–436.
- 23. Westers TM, Ireland R, Kern W, Alhan C, Balleisen JS, et al. (2012) Standardization of flow cytometry in myelodysplastic syndromes: a report from an international consortium and the European LeukemiaNet Working Group (2012) Leukemia. 26: 1730–1741.
- 24.
Shaffer LG, Slovak ML, Campbell LJ, eds. (2009) An International System for Human Cytogenetic Nomenclature: Recommendations of the International Standing Committee on Human Cytogenetic Nomenclature. Basel, Switzerland. Karger.
- 25. Oliveira GB, Pereira FG, Metze K, Lorand-Metze I (2001) Spontaneous apoptosis in chronic lymphocytic leukemia and its relationship to clinical and cell kinetic parameters. Cytometry 6: 329–335.
- 26. Bedin V, Adam RL, Sá BCS, Landman G, Metze K (2010) Fractal dimension is an independent prognostic factor for survival in melanoma. BMC Cancer 10: 260.
- 27. Lourenço GJ, Lorand-Metze I, Delamain MT, Miranda ECM, Kameo R, et al. (2010) Polymorphisms of glutathione S-transferase mu 1, theta 1 and pi 1 genes and prognosis in Hodgkin lymphoma. Leuk & Lymph 51: 2215–2221.
- 28. Ferro DP, Falconi MA, Adam RL, Ortega MM, Lima CSP, et al. (2011) Fractal Characteristics of May-Grünwald-Giemsa stained Chromatin are Independent Prognostic Factors for Survival in Multiple Myeloma. PLoS One 6(6): e20706.
- 29.
Altman DG, Royston P. (2000) What do we mean by validating a prognostic model? Statist Med 19;453–473.
- 30. Choodari-Oskooei B, Royston P, Parmar MKB (2010) A simulation study of predictive ability measures in a survival model I: explained variation measures. Statist Med 31: 2627–2643.
- 31. Sauerbrei W, Schumacher M (1992) A bootstrap resampling procedure for model building: application to the Cox regression model. Statist Med. 11: 2093–109.
- 32. Elston LB, Sueiro FAR, Cavalcanti JN, Metze K (2009) The Importance of the Mitotic Index as a Prognostic Factor for Survival of Canine Cutaneous Mast Cell Tumors: A Validation Study. Vet Pathol 46: 362–364.
- 33. Rybka MO, Cintra ML, de Souza EM, Metze K (2008) Density of dendritic cells around basal cell carcinoma is related to tumor size, anatomical site and stromal characteristics, and might be responsible for the response to topical therapy. International Journal of Dermatology 47: 1240–1244.
- 34. Adam RL, Silva RC, Pereira FG, Leite NJ, Lorand-Metze I, et al. (2006) The fractal dimension of nuclear chromatin as a prognostic factor in acute precursor B lymphoblastic leukemia. Cell Oncol. 28(1–2): 55–59.
- 35. Pinheiro RF, Metze K, Silva MR, Chauffaille ML (2009) The ambiguous role of interferon regulatory factor-1 (IRF-1) immunoexpression in myelodysplastic syndrome. Leuk Res. 33: 1308–1312.
- 36. Delamain MT, Metze K, Marques JF Jr, Reis AR, de Souza CA, et al. (2006) Optimization of CD34+ collection for autologous transplantation using the evolution of peripheral blood cell counts after mobilization with chemotherapy and G-CSF. Transfus Apher Sci. 34: 33–40.
- 37. Delamain MT, Marques JF Jr, de Souza CA, Lorand-Metze I, Metze K (2008) An algorithm based on peripheral CD34+ cells and hemoglobin concentration provides a better optimization of apheresis than the application of a fixed CD34 threshold. Transfusion. 48: 1133–1137.
- 38. Clark TG, Bradburn MJ, Love SB, Altman DG (2003) Survival Analysis Part IV: Further concepts and methods in survival analysis Br J Cancer. 89: 781–786.
- 39.
Smith B, Ryan MAK (2003) Survival analysis using Cox proportional hazards modeling for single and multiple event time data. Cary: SAS Institute, Inc. 254–28. Proceedings of the twenty-eighth annual SAS users group international conference.
- 40. Shao L-l, Zhang L, Hou Y, Yu S, Liu X-g, et al. (2012) Th22 cells as well as Th17 cells expand differentially in patients with early-stage and late-stage myelodysplastic syndrome. PLoS ONE 7(12): e51339.
- 41. Vido JR, Adam RL, Lorand-Metze I, Metze K (2011) Computerized texture analysis of atypical immature myeloid precursors in patients with myelodysplastic syndromes: an entity between blasts and promyelocytes. Diagn Pathol 6: 93.
- 42. Matarraz S, Teodosio C, Fernandez C, Albors M, Jara-Acevedo M, et al. (2012) The proliferation index of specific bone marrow cell compartments from myelodysplastic syndromes is associated with the diagnostic and patient outcome. PLOS ONE 7(8): e44321.
- 43. Ribeiro E, Lima CSP, Metze K, Lorand-Metze I (2004) Flow cytometric analysis of the expression of Fas/Fasl in bone marrow CD34+ cells in myelodysplastic syndromes: relation to disease progression. Leukemia & Lymphoma 45: 309–313.
- 44. Hellström-Lindberg E, Malcovati L (2008) Supportive care and use of hematopoietic growth factors in myelodysplastic syndromes. Semin Hematol 45: 14–22.
- 45. van de Loosdrecht AA, Westers TM, Westra AH, Drager AM, van der Velden V, et al. (2008) Identification of distinct prognostic subgroups in low- and intermediate-1-risk myelodysplastic syndromes by flow cytometry. Blood 111: 1067–1077.
- 46. Maftoun-Banankhah S, Maleki A, Karandikar NJ, Arbini AA, Fuda FS, et al. (2008) Multiparameter flow cytometric analysis reveals low percentage of bone marrow hematogones in myelodysplastic syndromes. Am J Clin Pathol 129: 300–308.
- 47. Wells DA, Benesch M, Loken MR, Vallejo C, Myerson D, et al. (2003) Myeloid and monocytic dispoiesis as determinated by flow cytometry scoring in myelodysplastic syndromes correlates with the IPSS and with outcome after hemopoietic stem cell transplantation. Blood 102: 394–405.
- 48. Gaipa G, Coustan-Smith E, Todisco E, Maglia O, Biondi A, et al. (2002) Characterization of CD34+, CD13+ CD33− Cells, a rare subset of immature human hematopoietic cells. Haematologica 87: 347–356.
- 49. Doehring LC, Heeger C, Aherrahrou Z, Kaczmarek PM, Erdmann J, et al. (2010) Myeloid CD34+CD13+ precursor cells transdifferentiate into chondrocyte-like cells in atherosclerotic intimal calcification Am J Pathol. 177: 473–480.
- 50. Metze K (2011) Pitfalls in the assessment of prognostic factors. Lancet Oncol 12: 1095–1096.
- 51. Metze K (2011) Dichotomizing continuous prognostic factors can cause paradoxical results in survival models. J Am Coll Surg 212: 132–134.
- 52. Metze K (2008) Dichotomization of continuous data - a pitfall in prognostic factor studies. Pathol Res Pract 204: 213–214.
- 53. Senent L, Arenillas L, Luño E, Ruiz JC, Sanz G, et al. (2013) Reproducibility of the World Health Organization 2008 criteria for myelodysplastic syndromes. Haematologica 98: 568–575.
- 54. Hiddemann W, Clarkson BD, Büchner TH, Melamed MR, Andreef M (1982) Bone marrow cell count per cubic millimeterbone marrow: a new parameter for quantitating therapy-induced cytoreduction in acute leukemia. Blood 59: 216–225.
- 55. Satoh C, Dan K, Yamashita T, Jo R, Tamura H, et al. (2008) Flow cytometric parameters with little interexaminer variability for diagnosing low-grade myelodysplastic syndromes. Leuk Res 32: 699–707.
- 56. Kern W, Haferlach C, Schnittinger S, Haferlach T (2010) Clinical utility of multiparameter flow cytometry in the diagnosis of 1012 patients with suspected myelodysplastic syndrome. Cancer 116: 4549–4563.