Impact Factor 3.845 | CiteScore 3.92
More on impact ›

Original Research ARTICLE

Front. Pharmacol., 07 October 2019 | https://doi.org/10.3389/fphar.2019.01155

Early and Accurate Prediction of Clinical Response to Methotrexate Treatment in Juvenile Idiopathic Arthritis Using Machine Learning

  • 1Department of Pharmacy, Guangzhou Women and Children’s Medical Center, Guangzhou Medical University, Guangzhou, China
  • 2Institute of Clinical Pharmacology, School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
  • 3Institute of Pediatrics, Guangzhou Women and Children’s Medical Center, Guangzhou Medical University, Guangzhou, China
  • 4Pediatric Allergy Immunology & Rheumatology Department, Guangzhou Women and Children’s Medical Center, Guangzhou Medical University, Guangzhou, China
  • 5Department of Medical, Guangzhou Women and Children’s Medical Center, Guangzhou Medical University, Guangzhou, China

Background and Aims: Accurately predicting the response to methotrexate (MTX) in juvenile idiopathic arthritis (JIA) patients before administration is the key point to improve the treatment outcome. However, no simple and reliable prediction model has been identified. Here, we aimed to develop and validate predictive models for the MTX response to JIA using machine learning based on electronic medical record (EMR) before and after administering MTX.

Materials and Methods: Data of 362 JIA patients with MTX mono-therapy were retrospectively collected from EMR between January 2008 and October 2018. DAS44/ESR-3 simplified standard was used to evaluate the MTX response. Extreme gradient boosting (XGBoost), support vector machine (SVM), random forest (RF), and logistic regression (LR) algorithms were applied to develop and validate models with 5-fold cross-validation on the randomly split training and test set. Data of 13 patients additionally collected were used for external validation.

Results: The XGBoost screened out the optimal 10 pre-administration features and 6 mix-variables. The XGBoost established the best model based on the 10 pre-administration variables. The performances were accuracy 91.78%, sensitivity 90.70%, specificity 93.33%, AUC 97.00%, respectively. Similarly, the XGBoost developed a better model based on the 6 mix-variables, whose performances were accuracy 94.52%, sensitivity 95.35%, specificity 93.33%, AUC 99.00%, respectively.

Conclusion: Based on common EMR data, we developed two MTX response predictive models with excellent performance in JIA using machine learning. These models can predict the MTX efficacy early and accurately, which provides powerful decision support for doctors to make or adjust therapeutic scheme before or after treatment.

Introduction

Methotrexate (MTX) is the first line treatment for the majority of patients with juvenile idiopathic arthritis (JIA). However, the efficacy of MTX varies greatly among individuals, with about 30 to 70% of JIA patients being effective (Ruperto et al., 2004; Foell et al., 2010). Patients who respond to MTX poorly are given biologicals alone or in co-treatment with MTX. Biologics can lead to more efficient disease control, but abuse of biologics can result in high costs and serious adverse reactions. Additionally, it usually takes 3–6 months before a decision is made as to MTX efficacy (Martini et al., 2019). Patients receiving “trial-and-error” therapy for such a long time may delay treatment, resulting in irreversible joint damage and even adverse reactions. Therefore, early identification of whether the patient is effective before starting MTX and then selection of appropriate therapy (MTX alone or combined with biologics) are of great significance for preventing disease progression. This means that it is very necessary to establish an efficacy prediction model before the onset of MTX in JIA.

Although MTX has been used to treat JIA for a long time, being able to predict who will respond to MTX is still very limited. To date, only Bulatovic et al. (2012) reported a predictive model for MTX response to JIA. However, the limitations of this model are as follows: the prediction accuracy was not high (the area under the curve, AUC, was only 72%); model variables contained controversial single nucleotide polymorphisms (SNPs), which required additional and expensive testing, thus limiting its widely available in clinical application. Moreover, this study only employed one traditional logistic regression algorithm, which is not applicable to the modeling of non-independent variables. In addition to this study, other studies on MTX response to JIA were only limited to discovering which indicators would affect the efficacy of MTX. But they did not provide a model for clinical application, so that it could not be easily applied in clinical practice (Hinks et al., 2011; Yanagimachi et al., 2011; Cobb et al., 2014; Zajc Avramovic et al., 2017).

Therefore, a simple, efficient and accurate MTX response prediction model is urgently needed to provide references for clinicians before treatment. In recent years, the predictive model developed by machine learning based on electronic medical record (EMR) data has played an excellent role in disease diagnosis, treatment, and prognosis. For example, in our previous work, we used a machine learning technique to acquire pediatric EMR and developed an auxiliary decision-making system for diseases diagnosis, which is comparable to that of human physicians (Liang et al., 2019); Machine learning is also used to predict the efficacy and prognosis of diseases in other diseases (Motwani et al., 2017; Browning et al., 2019). Similarly, in rheumatoid diseases, researchers used machine learning to establish disease diagnosis classifier, mortality prediction model and MTX related hepatotoxicity automatic recognizer basing on EMR data (Liao et al., 2010; Lin et al., 2015; Lezcano-Valverde et al., 2017). These results provide practical tools for the management of patients. However, currently, there are no reports about the prediction model of MTX response in JIA using machine learning only basing on EMRs.

