J Cancer 2020; 11(2):508-519. doi:10.7150/jca.30363

Research Paper

Evaluation of different scoring systems and gene mutations for the prognosis of myelodysplastic syndrome (MDS) in Chinese population

Meng-yi Du1, Min Xu1, Jun Deng1, Lin Liu1, Tao Guo1, Ling-hui Xia1, Yu Hu1,2* Corresponding address, Heng Mei1,2* Corresponding address

1. Institute of Hematology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430022, China.
2. Collaborative Innovation Center of Hematology, Huazhong University of Science and Technology, Wuhan, Hubei 430022, PR China.
* Heng Mei and Yu Hu contributed equally to this study.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Du My, Xu M, Deng J, Liu L, Guo T, Xia Lh, Hu Y, Mei H. Evaluation of different scoring systems and gene mutations for the prognosis of myelodysplastic syndrome (MDS) in Chinese population. J Cancer 2020; 11(2):508-519. doi:10.7150/jca.30363. Available from http://www.jcancer.org/v11p0508.htm

File import instruction

Abstract

MDS is a heterogeneous disease with diverse clinical manifestations, and an effective prognostic evaluation tool for MDS patients is needed. To achieve more accurate prognosis assessment for Chinese MDS patients, here we examined several scoring systems and explored the implications of gene mutations. The prognostic conditions were stratified against three different score systems (International Prognostic Scoring System (IPSS), WHO Prognostic Scoring System (WPSS), and Revised International Prognostic Scoring System (IPSS-R)) were retrospectively applied to 110 de novo MDS patients in study cohort in our hospital and the prognostic conditions were stratified respectively. IPSS-R out-performed the others, since it had less overlaps in survival curve, especially in the relatively low-risk group. Furthermore, genetic mutations were identified in 84 out of 110 patients and their association with overall survival (OS) were determined. Among them, sixty-three percent patients had at least one-point mutation, including thirty-five patients with normal karyotypes. The presence of TP53 mutations, but not TET2, DNMT3A or ASXL1 mutations was significantly correlated with shorter OS. A new model incorporating IPSS-R and TP53 mutations into survival analysis was proposed, and the prognostic value of this model was validated to be predominant in a 190-primary MDS patient independent cohort. Our data suggested that IPSS-R was more suitable for Chinese population. Attentions should be paid to the unfavourable mutations that might exert impact on the survival, especially in patients with relatively low risk.

Keywords: myelodysplastic syndrome, risk assessment, prognosis, gene mutations

Introduction

The myelodysplastic syndrome (MDS) represents a heterogeneous group of clonal hematopoietic disorders with diverse clinical manifestations. Approximately, 20%-30% of MDS patients may progress to acute myeloid leukemia [1, 2]. Clinically, accurate prediction of its prognosis is the cornerstone of guiding optimal MDS treatment. Prognostic models are useful tools for assessing risk of leukemic transformation, predicting life expectancy and making treatment plans [3]. Over the past two decades, three different scoring systems (IPSS[4], WPSS[5] and IPSS-R [6-8]) have been developed and successfully used worldwide in clinical practice for the evaluation of MDS prognosis[9]. Nonetheless, each system has its limitations and great divergences exists in various institutions when choosing a system for application [10]. All three systems cannot accurately assess cytogenetically normal patients at relatively low-risk[8, 11]. Moreover, MDS prognosis varies greatly with genetic profiles and races[12-15]. In this study, we aimed to identify the best scoring system for Chinese MDS patients by comparing the abovementioned ones. Importantly, we determined genetic mutations in cytogenetically normal patients to provide more accurate prognostic stratification.

Materials and methods

General information

A total of 300 patients, 110 for study cohort and 190 for validation cohort, with primary MDS were enrolled from the Institute of Hematology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China. Patients diagnosed with MDS were retrospectively collected from Jan. 2014 to Dec. 2018. This study has been approved by the Ethics Committees of Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China. Written informed consents were obtained from the patients or, in some special cases, from their relatives, and the study was conducted in strict accordance to the Helsinki Declaration (2013). The follow-up started from the date of enrolment to Jun.2017 for study cohort or Jun. 2019 for validation cohort or death and was conducted by telephoning patients or their attending doctors. In our series, no case was lost in the follow-up study. Patient characteristics were collected, including age, sex, initial peripheral blood data (absolute neutrophil count (ANC), hemoglobin (HB), platelet count (PLT)), initial percentage of bone marrow (BM) blasts, karyotype, WHO classification and survival status. Data concerning treatments received by patients during the follow-up period were recorded.

Diagnostic criteria

MDS diagnosis was made if (1) stable cytopenia ≥ 6 months (if stable cytopenias lasted for only 2 months, a specific karyotype or bi-lineage dysplasia were needed); (2) other potential disorders that are likely to trigger dysplasia and/or cytopenia were ruled out; (3) equal or greater than one of the three following requirements were satisfied: ① dysplasia (one or more of the 3 major bone marrow lineages≥10%); ② blast cell count ranging between 5% to 19%; ③ a specific MDS-associated karyotype [e.g., del(5q), del(20q), +8, or -7/del(7q)].

Furthermore, we used several other indicators that might help diagnose MDS, including aberrant immune-phenotypes, bone marrow abnormalities, or the presence of some molecular markers, signs or manifestations (i.e. abnormal CD34 antigen expression, fibrosis, dysplastic megakaryocytes, atypical localization of immature progenitors, and myeloid clonality)[16].

Methods

