- 1Department of Hematology, Shidong Hospital of Shanghai Yangpu District, Shanghai, China
- 2School of Computer Science and Technology, Tongji University, Shanghai, China
- 3Laboratory Diagnosis Department, Shanghai Kingmed Center for Clinical Laboratory Co., Ltd., Shanghai, China
Introduction: Acute myeloid leukemia (AML) is a hematological malignancy that requires accurate diagnosis and continuous monitoring to guide effective treatment. Flow cytometry is widely used because it enables the detection of minimal residual disease. However, current methods often rely on uniform marker panels, overlooking the heterogeneity that arises when different markers or staining protocols are used across patients. In addition, remission states are frequently neglected, despite their clinical importance for disease management and prognosis.
Methods: To address these challenges, we developed a machine learning–based classification framework that integrates heterogeneous flow cytometry data. A dataset comprising 53 markers was collected, and six different machine learning classifiers were trained to distinguish between AML, complete remission (AML-CR), and normal samples. Model performance was evaluated using accuracy, precision, recall, F1 score, and area under the ROC curve (AUC).
Results: Among the classifiers evaluated, the Random Forest model demonstrated the highest performance, achieving an accuracy of 94.92%, an F1-score of 94.13%, a precision of 94.58%, a recall of 93.74%, and an AUC of 94.83%. These results indicate that machine learning can effectively classify AML and remission states from heterogeneous flow cytometry data.
Discussion: This study highlights the value of machine learning in overcoming limitations of traditional flow cytometry analysis. By accommodating marker heterogeneity and incorporating remission states, the proposed framework provides a more robust and clinically relevant tool for AML diagnosis and monitoring. The findings suggest that machine learning models, particularly Random Forest, hold strong potential for improving precision in hematological diagnostics. The code for this study is publicly available at https://zenodo.org/records/15110287.
1 Introduction
Acute myeloid leukemia (AML) is a malignant clonal disease originating from the abnormal proliferation and differentiation of hematopoietic stem cells. It represents the most prevalent form of adult leukemia (1). Significant improvements in AML prognosis have been achieved through advancements in chemotherapy, targeted therapy, transplantation techniques, CAR-T therapy, and the ongoing refinement of supportive care (2–4). Flow cytometry, which utilizes specific antibodies to label surface antigens on leukemia cells, is capable of identifying and quantifying as few as 0.01% leukemia cells. It has emerged as a critical tool for the diagnosis and monitoring of AML, widely applied in the analysis of diverse cell populations (5, 6) and the assessment of minimal residual disease (MRD) to evaluate disease prognosis (7, 8). However, flow cytometry relies heavily on manual operation, which is associated with significant drawbacks, including time-consuming procedures, high subjectivity, and inconsistent results, potentially leading to missed diagnoses or misdiagnoses (9). Therefore, investigating the potential for automated diagnosis based solely on flow cytometry data and developing an intelligent diagnostic system that is automated, precise, and broadly applicable holds significant clinical value for early diagnosis and treatment.
In recent years, machine learning (ML) has achieved remarkable advancements in intelligent medical diagnosis, particularly in disease classification, prediction, and personalized treatment (10–12). ML excels at automatically learning from large-scale datasets and uncovering underlying patterns (13), offering unparalleled advantages over traditional methods, especially when processing complex and high-dimensional biomedical data such as flow cytometry data (14). However, existing models may not be directly applicable to real-world flow cytometry diagnostic scenarios due to the unstandardized nature of flow cytometry data, which often fails to meet the input requirements of these models. Specifically, in practical settings, the performance limitations of flow cytometry instruments necessitate the use of multiple panels to obtain comprehensive data for a single patient. For instance, the same marker may be labeled with different fluorescent dyes, leading to variations in the measured values. Such discrepancies are rarely encountered in publicly available standardized datasets, where each sample typically employs a consistent combination of markers and dyes (15). Furthermore, existing studies on flow cytometry datasets frequently overlook the analysis of patients with complete acute myeloid leukemia remission (AML-CR), which oversight limits the comprehensive understanding of patient data distribution, impairs the evaluation of treatment efficacy, and hinders the monitoring of disease relapse risk.
To address this issue, we collected a dataset from real-world diagnostic scenarios, encompassing samples with diverse combinations of antibodies and dyes. A key advantage of this dataset is its inclusion of AML-CR samples, enabling us to investigate variations in cellular populations across different disease stages. Subsequently, we calculate the statistical properties of each marker to standardize the samples into a consistent format. This step preserves the distribution information of markers while ensuring compatibility with the input requirements of ML models. Finally, ML algorithms are employed to automatically extract feature information from the standardized flow cytometry data and construct robust classification models.
Specifically, this study collected flow cytometry data from 59 AML patients, 34 AML-CR patients, and 101 bone marrow flow cytometry-normal (Norm) patients, encompassing the expression profiles of various cell surface markers. Utilizing multiple ML algorithms for feature extraction, selection, and modeling, we developed a diagnostic model capable of distinguishing among the three patient groups. The model underwent rigorous feasibility analysis, performance validation, and comprehensive evaluation. Extensive experimental results demonstrate the efficacy of the proposed method in AML diagnosis. Furthermore, we conducted additional analysis to assess the importance of markers within the model. On one hand, the model’s findings align with clinical knowledge, mutually reinforcing each other. On the other hand, the identification of potentially significant features may offer novel insights into disease mechanisms. This approach not only enhances our understanding of the immunological characteristics of AML but also equips clinicians with more scientific and efficient diagnostic tools.
2 Methods
2.1 Study population
This retrospective study analyzed flow cytometry data from patients treated at Shidong Hospital, Yangpu District, Shanghai, between January 2019 and October 2024. The study included samples from patients with acute myeloid leukemia (non-M3 type, AML), AML in complete remission (AML-CR), and those with normal bone marrow flow cytometry results (Norm). The normal group comprised patients with cytopenia or cytosis caused by non-neoplastic conditions, including nutritional anemia, immune thrombocytopenia, and primary thrombocythemia. Inclusion criteria: The study subjects are AML patients aged between 18 and 70 years. The diagnosis and classification of leukemia are based on the World Health Organization 5th Edition Classification of Hematologic and Lymphoid Tumors (16), with comprehensive evaluation considering clinical manifestations, morphology, cytogenetics, and molecular results (17). Complete remission of bone marrow after treatment is assessed according to the 4th Edition of Diagnostic and Efficacy Criteria for Hematologic Diseases, with a blast cell percentage of <5% defined as complete remission. Exclusion criteria: Patients with other hematologic disorders, severe infections, or other systemic diseases that may affect flow cytometry results, as well as cases with incomplete or obviously abnormal data.
Based on the above criteria, a total of 59 AML samples, 34 AML-CR samples, and 101 Norm samples are included in the study. The dataset comprised 53 distinct markers, as detailed in Table 1, which are utilized for subsequent feature engineering to extract and optimize classification features. The data are randomly split into training and testing sets at a ratio of 7:3 while maintaining the proportional distribution of each category in both sets. The complete workflow for sample screening and data processing is illustrated in Figure 1. The automatic diagnostic process of the flow cytometry-based model begins with data collection, where patient samples are acquired and raw flow cytometry data are generated. These data then undergo scanning, involving preprocessing and quality control to ensure accuracy and consistency. Next, during feature aggregation, relevant cellular features are extracted and combined to form comprehensive representations of each sample. The aggregated features are used in the training and inference phase, where the machine learning model is trained on labeled datasets and subsequently applied to new patient data for prediction. Finally, the model produces diagnostic results that assist clinicians in making decisions.