The purpose of this study is to develop simple, efficient and accurate models using machine learning for early predicting the efficacy of MTX in JIA based on integrating temporal features before and after starting MTX within three months.

Methods

Study Design and Population

We retrospectively collected the EMR data of children with JIA who visited Guangzhou Women and Children’s Medical Center from January 2008 to October 2018. Inclusion criteria were: (1) patients were new-onset and met the International League of Associations for Rheumatology criteria for JIA (Petty et al., 2004; Martini et al., 2019). (2) The onset age was 1–16 years old. (3) Patients received monotherapy with MTX for at least 3 months. (4) Co-treatment with non-steroidal anti-inflammatory drugs or corticosteroids were allowed. Exclusion criteria were: (1) combined therapy with other interfering drugs (e.g. biologic agents, sulfasalazine, etc.) within 3 months. (2) MTX therapy did not reach 3 months. (3) Serious missing of medical records. A total of 674 JIA children using MTX were screened out, but 362 patients were eventually included for modeling and validating. Furthermore, we continued to collect 13 JIA patients from November 2018 to January 2019 for external verification.

The study was performed according to the Helsinki declaration. Ethical approval was obtained from the ethics committee of this center (no. 2016021645). This study was a part of a large clinical trial (NCT81603203). All data were anonymous and no identifiable personal data of patients were available for the analysis. No additional informed consent was required.

Assessment of MTX Clinical Response

Weekly MTX was given to all patients by either oral or subcutaneous route at 10–15 mg/m2. Baseline disease activity was calculated before MTX treatment. Early response to MTX was evaluated at 3 months after using MTX. Since it is a retrospective study, it is difficult to collect subjective features such as patient/parent and physician’s global assessment of disease. Therefore, JADAS or ACRpedi (see Table 1 for the full name) scoring tools could not be applied to evaluate the response (Giannini et al., 1997; Consolaro et al., 2009). DAS44/ESR-3, a simplified standard related to the European League of Associations for Rheumatology criteria, was the most suitable choice for this retrospective study (Ranganath et al., 2007; Consolaro et al., 2009). The simplified formula of the disease activity is as follows: y=0.53938RAI+0.06465*SJC44+0.33ln(ESR)+0.224 (RAI, Ritchie articular index; SJC, Swollen joint count; ESR, erythrocyte sedimentation rate). The response was defined as a significant change of DAS44 scores from baseline to 3 months after starting MTX. Good response was defined as a significant decrease in DAS44 (>0.6), while a decrease of ≤0.6 was non-response.

TABLE 1
www.frontiersin.org

Table 1 The full name and abbreviation name of variables.

Clinical Variables

All data were collected from EMRs before administration of MTX (baseline) and within 3 months after administration (nearly 3 months). We collected a lot of variables, including: joint conditions [tender joint count (TJC), SJC, RAI, and JIA subtypes (oligoarticular, polyarticular, and other subtypes), joint imaging, etc.], the acute phase of inflammatory products [C-reactive protein (CRP), ESR], demographic data (age, gender, etc.), immune-related indicators [rheumatoid factor (RF), rheumatoid factor IgG (RF-IgG), antinuclear antibodies, anti-cyclic citrullinated peptide antibody, etc.], kidney function, liver function [total bilirubin (TBIL), direct bilirubin (DBIL), indirect bilirubin (IBIL), etc.], blood coagulation function [active partial thrombin time (APTT), prothrombin time (PT), thrombin time (TT), fibrinogen (FIB), etc.], blood routine testing, relevant lymphocytes (CD3+abs, CD3+CD4+, CD3+CD8+, etc.) and other data. See Table 1 for a list of all variables. Because some variables were seriously missing, they were not be used for modeling. All variables used for modeling are shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1 The variables used for modeling and their importance ranking (in order of median importance). The left part of the figure shows the variables used in the pre-administration variables model. The right part of the figure shows a mixture of variables before and after administration (variables collected within 3 months after administration with MTX). The shorter the transverse column (i.e. the smaller the value), the greater importance of the median ranking of the variable (see the top and left part).

Machine Learning

We collected two different sets of variables. Thus, two groups of models were established based on variables before the onset of MTX and mix-variables within 3 months after starting MTX respectively, for finding the best model. In the first group of models, referred to as pre-administration variables models (MTX-A), we included 46 variables (see the left part of Figure 1). In the second group of models, referred to as mix-variables models (MTX-B), we extended MTX-A by adding 32 new variables (see the right part of Figure 1). The main process can be divided into three steps: (1) data processing, (2) feature selection, (3) model generation and validation. Five-fold stratified cross-validation was used for assessing the performance and general error estimation of feature selection and model generation. Figure 2 shows the flowchart of this work. Machine learning techniques were implemented in Python 3 (Python 3.6.5) using the package Scikit-learn (Scikit-learn 0.19.1).