Three scoring systems (IPSS/WPSS/IPSS-R) were assessed in 110 study-cohort patients accordingly, in terms of mortality and median survival time. Then, the next genetic sequencing (NGS) tests were done in 84 out of 110 patients. Incidence of gene mutations and pattern of co-existing mutations were calculated, the associations of mutations with phenotypes and OS were determined. Moreover, the validation of a new prognostic model established by professor HOU and his colleges [17](HOU model) was taken, and an integrated scoring system based on our data comprehending IPSS-R and TP53 mutations (DU model) into survival analysis was proposed by us. The prognostic value of IPSS-R, HOU model and our model was then tested in the validation cohort of 190 primary MDS patients. Blood test, bone marrow cell morphology, bone marrow biopsy, cytogenetics and other tests were conducted in the patients. Patients in all groups received the same examinations and tests.

Grouping

According to the WHO (2008) classification [18], patients in our serious were classified into 5 categories in terms of myelogram and blood picture (Table 1). We didn't divide patients into groups according to WHO (2016) classification[19], because most of them were diagnosed before the year of 2016. Additionally, patients were grouped according to IPSS, WPSS and IPSS-R, by following the guidance of the 2017.V2 NCCN (Table 2). Relatively low-risk group referred to low-risk group and intermediate-1 group against IPSS; very low-risk group, low-risk group and intermediate risk group against WPSS; very low-risk group, low-risk group and intermediate risk group (IPSS-R intermediate patients may be managed as lower risk if their score is ≤3.5) against IPSS-R. The other subgroups against individual systems were designated as relatively high-risk group. All data used to calculate the risk stratification of the prognosis was the data at the time of diagnosis.

Gene sequencing

2-3 ml bone marrow was collected from 84 of 110 study-cohort and 190 validation-cohort MDS patients and then heparinized. Mononuclear cells were harvested from this bone marrow samples by using Ficoll solution. Then DNA was extracted from mononuclear cells using DNA extraction kit (Tiangen Biochemical Technology Co. Ltd, Beijing, China), and preserved at 2-8°C.

 Table 1 

Clinical characteristics of the 110 Chinese MDS patients in the study cohort, 84 patients that received next generation sequencing test, the groups of patients divided by TP53 mutation status, and the 190-patients validation cohort.

Study cohortNGS a patient in study cohortTP53-MUT
in study cohort
TP53-WILD
in study cohort
P value e
Validation cohort
Patients110841272-190
Age(year)49.5±15.850.2±1.8551.5±4.4850.0±2.030.7852.9±15.5
HB(g/l)69.5(58.0-92.8)93.2(67.5-110.5)62.5(51.0-73.0)68.5(59.0-96.8)0.0876.4(63.0-112.1)
PLT(*10^9/l)45.5(21.0-121.1)50.5(21.0-138.5)32.5(18.5-63.3)58(21.0-149.3)0.0849.2(38.1-130.2)
ANC(*10^9/l)1.12(0.63-2.13)1.20(0.64-2.43)0.73(0.27-1.28)1.44(0.72-2.51)0.021.45(0.78-2.77)
Classification b0.58
RCUD c970713
RARS
RCMD
2
26
2
17
012
16
8
43
RAEB-I282442049
RAEB-II383062466
MDS-U74139
5q-00002
All110841272190
Karyotype d0.001
Normal715335083
1 abnormal231621462
2 abnormal321111
≥3 abnormal13136734
IPSS0.001
Low650531
Intermediate-1493523350
Intermediate-2463653140
High985369
WPSS0.23
Very Low210117
Low870717
Intermediate271911839
High524273563
Very High14124845
MDS-U73039
IPSS-R0.002
Very Low440419
Low322202226
Intermediate342442049
High262021832
Very High14146864
OS11.0(4.0-30.0)9.0(3.0-25.5)2.5(1.3-4.8)11.5(5.0-33.8)<0.00112.0(5.0-27.0)
State0.01
live5946244112
die433192269
sAML87169

RCUD: refractory cytopenia with unilineage dysplasia; RARS: refractory anemia with ring sideroblast; RCMD: refractory cytopenia with multilineage dysplasia; RAEB-I: refractory anemia with excess blasts-I; RAEB-II: refractory anemia with excess blasts-II; MDS-u: myelodysplastic syndromes, unclassifiable; 5q-: MDS associated with isolated del (5q); sAML: secondary acute myelocytic leukemia; PB: peripheral blood; HB: hemoglobin; PLT: platelet count; ANC: absolute neutrophil count

a: 84 patients who were given target next gene sequencing (NGS) out of all 110 study cohort patients

b: all 110 Chinese MDS patients were classified by WHO (2008) criteria

c: RCUD includes refractory anemia (RA); refractory neutropenia (RN); refractory thrombocytopenia (RT)

d: the risk group of cytogenetics in IPSS and IPSS-R are different. For IPSS, Cytogenetics: Good = normal, -Y alone, del(5q) alone, del(20q) alone; Poor = complex (≥3 abnormalities) or chromosome 7 anomalies; Intermediate = other abnormalities. For IPSS-R, Cytogenetic risks: Very good = -Y, del(11q); Good = normal, del(5q), del(12p), del(20q), double including del(5q); Intermediate = del(7q), +8, +19, i(17q), any other single or double independent clones; Poor = -7, inv(3)/t(3q)/del(3q), double including -7/del(7q), complex: (3 abnormalities); Very poor = complex: >3 abnormalities. Before the prognostic scoring system, which could fit Chinese better, were identified, we divided our patients only by number of cytogenetic abnormalities, other than any of IPSS or IPSS-R.

e: Proportions of TP53-mut and TP53-wild patients were compared by utilizing Chi-square test. A P<0.05 was considered to be statistically significant.

