Predicting delayed methotrexate elimination in pediatric acute lymphoblastic leukemia patients: an innovative web-based machine learning tool developed through a multicenter, retrospective analysis

Background High-dose methotrexate (HD-MTX) is a potent chemotherapeutic agent used to treat pediatric acute lymphoblastic leukemia (ALL). HD-MTX is known for cause delayed elimination and drug-related adverse events. Therefore, close monitoring of delayed MTX elimination in ALL patients is essential. Objective This study aimed to identify the risk factors associated with delayed MTX elimination and to develop a predictive tool for its occurrence. Methods Patients who received MTX chemotherapy during hospitalization were selected for inclusion in our study. Univariate and least absolute shrinkage and selection operator (LASSO) methods were used to screen for relevant features. Then four machine learning (ML) algorithms were used to construct prediction model in different sampling method. Furthermore, the performance of the model was evaluated using several indicators. Finally, the optimal model was deployed on a web page to create a visual prediction tool. Results The study included 329 patients with delayed MTX elimination and 1400 patients without delayed MTX elimination who met the inclusion criteria. Univariate and LASSO regression analysis identified eleven predictors, including age, weight, creatinine, uric acid, total bilirubin, albumin, white blood cell count, hemoglobin, prothrombin time, immunological classification, and co-medication with omeprazole. The XGBoost algorithm with SMOTE exhibited AUROC of 0.897, AUPR of 0.729, sensitivity of 0.808, specificity of 0.847, outperforming the other models. And had AUROC of 0.788 in external validation. Conclusion The XGBoost algorithm provides superior performance in predicting the delayed elimination of MTX. We have created a prediction tool to assist medical professionals in predicting MTX metabolic delay. Supplementary Information The online version contains supplementary material available at 10.1186/s12911-023-02248-7.