FIGURE 2
www.frontiersin.org

Figure 2 The flowchart of model developing and validation. ET, extremely randomized trees; GBDT, gradient boosting decision tree; RF, random forest; XGBoost, extreme gradient boosting; SVM, support vector machine; LR, logistic regression.

Data Preprocessing

Some variables have been removed with >30% missing rates. In order to get a higher quality of the dataset, the missing values were filled with mean values of a group stratified by MTX response. For example, we used mean values of good response and non-response group to respectively fill the data of CRP in different outcome groups.

Feature Selection

Appropriate feature subsets were selected using ensemble models, including extremely randomized trees (ET), gradient boosting decision tree (GBDT), random forest (RF) and extreme gradient boosting (XGBoost). Firstly, data transformation was carried out on continuous variables to form four kinds of data: min-max normalization, z-score normalization, L2 normalization and original. Secondly, the above four algorithms were used to analyze the four forms of data and build 16 models using five-fold stratified cross-validation, so as to obtain the median importance ranking of variables in all models which is as the final features importance ranking. Finally, the next aim is to determine the feature set with the least variables but the highest predictive accuracy. The XGBoost algorithm was used to find the minimum-size list of features by forwarding feature selection, as follows: (1) beginning with the head of the ranked list of variables (the most important variable), XGBoost algorithm iteratively generates a new model by adding one variable at a time, and calculates its classification accuracy. (2) The list with the minimum size and optimum accuracy is therefore selected.

Model Generation and Validation

A cohort of 362 patients was randomly split into the training set and test set according to the ratio of 80:20. XGBoost, RF, support vector machine (SVM) and logistic regression (LR) algorithms were applied to develop classifiers respectively in our study. The classifiers were trained on the training set (n = 289), using the training set feature values (the minimum-size list of features) as input. Thus, each set of variables had 4 types of classifiers. Five-fold stratified cross-validation was used for internal validation. After training, classifiers were asked to predict the response of the test set (n = 73). For each set of variables, the four classifiers were compared with each other in terms of accuracy, and then the best classifier was selected as the final predictor. We further performed an external validation of the above classifiers with the subsequent collection of 13 patients.

Results

Patient Characteristics

A total cohort of 362 patients with JIA was included in developing models, and 13 patients were subsequently collected to external validate the best model. Table 2 describes the baseline characteristics of our study population. According to DAS44/ESR-3 simplified standard, 213 patients were rated as good response and 149 patients as non-response.

TABLE 2
www.frontiersin.org

Table 2 Baseline patient characteristics.

Feature Selection

Median importance ranking of all variables before using MTX and mix-variables before and after administering with MTX were shown in Figure 1 (left part and right part). The XGBoost algorithm was applied for selecting the minimum size and optimum accuracy features subset, and the results of this process were shown in Figure 3. We examined the predictive performance of the most prominent feature and identified the point at which there was no considerable gain in accuracy, sensitivity, and AUC, when adding the feature of the next highest ranking one to the model. The optimum values were obtained when these three measurements defined the most discriminative features. In the MTX-A predictors, the three measurements reached the optimum when 10 feature subsets were selected (see Figure 3A). The 10 selected significant variables are listed above the dotted red line in the left part of Figure 1. In the MTX-B predictors, the three measurements achieve maximum performance when 6 feature subsets were screened out (see Figure 3B). Variables above the dotted red line in the right part of Figure 1 are these 6 features. The degree of contribution of all the above selected variables to response and the formulas behind modeling were described in detail in the Supplementary.

FIGURE 3
www.frontiersin.org

Figure 3 The forward feature selection results for the model of pre-administration variables (the A part of figure 3) and the model of mix-variables before and after administration (the B part of figure 3). We examined the predictive performance of the most prominent feature and identified the point at which there was no considerable gain in accuracy, sensitivity, and area under the curve (AUC), when adding the feature of the next highest ranking one to the model. The optimum values were obtained when these three measurements defined the most discriminative features.

Model Performance and Comparison

Table 3 shows the classification accuracy results of the models which were evaluated using the test set. Of the MTX-A and MTX-B predictors, both the XGBoost models showed the best predictive performances. Therefore, the XGBoost models were selected as the final predictors. The performance of MTX-A XGBoost predictor was as follows: sensitivity 90.70% (95%CI: 82.0–99.4%), specificity 93.33% (95%CI: 84.4–100%), accuracy 91.78% (95%CI: 85.5–98.1%) and AUC 0.97. However, the MTX-B XGBoost predictor have a better performance, which achieves sensitivity 95.35% (95%CI: 89.0–100%), specificity 93.33% (95%CI: 84.4–100%), accuracy 94.35% (95%CI: 89.3–99.7%) and AUC 0.99. Figure 4 shows the mixed matrix results of each model in the MTX-A and MTX-B predictors of the test set.

TABLE 3
www.frontiersin.org

Table 3 The classification performance results of the models.

FIGURE 4
www.frontiersin.org