Figure 1. Automated diagnostic workflow for flow cytometry data analysis. For each patient, multiple LMD files are processed to calculate statistical features, which are then aggregated into a standardized format. The formatted data are subsequently input into a machine learning (ML) model for training and classification. LR, Logistic Regression; RF, Random Forest; MLP, Multilayer Perceptron; SVM, Support Vector Machine; XGB, Extreme Gradient Boosting; GBM, Gradient Boosting Machine.
2.2 Data selection
Fresh bone marrow samples (3 − 5mL), anticoagulated with heparin or EDTA, are collected and thoroughly mixed before storage at room temperature. The leukocyte count is determined using an automated hematology analyzer. Based on the leukocyte count, the sample is either diluted with PBS or concentrated by centrifugation at 1700 rpm to adjust the leukocyte concentration to 1 × 107/mL. Tubes are prepared according to the sample and the specific antibody panel. For membrane staining, a pre-prepared antibody cocktail is added to each tube based on the selected antibody combination. The sample is thoroughly mixed (at least 5 inversions), and the calculated volume of the diluted or concentrated sample is added to the bottom of the tube. After gentle mixing, the sample is incubated in the dark for 15–20 minutes to ensure optimal staining efficiency. Subsequently, red blood cell lysis buffer is added, followed by an additional 10-minute incubation in the dark until complete lysis is achieved. Centrifuge the sample at 1700 rpm for 5 minutes, then discard the supernatant. The pellet is resuspended in 2 mL of PBS, mixed thoroughly, and centrifuged again at 1700 rpm for 5 minutes. After discarding the supernatant, 600 µL of 1% paraformaldehyde fixation solution is added to resuspend the cells, which are then subjected to flow cytometry analysis. For intracellular staining, the above steps are followed according to the reagent manufacturer’s instructions. For detection of surface or cytoplasmic immunoglobulin light chains, the sample is washed three times with PBS before antibody addition. Data acquisition is performed on a Navios 10 COLORS/3 LASER flow cytometer, ensuring at least 5 × 105 events are collected per sample. The antibody panel included surface markers such as CD34, CD38, CD45, and CD117. All data are stored in LMD file format.
2.3 Data pre-processing and feature engineering
Each sample corresponds to tens of thousands of cells, and each cell carries multiple marker results. In clinical diagnosis, physicians often focus on the distribution patterns of these markers across different cell populations. Simply averaging the marker values for all cells may overlook important distribution characteristics. To address this, we incorporated additional statistical measures, including mean, standard deviation (std), median, skewness, and kurtosis, to capture the variability and asymmetry of marker distributions more comprehensively. These statistical measures effectively capture the distributional differences of markers across cell populations. The mean reflects the overall expression level of a marker, while the standard deviation quantifies variability among cells. The median reduces the influence of extreme values, skewness reveals distribution asymmetry, and kurtosis indicates the sharpness of the distribution or the presence of outliers. The formulas for calculating skewness and kurtosis are shown in Equations 1, 2, respectively.
where is the th data point, is the sample mean, is the sample standard deviation, and is the total number of samples. Therefore, the 53 marker values of all cells from each patient are transformed into a 265-dimensional feature vector, where each marker is represented by the five aforementioned distributional features. This transformation enables a more comprehensive representation of marker distribution across cells. The numerical distribution of the transformed data in the training and test sets is illustrated in Figure 2.