Predicting delayed methotrexate elimination in pediatric acute lymphoblastic leukemia patients: an innovative web-based machine learning tool developed through a multicenter, retrospective analysis Introduction Acute lymphoblastic leukemia (ALL) is a prevalent neoplasm in childhood.The incidence of ALL in children below 15 years of age is 0.004%, which accounts for about 35% of all cases of pediatric malignancies [1][2][3].Epidemiological studies of ALL indicate that the cumulative incidence is 1/2000 under the age of 15 [4].
Methotrexate (MTX) is a crucial antineoplastic agent in ALL therapy, which inhibits the synthesis of tumor cells by restraining dihydrofolate reductase.In clinical practice, high-dose methotrexate (HD-MTX) can significantly increase the blood drug concentration and permeate blood-brain and blood-testis barriers, so it is recommended as a common chemotherapy approach for ALL treatments.Although HD-MTX is deemed an effective ALL treatment, prolonged exposure to HD-MTX can cause hepatotoxicity, nephrotoxicity, and neurotoxicity [2,[5][6][7][8].A study from China revealed the rate of delayed MTX elimination was as high as 12.1%, which is a nonnegligible rate [9].A clinical trial has demonstrated that 2-12% of patients develop acute kidney injury (AKI) despite appropriate support during HD-MTX treatment [2,10].Furthermore, the severity of adverse reactions of MTX is linked to the concentration and duration of drug exposure.Since the liver and immune system of children are not yet fully developed, their tolerance and metabolic capacity to potential liver toxicity of drugs are inadequate.Therefore, children are more prone to delayed MTX elimination, which could affect their prognosis or lead to other adverse outcomes.Consequently, it is crucial to find ways to reduce the delayed elimination of MTX and the incidence of side effects.
To address the problem of delayed MTX elimination, the current approach is to monitor MTX concentration at 24 h, 48 h, and 72 h post-administration and to administer calcium leucovorin rescue agent and urine alkalization if necessary to accelerate MTX elimination.However, the risk of delayed elimination cannot be predicted based on patient's signs and data before medication.Therefore, early warning and timely intervention are crucial to effectively reduce the risk of delayed MTX elimination and prevent serious adverse drug reactions.
Artificial intelligence (AI) has been widely used in the medical field.In previous studies, machine learning (ML) was used to classify diseases and analyze the survival of prognosis [11,12].Researchers not only extracted disease features for building models, but also achieved high accuracy.This can reduce the fluctuation of patient incidence rate and save on medical costs.Therefore, it is necessary to apply ML to predict the metabolic delay of methotrexate.Researchers, such as Wang Yang [13], Yang Fan [14], and Min Zhang [7], have begun using ML to develop prediction models for delayed MTX elimination.However, previous studies have encountered various issues such as small sample sizes, inadequate representation, limited model construction methods, and insufficient comparability.Additionally, predictive indicators failed to fully consider patient clinical data and relevant clinical laboratory indicators.
This study aims to assess the potential correlation between premedication indicators and delayed elimination of MTX by integrating electronic medical data from multiple centers.Furthermore, a prediction model will be developed using ML methods and a web-based tool to offer an early warning for the delayed elimination of MTX in clinical settings.

Statistical analysis
The patients were randomly divided into a training set and a test set at a ratio of 7:3 using a random number table .The training set was utilized to select predictors and construct the prediction model, while the test set was used to evaluate the performance of the model.All statistical analyses were conducted in R for Windows (version 3.6.1,https://www.r-project.org/) and SPSS 25.0 (IBM Corporation, Armonk, NY, USA).The random forest algorithm was used to fill in missing values that were less than 30%.
Initially, the normality of continuous variables was assessed using the Shapiro-Wilk test.The t-test was utilized for normal data, while the Mann-Whitney test was used for non-normal data in the univariate analysis.Additionally, the Pearson chi-square test was used for categorical variables.The significant indicators selected by univariate analysis were further filtered using the least absolute shrinkage and selection operator (LASSO) regression method.To address the issue of imbalanced data sets, we conducted three different sampling methods on imbalanced datasets.Oversampling, under-sampling and Synthetic Minority Oversampling Technique (SMOTE) was employed to balance the data sets.MLbased prediction models were constructed using the predictors filtered by LASSO.In the model construction, four ML models were developed, including extreme gradient boosting (XGBoost), random forest classifier (RFC), adaptive boosting (AdaBoost), and light gradient boosting machine (LightGBM).The grid search algorithm was employed to determine the optimal parameters of the model.The area under the receiver operating characteristic curve (AUROC) and the area under the precisionrecall curve (AUPR) were used to evaluate the model performance.Additionally, SHapley Additive exPlanation (SHAP) was utilized to interpret the chosen model and complete SHAP visualization.The entire statistical analysis process is shown in Fig. 1.In the previous research on the prediction model of MTX delayed elimination, in addition to using ML, logistic regression was also used.So, we also build a logistic regression nomogram to compare its performance with optimal ML.Finally, we use an external validation set to ensure the generalization and consistency of the model.

Study population
In our research's dataset (1729 cases), there were 329 and 1400 cases with and without metabolic delays, respectively.After proportionally dividing the dataset with a ratio of 7:3, the training set (1210 cases) comprised 230 patients with metabolic delay and 980 patients without metabolic delay.In the test set (519 cases), 99 patients experienced metabolic delay and 420 patients did not.The external validation set includes 1090 data cases.

Feature selection and data preprocessing
Upon conducting the Shapiro-Wilk test, it was found that all variables were non-normally distributed (as per Additional Table 2).Consequently, we employed the Mann-Whitney U test to compare the continuous variables.Our analysis revealed that age, weight, Cr, UA, TBIL, ALB, ALT, PCV, WBC, HGB, LDH, and PT were statistically significant between the two groups.Furthermore, the Chi-square test indicated that immunological classification, ALL risk level, and co-medication with omeprazole displayed significant differences between the two groups (as illustrated in Table 1).We subsequently performed a LASSO regression analysis on the 15 significant predictors identified through univariate analysis.The paths of the coefficients with different log-transformed λ values in LASSO regression model was displayed in Fig. 2, which clearly demonstrates the significance of several variables, with the influence on delayed MTX elimination increasing as the line moves closer to zero.Moreover, the crossvalidation error plot of the LASSO regression model was depicted in Fig. 3. To create a more simplified model, we selected the top 11 variables that had the greatest impact on the outcome.Ultimately, the LASSO method identified eleven indicators, including age, weight, Cr, UA, TBIL, ALB, WBC, HGB, PT, immunological classification, and co-medication with omeprazole, which were used to develop our predictive models.

Model evaluation and interpretation
The variables selected previously were utilized as input variables to establish a prediction model for delayed MTX elimination, with the occurrence of delayed MTX elimination being designated as the outcome event (yes = 1, no = 0).Ultimately, a total of 230 patients with delayed MTX elimination and 980 patients without delayed MTX elimination were included in the training set to develop the predictive model.The test set was then used to validate the predictive ability of the established model.The performance of the delayed MTX elimination risk prediction models with different sampling methods are showed in Additional Table 3.
We chose the XGBoost model sampled by SMOTE as the optimal model for this study.The AUROC   4. The AUPR value is more sensitive to sample distribution, and the precision-recall (P-R) curve to showcase the model's precision and recall performance (Fig. 5).The AUROC value of XGBoost using SMOTE is 0.897(0.857-0.937)and it had an area under the P-R curve (AUPR) of 0.729.In addition, XGBoost sensitivity in SMOTE is 0.808.The higher the sensitivity, the better the model's ability to correctly identify delayed elimination, and the lower the missed diagnosis rate.The comparison process for selecting the optimal model can be found in Additional File 1.We apply the optimal model to predict external validation sets.It was found that AUROC = 0.788 (0.753-0.822) in external validation, indicating good discrimination ability.We apply the optimal model to predict external validation sets.We used the optimal model to predict the external validation set, and the model demonstrated good performance.Among them, AUROC = 0.788 (0.753-0.822),AUPR = 0.648, specificity = 0.813 (0.780-0.840), sensitivity = 0.680 (0.625-0.735).
As illustrated in Fig. 6, the summary graph of SHAP elucidates the prediction of all samples.The SHAP values of each sample's variable were plotted by scatter plot, and the relationship between SHAP values and outcomes was analyzed.In the XGBoost model, the SHAP summary plot ranked the importance of delayed MTX elimination variables as co-medication with omeprazole, Cr, UA, WBC, HGB, Age, HGB, ALB, immunological classification, weight, PT and TBIL.Additionally, a dependence plot was generated to assess the relationship between the variables and the predicted influence (Additional Figs.[1][2][3][4][5][6][7][8][9][10][11].The dependency graph lucidly portrays how individual variables affect the model's predictions. We constructed a Logistic regression nomogram using the 11 screened indicators.Figure 7 shows an example of using nomogram to predict MTX delayed elimination.The total score corresponds to the probability  value on the risk axis, and a higher total score indicates a higher risk of MTX delayed elimination.We evaluated the nomogram with an AUROC of 0.886(0.844-0.929)as shown in Additional Fig. 12.

Discussion
Several research studies have illustrated that prolonged elimination after administering HD-MTX to children with ALL may result in serious adverse effects, particularly in those with atypical renal function [2,[5][6][7][8]10].We formulated a risk assessment algorithm for predicting delayed MTX elimination based on pre-medication information.This can facilitate healthcare professionals in recognizing the possibility of delayed MTX elimination in children with ALL.
In this study, age, weight, Cr, UA, TBIL, ALB, WBC, HGB, PT, immunological classification, and concurrent use of omeprazole were recognized as risk factors for delayed MTX elimination.Most of these autonomous risk factors have been reported in preceding research [8,[13][14][15][17][18][19][20][21][22][23][24][25].For instance, Nakano T discovered that age, MTX dosage, and TBIL were independent risk factors for delayed MTX elimination [8].Xu's research revealed that scrutinizing serum Cr concentration can proficiently anticipate the delay of MTX elimination, and that patients with delayed metabolism have elevated serum Cr levels [22].A Japanese study found that serum UA levels were correlated with nephrotoxicity prompted by delayed MTX elimination [23].Another analysis indicated that MTX toxicity could be engendered by combining proton pump inhibitors (such as omeprazole), penicillin family antibiotics, and specific antimicrobial agents [24][25][26][27].We have retained most of the previous studies on influencing factors, while additionally incorporating FIB, PT, chloride in cerebrospinal fluid, and cerebrospinal fluid transparency.These parameters are easily obtainable in medical facilities, and multifarious possibilities of causing MTX metabolism delay are exhaustively contemplated.For instance, HD-MTX therapy will prolong thrombin time and diminish FIB [28].Additionally, distinct dosages of MTX exhibit notable drug concentrations in serum and cerebrospinal fluid [29].The predictors WBC, HGB, and PT are seldom mentioned in preceding studies and require further validation.
Recently, ML techniques have garnered increasing attention in clinical research and emerged as a powerful instrument for addressing numerous healthcare problems [30][31][32].In this investigation, we compared the performance of different ML models in different sampling methods for imbalanced data.Among these models' evaluation, we found that the XGBoost in SMOTE and Light-GBM in oversampling were comparable in performance.However, XGBoost demonstrated the better AUPR value and sensitivity.Nitesh Chawla et al. described that smote works by selecting the nearest instances in the feature space, drawing a line between the instances in the feature space, and drawing a new sample along a point of the line [33].Consequently, we ultimately opted for XGBoost in SMOTE to construct the final prediction model.XGBoost is extensively utilized by data scientists and delivers the most cutting-edge outcomes on a plethora of issues.For instance, XGBoost forestalls overfitting and has the ability to handle voluminous data [34].Luu Ho Thanh Lam et al. selected XGBoost as the optimal model after SMOTE, to classify the molecular subtypes of lowgrade glioma [35].Nwanosike EM et al. evaluated the advancements of ML algorithms in clinical applications, and the XGBoost algorithm exhibited the highest potential for clinical implementation [36].We have also implemented the optimal prediction model on the web page to provide a reliable tool for clinical medical professionals and researchers.The web page address is https://cqmugj.shinyapps.io/mtx_jc/.
We constructed a nomogram, which was commonly used in previous studies to predict MTX delayed elimination.We found that the AUROC value using the nomogram was smaller than that of the optimal model (XGBoost).On the other hand, nomogram is a non-parametric model that requires the total score to obtain the probability.And it can't automatically calculate the result, which is a bit inconvenient compared to ML.In addition, the model's AUROC and specificity after external validation indicated that it had good discrimination and a low misdiagnosis rate.And the result also reflected the transportability and generalization ability of the model.On the other hand, it indicates that the model has good consistency in different time periods compared to the model development queue.
The current research is mainly to accurately diagnose the adverse reaction or MTX delayed elimination by using the post medication test index of methotrexate combined with ML.We summarize some similar studies and draw a Table 2.For example, Hu et al. created an ML-based model for predicting low-dose MTX-related hepatotoxicity with an AUC of 0.97 but only accuracy of 0.64 [37].Zhan et al. employed an artificial intelligence algorithm to forecast neutropenia and fever caused by high-dose MTX in children with B-cell ALL, with an AUC of 0.870 [38].The performance of our model is similar to that of Zhan M et al. [7], but inferior to Schmidt, Fig. 7 A constructed nomogram for prediction of delayed MTX elimination in Pediatric ALL Patients D [13].In addition, we summarized some researches on analysis of MTX delayed elimination factors in recent years (see Additional Table 1).However, few studies have integrated the identified risk factors and applied them directly to the prediction of delayed MTX elimination.Zhan M et al. used hematocrit, risk classification, dose, SLC19A1 rs2838958, and sex indicators to develop a prediction model for delayed elimination of MTX.The highest AUC of the model was 0.807 (95% CI, 0.724-0.889)[7].They used fewer variables and included genetic factors to build a prediction model with better performance.However, our predictors are easily obtainable and it is of great value in identifying MTX metabolic delay.
Nonetheless, the study has certain limitations.Firstly, the incidence, treatment, and individual differences in ALL across different regions may hinder the applicability of the model.Secondly, some variables, such as the genetic characteristics of the affected children and their living environment, have not been included.Thirdly, our study was retrospective research, the examination of some cases was done with inadequate equipment and training, and some indicators with missing values greater than 30% (e.g.urine volume) were not included in the model.Finally, the generalization ability of the model should be further confirmed through multi-center external validation in future studies.

Conclusions
In summary, this investigation illustrates that factor such as age, body weight, creatinine, uric acid, total bilirubin, albumin, white blood cell count, hemoglobin, prothrombin time, cellular morphological classification, and concomitant use of omeprazole could be served as predictors for delayed MTX elimination.Through the application of XGBoost after SMOTE, delayed MTX elimination can be effectively identified in children diagnosed with ALL.Our predictive model provides a reliable means for monitoring the metabolic delay of MTX, even in the absence of MTX plasma concentration monitoring.By utilizing this tool, medical professionals can take timely targeted measures to prevent the occurrence of MTX-related adverse drug events.

Fig. 2
Fig. 2 Coefficient regression graphThe horizontal coordinate is the magnitude of the λ value in the LASSO regression model.As the λ value changes, the later the coefficient is compressed to zero the more influential the variable is.The graphs show that age, TBIL, and Immunological Classification are highly significant

Fig. 4 Fig. 3 Fig. 6 Fig. 5
Fig. 4 ROC curve of 4 ML models for predicting MTX delayed elimination in the testing set

Table 1
Characteristics of patients with and without delayed MTX elimination performance of the delayed MTX elimination risk prediction model with SMOTE is illustrated in Fig.

Table 2
Summary table of machine learning applied to MTX delayed elimination or Adverse reactions