Figure 4 The mixed matrix results of each model in the MTX-A and MTX-B predictors of the test set. For example, when the true lables were NR(non-response) and predicted lables were NR, it indicated that the number of NR was correctly predicted. However, when true lables were NR, and predicted labels were GR (good response), it indicated the number of NR incorrectly predicted to GR. As can be seen from the figure, the predicted values of model A (XGBoost model) and model E (XGBoost model) are the closest to the real values, indicating the best prediction performance. The numbers in the pink grids represent the number of cases that were accurately predicted.

External Verification and Clinical Application

Data of 13 patients newly collected were applied to externally validate the two XGBoost predictors. The performances of MTX-A XGBoost predictor were sensitivity 100.0%, specificity 86.1% and accuracy 90.0%, respectively. But the MTX-B XGBoost predictor got better performance results, with 100% sensitivity, specificity, and accuracy. Furthermore, we applied these predictors to two randomly selected patients to predict their outcomes. We input their clinical variables into the two predictors. Both models produced correct predictive outcomes (see Table 4).

TABLE 4
www.frontiersin.org

Table 4 The application of XGBoost predictors for clinical patients to predict their response.

Discussion

We developed and validated two prediction models for MTX response in a large JIA cohort according to EMR data using machine learning. Models developed by XGBoost showed the best performance. CRP, CD3+Abs, RF-IgG, TJC, DBIL, IBIL, APTT, PT, TT, and FIB were variables screened out by the pre-administration variables model; and the mix-variables model filtered out features as follows: CRP/3m, CD3+CD4+/3m, CD3+CD8+/3m, RF-IgG/3m, TBIL/3m, and FIB, which were collected before and after administering MTX. The pre-administration variables model could accurately identify ninety-seven percent of patients whether responded to MTX (AUC 97%); 99% of patients were distinguished by the mix-variables model (AUC 99%).

To our knowledge, this is the first article to establish an accurate and easy-to-use predictive model for the efficacy of MTX in JIA based solely on EMRs using machine learning. Bulatovic et al. (2012) reported the first and so far the only one efficacy prediction model, which including variables of SNPs and ESR. The identified accuracy, sensitivity, specificity, positive predictive value, negative predictive value were 72, 78, 49, 83, and 41%, respectively. These results were all lower than ours. SNPs are important in precision medicine. However, a GWAS study found no direct correlation between SNPs related to the pathways of MTX and the efficacy in 759 JIA patients (Cobb et al., 2014). This suggests that SNPs may not necessary in revealing the outcome of JIA (Albers et al., 2009; Roszkiewicz and Smolewska, 2017). Moreover, the expression and activity of enzymes in children may be affected by growth, and the genotype may not directly reflect the phenotype (Zhao et al., 2010; Roszkiewicz and Smolewska, 2017). In addition, genotype detection increases the cost and time, which is not as convenient, efficient and cheap as conventional detection. These suggest that clinical phenotypic variables may be more important in influencing outcomes (Roszkiewicz and Smolewska, 2017). Our results verified the above view.

Of course, the low prediction accuracy of Bulatovic’s study may be because of using only the traditional logistic regression, which may not be the optimal method. Machine learning is hot in recent years, which has gained remarkable achievements in biomedicine (Austin et al., 2013; Lin et al., 2015; Motwani et al., 2017; Browning et al., 2019; Liang et al., 2019). Specifically, machine-learning approaches may offer advantages over conventional techniques. The advantages of machine learning over traditional modeling methods are as follows: (1) machine learning can deal with more complex, high-dimensional and interactive variables, but the latter has limited ability to fix those problems. (2) Traditional modeling has poor generalization, though the former can model with strong generalization and better accuracy (Kruppa et al., 2012; Lee et al., 2018). Therefore, several advanced machine learning methods including SVM, RF, and XGBoost as well as LR were applied to model and compare the results in our study. Of our results, the performances of XGBoost models were the best, and the results of LR models were mostly poor. Further, the performances of those four modeling methods were all better than those of Bulatovic’s study. These results also confirmed that XGBoost could effectively avoid overfitting and improve prediction performance. LR, as a kind of traditional analysis, may appear low-fitting (Lee et al., 2006). Therefore, when dealing with similar classification problems, it may be more appropriate to select advanced algorithms.