Figure 2. A detailed visualization of the dataset is presented. The left panel illustrates the distribution of markers across different samples, while white indicates the presence of corresponding markers and black indicates their absence. The right panel displays the distribution of features extracted from different samples. Here, Markers refers to the 53 cell surface markers analyzed, while Features represents the 265 statistical features derived from the 53 markers.
2.4 Model establishment and evaluation
The flow cytometry data is randomly split into training and testing sets in a 7:3 ratio. To ensure a comprehensive evaluation, this study employs six widely used ML algorithms, encompassing various classical approaches. These include linear models [Logistic Regression, LR (18)], ensemble methods [Random Forest, RF (19); Extreme Gradient Boosting, XGBoost (20); Gradient Boosting Machine, GBM (21)], neural networks [Multilayer Perceptron, MLP (22)], and support vector machines [Support Vector Machine, SVM (23)]. Each algorithm represents a distinct learning paradigm: linear models effectively capture linear relationships and are suitable when interpretability and simplicity are important, ensemble methods enhance predictive performance by aggregating multiple weak learners and work well for complex, nonlinear, and noisy data, neural networks excel at modeling complex nonlinear patterns and are particularly effective with large datasets and intricate feature interactions, support vector machines are well-suited for high-dimensional classification tasks, especially when the classes are separable with clear margins. All models are implemented using scikit-learn==1.6.1 in Python 3.10. In the training set, all flow cytometry data are divided into five parts, and five-fold cross-validation is applied. This strategy helps assess the stability and generalization ability of the model by cycling through each subset as the validation set. During cross-validation, grid search is used to optimize the hyperparameters. The specific parameter search space for each method is shown in Table 2. The evaluation metrics include accuracy, F1-score, precision, recall, and area under the curve (AUC). These metrics collectively reflect the model’s performance in the classification task, with particular significance for F1-score and AUC when handling imbalanced data. In addition, to analyze the contribution of features to the prediction, shap==0.46.0 (24) is used to estimate feature importance. Based on the concept of Shapley values, SHAP assigns an importance value to each feature, thereby helping to explain the model’s decision-making process and enhancing the model’s interpretability and reliability.
3 Results
3.1 Clinical characteristics of patients
We enrolled 59 patients with acute myeloid leukemia (AML, non-M3), 34 patients who achieved complete bone marrow remission after treatment (AML-CR), and 101 individuals with normal bone marrow flow cytometry results (Norm). The dataset is randomly divided into a training set (n=135) and a test set (n=59) in a 7:3 ratio. The collected variables encompassed demographic characteristics (age and sex), routine blood parameters (white blood cell count, hemoglobin level, and platelet count), the proportion of bone marrow blasts/immunized cells, and 53 commonly used markers from flow cytometry data. The baseline characteristics of the participants are summarized in Table 3. The mean age of patients in the training set is 66.59 years (66.59 ± 14.33), while that in the test set is 66.66 years (66.66 ± 13.81), with no statistically significant difference between the two groups (P > 0.05). Chi-square analysis revealed no significant differences (P > 0.05) in any of the examined variables between the training and test sets, indicating similar distribution patterns across all factors. These results demonstrate that both the training and test sets are well-balanced and appropriate for subsequent predictive analysis.
3.2 Model performance
Table 4 and Figure 3 present the performance of different models in the classification task. Overall, both RF and XGB achieved superior results across all metrics, with F1-score, precision, and recall reaching 0.9413, 0.9458, and 0.9374, respectively, and an accuracy of 0.9492. These results indicate strong generalization capabilities for these two models in the classification task. LR exhibited a high AUC value (0.9705); however, its F1-score (0.9175) and precision (0.9098) are slightly lower than those of RF and XGB, suggesting some degree of misclassification. Both MLP and SVM demonstrated identical classification performance, with an accuracy of 0.9153 and an F1 score of 0.8932. Notably, SVM achieved the highest AUC value (0.9741), although its other metrics are lower than those of RF and XGB. GBM performed slightly worse than the other models, with the lowest accuracy (0.8983) and F1 score (0.8806). Although its AUC value reached 0.9491, its precision (0.8752) and recall (0.8974) exhibited a certain gap, possibly due to the model’s limited ability to distinguish between specific classes. Overall, RF and XGB demonstrated robust performance across all metrics, making them the most suitable candidates for this task. Meanwhile, SVM, with the highest AUC value, may offer advantages in certain application scenarios. Additionally, confusion matrix analysis revealed that the misclassified samples are evenly distributed across different categories, indicating that these models maintained a balanced classification error across classes.