A total of 25 well-known genes from the MDS genomes (Figure 2.) were analysed by targeted deep sequencing. Genomic DNA from each sample was sequenced. Briefly, KAPA hyper library kit and corresponding adaptors were used to construct library, and then Agencourt AMPure XP reagent (Beckman Coulter, Inc.) was used to purify magnetic beads. The sample was mixed with MiSeq Reagent Kit V3 (600 cycles)/MiSeq Reagent Kit V2 (500 cycles) or FC-420-1004 MiniSeq Mid Output Kit (300 cycles)/Miniseq High Output Kit (300 cycles) (Yuanqi Biopharmaceutical Technology Co. Ltd., Shanghai, China), and then sequenced on an illumina MiSeq or illumina MiniSeq sequencer (Illumina, California, American), respectively. Read pairs were aligned to Refseq hg19 (downloaded from UCSC Genome Browser) by Burrows-Wheeler Aligner (BWA) version 0.7.13-r1126. Samtools version 1.3 was used to generate chromosomal coordinate-sorted bam files. The mean depth of each sample was 1500x-2000x, with an average >90% of the target sequence which covered sufficiently deep for variant calling. The designed sequencing primers covered all mutational hot spots of the candidate genes. Samtools mpileup were applied to call single nucleotide variants (SNVs) and indels. Homemade pipeline was used to filter SNVs and indels detected by the above software, excluding: 1) mutations reported with low confidence; 2) mutations reported in 1000 Genomes (database of single nucleotide polymorphism (dbSNP) 138) as common single nucleotide polymorphisms (SNPs) and not included in COSMIC version v83. SNVs and indels were annotated using Annovar Software. All the somatic functional mutations, including nonsynonymous SNVs, frameshift or in-frame indels, stop-gain and stop-loss were obtained. Visual inspection was used to exclude potential false positive results. Finally, Integrative Genomics Viewer (http://www.igv.org/) software was used for data visualization and verification.

 Table 2 

The mortality and median overall survival of different groups stratified according to IPSS/ WPSS/ IPSS-R/HOU model/DU model respectively

THE OUTCOME OF DIFFERENT PROGNOSTIC SCORING SYSTEM
SYSTEMNO.RATEFATALMORTALI TYMEDIAN OS(m)95% CITarone-WareP
IPSS (STUDY COHORT)LLUL
low60.0520.335515.095.030.847<0.001
int-1490.45120.242011.828.2
int-2460.42220.4863.28.9
high90.0870.78204.9
all11043
WPSS (STUDU COHORT)
very low20.0210.50340.085.431.584<0.001
low80.0710.133325.9156.1
int270.2580.30721.449.3
high520.47230.4439.521.0
very high140.1390.6421.98.8
MDS-U a70.0610.14200.050.8
all11043
IPSS-R (STUDY COHORT)
very low40.0410.256012.0108.055.854<0.001
low320.2990.28309.250.8
int340.31130.38116.715.3
high260.24120.4652.57.5
very high140.1380.5720.83.2
all11043
IPSS-R (84 NGS of 110 STUDY COHORT)
very low40.0500.006012.0108.047.489<0.001
low220.2620.093512.058.0
int240.2980.3384.211.8
high200.24110.5552.87.2
very high140.17100.7120.83.2
all8431
HOU Model (84 NGS of 110 STUDY COHORT)
low210.2550.242418.130.010.3380.016
int70.0830.4312027.4
high260.31100.3883.013.0
very high300.36130.4350.79.3
all8431
DU Model (84 NGS of 110 STUDY COHORT)
very low40.0500.006012.0108.052.256<0.001
low220.2620.093512.058.0
int200.2460.30116.615.3
high240.29130.5441.96.1
very high140.17100.7120.83.2
all8431
IPSS-R (VALIDATION COHORT)
very low190.1030.166853.882.281.437<0.001
low260.1460.233321.844.2
int490.26160.33149.418.6
high320.17190.59117.814.2
very high640.34340.5364.47.6
all19078
HOU Model (VALIDATION COHORT)
low30.0010.33230.051.839.864<0.001
int510.27140.273312.054.0
high620.33250.40126.917.1
very high740.39380.5175.38.7
all19078
DU Model (VALIDATION COHORT)
very low190.1030.166853.882.287.109<0.001
low250.1360.243522.048.0
int470.25150.321612.619.4
high850.45440.5275.38.7
very high140.07100.7130.06.7
all19078

IPSS: International Prognostic Scoring System; WPSS: WHO Prognostic Scoring System; IPSS-R: Revised International Prognostic Scoring System; NGS patients: patients that received next generation sequencing test; HOU model: the prognostic model for MDS established by HOU and his colleague; DU model: the prognostic model for MDS established based on our cohort; LL: low-limit; UL: up-limit

a: the classification of MDS-U can't be stratified by WPSS

Treatments

The treatment for relatively low-risk group of patients primarily aimed to optimize their hematopoietic function and improve the quality of life. Patients with HB <60 g / L or severe anemia were infused with suspended erythrocytes. Patients with PLT <20*10^9 / L or suffering from active bleeding were given hemostatic drugs or transfused with platelets. Patients with isolated del (5q-) were administered thalidomide. Severe agranulocytosis was treated with colony-stimulating factor and iron, among others. The treatment for high-risk MDS was aimed to delay disease progression and prolong survival. Patients above 65 or in poor conditions were given demethylating drugs (decitabine). Patients under 65 in relatively stable conditions were treated with decitabine in combination with CAG regimen (aclarubicin, cytosine arabinoside, and granulocyte colony stimulating factor).

Statistical analysis

The comparison of median survival time was performed using the Tarone-Ware test. Survival analysis started from the date of diagnosis to the date of death or last contact. Survival curves were prepared by employing the Kaplan-Meier method. Spearman correlation analysis was used to ascertain how phenotypes correlated with mutations. Proportions of mutational co-existence were compared by utilizing Chi-square test. Cox proportional hazard regression model was used for univariable and multivariable analysis. A P<0.05 was considered to be statistically significant. SPSS 22.0 software package (version 22; IBM Corporation, Armonk NY) was used for statistical analysis.

Results

Follow-up results of study-cohort

The series included 60 males and 50 females, with a male-to-female ratio at 1.2:1. The median age of these patients was 49.5 years old (Standard Deviation (SD) =15.8), ranging from 13 to 89. All the patients were followed up for an average of 11 (4 ~ 30) months. Forty-three patients died, with a mortality of 39.0%. Eight cases converted into sAML, with the conversion rate at 7.3% (Table 1). Furthermore, we evaluated the outcome of different prognostic scorning system. We found that risk stratification using IPSS-R showed the best prognostic values as the worse the prognosis the higher the mortality in patients. In contrast, IPSS and WPSS were somehow deficient in predicting prognosis because high mortality was observed in the lowest risk group using these two systems respectively (Table 2).

Survival analysis

We next performed the survival analysis. We found that the median survival time was reduced progressively with the increased risk stratified according to IPSS, WPSS and IPSS-R, respectively (Table 2). The differences in the median survival time among strata within the system were of statistical significance. However, the median OS of the relatively low risk groups in all three scoring systems were significantly shorter than the referenced data in NCCN Guideline[1]. For example, the median OS of low-risk group in IPSS-R was 63 months in NCCN Guideline and 35 months in our study. Moreover, the survival curve showed that the differences in prognosis among relatively low-risk groups according to WPSS were not evident (Figure 1A). The possible explanation is that there were many intersections between relatively low risk groups, only the very-high risk group did not cross with others, making it difficult to differentiate the prognosis in the low-risk groups. As for IPSS, there were fewer intersections among subgroups, but crossover still existed between low risk group and intermediate-1 risk group (Figure 1B). Besides, early death existed in the lowest risk group of both IPSS and WPSS. In contrast, IPSS-R had least intersections with clear survival curves for each group, which rendered it easy to judge the prognostic differences among different risk groups (Figure 1C).

Incidence of gene mutations and pattern of co-existing mutations

To help improve the accuracy of risk stratification for MDS patients, 84 of all 110 study-cohort patients were subjected to gene sequencing. Fifty-three (63%) patients were found to have gene mutations, with 32 (29%) patients having merely one mutation, 15 (13.6%) having two mutations, 6 (5.5%) having no less than 3 mutations, and one patient having 9 gene mutations. The most frequently mutated genes were TET2, TP53, DNMT3A, ASXL1 and RUNX1 (Figure 2). We next analysis the distribution of the mutations and found that 82.8% were distributed in RCMD, RAEB-I and RAEB-II patients. The higher the risk of classification, the higher the quantity of TP53 mutations (Figure 3A). Moreover, TET2, TP53, DNMT3A and ASXL1 gene mutations occurred in high frequency in all three groups of patients (Figure 3B/C/D).

TET2, ASXL1 mutations were more likely to co-exist with other mutations, occurring at a ratio of 64.7% (χ2=20.516, P<0.0001) and 66.7% (χ2=19.329, P<0.0001), respectively. ASXL1 mutations tended to occur in patients with more mutations. In fact, 44.4% patients carrying ASXL1 mutation had 3 or more gene mutations (P<0.05). In contrast, with TP53 mutations the frequency of co-existence was 33.3% (χ2=12.726, P=0.013).

 Figure 1 

Kaplan-Meier curves for overall survival of 110 study cohort patients in our study with survival data. (A) Survival of patients stratify according to WPSS. (B) Survival of patients stratify according to IPSS. (C) Survival of patients stratify according to IPSS-R.

J Cancer Image (Click on the image to enlarge.)
 Figure 2 

Mutation spectrum of both study cohort and validation cohort. 53 out of 84 study-cohort and 124 out of 190 validation-cohort patients had at least one mutation of the 25 MDS-associated genes. Each column represents an individual patient sample, and each coloured cell represents a mutation of the gene.

J Cancer Image (Click on the image to enlarge.)
 Figure 3 

Gene mutation status in 84 study cohort patients that received next generation sequencing test. (A) Mutation status of five prime genes in 84 gene sequencing patients. (B) Mutation status of all MDS-related genes in RCMD patients. (C) Mutation status of all MDS-related genes in RAEB-I patients. (D) Mutation status of all MDS-related genes in RAEB-II patients.

J Cancer Image (Click on the image to enlarge.)

Association between phenotypes and frequent mutations

Next, we investigated the association between phenotypes and frequent mutations. Indeed, correlations were found between mutations and phenotypes, such as age, HB/PLT/ANC level, karyotype, classification and IPSS/WPSS/IPSS-R risk level. The results showed that TP53 was associated with thrombocytopenia (r= -0.25, P=0.02), complex karyotype (r=0.38, P<0.001) and higher IPSS/IPSS-R risk level (r=0.35, P=0.001, r=0.33, P=0.002). TET2 mutation was associated with older age (r=0.23, P=0.035) and ASXL1 mutation with lower ANC level (r= -0.21, P=0.049).

Associations between gene mutations and overall survival (OS)

The associations between gene mutations and OS was evaluated next. In univariate COX analysis, comparison between patients harboring TP53 mutations (14.3%) and those with wild-types (wt) showed that TP53 mutations were associated with a shorter overall survival (HR=5.67 [2.58-12.45], P<0.001; median OS, 2.5 months for patients with TP53 mutations vs. 11.5 months for those with wt). There was no statistical difference in the survival between patients who harbored mutated genes other than TP53 and the wild type patients (Figure 4A). Compared to wt patients, patients with TET2-mutation (20.3% of the patients; HR=0.68 [0.30-1.56], P=0.37), patients with DNMT3A-mutation (10.7% of the patients; HR=0.77 [0.27-2.17], P=0.62), or patients with ASXL1-mutation (10.7% of the patients; HR=0.54 [0.17-1.76], P=0.31) had no significant association with poor prognosis (Figure 4B/4C/4D). Taking both IPSS-R and TP53 mutation into multivariate analysis, the results showed that TP53 was significantly associated with overall survival (HR=5.25 [2.25- 12.25], P<0.001) regardless of the risk stratification. (Table 3) Addition of age to this model did not affect the overall results.

In addition, a new risk model was evaluated in 84 NGS-patients in our study. Hou et al[17] developed the risk model that incorporating the weighted coefficients of clinical and genetic factors (age ×0.025-IPSS-R lower risk group ×1.184 +CBL × 0.829+IDH2 ×0.829 +DNMT3A ×0.452 +ASXL1 × 0.442+TP53 ×2.254) for Chinese patients in 2018. Four risk groups were proposed: low (score <-0.5), intermediate (score -0.5~0.5), high (score 0.51~1.5) and very high (score >1.5). In our cohort, the median OS was 24, 12, 8 and 5 months for low, intermediate, high, and very high subgroups, respectively. The survival curve showed more intersections than IPSS-R (Figure 5A/B), and the mortality rate did not increase with risk level parallel.

 Table 3 

Univariate and multivariate analysis (Cox regression) for the overall survival in 84 study-cohort MDS patients that received next generation sequencing test

Univariate COX analysisPHR95 CI
TP53<0.0015.672.58-12.45
TET20.370.680.30-1.56
DNMT3A0.620.770.27-2.17
ASXL10.310.540.17-1.76
RUNX10.210.280.04-2.02
IPSS-R<0.0017.112.42-20.91
Multivariate COX analysisPBHR95 CI
All patientsTP53<0.0011.6585.252.25- 12.25
IPSS-R0.0011.8536.382.15- 18.93
Relatively low-risk patientsTP530.0201.6565.241.30-21.09
IPSS-R<0.0011.6585.252.06- 13.38
Relatively high-risk patientsTP530.0011.223.381.62- 7.07
IPSS-R0.0031.203.321.49- 7.38
 Figure 4 

Kaplan-Meier curves for overall survival in the 84 study cohort patients that received next generation sequencing test. (A) Survival of patients with no mutation, with TP53 mutations and with the other mutations. (B) Survival of patients with no mutation, with TET2 mutations and with the other mutations. (C) Survival of patients with no mutation, with DNMT3A mutations and with the other mutations. (D) Survival of patients with no mutation, with ASXL1 mutations and with the other mutations.

J Cancer Image (Click on the image to enlarge.)

To improve the risk stratification for MDS patients, an integrated scoring system based on our data incorporating IPSS-R and TP53 mutations into survival analysis was proposed. This model was developed incorporating the weighted coefficients of these two factors: TP53×1.658 + IPSS-R risk level×1.853. Five risk groups were proposed: very low (score < 2; n=4), low (score 2~5; n=21), intermediate (score 5.1~7; n=21), high (score 7.1~10; n =24), very high (score > 10; n=14). The median OS was 60, 35, 11, 4 and 2 months for very low, low, intermediate, high, and very high subgroups, respectively. This clinically relevant integrated scoring system divided the MDS patients into five groups with different clinical outcomes (P< 0.001) (Figure 5C). The median OS in relatively low-risk group of our model was longer than the same risk group of IPSS-R. Interestingly, patients who were in relatively low risk group according to IPSS-R but in relatively high-risk groups according to our new model were all TP53 mutated and cytogenetically normal.

To verify the prognostic value of IPSS-R, HOU model and our model, a validation cohort containing 190 next gene sequencing tested patients was then established. The validation cohort was made up by 104 males and 84 females with a median age of 52.9±15.5 years. Similarly, the patients concentrated on the classification of RCMD, RAEB-I and RAEB-II, and the overall mortality rate is also about 50%. Besides, the frequency of TP53/ TET2/ DNMT3A/ ASXL1/ RUNX1 gene mutation was high and the co-existing spectrum with TP53 mutation was narrow. Then, the three models (IPSS-R/ HOU model/ DU model) were assessed by these 190 patients accordingly. Our model was outperformed than IPSS-R and HOU model among OS, Kaplan-Meier curves and mortality rate. For IPSS-R, there was smaller discrepancy of OS between different risk stratified subgroups, and un-reasonable high mortality in high risk group (Table 2) than ours. For HOU model, a few patients were divided into low risk group (Table 2), which indicated that it could not figure out superior prognosis patients pretty well, let along more precise risk stratification of every subgroup for there was more crossover even than IPSS-R in Kaplan-Meier curves for overall survivals (Figure 5D/E/F).

Discussion

This study showed that IPSS-R could better predict the prognosis in terms of mortality and overall survival of MDS patients, especially relatively low-risk MDS patients. Previous studies demonstrated that IPSS-R out-performed IPSS and WPSS in the prediction of transformation to acute leukemia and progression free survival (PFS) [20, 21]. The better performance of IPSS-R might be attributed to finer stratification of Chromosomal karyotypes and quantification of peripheral blood loss. Initially, IPSS-R was used only for the establishment of primary patients before treatment, but not employed as a tool for therapy-related or prognosis monitoring of MDS[6, 22, 23]. However, studies found that IPSS-R worked well in treatment-related and post-treatment risk stratification, as well as dynamic prognosis monitoring in patients received lenalidomide or demethylating agents[23-26].

 Figure 5 

Kaplan-Meier curves for overall survival in MDS patients that received next generation sequencing test with survival data. Survival of 84 study-cohort patients: (A) according to IPSS-R. (B) according to HOU model. (C) according to DU model; Survival of 190 validation-cohort patients: (D) according to IPSS-R. (E) according to HOU model. (F) according to DU model.

J Cancer Image (Click on the image to enlarge.)

Nevertheless, can IPSS-R alone be used in clinical practice? Indeed, IPSS-R still has its limitations. For example, WPSS has been proved to be effective in the prediction of OS, PFS, cumulative incidence of relapse (CIR) and cumulative incidence of non-relapse mortality (CINRM) in the post-transplantation patients [27]. However, IPSS-R was only used in pre-transplantation assessment [28]. In addition, for relatively low risk, so far researchers haven't agreed on whether IPSS-R or Low Risk Prognostic Scoring System (LR-PSS) [29] can better predict the prognosis. Previous studies showed that the two scoring systems were not satisfactory in identifying unfavourable OS, and LR-PSS was comparatively more sensitive than IPSS-R[8, 11]. In our study, the median OS of relatively low risk group was much shorter than the data provided by NCCN Guideline. Because patients with relatively low risk were given symptomatic and supportive treatment rather than chemotherapy or HMAs, the effect of treatment level on outcome should be ruled out. Therefore, there might be patients with poor prognosis mixed into relatively low risk groups.

With the development of gene sequencing technologies, more and more gene mutations were identified to be associated with MDS. Several studies revealed that the mutation spectra of MDS and some of the mutations were intimately related to the prognosis of MDS [30, 31]. In this study, TP53 gene mutation was found to be associated with thrombocytopenia, complex karyotypes and reduced OS. It has been generally accepted that TP53 mutation is an independent factor for poor prognosis[32-34] in MDS patients. Moreover, we assessed the prognostic impact of TP53 mutations on IPSS-R in MDS patients. So far there have been few studies integrating molecular data into IPSS-R in MDS patients[17, 30, 33]. Haferlach et al.[30] utilized a combination of conventional factors (age, gender, and IPSS-R) and mutations in 14 genes (including TP53) as a novel prognostic model which had not been published in detail. The patients were stratified into four risk groups according to this model, which was demonstrated to be more accurate than IPSS-R in validation cohort. In a study of Nazha et al.[33], incorporation of three mutations (EZH2, SF3B1, and TP53) into IPSS-R could improve the predictive power in 508 USA patients with primary and secondary MDS. Recently, Hou et al.[17] developed an integrated risk-stratification model incorporating conventional risk factors (age and IPSS-R) and 5 gene mutations (CBL, IDH2, DNMT3A, ASXL1, TP53), which might be suitable for Chinese patients. Actually, the implication of mutations on clinical outcomes in Asian MDS patients may be different from that in western patients, due to difference in disease natures of MDS and racial back-ground between these two populations[13-15]. We next evaluated if the model established by HOU (HOU model) was suitable for MDS patients in our cohort. However, the results showed that patients could be better stratified by IPSS-R than the HOU model in aspect of survival analysis, median OS and mortality rate. Therefore, to improve prognostic prediction in Chinese MDS patients, we here proposed a new prognostic model incorporating TP53 mutations and IPSS-R based on our data. Our model can distinguish poor outcome patients with normal karyotype from real low-risk patients, thus indicating better prognostic power than IPSS-R, which have been verified in both study and validation group.

Apart from TP53, impact of most MDS-related gene mutations on the MDS prognosis remain controversial. One important reason is ethnic variation. For example, ASXL1 mutation did not predict outcome but only a trend of adverse prognosis of patients with MDS in our cohort and other Chinese population, in contrast, ASXL1 mutation was related to poor prognosis in European patients [12-15]. Among the mutated genes, TET2, SRSF2 and SETBP1 have been subjected to META analysis [35-37] for no unanimous effect of them on prognosis reached. Moreover, co-mutations of genes further complicate the analysis and the interpretation of their influence on MDS prognosis. Our present study showed that TET2, ASXL1 and DNMT3A mutations were common in MDS patients (with mutation ratios 20.3%, 10.7% and 10.7%, respectively) and they all had co-mutations, but they were not found to be correlated with MDS prognosis. Patients with TET2 or DNMT3A mutations showed a higher OS rate at the beginning of the disease course when compared with other groups, probably because a large part of our patients received decitabine as main therapy. Of note, some studies reported that patients with TET2 and DNMT3A mutations had higher rates of response to HMAs and those with ASXL1 mutation had longer OS [38-40]. Nonetheless, the impact of these mutations on the prognosis was moderate. Our further study will aim at identifying mutations that affect the efficacy of HMAs and prognosis of MDS and understanding the impacts they exert on the efficacy and prognosis.

The scoring systems of prognosis are designed to provide guides or information for better treatment. MDS patients with TET2 mutation could predict better response and a trend for longer OS to hypomethylating agents [38]. However, the OS showed completely no difference between those with and without TET2 mutation when treatment was not taken into consideration[35]. Therefore, the effect of gene mutation is varied due to non-uniform therapeutic strategies. Until now, most meta-analyses didn't take into consideration the impact of treatment protocols but only analysed the general effect of gene mutations on prognosis, which is an important reason why a consensus has not been reached. Future studies should, on the basis of finer stratification, focus on the relationship between gene mutations (and other factors) and different treatment protocols.

The limitation of our study is that the total number of patients were not that variety, and all patients involved in were from one center but not from multicenter. Our ongoing work also aims to explore the outcome of gene mutation on outcome when treated with specific treatment.

In summary, IPSS-R is a scoring system that well fits Chinese MDS patients. A population-specific scoring system that integrates the second-generation gene sequencing technologies can help achieve more accurate prognostic assessment and thereby provide guidance for more effective treatment for Chinese MDS patients.

Acknowledgements

We thank the Laboratory Department, the Medical Records Department, the Yuanqi Biopharmaceutical Technology Co. Ltd. (Shanghai, China) and all the clinical departments which provided us clinical and laboratory information of patients. The study was funded and supported by the National Natural Science Foundation of China (Grant No. 81570116; No. 81873434).

Competing Interests

The authors have declared that no competing interest exists.

References

1. Greenberg PL, Stone RM, Al-Kali A, Barta SK, Bejar R, Bennett JM. et al. Myelodysplastic Syndromes, Version 2.2017, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw. 2017;15:60-87

2. Ma X. Epidemiology of myelodysplastic syndromes. Am J Med. 2012;125:S2-5

3. Cazzola M, Malcovati L. Prognostic classification and risk assessment in myelodysplastic syndromes. Hematol Oncol Clin North Am. 2010;24:459-68

4. Greenberg P, Cox C, LeBeau MM, Fenaux P, Morel P, Sanz G. et al. International scoring system for evaluating prognosis in myelodysplastic syndromes. Blood. 1997;89:2079-88

5. Malcovati L, Della PMG, Strupp C, Ambaglio I, Kuendgen A, Nachtkamp K. et al. 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. 2011;96:1433-40

6. Greenberg PL, Tuechler H, Schanz J, Sanz G, Garcia-Manero G, Sole F. et al. Revised international prognostic scoring system for myelodysplastic syndromes. Blood. 2012;120:2454-65

7. van Spronsen MF, Ossenkoppele GJ, Holman R, van de Loosdrecht AA. Improved risk stratification by the integration of the revised international prognostic scoring system with the myelodysplastic syndromes comorbidity index. Eur J Cancer. 2014;50:3198-205

8. Zeidan AM, Sekeres MA, Wang X-F, Al Ali N, Garcia-Manero G, Steensma DP. et al. Comparing the prognostic value of risk stratifying models for patients with lower-risk myelodysplastic syndromes: Is one model better?. American Journal of Hematology. 2015;90:1036-40

9. Zeidan AM, Komrokji RS. There's risk, and then there's risk: The latest clinical prognostic risk stratification models in myelodysplastic syndromes. Curr Hematol Malig Rep. 2013;8:351-60

10. Germing U, Kundgen A. Prognostic scoring systems in MDS. Leuk Res. 2012;36:1463-9

11. Valcarcel D, Sanz G, Ortega M, Nomdedeu B, Luno E, Diez-Campelo M. et al. Use of newer prognostic indices for patients with myelodysplastic syndromes in the low and intermediate-1 risk categories: a population-based study. Lancet Haematol. 2015;2:e260-6

12. Chen TC, Hou HA, Chou WC, Tang JL, Kuo YY, Chen CY. et al. Dynamics of ASXL1 mutation and other associated genetic alterations during disease progression in patients with primary myelodysplastic syndrome. Blood Cancer Journal. 2014;4:e177

13. Chen B, Zhao WL, Jin J, Xue YQ, Cheng X, Chen XT. et al. Clinical and cytogenetic features of 508 Chinese patients with myelodysplastic syndrome and comparison with those in Western countries. Leukemia. 2005;19:767-75

14. Lee JH, Lee JH, Shin YR, Lee JS, Kim WK, Chi HS. et al. Application of different prognostic scoring systems and comparison of the FAB and WHO classifications in Korean patients with myelodysplastic syndrome. Leukemia. 2003;17:305

15. Yang YT, Hou HA, Liu CY, Lin CC, Chou WC, Lee FY. et al. IPSS-R in 555 Taiwanese patients with primary MDS: Integration of monosomal karyotype can better risk-stratify the patients. AM J HEMATOL. 2014;89:E142-E9

16. Valent P, Horny HP, Bennett JM, Fonatsch C, Germing U, Greenberg P. et al. Definitions and standards in the diagnosis and treatment of the myelodysplastic syndromes: Consensus statements and report from a working conference. Leuk Res. 2007;31:727-36

17. Hou HA, Tsai CH, Lin CC, Chou WC, Kuo YY, Liu CY. et al. Incorporation of mutations in five genes in the revised International Prognostic Scoring System can improve risk stratification in the patients with myelodysplastic syndrome. Blood Cancer J. 2018;8:39

18. Vardiman JW, Thiele J, Arber DA, Brunning RD, Borowitz MJ, Porwit A. et al. The 2008 revision of the World Health Organization (WHO) classification of myeloid neoplasms and acute leukemia: rationale and important changes. Blood. 2009;114:937-51

19. Arber DA, Orazi A, Hasserjian R, Thiele J, Borowitz MJ, Beau MML. et al. The 2016 revision to the World Health Organization classification of myeloid neoplasms and acute leukemia. Blood. 2016;127:2391

20. de Swart L, Smith A, Johnston TW, Haase D, Droste J, Fenaux P. et al. Validation of the revised international prognostic scoring system (IPSS-R) in patients with lower-risk myelodysplastic syndromes: a report from the prospective European LeukaemiaNet MDS (EUMDS) registry. Br J Haematol. 2015;170:372-83

21. Della PM, Tuechler H, Malcovati L, Schanz J, Sanz G, Garcia-Manero G. et al. Validation of WHO classification-based Prognostic Scoring System (WPSS) for myelodysplastic syndromes and comparison with the revised International Prognostic Scoring System (IPSS-R). A study of the International Working Group for Prognosis in Myelodysplasia (IWG-PM). Leukemia. 2015;29:1502-13

22. Killick SB, Carter C, Culligan D, Dalley C, Das-Gupta E, Drummond M. et al. Guidelines for the diagnosis and management of adult myelodysplastic syndromes. Br J Haematol. 2014;164:503-25

23. Pfeilstocker M, Tuchler H, Schonmetzler A, Nosslinger T, Pittermann E. Time changes in predictive power of established and recently proposed clinical, cytogenetical and comorbidity scores for Myelodysplastic Syndromes. Leuk Res. 2012;36:132-9

24. Lamarque M, Raynaud S, Itzykson R, Thepot S, Quesnel B, Dreyfus F. et al. The revised IPSS is a powerful tool to evaluate the outcome of MDS patients treated with azacitidine: the GFM experience. Blood. 2012;120:5084-5

25. Mishra A, Corrales-Yepez M, Ali NA, Kharfan-Dabaja M, Padron E, Zhang L. et al. Validation of the revised International Prognostic Scoring System in treated patients with myelodysplastic syndromes. Am J Hematol. 2013;88:566-70

26. Zeidan AM, Al AN, Barnard J, Padron E, Lancet JE, Sekeres MA. et al. Comparison of clinical outcomes and prognostic utility of risk stratification tools in patients with therapy-related vs de novo myelodysplastic syndromes: a report on behalf of the MDS Clinical Research Consortium. Leukemia. 2017;31:1391-7

27. Ma L, Hao S, Diong C, Goh YT, Gopalakrishnan S, Ho A. et al. WPSS is a strong prognostic indicator for clinical outcome of allogeneic transplant for myelodysplastic syndrome in Southeast Asian patients. Ann Hematol. 2015;94:761-9

28. Alzahrani M, Power M, Abou MY, Barnett M, Broady R, Forrest D. et al. Improving Revised International Prognostic Scoring System Pre-Allogeneic Stem Cell Transplantation Does Not Translate Into Better Post-Transplantation Outcomes for Patients with Myelodysplastic Syndromes: A Single-Center Experience. Biol Blood Marrow Transplant. 2018;24:1209-15

29. Garcia-Manero G, Shan J, Faderl S, Cortes J, Ravandi F, Borthakur G. et al. A prognostic score for patients with lower risk myelodysplastic syndrome. Leukemia. 2008;22:538-43

30. Haferlach T, Nagata Y, Grossmann V, Okuno Y, Bacher U, Nagae G. et al. Landscape of genetic lesions in 944 patients with myelodysplastic syndromes. Leukemia. 2014;28:241-7

31. Papaemmanuil E, Gerstung M, Malcovati L, Tauro S, Gundem G, Van Loo P. et al. Clinical and biological implications of driver mutations in myelodysplastic syndromes. Blood. 2013;122:3616-27

32. Bejar R, Stevenson KE, Caughey B, Lindsley RC, Mar BG, Stojanov P. et al. Somatic mutations predict poor outcome in patients with myelodysplastic syndrome after hematopoietic stem-cell transplantation. J Clin Oncol. 2014;32:2691-8

33. Nazha A, Narkhede M, Radivoyevitch T, Seastone DJ, Patel BJ, Gerds AT. et al. Incorporation of molecular data into the Revised International Prognostic Scoring System in treated patients with myelodysplastic syndromes. Leukemia. 2016;30:2214-20

34. Sallman DA, Komrokji R, Vaupel C, Cluzeau T, Geyer SM, McGraw KL. et al. Impact of TP53 mutation variant allele frequency on phenotype and outcomes in myelodysplastic syndromes. Leukemia. 2016;30:666-73

35. Lin Y, Lin Z, Cheng K, Fang Z, Li Z, Luo Y. et al. Prognostic role of TET2 deficiency in myelodysplastic syndromes: A meta-analysis. Oncotarget. 2017;8:43295-305

36. Zheng X, Zhan Z, Naren D, Li J, Yan T, Gong Y. Prognostic value of SRSF2 mutations in patients with de novo myelodysplastic syndromes: A meta-analysis. PLoS One. 2017;12:e0185053

37. Shou LH, Cao D, Dong XH, Fang Q, Wu Y, Zhang Y. et al. Prognostic significance of SETBP1 mutations in myelodysplastic syndromes, chronic myelomonocytic leukemia, and chronic neutrophilic leukemia: A meta-analysis. PLoS One. 2017;12:e0171608

38. Bejar R, Lord A, Stevenson K, Bar-Natan M, Perez-Ladaga A, Zaneveld J. et al. TET2 mutations predict response to hypomethylating agents in myelodysplastic syndrome patients. Blood. 2014;124:2705-12

39. Itzykson R, Kosmider O, Cluzeau T, Mansat-De Mas V, Dreyfus F, Beyne-Rauzy O. et al. Impact of TET2 mutations on response rate to azacitidine in myelodysplastic syndromes and low blast count acute myeloid leukemias. Leukemia. 2011;25:1147-52

40. Traina F, Visconte V, Elson P, Tabarroki A, Jankowska AM, Hasrouni E. et al. Impact of molecular mutations on treatment response to DNMT inhibitors in myelodysplasia and related neoplasms. Leukemia. 2014;28:78-87

Author contact

Corresponding address Corresponding authors: Heng Mei, MD, Ph. D, Institute of Hematology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, No.1277 Jiefang Avenue, Wuhan 430022, China. Phone: 86-027-85726387, Fax: +862785726387, E-mail: hmeiedu.cn. Yu Hu, PhD, Institute of Hematology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430022, China; Collaborative Innovation Center of Hematology, Huazhong University of Science and Technology, Wuhan, Hubei 430022, PR China Email: dr_huyucom.


Received 2018-10-2
Accepted 2019-9-14
Published 2020-1-1