In addition to the above Bulatovic’s study, other literature also explored which variables were associated with the efficacy of MTX (Yanagimachi et al., 2011; Cobb et al., 2014; Zajc Avramovic et al., 2017). But they did not establish a model, which could not be applied easily in clinical. For instance, some studies showed SNPs were associated with MTX efficacy (Hinks et al., 2011; de Rotte et al., 2012; Zajc Avramovic et al., 2017) clinical variables like TJC, ESR, CRP, etc. were also correlated with MTX response (Albers et al., 2009; Yanagimachi et al., 2011; Cobb et al., 2014; Franova et al., 2016). Variables considered in our study (n = 78) were more than and mostly different from those of the reported studies (Albers et al., 2009; Yanagimachi et al., 2011; Cobb et al., 2014; Franova et al., 2016). TJC, CRP and RF-IgG screened out by our study were reported to associate with disease activity, prognosis, efficacy of rheumatoid arthritis or JIA (Zborovskii et al., 1999; Cabral et al., 2005; Jaskowski et al., 2010). CD3+, CD4+, and CD8+ are the most important T lymphocytes, which are widely distributed in the joint synovial membrane and fluid and play an important role in the pathogenesis and classification of JIA (Finnegan et al., 2011; Hui et al., 2012; Van Nieuwenhove et al., 2019). One mechanism of MTX action is to regulate or inhibit T cell immunity to achieve the therapeutic effect (Johnston et al., 2005). In this study, these cells were contributed to MTX outcome, which is consistent with the MTX effect mechanism and similar to some reports (Isaacs, 2007; Bulatovic Calasan et al., 2015). Liver function variables (TBIL, DBIL, and IBIL) were also screened out. It is well-known that MTX is metabolized by the liver, so its pharmacokinetics will be affected by liver function, which will then affect pharmacodynamics. Coagulation markers (APTT, PT, TT, and FIB) also contributed to MTX efficacy in our study. We know ESR is an important variable in calculating disease activity and MTX efficacy. While ESR depends mainly on plasma FIB (Bedell and Bush, 1985). Additionally, FIB affects the adhesion, spread, proliferation of endothelial cells, and the repair of joint synovial tissue. So FIB may have an indirect contribution to MTX treatment (Carney et al., 1992; Cid et al., 1993; Grober et al., 1993). As for the degree of contribution of all the above selected variables to the outcome (efficacy), we described details in the Supplementary. It can be seen from the Supplementary that CRP was the most significant variable for the pre-administration model, while for the mix-variables model, RF-IgG/3m was the most important.

Additionally, from our results, the mix-variables model was better than the pre-administration variables model. This indicated that temporal features after administration but before evaluating efficacy also had an important influence on the treatment outcome. Just as variables before and during pregnancy can both have an impact on pregnancy outcome. This suggests that in addition to considering the influence of pre-administration variables, the efficacy should be evaluated in combination with post-administration features. The limits of this study were the relatively small sample size, insufficient representativeness of externally verified samples and retrospective research. The next step is to conduct prospective studies to model and validate.

Conclusions

In summary, for the first time, based on EMR we used advanced machine learning to establish two early predictive models for the MTX efficacy in JIA, including pre-administration variables model and mix-variables model. The latter model performed better. The variables screened by the models were closely related to the pathogenesis of diseases, pharmacokinetics and pharmacodynamics of MTX, and could be fully explained. Interestingly, the coagulation indicators filtered out by our models may indicate the new pathogenesis of JIA and the unelucidated mechanism of MTX. This model is simple, efficient and accurate, and can be easily generalized by clinicians and pharmacists to make early treatment decisions to patients of different ethnic groups.

Data Availability Statement

The datasets generated for this study are available on request to the Corresponding Author.

Ethics Statement

The studies involving human participants were reviewed and approved by drug clinical research ethics committee of Guangzhou women and children’s medical center. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. Written informed consent was Obtained From the Minor(S)’ Legal Guardian/Next of kin for the publication of any potentially identifiable images or data included in this article.

Author Contributions

Conceptualization: XM, YH, HL (16th author) and HZ. Data curation: XM, HWL, YC, FH, SZ and HL (13th author). Formal analysis: XC, HXL, LP, HL (16th author) and HZ. Funding acquisition: XM and MH. Investigation: XM, HWL, YC and FH. Methodology: XM, XC, JL, HXL and LP. Project administration: XM, MH, YH, HL (16th author) and HZ. Resources: HWL, SZ, PZ, YX and HZ. Software: XC and HXL. Supervision: XM. Validation: PZ and YX. Writing — original draft: XM and XC. Writing — review and editing: XM, XC, HWL, JL, FZ, YH, HYL and HZ.

Funding

This research was supported by grants from the National Natural Science Foundation of China (grant no. 81603203), the National Key Research and Development Program (grant no. 2017YFC0909303), the National Key Research and Development Program (Grant No. 2016YFC0905001), Health Commission of Guangdong Province (Grant No. A2016400), Guangdong Pharmaceutical Association Program (Grant No. 2015FS10 and 2015SW05), Guangzhou Institute of Pediatrics/Guangzhou Women and Children’s Medical Center (Grant No. YIP-2018-020).

Conflict of Interest

The 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.

Acknowledgments

We thank the physicians and nurses from the division of Immunology for their cooperation. Thanks for the help of Bing Zhu and Yinghua Li from the central lab, Xiaoqiong Gu from clinical biological resource bank and Xiaojun Cao from the department of science and education.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.01155/full#supplementary-material

References

Albers, H. M., Wessels, J. A., van der Straaten, R. J., Brinkman, D. M., Suijlekom-Smit, L. W., Kamphuis, S. S., et al. (2009). Time to treatment as an important factor for the response to methotrexate in juvenile idiopathic arthritis. Arthritis Rheum. 61, 46–51. doi: 10.1002/art.24087