Figure 3. The performance analysis of different models. The first row shows the receiver operating characteristic (ROC) curve, the second row shows the precision-recall (PR) curve, and the third row shows the confusion matrix. LR, linear regression; RF, random forest; MLP, multi-layer perceptron; SVM, support vector machine; XGB, extreme gradient boosting; GBM, gradient boosting machine.
3.3 Feature importance analysis
Figure 4 displays the top 10 most important features in the decision-making process for each model. Notably, CD117 and CD34 are identified as the most influential markers for distinguishing AML in the GBM, RF, and XGB models. These findings resonate with clinical practice, where both CD117 and CD34 are extensively utilized for AML differentiation and classification. This congruence between model predictions and clinical practices underscores the potential of data-driven approaches in medical diagnostics and highlights the pivotal roles of CD117 and CD34 in AML diagnosis. HLA-DR significantly affects the classification of AML and AML-CR in the GBM, RF, and LR models, while CD45 plays a crucial role in the GBM and SVM models, suggesting its potential contribution to the regulation of immunophenotypic heterogeneity in AML. Lower expression of HLA-DR is commonly observed in immature leukemia cells, particularly in M0/M1 subtypes, indicating a differentiation block that may facilitate immune evasion or serve as a marker for specific stages of differentiation. Furthermore, CD45 expression varies throughout the stages of myeloid differentiation, with its heterogeneity, in conjunction with differential expression of HLA-DR, contributing to the formation of different immune subtypes, thus indicating a differentiation blockade and the coexistence of multiple stages of differentiation. Additionally, the markers in the Norm group exhibit lower specificity, raising concerns about the risk of overfitting in some models due to the smaller phenotypic variability observed in normal samples. These insights suggest that while the models generally perform well, caution is warranted regarding the risk of overfitting, especially when distinguishing between normal and diseased states.

