Automated analysis of acute myeloid leukemia minimal residual disease using a support vector machine
Metrics: PDF 1285 views | HTML 1324 views | ?
Wanmao Ni1, Beili Hu2, Cuiping Zheng3, Yin Tong4, Lei Wang1, Qing-qing Li1, Xiangmin Tong1, Yong Han1
1Key Laboratory of Tumor Molecular Diagnosis and Individualized Medicine of Zhejiang Province, Zhejiang Provincial People’s Hospital, Hangzhou, Zhejiang, P. R. China
2Medical of College, Zhejiang University, Hangzhou, Zhejiang, P. R. China
3The Center Hospital of Wenzhou, Wenzhou, Zhejiang, P. R. China
4Department of Hematology, Shanghai General Hospital, Shanghai, P. R. China
Xiangmin Tong, email: firstname.lastname@example.org
Yong Han, email: email@example.com
Keywords: flow cytometry, immunophenotyping, support vector machine, acute myeloid leukemia, minimal residual disease
Received: July 04, 2016 Accepted: September 29, 2016 Published: October 04, 2016
We investigated the ability of support vector machines (SVM) to analyze minimal residual disease (MRD) in flow cytometry data from patients with acute myeloid leukemia (AML) automatically, objectively and standardly. The initial disease data and MRD review data in the form of 159 flow cytometry standard 3.0 files from 36 CD7-positive AML patients in whom MRD was detected more than once were exported. SVM was used for training with setting the initial disease data to 1 as the flag and setting 15 healthy persons to set 0 as the flag. Based on the two training groups, parameters were optimized, and a predictive model was built to analyze MRD data from each patient. The automated analysis results from the SVM model were compared to those obtained through conventional analysis to determine reliability. Automated analysis results based on the model did not differ from and were correlated with results obtained through conventional analysis (correlation coefficient c = 0.986, P > 0.05). Thus the SVM model could potentially be used to analyze flow cytometry-based AML MRD data automatically, objectively, and in a standardized manner.
Acute myeloid leukemia (AML) is the most common form of leukemia in adults. Moreover, it has the lowest survival rate of any type of leukemia . During AML treatment, detection of minimal residual disease (MRD) using real-time quantitative polymerase chain reaction (RQ-PCR) and flow cytometry (FCM) provide a powerful basis for adjusting the diagnosis and treatment [2–8]. The presence of MRD is also strongly associated with risk of relapse and prognosis [9, 10]. RQ-PCR-based MRD detection is highly sensitive, but it depends on AML patients expressing specific fusion genes (e.g. AML1/ETO) or mutant genes (e.g. c-kit mutant), or overexpressing certain genes (e.g. Wilm’s oncogene) [11–14]. Consequently, detection is limited. FCM can detect leukemia associated immunophenotyping (LAIP) in 70%–75% of AML patients’ at an initial stage. The sensitivity of FCM is about 10−4 for MRD in AML whose immunophenotyping is obviously different from the normal cells, whereas the sensitivity is only 10−3 when immunophenotyping is partially overlapped with normal cells [15, 16].
Conventional FCM-based MRD detection has several disadvantages. (1) It requires LAIP, so more antibodies are consumed, though the detection schemes can be standardized . (2) Because it requires the analysts to have enough experience with the antigen differentiation regularity to rule out reactive hyperplasia of normal cells, different analysts may draw different conclusions from the same data, making the results subjective. (3) It mainly uses a scatter diagram in two-dimensional space, so multidimensional FCM data are not completely utilized, and multidimensional space recognition cannot be carried out.
Support vector machines (SVMs) are supervised learning algorithms used with neural networks, and have been widely used for analysis of DNA microarray data [17, 18] and facial recognition [19, 20], among other applications. An SVM can simultaneously learn all of the features derived from definitively classified training data in multidimensional space, build a recognition model for these data, and then recognize and classify unknown data using this model. With multidimensional recognition, the SVM model has greater accuracy and precision than one- or two-dimensional recognition. In recent years, SVMs have been gradually applied to FCM data analysis [21–27] – e.g. it can distinguish malignant lymphocytes from benign lymphocytes. Toedling et al. employed an SVM to detect residual disease in acute B lymphocytic leukemia with 99.78% sensitivity and 98.87% specificity . The small and rapid LIBSVM library developed by Chang et al. is a mature and perfect class library for all SVM algorithms .
In the present study, the LIBSVM was introduced into the data analysis to overcome the disadvantages of FCM-based MRD detection. By combining SVM multidimensional training and recognition characteristics with the multidimensional advantages of FCM, MRD can be quantified more precisely, objectively and economically.
In this study, 159 AML data sets were selected from 36 patients, including 22 men and 14 women with an average age 42 years (range, 17 to 72 years). The distribution of AML subtypes included 10 M0, 3 M1, 15 M2, 6 M5, and 2 M6. The M3, M4, M7 and other subtypes seldom expressed CD7, so few cases were selected.
The effect of training cell number
Combing the data from the patients’ initial immunophenotyping with the data from the 15 healthy individuals generated a training data file with more than 2 × 105 events. However, because with LIBSVM the length of the training time was proportional to size of the training data file, using the complete training data set it took 24–36 h to finish the parameter optimization and model building with low efficiency. Therefore, to optimize parameters and model building more quickly, 104 events were selected from the training data through stratified random sampling, which only took 5–10 min. In the patient group, the residual leukemic cell fractions in 10 AML MRD data sets were analyzed using the SVM. Notably, sets containing the 104 sampled events did not significantly differ from the complete data sets (Paired t-test P = 0.0792) (Table 1 and Figure 1). Therefore, to improve the training efficiency, 104 events selected through stratified random sampling were applied for training and model building.
Table 1: Comparison between different event numbers for training and calculating AML MRD
2 × 105 events
Figure 1: Comparison between different event numbers for training and calculating AML MRD. (104 events vs. 2 × 105 events, P > 0.05)
Optimization of training parameters C and γ
The grid.py script in LIBSVM was used to optimize parameters C and γ for each patient, based on the 104 sampled event groups. Through optimization training for the 36 patients, 36 groups of C = 0.50–32768 (median, 8.00), γ = 0.13–8.00 (median, 8.00), and a corresponding optimized cross-validation accuracy of 98.25%–99.95% (median, 99.46%) were obtained, which meant that the accuracy of the individual-specific predictive models for the 36 patients could reach 98.25%–d99.95% when they were used in the MRD SVM analysis.
Correlation between SVM group and manual group
The 36 predictive models were used in the automated SVM analysis of MRD data, and the MRD data was also analyzed conventionally, yielding 159 groups of paired data. The MRD cell fraction was determined be 0.006%–82.180% using conventional analysis and 0.006%–77.200% using the automated analysis. The correlation coefficient was 0.986, and a two-tailed paired t-tests showed the correlation to be significant (P < 0.05) and without a significant statistical difference between the results (P = 0.134). Using the Bland-Altman comparison method, only 11 of the 159 data pairs were out of 95% limits of agreement (Figure 2). From the scatter diagrams, it was apparent that the distribution of leukemia cells determined using SVM analysis was similar to that obtained using manual analysis (Figure 3).
Figure 2: Bland-Altman comparison of SVM and manual analysis results. Of the 159 pairs of data, only 11 were outside the 95% limits of agreement, which was from −4.4 to 3.9.
Figure 3: Comparison of the automatic SVM and manual analyses of typical AML patient results. For clarity, each scatter diagram shows 104 events. (A) The leukemic cell fraction was 24.672%, according to the SVM predictive model building of this MRD. The leukemia cell events are in red and the normal cells are in blue. (B) According to the manual analysis, the leukemia cell fraction was 24.466%, based on initial immunophenotyping of the patients. The gate was set by each step, and the MRD ratio was calculated as “P1 and P2 and P3”. The scatter diagrams were CD7/CD117, SSC/CD45, and CD117/HLA-DR, from top to bottom.
An SVM is a supervised learning algorithm that can learn the characteristics of known objects in multiple dimensions then build predictive models with which to classify data of unknown classification . LIBSVM is an excellent, easy and mature library. For FCM data, SVM has the advantage of being capable of multidimensional analysis, especially for 4 or more colors, and avoids the artificial misjudgment and experience requirements of an analyst. In this study, when LIBSVM was applied for automated MRD analysis, the results did not significantly differ from those obtained conventionally, and the results of the analysis could be displayed in different colors on scatter diagrams .
While the SVM is learning the known classification data and model building, parameter optimization can affect the ability of the model to accurately analyze unclassified data. C and γ are both important parameters for optimization . C is the penalty coefficient, which controls the model’s ability to generalize. If C is too large or too small, the ability of values to float will be poor . Parameter γ controls the degrees of freedom in the nonlinear model – i.e. the number of support vectors. Only when C and γ are optimum does the model have the highest prediction accuracy (CV rate). Each patient expressed different levels of CD7, CD117, CD45, and HLA-DR, so the distributions differed in the multidimensional space. Consequently, when leukemia cells were mixed with normal cells to form the training data file, the dividing plane between the leukemic and normal cells differed in the multidimensional space. Thus, finding the corresponding optimal dividing plane for the patient, which would identify the optimal parameters of the individual-specific C and γ, was the key to establishing the SVM model . The grid.py script in the LIBSVM class library was therefore used to optimize the parameters. However, all the data required cross-validation 5 times, and C and γ was tested step by step, so the optimization time and data volume were closely related. We found that if parameter optimization and model building were done using the complete training data (derived from patients with initial leukemia cells and 15 cases of normal data, including more than 20 million events), it would take 24–36 h to complete, which is not suitable for practical application. We therefore used a stratified random method to extract 104 data as training data. Comparison with the automatic analysis results obtained with MRD models using two different size data sets revealed no significant difference between the models (P > 0.05).
When the optimized parameters were used to establish a predictive model, the results of the SVM automated MRD analysis did not differ from those obtained using the manual method (P > 0.05), with a correlation coefficient of 0.986. By comparing the results obtained using the SVM with those obtained manually from scatter diagrams, we found that the SVM properly studied the distribution characteristics of the patient’s initial leukemia cells in multiple dimensions and accurately identified the residual leukemia cells.
In conclusion, automated SVM and manual analyses showed good consistency, which could reduce the experience requirements for MRD analysis. However, the sample size of AML patients involved was not large enough, and additional data from CD7-positive AML patients is needed. Moreover, AML patients exhibiting other expression patterns (e.g. CD19 expression without CD7 expression or abnormal CD33 and CD13 expression patterns without CD7 expression) should also be analyzed using an SVM in order to determine the scope of its utility. In short, the introduction of LIBSVM and other SVM algorithm libraries can be completely applied to the multidimensional features of FCM data, which would make interpretation of FCM data more objective and would provide a basis for achieving automated FCM data analysis.
MATERIALS AND METHODS
Immunophenotypes vary among AML patients, and there is no consensus on the LAIP for MRD detection. We studied one specific LAIP containing CD7, CD117, CD45 and HLA-DR. Information with the features listed below for patients treated from 2010 to 2012 were selected from the FCM database in our division. (1) Initial immunophenotyping schemes contained CD7, CD117, CD45 and HLA-DR, and these four antigens were detected in the same tube. (2) Leukemia residual disease lesion detection was carried out no less than twice. (3) The specimens were bone marrow. (4) The cell counts in the data file were at least 105. Using these criteria, 159 data files from 36 patients were selected. An individual-specific predictive model was built for each patient.
Reagents and instruments
CD7 (Becton Dickinson Biosciences, USA) was labeled with fluorescein isothiocyanate (FITC) using clone No. 4H9. CD117 (Beckman Coulter, Inc. USA) was labeled with phycoerythrin (PE) using clone No. 95C3. HLA-DR (Becton Dickinson Biosciences, USA) was labeled with allophycocyamin (APC) using clone No. G46-6. CD45 (Becton Dickinson Biosciences, USA) was labeled with phycoerythrin cyanin 7 (PE-cy7) using clone No. J33.
A BD FACS Calibur equipped with two lasers and four colors was applied for the MRD study. LIBSVM version 3.16 was used for model building, parameter optimization and data prediction for automated analysis. ACEA NovoExpressTM, which was developed by ACEA Bioscience Inc. to have a humanization design interface and functions, was used for manual analysis.
Automated MRD analysis using LIBSVM
The individual model for each patient was built through the processes of training data file derivation, parameter optimization, and individual model building.
ACEA NovoExpressTM was used to open the flow cytometry standard (FCS) data files for patients’ initial immunophenotyping. The leukemia cells at diagnosis were selected and set as P1, after which P1 were exported in the comma separated value (CSV) file format. Thereafter, CSV files were imported into MATLAB R2011a and saved as matrices in which the columns were fluorescence, rows were cells, and the matrix values were fluorescence intensities. The matrix data were written into training data files named [patientname.train.txt] in the form required for LIBSVM, and the data flag was set as 1. FCS data from 15 healthy individuals were read in the same way, and the data were added to the aforementioned [patientname.train.txt] file. That data flag was set as 0. To normalize the data and speed up model building, ‘svm-scale.exe’ was used to adjust the data in the [patientname.train.txt] file following the LIBSVM requirements so as to construct the [patientname.train.scale] file. At the same time, the [patientname.train.range] file was formed to enable adjustment of the data in the prediction stage.
After the data were ready, the LIBSVM default kernel (radial basis function) was employed to build individual-specific models. In addition, using the penalty coefficient C and core parameter γ, model prediction was optimized to form the optimized model. For this purpose, grid.py in LIBSVM was called to carry out the cross-validation of the training data file formed above. The default parameter values were set through 5-time cross-validation with the command ‘grid.py –svmtrain “svmtrain.exe” -log2c -gnuplot “gnuplot.exe” “[patientname.train.scale]”’. The optimized C and γ were obtained to build the patients’ individual-specific training model with the command ‘svm-train.exe–c–g “[patientname.train.scale]” “[patientname.train.model]”’.
For the automated MRD analysis, the original data for a patient’s MRD were processed in the standard LIBSVM form and a data flag was set to 1, as mentioned above. To recognize MRD data with different detection times, the file was named [patientname+detectiontime.test], and the aforementioned [patientname.train.range] file was used to scale [patientname+detectiontime.test] into [patientname+detectiontime.test.scale]. The svm-predict.exe program in the LIBSVM software package was then combined with the individual patient’s specific training model to automatically analyze the [patientname+detectiontime.test.scale] file and predict the accuracy of flag 1; that is, the accuracy of the patient’s residual leukemic cell fraction. Finally, the predicted leukemic cells (in red) and normal cells (in blue) were plotted using the “scatter” function in MATLAB.
Conventional manual analysis of MRD
The three scatter diagrams CD7/CD117, SSC/CD45 and FSC/HLA-DR were plotted for patients’ initial immunophenotyping. The P1, P2 and P3 gates were set for leukemia cells in the scatter diagrams. Finally, the leukemic cell fraction was obtained through the logical combination “P1 AND P2 AND P3”.
The data groups in the automated SVM analysis were compared to the corresponding data in the conventional manual analysis, one by one, using two-tailed paired t tests. Values of P < 0.05 were considered significant.
Funding: This work was supported by the National Science Foundation of China (81570198), the National Science Foundation of Zhejiang Province (LY14H160001, 2014C33159), and the TCM Foundation of Zhejiang Province (2016ZZ007, 2016ZB018).
CONFLICTS OF INTEREST
Wanmao Ni, Yin Tong and Yong Han performed the research, Beili Hu, Xiangmin Tong, Yong Han designed the research study, Wanmao Ni, Xiangmin Tong, Yin Tong and Lei Wang analyzed the data, Wanmao Ni, Beili Hu, Cuiping Zheng, Yin Tong, Lei Wang, Qingqing Li, Xiangmin Tong, and Yong Han wrote and revised the paper.
1. Deschler B, Lubbert M. Acute myeloid leukemia: epidemiology and etiology. Cancer. 2006; 107:2099–2107.
2. Campana D. Role of minimal residual disease evaluation in leukemia therapy. Curr Hematol Malig Rep. 2008; 3:155–160.
3. Rossi G, Minervini MM, Carella AM, de Waure C, di Nardo F, Melillo L, D’Arena G, Zini G, Cascavilla N. Comparison between multiparameter flow cytometry and WT1-RNA quantification in monitoring minimal residual disease in acute myeloid leukemia without specific molecular targets. Leuk Res. 2012; 36:401–406.
4. DiNardo CD, Luger SM. Beyond morphology: minimal residual disease detection in acute myeloid leukemia. Curr Opin Hematol. 2012; 19:82–88.
5. Kern W, Haferlach C, Haferlach T, Schnittger S. Monitoring of minimal residual disease in acute myeloid leukemia. Cancer. 2008; 112:4–16.
6. Freeman SD, Jovanovic JV, Grimwade D. Development of minimal residual disease-directed therapy in acute myeloid leukemia. Semin Oncol. 2008; 35:388–400.
7. van der Velden VH, Boeckx N, van Wering ER, van Dongen JJ. Detection of minimal residual disease in acute leukemia. J Biol Regul Homeost Agents. 2004; 18:146–154.
8. Hokland P, Ommen HB, Nyvold CG, Roug AS. Sensitivity of minimal residual disease in acute myeloid leukaemia in first remission--methodologies in relation to their clinical situation. Br J Haematol. 2012; 158:569–580.
9. Xu XJ, Feng JH, Tang YM, Shen HQ, Song H, Yang SL, Shi SW, Xu WQ. Prognostic significance of flow cytometric minimal residual disease assessment after the first induction course in Chinese childhood acute myeloid leukemia. Leuk Res. 2013; 37:134–138.
10. Loken MR, Alonzo TA, Pardo L, Gerbing RB, Raimondi SC, Hirsch BA, Ho PA, Franklin J, Cooper TM, Gamis AS, Meshinchi S. Residual disease detected by multidimensional flow cytometry signifies high relapse risk in patients with de novo acute myeloid leukemia: a report from Children’s Oncology Group. Blood. 2012; 120:1581–1588.
11. Moon HW, Shin S, Kim HY, Kim YR, Cho HI, Yoon SS, Park S, Kim BK, Chun H, Kim HC, Park CJ, Min YH, Lee DS. Therapeutic use of granulocyte-colony stimulating factor could conceal residual malignant cells in patients with AML1/ETO+ acute myelogenous leukemia. Leukemia. 2006; 20:1408–1413.
12. Bergmann L, Miething C, Maurer U, Brieger J, Karakas T, Weidmann E, Hoelzer D. High levels of Wilms’ tumor gene (wt1) mRNA in acute myeloid leukemias are associated with a worse long-term outcome. Blood. 1997; 90:1217–1225.
13. Inaba H, Coustan-Smith E, Cao X, Pounds SB, Shurtleff SA, Wang KY, Raimondi SC, Onciu M, Jacobsen J, Ribeiro RC, Dahl GV, Bowman WP, Taub JW, et al. Comparative analysis of different approaches to measure treatment response in acute myeloid leukemia. J Clin Oncol. 2012; 30:3625–3632.
14. Kotliar N, Koziner B. [Use of polymerase chain reaction (PCR) in oncohematology. Detection of minimal residual disease]. Medicina (B Aires). 1995; 55:159–166.
15. Campana D, Coustan-Smith E. Detection of minimal residual disease in acute leukemia by flow cytometry. Cytometry. 1999; 38:139–152.
16. Campana D, Coustan-Smith E. Minimal residual disease studies by flow cytometry in acute leukemia. Acta Haematol. 2004; 112:8–15.
17. Furey TS, Cristianini N, Duffy N, Bednarski DW, Schummer M, Haussler D. Support vector machine classification and validation of cancer tissue samples using microarray expression data. Bioinformatics. 2000; 16:906–914.
18. Brown MP, Grundy WN, Lin D, Cristianini N, Sugnet CW, Furey TS, Ares M, Haussler D. Knowledge-based analysis of microarray gene expression data by using support vector machines. Proc Natl Acad Sci USA. 2000; 97:262–267.
19. Osuna E, Freund R, Girosit F. Training support vector machines: an application to face detection. Computer Vision and Pattern Recognition, 1997 Proceedings, 1997 IEEE Computer Society Conference on: IEEE), (1997) pp. 130–136.
20. Guo G, Li SZ, Chan K. . Face recognition by support vector machines. Automatic Face and Gesture Recognition, 2000 Proceedings Fourth IEEE International Conference on: IEEE), (2000) pp. 196–201.
21. Bashashati A, Lo K, Gottardo R, Gascoyne RD, Weng A, Brinkman R. A pipeline for automated analysis of flow cytometry data: preliminary results on lymphoma sub-type diagnosis. Conf Proc IEEE Eng Med Biol Soc. 2009; 2009:4945–4948.
22. Baudry JP, Raftery AE, Celeux G, Lo K, Gottardo R. Combining Mixture Components for Clustering. J Comput Graph Stat. 2010; 9:332–353.
23. Boedigheimer MJ, Ferbas J. Mixture modeling approach to flow cytometry data. Cytometry A. 2008; 73:421–429.
24. Finak G, Bashashati A, Brinkman R, Gottardo R. Merging mixture components for cell population identification in flow cytometry. Adv Bioinformatics. 2009:247646.
25. Lo K, Brinkman RR, Gottardo R. Automated gating of flow cytometry data via robust model-based clustering. Cytometry A. 2008; 73:321–332.
26. Rajwa B, Venkatapathi M, Ragheb K, Banada PP, Hirleman ED, Lary T, Robinson JP. Automated classification of bacterial particles in flow by multiangle scatter measurement and support vector machine classifier. Cytometry A. 2008; 73:369–379.
27. Fiser K, Sieger T, Schumich A, Wood B, Irving J, Mejstrikova E, Dworzak MN. Detection and monitoring of normal and leukemic cell populations with hierarchical clustering of flow cytometry data. Cytometry A. 2012; 81:25–34.
28. Toedling J, Rhein P, Ratei R, Karawajew L, Spang R. Automated in-silico detection of cell populations in flow cytometry readouts and its application to leukemia disease monitoring. BMC Bioinformatics. 2006; 7:282.
29. Chang C-C, Lin C-J. LIBSVM: A library for support vector machines. ACM Trans Intell Syst Technol. 2011; 2:1–27.
30. Statnikov A, Wang L, Aliferis CF. A comprehensive comparison of random forests and support vector machines for microarray-based cancer classification. BMC Bioinformatics. 2008; 9:319.
31. Dong B, Cao C, Lee SE. Applying support vector machines to predict building energy consumption in tropical region. Energy Build. 2005; 37:545–553.
32. Huang C-L, Liao H-C, Chen M-C. Prediction model building and feature selection with support vector machines in breast cancer diagnosis. Expert Syst Appl. 2008; 34:578–587.
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.