PubMed Abstract | CrossRef Full Text | Google Scholar

Austin, P. C., Tu, J. V., Ho, J. E., Levy, D., Lee, D. S. (2013). Using methods from the data-mining and machine-learning literature for disease classification and prediction: a case study examining classification of heart failure subtypes. J. Clin. Epidemiol. 66, 398–407. doi: 10.1016/j.jclinepi.2012.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bedell, S. E., Bush, B. T. (1985). Erythrocyte sedimentation rate. From folklore to facts. Am. J. Med. 78, 1001–1009. doi: 10.1016/0002-9343(85)90224-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Browning, M., Kingslake, J., Dourish, C. T., Goodwin, G. M., Harmer, C. J., Dawson., G. R. (2019). Predicting treatment response to antidepressant medication using early changes in emotional processing. Eur. Neuropsychopharmacol. 29, 66–75. doi: 10.1016/j.euroneuro.2018.11.1102

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulatovic Calasan, M., Vastert, S. J., Scholman, R. C., Verweij, F., Klein, M., Wulffraat, N. M., et al. (2015). Methotrexate treatment affects effector but not regulatory T cells in juvenile idiopathic arthritis. Rheumatology (Oxford). 54, 1724–1734. doi: 10.1093/rheumatology/kev101

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulatovic, M., Heijstek, M. W., Van Dijkhuizen, E. H., Wulffraat, N. M., Pluijm, S. M., de Jonge., R. (2012). Prediction of clinical non-response to methotrexate treatment in juvenile idiopathic arthritis. Ann. Rheum. Dis. 71, 1484–1489. doi: 10.1136/annrheumdis-2011-200942

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabral, D., Katz, J. N., Weinblatt, M. E., Ting, G., Avorn, J., Solomon., D. H. (2005). Development and assessment of indicators of rheumatoid arthritis severity: results of a Delphi panel. Arthritis Rheum. 53, 61–66. doi: 10.1002/art.20925

PubMed Abstract | CrossRef Full Text | Google Scholar

Carney, D. H., Mann, R., Redin, W. R., Pernia, S. D., Berry, D., Heggers, J. P., et al. (1992). Enhancement of incisional wound healing and neovascularization in normal rats by thrombin and synthetic thrombin receptor-activating peptides. J. Clin. Invest. 89, 1469–1477. doi: 10.1172/JCI115737

PubMed Abstract | CrossRef Full Text | Google Scholar

Cid, M. C., Grant, D. S., Hoffman, G. S., Auerbach, R., Fauci, A. S., Kleinman., H. K. (1993). Identification of haptoglobin as an angiogenic factor in sera from patients with systemic vasculitis. J. Clin. Invest. 91, 977–985. doi: 10.1172/JCI116319

PubMed Abstract | CrossRef Full Text | Google Scholar

Cobb, J., Cule, E., Moncrieffe, H., Hinks, A., Ursu, S., Patrick, F., et al. (2014). Genome-wide data reveal novel genes for methotrexate response in a large cohort of juvenile idiopathic arthritis cases. Pharmacogenomics J. 14, 356–364. doi: 10.1038/tpj.2014.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Consolaro, A., Ruperto, N., Bazso, A., Pistorio, A., Magni-Manzoni, S., Filocamo, G., et al. (2009). Development and validation of a composite disease activity score for juvenile idiopathic arthritis. Arthritis Rheum. 61, 658–666. doi: 10.1002/art.24516

PubMed Abstract | CrossRef Full Text | Google Scholar

de Rotte, M. C., Bulatovic, M., Heijstek, M. W., Jansen, G., Heil, S. G., van Schaik, R. H., et al. (2012). ABCB1 and ABCC3 gene polymorphisms are associated with first-year response to methotrexate in juvenile idiopathic arthritis. J. Rheumatol. 39, 2032–2040. doi: 10.3899/jrheum.111593

PubMed Abstract | CrossRef Full Text | Google Scholar

Finnegan, S., Clarke, S., Gibson, D., McAllister, C., Rooney., M. (2011). Synovial membrane immunohistology in early untreated juvenile idiopathic arthritis: differences between clinical subgroups. Ann. Rheum. Dis. 70, 1842–1850. doi: 10.1136/ard.2010.148635

PubMed Abstract | CrossRef Full Text | Google Scholar

Foell, D., Wulffraat, N., Wedderburn, L. R., Wittkowski, H., Frosch, M., Gerss, J., et al. (2010). Methotrexate withdrawal at 6 vs 12 months in juvenile idiopathic arthritis in remission: a randomized clinical trial. JAMA 303, 1266–1273. doi: 10.1001/jama.2010.375

PubMed Abstract | CrossRef Full Text | Google Scholar