Figure 4. The feature importance analysis of different models based on SHAP. Each row represents the SHAP feature importance analysis results of a method on acute myeloid leukemia (AML), acute myeloid leukemia complete remission (AML-CR), and normal (Norm) classes, respectively. LR, linear regression; RF, random forest; MLP, multi-layer perceptron; SVM, support vector machine; XGB, extreme gradient boosting; GBM, gradient boosting machine.
4 Discussion
Acute myeloid leukemia (AML) is a highly aggressive and heterogeneous hematopoietic malignancy, representing 15-20% of all leukemia diagnoses worldwide. The global annual incidence is estimated at 1.0–1.5 per 100,000 population, corresponding to approximately 40,000–50,000 new cases each year, with around 80% occurring in adults (aged ≥ 18 years) (25). Despite significant advances in therapeutic modalities, including chemotherapy, allogeneic stem cell transplantation, and novel targeted/immunotherapeutic agents, the overall prognosis of AML remains unsatisfactory. Patients with high-risk features (e.g., advanced age [≥60 years], TP53 mutations, or complex karyotypes) exhibit particularly dismal outcomes, with 5-year overall survival rates as low as 5-15% (4). The therapeutic landscape remains particularly constrained for patients with relapsed/refractory (R/R) AML. Disease progression is frequently complicated by life-threatening cytopenias (including transfusion-dependent anemia), opportunistic infections, and progressive multiorgan dysfunction - clinical manifestations directly attributable to the intrinsically aggressive biology of leukemic cells. This underscores the critical unmet need for advanced diagnostic modalities to guide precision therapeutic strategies. Flow cytometry-based immunophenotyping analysis offers critical support for the early detection of disease progression and the development of personalized treatment strategies by dynamically monitoring cell characteristics.
Conventional flow cytometry methods for manual detection typically rely on gating and clustering techniques to categorize cells into multiple subpopulations, performing multi-parameter analysis of surface and intracellular markers. These methods are widely used for the diagnosis and classification of hematologic malignancies (26). Conventional flow cytometry data analysis remains labor-intensive and subject to inter-operator variability due to its reliance on manual gating expertise. This inherent limitation has spurred the rapid adoption of machine learning algorithms in hematological diagnostics, enabling automated analysis of cellular morphology, immunophenotypic patterns, and histopathological features with enhanced reproducibility (27–29). Among them, Beni et al. (30) introduces a multi-cell classification benchmark dataset; Hu et al. (31) uses deep convolutional neural networks for cytomegalovirus classification based on flow cytometry data; Li et al. (32) transforms SW-480 epithelial cancer cell flow cytometry data into images and uses convolutional neural networks for classification. Although these methods have positively contributed to the improvement of flow cytometry diagnostics, their practical implementation continues to encounter substantial challenges: (1) the inherently high-dimensional nature of patient-level FCM data; (2) substantial inter-sample variability introduced by both biological heterogeneity (e.g., treatment response status) and technical factors (e.g., instrument configuration); and (3) the critical knowledge gap regarding immunophenotypic patterns during remission phases. Addressing these challenges, we developed a novel cell-level data integration framework using real-world clinical FCM datasets. Notably, our study represents the first systematic incorporation of remission phase AML samples into classification models, thereby establishing a much-needed benchmark for treatment response monitoring.
This study pioneers the inclusion of AML patients in complete remission (AML-CR) within a classification system and has successfully developed an artificial intelligence model that accurately differentiates among healthy individuals, AML patients, and AML-CR. The model demonstrated > 90% accuracy in all baseline tests, confirming its validity. This high accuracy further indicates that AML-CR exhibits significantly distinct characteristics compared to the other two groups, offering new perspectives for clinical diagnosis. Notably, our machine learning model shows substantial speed advantages over both manual analysis and deep learning methods. In bone marrow flow cytometry testing, the entire process from sample processing and staining to data acquisition typically requires several to over ten minutes (33). Data analysis requires professionals to manually gate and analyze antigen expression patterns while integrating clinical background for interpretation. Flow cytometry specialists at Shanghai KingMed Diagnostics report that the analytical duration varies significantly (15 minutes to several hours) depending on sample complexity, clinical requirements, and operator experience. Our diagnostic approach completes single-patient data analysis in< 1 second, demonstrating two key advantages: (1) a 100-1000× improvement in processing speed compared to conventional methods, and (2) a substantial reduction in technologist workload. Clinical implementation of this method enables real-time assessment of disease status and treatment response, supporting timely therapeutic decision-making. By applying the proposed method to actual flow cytometry data analysis, physicians can evaluate patients’ disease status and treatment response more rapidly and accurately, thereby developing more personalized treatment plans. AI-assisted flow cytometry analysis is expected to play an important role in primary medical institutions lacking flow cytometry diagnostic specialists, helping more patients benefit.
Multimodal SHAP analysis demonstrates that the key markers of the AML group (CD117, CD34, HLA-DR) exhibit strong concordance with established clinical diagnostic criteria, thereby validating their pivotal role as core immunophenotypic markers for leukemia cell identification and classification. We identified significant differences in the expression patterns of specific cell surface markers among de novo patients with AML, healthy individuals, and post-treatment AML patients who achieved complete remission. Notably, CD34 and CD117 expression levels are significantly higher in AML patients compared to both healthy individuals and remission-phase patients, whereas CD45 expression is comparatively reduced. These findings suggest a potential mechanistic link between aberrant marker expression and AML pathogenesis, therapeutic efficacy, and relapse, advancing our understanding of the disease’s biology. Furthermore, flow cytometric profiling of AML patients in complete remission facilitates the identification of therapy-responsive immunological markers and provides early warning signs for potential disease recurrence. Notably, the markers identified in the normal group demonstrated relatively low specificity, potentially reflecting physiological variations in immune homeostasis. To improve model robustness, adversarial training or sample size expansion approaches should be considered to enhance interference resistance.
Notwithstanding the meaningful contributions of this work, certain limitations merit consideration. Chief among these is the restricted generalizability inherent to single-center studies with limited sample sizes. Second, the intrinsic heterogeneity and technical variability in flow cytometry data may introduce measurement noise and analytical interference, which could adversely affect the model’s predictive accuracy. Additionally, the biological implications of specific cell surface markers warrant further investigation. Future studies should employ expanded cohorts incorporating diverse AML subtypes and treatment phases to optimize and validate the machine learning model’s performance. Furthermore, integrating genomic sequencing data and other multi-omics information to develop a multi-omics fusion model, alongside multimodal imaging data and clinical information, could provide a more comprehensive AML diagnostic and prognostic assessment tool, thereby advancing the precision and scientific rigor of clinical decision-making. The current study is limited to non-M3 AML patients. Future research should extend the model’s applicability to additional AML subtypes, including acute promyelocytic leukemia (M3) and other rare variants, to enhance its clinical utility. Looking forward, translational application of this model to the broader spectrum of hematologic neoplasms, such as acute lymphoblastic leukemia, may establish new frameworks for AI-powered.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by Ethics Committee of Shidong Hospital, Yangpu District, Shanghai, China. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants' legal guardians/next of kin because this was a retrospective study that utilizes medical records obtained from previous clinical diagnosis and treatment.
Author contributions
RB: Data curation, Writing – review & editing, Writing – original draft. MF: Writing – review & editing, Methodology, Writing – original draft, Data curation. MW: Writing – original draft, Writing – review & editing. YL: Writing – review & editing, Writing – original draft. LH: Writing – review & editing, Writing – original draft. YY: Conceptualization, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Conflict of interest
Author YL was employed by Shanghai Kingmed Center for Clinical Laboratory Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Short NJ, Rytting ME, and Cortes JE. Acute myeloid leukaemia. Lancet. (2018) 392:593–606. doi: 10.1016/S0140-6736(18)31041-9
2. DiNardo CD, Erba HP, Freeman SD, and Wei AH. Acute myeloid leukaemia. Lancet. (2023) 401:2073–86. doi: 10.1016/S0140-6736(23)00108-3
3. Kantarjian H, Kadia T, DiNardo C, Daver N, Borthakur G, Jabbour E, et al. Acute myeloid leukemia: current progress and future directions. Blood Cancer J. (2012) 11:41. doi: 10.1038/s41408-021-00425-3
4. Döhner H, Wei AH, Appelbaum FR, Craddock C, DiNardo CD, Dombret H, et al. Diagnosis and management of AML in adults: 2022 recommendations from an international expert panel on behalf of the ELN. Blood J Am Soc Hematol. (2022) 140:1345–77. doi: 10.1182/blood.2022016867
5. Haferlach T and Schmidts I. The power and potential of integrated diagnostics in acute myeloid leukaemia. Br J Haematol. (2020) 188:36–48. doi: 10.1111/bjh.16360
6. McKinnon KM. Flow cytometry: an overview. Curr Protoc Immunol. (2018) 120:5–1. doi: 10.1002/cpim.40
7. McCarthy N, Gui G, Dumezy F, Roumier C, Andrew G, Green S, et al. Pre-emptive detection and evolution of relapse in acute myeloid leukemia by flow cytometric measurable residual disease surveillance. Leukemia. (2024) 38:1667–73. doi: 10.1038/s41375-024-02300-z
8. Buldini B, Maurer-Granofszky M, Varotto E, and Dworzak MN. Flow-cytometric monitoring of minimal residual disease in pediatric patients with acute myeloid leukemia: recent advances and future strategies. Front Pediatrics. (2019) 7:412. doi: 10.3389/fped.2019.00412
9. Riva G, Nasillo V, Ottomano AM, Bergonzini G, Paolini A, Forghieri F, et al. Multiparametric flow cytometry for MRD monitoring in hematologic Malignancies: clinical applications and new challenges. Cancers. (2021) 13:4582. doi: 10.3390/cancers13184582
10. Rotem O and Zaritsky A. Visual interpretability of bioimaging deep learning models. Nat Methods. (2024) 21:1394–7. doi: 10.1038/s41592-024-02322-6
11. Shehab M, Abualigah L, Shambour Q, Abu-Hashem MA, Shambour MKY, Alsalibi AI, et al. Machine learning in medical applications: A review of state-of-the-art methods. Comput Biol Med. (2022) 145:105458. doi: 10.1016/j.compbiomed.2022.105458
12. Song X, Liu X, Liu F, and Wang C. Comparison of machine learning and logistic regression models in predicting acute kidney injury: A systematic review and meta-analysis. Int J Med informatics. (2021) 151:104484. doi: 10.1016/j.ijmedinf.2021.104484
13. Szałata A, Hrovatin K, Becker S, Tejada-Lapuerta A, Cui H, Wang B, et al. Transformers in single-cell omics: a review and new perspectives. Nat Methods. (2024) 21:1430–43. doi: 10.1038/s41592-024-02353-z
14. Rosenberg CA, Rodrigues MA, Bill M, and Ludvigsen M. Comparative analysis of feature-based ML and CNN for binucleated erythroblast quantification in myelodysplastic syndrome patients using imaging flow cytometry data. Sci Rep. (2024) 14:9349. doi: 10.1038/s41598-024-59875-x
15. Huys EH, Hobo W, and Preijers FW. OMIP-081: A new 21-monoclonal antibody 10-color panel for diagnostic polychromatic immunophenotyping. Cytometry Part A. (2022) 101:117–21. doi: 10.1002/cyto.a.24511
16. Alaggio R, Amado C, Anagnostopoulos I, Attygalle AD, Araujo IBdO, Berti E, et al. The 5th edition of the World Health Organization classification of haematolymphoid tumours: lymphoid neoplasms. Leukemia. (2022) 36:1720–48. doi: 10.1038/s41375-022-01620-2
17. Khoury JD, Solary E, Abla O, Akkari Y, Alaggio R, Apperley JF, et al. The 5th edition of the World Health Organization classification of haematolymphoid tumours: myeloid and histiocytic/dendritic neoplasms. Leukemia. (2022) 36:1703–19. doi: 10.1038/s41375-022-01613-1
18. Cox DR. The regression analysis of binary sequences. J R Stat Soc Ser B: Stat Methodol. (1958) 20:215–32. doi: 10.1111/j.2517-6161.1958.tb00292.x
20. Chen T, He T, Benesty M, Khotilovich V, Tang Y, Cho H, et al. Xgboost: extreme gradient boosting. R package version, Vol. 1. (2015). pp. 1–4.
21. Natekin A and Knoll A. Gradient boosting machines, a tutorial. Front Neurorobotics. (2013) 7:21. doi: 10.3389/fnbot.2013.00021
22. Popescu MC, Balas VE, Perescu-Popescu L, and Mastorakis N. Multilayer perceptron and neural networks. WSEAS Trans Circuits Systems. (2009) 8:579–88. doi: 10.5555/1639537.1639542
23. Hearst MA, Dumais ST, Osuna E, Platt J, and Scholkopf B. Support vector machines. IEEE Intelligent Syst Their Appl. (1998) 13:18–28. doi: 10.1109/5254.708428
24. Lundberg SM and Lee SI. A unified approach to interpreting model predictions. Adv Neural Inf Process Systems. (2017) 30:4768–77. doi: 10.5555/3295222.3295230
25. Jani CT, Ahmed A, Singh H, Mouchati C, Al Omari O, Bhatt PS, et al. Burden of AML, 1990-2019: Estimates from the global burden of disease study. JCO Global Oncol. (2023) 9. doi: 10.1200/GO.23.00229
26. Del Principe MI, De Bellis E, Gurnari C, Buzzati E, Savi A, Consalvo MAI, et al. Applications and efficiency of flow cytometry for leukemia diagnostics. Expert Rev Mol Diagnostics. (2019) 19:1089–97. doi: 10.1080/14737159.2019.1691918
27. Ghete T, Kock F, Pontones M, Pfrang D, Westphal M, Höfener H, et al. Models for the marrow: A comprehensive review of AI-based cell classification methods and Malignancy detection in bone marrow aspirate smears. HemaSphere. (2024) 8:70048. doi: 10.1002/hem3.70048
28. Dehkharghanian T, Mu Y, Tizhoosh HR, and Campbell CJ. Applied machine learning in hematopathology. Int J Lab Hematol. (2023) 45:87–94. doi: 10.1111/ijlh.14110
29. Lewis JE and Pozdnyakova O. Digital assessment of peripheral blood and bone marrow aspirate smears. Int J Lab Hematol. (2023) 45:50–8. doi: 10.1111/ijlh.14082
30. Bini L, Mojarrad FN, Liarou M, Matthes T, and Marchand-Maillet S. FlowCyt: A comparative study of deep learning approaches for multi-class classification in flow cytometry benchmarking. arXiv. (2024). Available online at: https://arxiv.org/abs/2403.00024.
31. Hu Z, Tang A, Singh J, Bhattacharya S, and Butte AJ. A robust and interpretable end-to-end deep learning model for cytometry data. Proc Natl Acad Sci. (2020) 117:21373–80. doi: 10.1073/pnas.2003026117
32. Li Y, Mahjoubfar A, Chen CL, Niazi KR, Pei L, and Jalali B. Deep cytometry: deep learning with real-time inference in cell sorting and flow cytometry. Sci Rep. (2019) 9:11088. doi: 10.1038/s41598-019-47193-6
Keywords: acute myeloid leukemia, flow cytometry, machine learning, feature importance analysis, diagnostic model
Citation: Bao R, Feng M, Wang M, Liu Y, Hu L and Yao Y (2025) Detection of acute myeloid leukemia and remission states using heterogeneous flow cytometry data. Front. Oncol. 15:1638074. doi: 10.3389/fonc.2025.1638074
Received: 30 May 2025; Accepted: 09 September 2025;
Published: 30 September 2025.
Edited by:
Sharon R. Pine, University of Colorado Anschutz Medical Campus, United StatesReviewed by:
Paola Pacelli, University of Siena, ItalyZesong Yang, First Affiliated Hospital of Chongqing Medical University, China
Copyright © 2025 Bao, Feng, Wang, Liu, Hu and Yao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yonghua Yao, c2R5eXh5azJAMTYzLmNvbQ==; Liang Hu, cmFpbm1pbGtAZ21haWwuY29t