Franova, J., Fingerhutova, S., Kobrova, K., Srp, R., Nemcova, D., Hoza, J., et al. (2016). Methotrexate efficacy, but not its intolerance, is associated with the dose and route of administration. Pediatr. Rheumatol. Online J. 14, 36. doi: 10.1186/s12969-016-0099-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Giannini, E. H., Ruperto, N., Ravelli, A., Lovell, D. J., Felson, D. T., Martini, A. (1997). Preliminary definition of improvement in juvenile arthritis. Arthritis Rheum. 40, 1202–1209. doi: 10.1002/1529-0131(199707)40:7<1202::AID-ART3>3.0.CO;2-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Grober, J. S., Bowen, B. L., Ebling, H., Athey, B., Thompson, C. B., Fox, D. A., et al. (1993). Monocyte-endothelial adhesion in chronic rheumatoid arthritis. In situ detection of selectin and integrin-dependent interactions. J. Clin. Invest. 91, 2609–2619. doi: 10.1172/JCI116500

PubMed Abstract | CrossRef Full Text | Google Scholar

Hinks, A., Moncrieffe, H., Martin, P., Ursu, S., Lal, S., Kassoumeri, L., et al. (2011). Association of the 5-aminoimidazole-4-carboxamide ribonucleotide transformylase gene with response to methotrexate in juvenile idiopathic arthritis. Ann. Rheum. Dis. 70, 1395–1400. doi: 10.1136/ard.2010.146191

PubMed Abstract | CrossRef Full Text | Google Scholar

Hui, A. Y., McCarty, W. J., Masuda, K., Firestein, G. S., Sah, R. L. (2012). A systems biology approach to synovial joint lubrication in health, injury, and disease. Wiley Interdiscip. Rev. Syst. Biol. Med. 4, 15–37. doi: 10.1002/wsbm.157

PubMed Abstract | CrossRef Full Text | Google Scholar

Isaacs, J. D. (2007). T cell immunomodulation—the Holy Grail of therapeutic tolerance. Curr. Opin. Pharmacol. 7, 418–425. doi: 10.1016/j.coph.2007.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaskowski, T. D., Hill, H. R., Russo, K. L., Lakos, G., Szekanecz, Z., Teodorescu, M. (2010). Relationship between rheumatoid factor isotypes and IgG anti-cyclic citrullinated peptide antibodies. J. Rheumatol. 37, 1582–1588. doi: 10.3899/jrheum.091236

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnston, A., Gudjonsson, J. E., Sigmundsdottir, H., Ludviksson, B. R., Valdimarsson, H. (2005). The anti-inflammatory action of methotrexate is not mediated by lymphocyte apoptosis, but by the suppression of activation and adhesion molecules. Clin. Immunol. 114, 154–163. doi: 10.1016/j.clim.2004.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kruppa, J., Ziegler, A., Konig, I. R. (2012). Risk estimation and risk prediction using machine-learning methods. Hum. Genet. 131, 1639–1654. doi: 10.1007/s00439-012-1194-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Tu, Z., Deng, M., Sun, F., Chen, T. (2006). Diffusion kernel-based logistic regression models for protein function prediction. OMICS 10, 40–55. doi: 10.1089/omi.2006.10.40

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H. C., Yoon, S. B., Yang, S. M., Kim, W. H., Ryu, H. G., Jung, C. W., et al. (2018). Prediction of acute kidney injury after liver transplantation: machine learning approaches vs. logistic regression model. J. Clin. Med. 7, E428. doi: 10.3390/jcm7110428

PubMed Abstract | CrossRef Full Text | Google Scholar

Lezcano-Valverde, J. M., Salazar, F., Leon, L., Toledano, E., Jover, J. A., Fernandez-Gutierrez, B., et al. (2017). Development and validation of a multivariate predictive model for rheumatoid arthritis mortality using a machine learning approach. Sci. Rep. 7, 10189. doi: 10.1038/s41598-017-10558-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, H., Tsui, B. Y., Ni, H., Valentim, C. C. S., Baxter, S. L., Liu, G., et al. (2019). Evaluation and accurate diagnoses of pediatric diseases using artificial intelligence. Nat. Med. 25, 433–438. doi: 10.1038/s41591-018-0335-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, K. P., Cai, T., Gainer, V., Goryachev, S., Zeng-treitler, Q. Raychaudhuri, S., et al. (2010). Electronic medical records for discovery research in rheumatoid arthritis. Arthritis Care Res. (Hoboken). 62, 1120–1127. doi: 10.1002/acr.20184

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C., Karlson, E. W., Dligach, D., Ramirez, M. P., Miller, T. A., Mo, H., et al. (2015). Automatic identification of methotrexate-induced liver toxicity in patients with rheumatoid arthritis from the electronic medical record. J. Am. Med. Inform. Assoc. 22, e151–e161. doi: 10.1136/amiajnl-2014-002642

PubMed Abstract | CrossRef Full Text | Google Scholar

Martini, A., Ravelli, A., Avcin, T., Beresford, M. W., Burgos-Vargas, R., Cuttica, R., et al. (2019). Toward new classification criteria for juvenile idiopathic arthritis: first steps, pediatric rheumatology international trials organization international consensus. J. Rheumatol. 46, 190–197. doi: 10.3899/jrheum.180168

PubMed Abstract | CrossRef Full Text | Google Scholar

Motwani, M., Dey, D., Berman, D. S., Germano, G., Achenbach, S., Al-Mallah, M. H., et al. (2017). Machine learning for prediction of all-cause mortality in patients with suspected coronary artery disease: a 5-year multicentre prospective registry analysis. Eur. Heart J. 38, 500–507. doi: 10.1093/eurheartj/ehw188

PubMed Abstract | CrossRef Full Text | Google Scholar

Petty, R. E., Southwood, T. R., Manners, P., Baum, J., Glass, D. N., Goldenberg, J., et al. (2004). International League of Associations for Rheumatology classification of juvenile idiopathic arthritis: second revision, Edmonton, 2001. J. Rheumatol. 31, 390–392.

PubMed Abstract | Google Scholar

Ranganath, V. K., Yoon, J., Khanna, D., Park, G. S., Furst, D. E., Elashoff, D. A., et al. (2007). Comparison of composite measures of disease activity in an early seropositive rheumatoid arthritis cohort. Ann. Rheum. Dis. 66, 1633–1640. doi: 10.1136/ard.2006.065839

PubMed Abstract | CrossRef Full Text | Google Scholar

Roszkiewicz, J., Smolewska, E. (2017). In the pursuit of methotrexate treatment response biomarker in juvenile idiopathic arthritis—Are we getting closer to personalised medicine? Curr. Rheumatol. Rep. 19, 19. doi: 10.1007/s11926-017-0646-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruperto, N., Murray, K. J., Gerloni, V., Wulffraat, N., de Oliveira, S. K., Falcini, F., et al. (2004). A randomized trial of parenteral methotrexate comparing an intermediate dose with a higher dose in children with juvenile idiopathic arthritis who failed to respond to standard doses of methotrexate. Arthritis Rheum. 50, 2191–2201. doi: 10.1002/art.20288

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Nieuwenhove, E., Lagou, V., Van Eyck, L., Dooley, J., Bodenhofer, U., Roca, C., et al. (2019). Machine learning identifies an immunological pattern associated with multiple juvenile idiopathic arthritis subtypes. Ann. Rheum. Dis. doi: 10.1136/annrheumdis-2018-214354

CrossRef Full Text | Google Scholar

Yanagimachi, M., Naruto, T., Hara, T., Kikuchi, M., Hara, R., Miyamae, T., et al. (2011). Influence of polymorphisms within the methotrexate pathway genes on the toxicity and efficacy of methotrexate in patients with juvenile idiopathic arthritis. Br. J. Clin. Pharmacol. 71, 237–243. doi: 10.1111/j.1365-2125.2010.03814.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zajc Avramovic, M., Dolzan, V., Toplak, N., Accetto, M., Lusa, L., Avcin., T. (2017). Relationship Between Polymorphisms in Methotrexate Pathway Genes and Outcome of Methotrexate Treatment in a Cohort of 119 Patients with Juvenile Idiopathic Arthritis. J. Rheumatol. 44, 1216–1223. doi: 10.3899/jrheum.160950

PubMed Abstract | CrossRef Full Text | Google Scholar

Zborovskii, A. B., Sivordova, L. E., Derevianko, L. I., Zborovskaia, I. A., Zavodovskii, B. V. (1999). [Efficacy of D-penicillamine and methotrexate in the treatment of rheumatoid arthritis in relation to levels of circulating rheumatoid factors of different classes]. Ter. Arkh. 71, 60–63.

Google Scholar

Zhao, W., Fakhoury, M., Jacqz-Aigrain, E. (2010). Developmental pharmacogenetics of immunosuppressants in pediatric organ transplantation. Ther. Drug Monit. 32, 688–699 doi: 10.1097/FTD.0b013e3181f6502d

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: methotrexate, juvenile idiopathic arthritis, prediction model, machine learning, clinical response

Citation: Mo X, Chen X, Li H, Li J, Zeng F, Chen Y, He F, Zhang S, Li H, Pan L, Zeng P, Xie Y, Li H, Huang M, He Y, Liang H and Zeng H (2019) Early and Accurate Prediction of Clinical Response to Methotrexate Treatment in Juvenile Idiopathic Arthritis Using Machine Learning. Front. Pharmacol. 10:1155. doi: 10.3389/fphar.2019.01155

Received: 28 June 2019; Accepted: 09 September 2019;
Published: 07 October 2019.

Edited by:

James Cheng-Chung Wei, Chung Shan Medical University, Taiwan

Reviewed by:

Song Guo Zheng, The Ohio State University, United States
Adina Turcu-Stiolica, University of Medicine and Pharmacy of Craiova, Romania

Copyright © 2019 Mo, Chen, Li, Li, Zeng, Chen, He, Zhang, Li, Pan, Zeng, Xie, Li, Huang, He, Liang and Zeng. 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: Yanling He, 734242355@qq.com; Huiying Liang, lianghuiying@hotmail.com; Huasong Zeng, huasongxuqing@163.com

These authors have contributed equally to this work