Skip to main content
  • Research article
  • Open access
  • Published:

Baseline metabolic profiles of early rheumatoid arthritis patients achieving sustained drug-free remission after initiating treat-to-target tocilizumab, methotrexate, or the combination: insights from systems biology



We previously identified, in newly diagnosed rheumatoid arthritis (RA) patients, networks of co-expressed genes and proteomic biomarkers associated with achieving sustained drug-free remission (sDFR) after treatment with tocilizumab- or methotrexate-based strategies. The aim of this study was to identify, within the same patients, metabolic pathways important for achieving sDFR and to subsequently study the complex interactions between different components of the biological system and how these interactions might affect the therapeutic response in early RA.


Serum samples were analyzed of 60 patients who participated in the U-Act-Early trial ( number NCT01034137) and initiated treatment with methotrexate, tocilizumab, or the combination and who were thereafter able to achieve sDFR (n = 37); as controls, patients were selected who never achieved a drug-free status (n = 23). Metabolomic measurements were performed using mass spectrometry on oxidative stress, amine, and oxylipin platforms covering various compounds. Partial least square discriminant analyses (PLSDA) were performed to identify, per strategy arm, relevant metabolites of which the biological pathways were studied. In addition, integrative analyses were performed correlating the previously identified transcripts and proteins with the relevant metabolites.


In the tocilizumab plus methotrexate, tocilizumab, and methotrexate strategy, respectively, 19, 13, and 12 relevant metabolites were found, which were subsequently used for pathway analyses. The most significant pathway in the tocilizumab plus methotrexate strategy was “histidine metabolism” (p < 0.001); in the tocilizumab strategy it was “arachidonic acid metabolism” (p = 0.018); and in the methotrexate strategy it was “arginine and proline metabolism” (p = 0.022). These pathways have treatment-specific drug interactions with metabolites affecting either the signaling of interleukin-6, which is inhibited by tocilizumab, or affecting protein synthesis from amino acids, which is inhibited by methotrexate.


In early RA patients treated-to-target with a tocilizumab- or methotrexate-based strategy, several metabolites were found to be associated with achieving sDFR. In line with our previous observations, by analyzing relevant transcripts and proteins within the same patients, the metabolic profiles were found to be different between the strategy arms. Our metabolic analysis further supports the hypothesis that achieving sDFR is not only dependent on predisposing biomarkers, but also on the specific treatment that has been initiated.

Trial registration, NCT01034137. Registered on January 2010


Rheumatoid arthritis (RA) is a systemic disease characterized by inflammation and damage of the affected joints; although the cause is not known, both genetic as well as environmental factors are reported to be associated with the condition [1,2,3,4]. Initiating treatment, aiming for sustained remission or low disease activity early in the course of the disease, is important to preserve physical function and improve long-term prognosis [5,6,7]. Biological disease-modifying anti-rheumatic drugs (DMARDs) are to date mainly used as the second line of therapy in the management of early RA, although several studies showed their superior efficacy over traditional DMARDs in reducing disease activity and halting joint damage [8, 9]. Starting biological therapy in newly diagnosed RA patients as standard care, however, still remains highly controversial considering costs and unnecessary exposure to adverse events as a proportion of patients will be over-treated when using such an approach. Therefore, predictors are not only needed for treatment response to the currently recommended conventional DMARDs, but also for identifying patients for whom it would be favorable to initiate, as first therapy, a step-down biological-based strategy (i.e. tapering and finally discontinuing treatment) as achieving remission in the early stage of the disease improves the long-term clinical outcome.

Recent developments in “omics” technologies—such as genomics, transcriptomics, proteomics, and metabolomics—made it feasible to measure a broader spectrum of disease biomarkers for prediction of disease progression and development of personalized treatment strategies in RA [10]. Metabolomics is the non-targeted study of small-molecule metabolites and has become of increased interest in recent years due to the development and accessibility of new high-throughput technologies, including nuclear magnetic resonance spectroscopy and mass spectrometry (MS) [11]. Metabolites provide, under a given set of conditions, detailed information on cellular processes that are indicative for the disease state and are considered as the final downstream product of gene expression [12]. Especially in RA, metabolites are of particular interest as widespread cytokine-mediated inflammatory processes alter the cellular metabolism, when macrophages and lymphocytes become activated [13]. The role of these compounds in biomarker discovery has also been demonstrated previously, suggesting that metabolic analysis is potentially valuable in identifying markers for treatment response in patients with RA [14,15,16,17].

The aim of this study was to identify relevant metabolites and important metabolomic pathways associated with achieving sustained drug-free remission (sDFR) after a treat-to-target tocilizumab- or methotrexate-based strategy initiated in DMARD-naïve early RA patients. We previously identified, within the same patients, networks of co-expressed genes [18] and several inflammatory proteins [19] associated with sDFR and now, in the present study, by revealing metabolic biomarkers, exploring the systems biology of these patients in more detail by also integrating the findings of our previous studies.


Patient selection

From the U-Act-Early strategy trial, patients were selected who achieved sDFR, defined as being drug-free for ≥ 3 months until end of the study, after initiating tocilizumab, step-up methotrexate, or tocilizumab plus methotrexate therapy. As controls, we selected patients who never achieved a drug-free status at any time point during the study. A detailed description of the study design has been reported previously [20]. Briefly, DMARD-naïve patients with early RA were randomized (1:1:1) to one of the three strategy arms and treated to the target of sustained remission, defined as disease activity score assessing 28 joints (DAS28) < 2.6 with ≤ 4 swollen joints for ≥ 24 weeks. Tocilizumab was administered intravenously every four weeks at a dose of 8 mg/kg with a maximum of 800 mg. Methotrexate (oral) was given every week with a starting dose of 10 mg and was increased to 30 mg (or maximum tolerable dose) with steps of 5 mg every four weeks until remission was reached. If the treatment target was achieved, medication was tapered stepwise and finally discontinued, if remission persisted.

Metabolomic platforms

Baseline serum samples were measured on oxidative stress, amines, and oxylipins MS platforms, which have been applied previously and are validated [21,22,23]. The oxidative stress platform covers various isoprostane classes, signaling lipids from the sphingosine and sphinganine classes and their phosphorylated forms, as well as three classes of lysophosphatidic acids: lysophosphatidic acids, alkyl-lysophosphatidic acids, and cyclic-lysophosphatidic acids (all in the chain length species range of C14–C22). The amine platform covers amino acids and biogenic amines and the oxylipin platform covers classical and non-classical eicosanoids from different polyunsaturated fatty acids. In total, 263 metabolites were measured for each sample on the different platforms: 57 signaling lipid mediators, 128 oxylipins, and 78 amines. Serum samples were thawed on ice and vortexed before preparation procedures following inner standard protocols [22, 23]; extra samples were pooled for internal quality control (QC). For each platform, QC samples were added and the relative standard deviation (RSD) per metabolite was calculated. Only those metabolites that complied with the acceptance criteria (RSD < 15% for amines; RSD < 30% for oxidative stress and oxylipins) were selected for further data analysis. Additional information about the metabolite profiling on the different platforms is provided in Additional file 1. All metabolite analyses were performed by the Biomedical Metabolomics Facility Leiden department at Leiden University.

Data pre-processing

For peak determination and integration, signaling lipid mediators profiled by the oxidative stress platform were pre-processed by LabSolutions (Shimadzu, Version 5.65); peak-picking of oxylipins was performed with Agilent MassHunter Quantitative Analysis software (Agilent, Version B.05.01) and amines with MultiQuant Software for Quantitative Analysis (AB SCIEX, Version 3.0.2). For all metabolites, raw data correction was accomplished using selected internal standards by calculating the ratio of peak area of the target compound to the peak area of assigned internal standard from which a response ratio for each analyte was obtained. QC samples were used for evaluating the quality of the targeted compounds according to the in-house written protocol and the data were hereafter ready to be used for statistical analyses.

Statistical analyses

Baseline clinical characteristics are described as mean (standard deviation [SD]), median (interquartile range [IQR]) or as proportions (%); between-group differences (sDFR versus controls) were tested within each strategy arm using independent t test, Mann–Whitney U test, or Pearson χ2 test, respectively. A linear mixed model with a random intercept and baseline DAS28, week of visit, and group (sDFR versus controls) as fixed effects was built to evaluate, within the strategy arms, differences in disease activity over time. As metabolite concentrations are influenced by a variety of factors, we performed principal component analyses (PCA) to identify possible confounders. The following parameters were considered: age; body mass index, gender, ethnicity, disease duration, smoking, alcohol consumption, seropositivity for rheumatoid factor (RF) or cyclic citrullinated peptide (CCP), erythrocyte sedimentation rate (ESR), and C-reactive protein (CRP). Thereafter, supervised partial least square discriminant analyses (PLSDA) were performed for each class (lipids, amines, and oxylipins) to identify relevant metabolites within each strategy arm. Several multivariate discrimination techniques currently exist but the main advantage of PLSDA is the handling of collinearity and noisy data (i.e. more observations than samples), both common in metabolomics experiments [24]. Data were first normalized (natural log-transformed) and then standardized (z-score) to ensure that all metabolite scores are comparable by giving them equal weight. The variable importance on projection (VIP) was used for metabolite selection; this measure accumulates the importance of each variable, whereas a higher VIP score shows that it is more relevant to predict the outcome [25]. Metabolites with VIP > 1 in the first component were considered important as the squared sum of all VIP values is equal to “1”, i.e. the average VIP. Thereafter, the Mann–Whitney U test (sDFR versus controls) was performed within these selected metabolites to identify those who are most relevant (p < 0.10) within each strategy arm (not corrected for multiple testing), which were subsequently used for pathway analysis in the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Pathways were considered relevant when p ≤ 0.05. In addition, integrative analyses between the previously identified transcripts [18] and proteins [19] and relevant metabolites were performed by calculating statistically significant (p ≤ 0.05) transcript–protein and protein–metabolite correlations (Pearson correlation coefficients [PCC]). These network analyses were visualized using VisANT 5.0 software [26, 27]. All other analyses were performed using the web-based tool MetaboAnalyst version 4.0 [28] and R version 3.4.3 (R Foundation for Statistical Computing, Vienna, Austria).


Serum samples were analyzed for 60 patients (tocilizumab plus methotrexate arm: n = 14 sDFR, n = 5 controls; tocilizumab arm: n = 13 sDFR, n = 11 controls; methotrexate arm: n = 10 sDFR, n = 7 controls) and their baseline characteristics are summarized in Table 1. The mean (SD) age of all patients was 53 (14) years with a median (IQR) symptom duration of 23 (18–40) days; the majority was seropositive for rheumatoid factor (60%) or cyclic citrullinated peptide (60%). At baseline, the mean (SD) DAS28 of these patients was 4.9 (1.1) with a median (IQR) ESR of 20 (11–32) mm/1thh and median (IQR) CRP of 9 (3–18) mg/L. No significant differences were noted in clinical characteristics at baseline between the groups (sDFR versus controls) within the strategy arms (p ≥ 0.07). The mean (standard error) DAS28 scores over time of the sDFR and control groups are shown in Fig. 1. In the longitudinal analysis, significant lower DAS28 scores were found in the sDFR group within the tocilizumab plus methotrexate (mean – 1.18, 95% confidence interval [CI] – 0.87, – 1.50; p < 0.001), tocilizumab (mean – 0.87, 95% CI – 0.62, – 1.12; p < 0.001), and methotrexate (mean − 0.43, 95% CI -0.16, − 0.69; p = 0.009) arms, when compared to the control group.

Table 1 Baseline characteristics of the patients included in the analyses
Fig. 1
figure 1

Mean (SE) DAS28-ESR over time in those achieving sDFR (continuous line) vs controls (dotted line) within the three strategy arms. DAS28 disease activity score assessing 28 joints, ESR erythrocyte sedimentation rate, MTX methotrexate, sDFR sustained drug-free remission, SE standard error, TCZ tocilizumab

Metabolite biomarkers for sDFR

PCA revealed no clear confounders (data not shown) and therefore metabolite concentrations were not corrected for clinical characteristics in further analyses. PLSDA identified 35 metabolites (15 signaling lipid mediators, 14 oxylipins, 6 amines) in the tocilizumab plus methotrexate arm, 33 metabolites (9 signaling lipid mediators, 15 oxylipins, 9 amines) in the tocilizumab arm, and 33 metabolites (11 signaling lipid mediators, 13 oxylipins, 9 amines) in the methotrexate arm. Of these metabolites, 19, 13, and 12, respectively, were subsequently selected for further pathway analyses (Table 2). When comparing the metabolites between the strategy arms, one (LPA c16:1) showed overlap in the tocilizumab plus methotrexate and tocilizumab arm; two (12,13-DiHODE and LPA c16:0) in the tocilizumab plus methotrexate and methotrexate strategy, and one (3-Methylhistidine) in the tocilizumab and methotrexate strategy.

Table 2 Identified metabolites associated with achieving sustained drug-free remission in the three strategy arms

Pathway analyses

Pathway overviews for each strategy arm are shown in Additional file 2. The three most relevant KEGG pathways in the tocilizumab plus methotrexate arm were “histidine metabolism” (p < 0.001), “sphingolipid metabolism” (p = 0.004), and “arachidonic acid metabolism” (p = 0.037). Within the “histidine metabolism” pathway, a significant lower concentration of histamine (p = 0.002) was found in the sDFR group when compared to controls (Fig. 2). In the “arachidonic acid metabolism” pathway, production of PGD2 and 12,13-DiHODE, 9,10,13-TriHOME, and 9,12,13-TriHOME through lipoxygenase was higher in those who achieved sDFR. Within the “sphingolipid metabolism” pathway, significantly lower levels of sphinganine (Spha c18:0, p = 0.007) and sphingosine (Sph c18:1, p = 0.012) were found in the sDFR (versus controls) group, which are both related to ceramide generation. In the tocilizumab arm, three most relevant pathways were identified: “arachidonic acid metabolism” (p = 0.018); “lysine degradation” (p = 0.023); and “cysteine and methionine metabolism” (p = 0.030, Fig. 2). In the “arachidonic acid metabolism” pathway, higher concentrations of prostaglandin E2 (PGE2), 8-isoprostaglandin E2 (8-iso-PGE2), prostaglandin A2 (PGA2), 8-isoprostaglandin A2 (8-iso-PGA2), 8,9-DiHETrE, and 5,6-DiHETrE (p ≤ 0.08) were found in the sDFR group, compared to controls, implying a more active role of prostaglandins and isoprostanes, which are critical signaling molecules in various inflammatory disease including RA [29,30,31]. Furthermore, in the “lysine degradation” pathway, significantly higher concentrations were found in the sDFR group (versus controls) of L-pipecolic acid (p = 0.026) and in the “cysteine and methionine metabolism” pathway, slightly but not statistically significantly lower concentrations were found in the sDFR group for cystathionine (p = 0.052) and homocysteine (p = 0.09). Most metabolites in the methotrexate arm were associated with the “arginine and proline” and “histidine metabolism” pathways (p = 0.022 and p = 0.025, respectively). When compared to controls, lower oxylipin levels of 14,15-DiHETE, 19,20-DiHDPA, and 12,13-DiHODE were found in the sDFR group (p ≤ 0.10), which indicates fewer active cytochromes P450 (CYP450) and lipoxygenase-based fatty acids metabolism, while the observed increased L-proline (p = 0.040) and L-arginine (p = 0.08) levels suggest a more active role of this pathway (Fig. 2). In the “histidine metabolism” pathway, decreased levels of 3-methylhistidine (p = 0.06) and higher levels of anserine (p = 0.06) were observed in the sDFR group (versus controls). Other changes of amine levels included higher concentrations of lysine (p = 0.032) and lower concentrations of 5-hydroxy-tryptophan (p = 0.08). Another significant pathway in the methotrexate arm was “aminoacyl-tRNA biosynthesis” (p = 0.027).

Fig. 2
figure 2

Pathway analysis within the identified metabolites in the (a) tocilizumab plus methotrexate, (b) tocilizumab, and (c) methotrexate strategy arms. Metabolites depicted in red nodes have, on average, lower concentration in the sDFR group compared to controls; those depicted in green nodes have a higher concentration. *p ≤ 0.10, **p ≤ 0.05, ***p ≤ 0.01. Metabolites not included in the top three most relevant pathways are not displayed

Systems biology: from transcripts to proteins and metabolites

Figure 3 shows significant transcript–protein and protein–metabolite correlations within the three strategy arms. The average PCC between the biomarkers in the tocilizumab plus methotrexate arm was 0.54; these were 0.48 and 0.57 in the tocilizumab and methotrexate arm, respectively. Biomarkers showing > 10 connections within the networks were considered most important (i.e. signature biomarker) as they show the highest connectivity and therefore contribute most to the pathway analyses. The signature biomarkers in the tocilizumab plus methotrexate arm were the protein chemokine (C-C motif) ligand 5 (CCL5), as it was significantly correlated to seven metabolites and four transcripts, and the protein tumor necrosis factor receptor 1 (TNF-R1), being significantly correlated to three metabolites and 14 transcripts. In the tocilizumab arm, the proteins TIMP metallopeptidase inhibitor 1 (TIMP-1, 13 with transcripts) and thyroid peroxidase (TPO, 12 with transcripts) showed the highest number of correlations and in the methotrexate arm the protein granulocyte colony-stimulating factor (G-CSF, 4 with metabolites; 26 with transcripts).

Fig. 3
figure 3

Network correlation between transcriptomic (gray nodes), proteomic (orange nodes), and metabolomic (blue nodes) biomarkers in the (a) tocilizumab plus methotrexate, (b) tocilizumab, and (c) methotrexate strategy arms. Only significant transcriptomic–proteomic and proteomic–metabolomic correlations are displayed


We identified several small-molecule metabolites, by using high-throughput MS, associated with achieving sDFR after treatment with tocilizumab- or methotrexate-based strategies in newly diagnosed RA patients. In line with our previous observations, by measuring transcripts and proteins within the same patients, different metabolic profiles were found between the treatment strategies, further supporting the hypothesis that achieving sDFR is likely dependent on pre-treatment concentrations of specific biomarkers as no differences in clinical characteristics could be found. Although we did find different metabolic pathways between the treatment strategies when using the identified metabolites, the pathways within each strategy arm were found to be specific for the respective treatment, which shows the possibility of selecting biomarkers for prediction of a good treatment-specific response.

An important metabolic pathway within the tocilizumab plus methotrexate strategy was “sphingolipid metabolism,” in which ceramide synthases is closely associated with the consumption of both sphingosine, a lipid signaling molecule stimulating several cellular processes important in RA, such as cell growth, differentiation, and migration [32], and its derivative sphinganine [33]. These metabolites have been found to induce IL-6 production and may thus influence the treatment response to tocilizumab, a humanized monoclonal antibody against the IL-6 receptor. In the present analysis, lower baseline levels of sphingosine and sphinganine were found in the sDFR group of the tocilizumab plus methotrexate strategy arm when compared to controls although the opposite was expected (i.e. less inhibition resulting in higher disease activity); a possible explanation for these findings could be that patients with higher levels of IL-6 are more likely to respond to tocilizumab. To test this hypothesis, we evaluated within the same patients the absolute IL-6 concentrations in baseline (i.e. pre-treatment) serum, as was measured previously using Luminex® [19], but found no statistically significant difference between those who would later achieve sDFR and the controls (mean 65 versus 70 pg/mL, respectively). This result confirms the findings of others showing that baseline levels of IL-6 are not predictive of clinical outcomes of tocilizumab treatment [34]. Further studies are required to elucidate on how metabolite levels are altered and how these pathophysiological changes eventually affect the response to specific therapies. Ceramide synthase 6 (CerS6), an enzyme also important in sphingolipid biosynthesis, is mediated by methotrexate [35], which might indicate a role of sphingolipids in the treatment response to a tocilizumab plus methotrexate-based strategy in early RA.

In the tocilizumab strategy arm, several involved metabolic pathways were found to be associated with signaling of the IL-6 protein whereas increased levels of prostaglandins and isoprostanes were observed in the sDFR group when compared to controls. These metabolites are involved in the metabolism of arachidonic acid, a polyunsaturated fatty acid reported to be a key intermediate promoting inflammation. Prostaglandins are known to have a stimulating effect on IL-6 and could thus influence the response to tocilizumab therapy [36]. Other important metabolites involved in IL-6 signaling are L-pipecolic acid, which is produced during the degradation of lysine, and cysteine, which is synthesized from serine. Both lysine and serine are essential amino acids in humans as they are being used in the biosynthesis of proteins whereas lysine has been reported to downregulate the release of IL-6 [37].

In the methotrexate strategy, important metabolites were proline and arginine for which higher concentrations were observed in those achieving sDFR when compared to controls. Arginine, apart from its role in protein synthesis, serves as the precursor of proline and glutamate; glutamates are involved in the generation of glutamine, which stimulates purine and pyrimidine formation that is required for cell proliferation. These organic compounds thus seem to play a crucial role in the direct treatment response to methotrexate as this drug antagonizes folic acid [38], which is important for purine and pyrimidine formation, and inhibits via this pathway the synthesis of nucleic acids and subsequently protein synthesis [39]. One of the other important amino acids for protein synthesis is histidine, of which the post-translational modified product, 3-methylhistidine, was found to have a decreased metabolism in the present study in those achieving sDFR in the methotrexate strategy. Histidine is also involved in the local immune response as it is a precursor of histamine, a compound that induces permeability of capillaries allowing, for example, leukocytes and pro-inflammatory molecules to elicit an immune response. These amino acids thus might directly affect the treatment response to methotrexate therapy.

There are some limitations to this study. First, although not uncommon in metabolomic studies, the number of samples measured was relatively small, enhancing the likelihood of false-negative findings (i.e. type II error). To minimize this risk, we used for detecting relevant metabolites analyses suitable for handling such datasets consisting of more markers than samples. Second, before the samples were measured, serum was pre-processed according to usual guidelines, which differed for each platform, and compounds were, after analyses, corrected using internal standards. These variety of factors could potentially impair replicating or external validation of findings. Third, for defining remission, and thus tapering medication, we used DAS28 criteria which is highly dependent of acute phase response and might not always reflect an inflammation-free state of the patient. To minimize this risk, no more than four swollen joints were allowed during the remission period, enhancing the likelihood that patients were also clinically in remission. Nevertheless, American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) Boolean-based remission criteria might have been a more reliable tool for assessing disease activity as it is more stringent in assessing inflammation [40].


We have identified several relevant metabolites in baseline serum related to achieving sDFR after treatment with tocilizumab- or methotrexate-based strategies in DMARD-naïve RA patients. In line with our previous work on the analyses of transcripts and proteins, performed within the same patients, the identified metabolic pathways were shown to be specific for the treatment that was initiated. These results might provide further insight into the role of predisposing biomarkers for eventually achieving sDFR in early RA. Signature metabolite biomarkers have been identified which could potentially serve as key prognostic factors for developing personalized care but need to be validated in large replication studies. Further studies are also warranted to elucidate on the drug metabolism in those patients with refractory disease to, by initiating other therapies based on pharmacogenomics, achieve better treatment outcomes.



8-isoprostaglandin E2


American College of Rheumatology


Chemokine (c-c motif) ligand 5


Cyclic citrullinated peptide


Ceramide synthase 6


Confidence interval


C-reactive protein


Cytochromes P450


Disease activity score assessing 28 joints


Disease-modifying anti-rheumatic drugs


Erythrocyte sedimentation rate


European League Against Rheumatism


Granulocyte colony-stimulating factor




Interquartile range


Kyoto Encyclopedia of Genes and Genomes


Mass spectrometry


Principal component analyses


Prostaglandin A2


Prostaglandin E2


Partial least square discriminant analyses


Quality control


Rheumatoid arthritis


Rheumatoid factor


Standard deviation


Sustained drug-free remission


TIMP metallopeptidase inhibitor 1


Tumour necrosis factor receptor 1


Thyroid peroxidase


Variable importance on projection


  1. Pedersen M, Jacobsen S, Klarlund M, Pedersen BV, Wiik A, Wohlfahrt J, Frisch M. Environmental risk factors differ between rheumatoid arthritis with and without auto-antibodies against cyclic citrullinated peptides. Arthritis Res Ther. 2006;8(4):R133.

    Article  Google Scholar 

  2. Viatte S, Plant D, Raychaudhuri S. Genetics and epigenetics of rheumatoid arthritis. Nat Rev Rheumatol. 2013;9(3):141–53.

    Article  CAS  Google Scholar 

  3. Stolt P, Bengtsson C, Nordmark B, Lindblad S, Lundberg I, Klareskog L, Alfredsson L, EIRA study group. Quantification of the influence of cigarette smoking on rheumatoid arthritis: results from a population based case-control study, using incident cases. Ann Rheum Dis. 2003;62(9):835–41.

    Article  CAS  Google Scholar 

  4. Liao KP, Alfredsson L, Karlson EW. Environmental influences on risk for rheumatoid arthritis. Curr Opin Rheumatol. 2009;21(3):279–83.

    Article  Google Scholar 

  5. van Nies JA, Tsonaka R, Gaujoux-Viala C, Fautrel B, van der Helm-van Mil AH. Evaluating relationships between symptom duration and persistence of rheumatoid arthritis: does a window of opportunity exist? Results on the Leiden early arthritis clinic and ESPOIR cohorts. Ann Rheum Dis. 2015;74(5):806–12.

    Article  CAS  Google Scholar 

  6. Smolen JS, Landewe R, Bijlsma J, Burmester G, Chatzidionysiou K, Dougados M, Nam J, Ramiro S, Voshaar M, van Vollenhoven R, et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2016 update. Ann Rheum Dis. 2017;76(6):960–77.

    Article  Google Scholar 

  7. van der Linden MP, le Cessie S, Raza K, van der Woude D, Knevel R, Huizinga TW, van der Helm-van Mil AH. Long-term impact of delay in assessment of patients with early arthritis. Arthritis Rheum. 2010;62(12):3537–46.

    Article  Google Scholar 

  8. Alfonso-Cristancho R, Armstrong N, Arjunji R, Riemsma R, Worthy G, Ganguly R, Kleijnen J. Comparative effectiveness of biologics for the management of rheumatoid arthritis: systematic review and network meta-analysis. Clin Rheumatol. 2017;36(1):25–34.

    Article  Google Scholar 

  9. Hazlewood GS, Barnabe C, Tomlinson G, Marshall D, Devoe DJ, Bombardier C. Methotrexate monotherapy and methotrexate combination therapy with traditional and biologic disease modifying anti-rheumatic drugs for rheumatoid arthritis: A network meta-analysis. Cochrane Database Syst Rev. 2016;8:CD010227.

    Google Scholar 

  10. Priori R, Scrivo R, Brandt J, Valerio M, Casadei L, Valesini G, Manetti C. Metabolomics in rheumatic diseases: the potential of an emerging methodology for improved patient diagnosis, prognosis, and treatment efficacy. Autoimmun Rev. 2013;12(10):1022–30.

    Article  CAS  Google Scholar 

  11. Kang J, Zhu L, Lu J, Zhang X. Application of metabolomics in autoimmune diseases: insight into biomarkers and pathology. J Neuroimmunol. 2015;279:25–32.

    Article  CAS  Google Scholar 

  12. Morgan RP, Louise KC. Omic technologies: genomics, transcriptomics, proteomics and metabolomics. Obstet Gynaecol. 2011;13(3):189–95.

    Google Scholar 

  13. Chimenti MS, Triggianese P, Conigliaro P, Candi E, Melino G, Perricone R. The interplay between inflammation and metabolism in rheumatoid arthritis. Cell Death Dis. 2015;6:e1887.

    Article  CAS  Google Scholar 

  14. Cuppen BV, Fu J, van Wietmarschen HA, Harms AC, Koval S, Marijnissen AC, Peeters JJ, Bijlsma JW, Tekstra J, van Laar JM, et al. Exploring the inflammatory metabolomic profile to predict response to TNF-alpha inhibitors in rheumatoid arthritis. PLoS One. 2016;11(9):e0163087.

    Article  Google Scholar 

  15. Priori R, Casadei L, Valerio M, Scrivo R, Valesini G, Manetti C. 1H-NMR-based metabolomic study for identifying serum profiles associated with the response to etanercept in patients with rheumatoid arthritis. PLoS One. 2015;10(11):e0138537.

    Article  Google Scholar 

  16. Chimenti MS, Tucci P, Candi E, Perricone R, Melino G, Willis AE. Metabolic profiling of human CD4+ cells following treatment with methotrexate and anti-TNF-alpha infliximab. Cell Cycle. 2013;12(18):3025–36.

    Article  CAS  Google Scholar 

  17. Kapoor SR, Filer A, Fitzpatrick MA, Fisher BA, Taylor PC, Buckley CD, McInnes IB, Raza K, Young SP. Metabolic profiling predicts response to anti-tumor necrosis factor alpha therapy in patients with rheumatoid arthritis. Arthritis Rheum. 2013;65(6):1448–56.

    Article  CAS  Google Scholar 

  18. Teitsma XM, Jacobs JWG, Mokry M, Borm MEA, Petho-Schramm A, van Laar JM, Bijlsma JWJ, Lafeber FPJ. Identification of differential co-expressed gene networks in early rheumatoid arthritis achieving sustained drug-free remission after treatment with a tocilizumab-based or methotrexate-based strategy. Arthritis Res Ther. 2017;19(1):170.

    Article  Google Scholar 

  19. Teitsma XM, Jacobs JWG, Concepcion AN, Petho-Schramm A, Borm MEA, van Laar JM, Bijlsma JWJ, Lafeber F. Explorative analyses of protein biomarkers in patients with early rheumatoid arthritis achieving sustained drug-free remission after treatment with tocilizumab- or methotrexate-based strategies: from transcriptomics to proteomics. Clin Exp Rheumatol. 2018; [Epub ahead of print].

  20. Bijlsma JW, Welsing PM, Woodworth TG, Middelink LM, Petho-Schramm A, Bernasconi C, Borm ME, Wortel CH, ter Borg EJ, Jahangier ZN, et al. Early rheumatoid arthritis treated with tocilizumab, methotrexate, or their combination (U-Act-Early): a multicentre, randomised, double-blind, double-dummy, strategy trial. Lancet. 2016;388(10042):343–55.

    Article  CAS  Google Scholar 

  21. Schoeman JC, Harms AC, van Weeghel M, Berger R, Vreeken RJ, Hankemeier T. Development and application of a UHPLC-MS/MS metabolomics based comprehensive systemic and tissue-specific screening method for inflammatory, oxidative and nitrosative stress. Anal Bioanal Chem. 2018;410(10):2551–68.

    Article  CAS  Google Scholar 

  22. Strassburg K, Huijbrechts AM, Kortekaas KA, Lindeman JH, Pedersen TL, Dane A, Berger R, Brenkman A, Hankemeier T, van Duynhoven J, et al. Quantitative profiling of oxylipins through comprehensive LC-MS/MS analysis: application in cardiac surgery. Anal Bioanal Chem. 2012;404(5):1413–26.

    Article  CAS  Google Scholar 

  23. Noga MJ, Dane A, Shi S, Attali A, van Aken H, Suidgeest E, Tuinstra T, Muilwijk B, Coulier L, Luider T, et al. Metabolomics of cerebrospinal fluid reveals changes in the central nervous system metabolism in a rat model of multiple sclerosis. Metabolomics. 2012;8(2):253–63.

    Article  CAS  Google Scholar 

  24. Gromski PS, Muhamadali H, Ellis DI, Xu Y, Correa E, Turner ML, Goodacre R. A tutorial review: Metabolomics and partial least squares-discriminant analysis--a marriage of convenience or a shotgun wedding. Anal Chim Acta. 2015;879:10–23.

    Article  CAS  Google Scholar 

  25. Mehmood T, Hovde Liland K, Snipen L, Saebo S. A review of variable selection methods in Partial Least Squares Regression. Chemom Intell Lab Syst. 2012;118:62–9.

    Article  CAS  Google Scholar 

  26. Hu Z, Mellor J, Wu J, DeLisi C. VisANT: an online visualization and analysis tool for biological interaction data. BMC Bioinformatics. 2004;5:17.

    Article  Google Scholar 

  27. Hu Z, Snitkin ES, DeLisi C. VisANT: an integrative framework for networks in systems biology. Brief Bioinform. 2008;9(4):317–25.

    Article  CAS  Google Scholar 

  28. Chong J, Soufan O, Li C, Caraus I, Li S, Bourque G, Wishart DS, Xia J. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46(W1):W486–94.

    Article  Google Scholar 

  29. Basu S. Bioactive eicosanoids: role of prostaglandin F(2alpha) and F(2)-isoprostanes in inflammation and oxidative stress related pathology. Mol Cell. 2010;30(5):383–91.

    Article  CAS  Google Scholar 

  30. Hata AN, Breyer RM. Pharmacology and signaling of prostaglandin receptors: multiple roles in inflammation and immune modulation. Pharmacol Ther. 2004;103(2):147–66.

    Article  CAS  Google Scholar 

  31. Fattahi MJ, Mirshafiey A. Prostaglandins and rheumatoid arthritis. Arthritis. 2012;2012:239310.

    Article  Google Scholar 

  32. Hu PF, Chen Y, Cai PF, Jiang LF, Wu LD. Sphingosine-1-phosphate: a potential therapeutic target for rheumatoid arthritis. Mol Biol Rep. 2011;38(6):4225–30.

    Article  CAS  Google Scholar 

  33. Mullen TD, Hannun YA, Obeid LM. Ceramide synthases at the centre of sphingolipid metabolism and biology. Biochem J. 2012;441(3):789–802.

    Article  CAS  Google Scholar 

  34. Wang J, Devenport J, Low JM, Yu D, Hitraya E. Relationship between baseline and early changes in C-reactive protein and interleukin-6 levels and clinical response to tocilizumab in rheumatoid arthritis. Arthritis Care Res. 2016;68(6):882–5.

    Article  CAS  Google Scholar 

  35. Fekry B, Esmaeilniakooshkghazi A, Krupenko SA, Krupenko NI. Ceramide synthase 6 is a novel target of methotrexate mediating its antiproliferative effect in a p53-dependent manner. PLoS One. 2016;11(1):e0146618.

    Article  Google Scholar 

  36. Li P, Shan JX, Chen XH, Zhang D, Su LP, Huang XY, Yu BQ, Zhi QM, Li CL, Wang YQ, et al. Epigenetic silencing of microRNA-149 in cancer-associated fibroblasts mediates prostaglandin E2/interleukin-6 signaling in the tumor microenvironment. Cell Res. 2015;25(5):588–603.

    Article  CAS  Google Scholar 

  37. Al-Malki AL. Suppression of acute pancreatitis by L-lysine in mice. BMC Complement Altern Med. 2015;15:193.

    Article  Google Scholar 

  38. Wessels JA, Huizinga TW, Guchelaar HJ. Recent insights in the pharmacological actions of methotrexate in the treatment of rheumatoid arthritis. Rheumatology (Oxford). 2008;47(3):249–55.

    Article  CAS  Google Scholar 

  39. Guma M, Tiziani S, Firestein GS. Metabolomics in rheumatic diseases: desperately seeking biomarkers. Nat Rev Rheumatol. 2016;12(5):269–81.

    Article  CAS  Google Scholar 

  40. Bykerk VP, Massarotti EM. The new ACR/EULAR remission criteria: rationale for developing new criteria for remission. Rheumatology (Oxford). 2012;51(Suppl 6):vi16–20.

    Google Scholar 

Download references


The authors thank all participating institutions and personnel who were involved in the study, particularly (research) nurses, physicians’ assistants, and rheumatologists. Lastly, the authors thank the participating patients for their kind willingness to participate.


The U-Act-Early trial was funded by Roche Nederland BV.

Availability of data and materials

The generated dataset is available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



XMT: study design, data collection, data interpretation, and writing. WY: study design, data collection, data interpretation, and writing. JWGJ: study design, data interpretation, and writing. AP-S: data interpretation and writing. MEAB: data interpretation and writing. ACH: study design, data interpretation, and writing. TH: study design, data interpretation, and writing. JMVL: data interpretation and writing. JWJB: study design, data interpretation, and writing. FPJGL: study design, data interpretation, and writing. All authors read and approved the manuscript.

Corresponding author

Correspondence to Xavier M Teitsma.

Ethics declarations

Ethics approval and consent to participate

The medical ethics research committee of the University Medical Center Utrecht approved the study for all participating hospitals. All patients signed informed consent before study entry.

Consent for publication

Not applicable.

Competing interests

The department of the authors who included patients (JWGJ and JWJB) in the U-Act-Early trial received reimbursements from Roche Nederland BV. JWJB reported grants and fees from Roche, AbbVie, Bristol-Myers Squibb, Merck Sharp & Dohme, Pfizer, and UCB. JMvL received fees from Arthrogen, MSD, Pfizer, Eli Lilly, and BMS and research grants from Astra Zeneca, Roche-Genentech. FPJGL reports grants from Roche. AP-S is an employee of F Hoffmann-La Roche and MEAB is an employee of Roche Nederland BV. XMT, WY, ACH, and TH declare no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Additional information regarding the metabolite profiling on the three platforms. (DOCX 15 kb)

Additional file 2:

Overview of the pathway analysis in the (a) tocilizumab plus methotrexate, (b) tocilizumab, and (c) methotrexate strategy arms. The top three most relevant pathways in the tocilizumab plus methotrexate arm were (1) “histidine metabolism,” (2) “sphingolipid metabolism,” and (3) “arachidonic acid metabolism;” in the tocilizumab arm, these were (1) “arachidonic acid metabolism,” (2) “lysine degradation,” and (3) “cysteine and methionine metabolism;” in the methotrexate arm, these were (1) “arginine and proline metabolism,” (2) “histidine metabolism,” and (3) “aminocyl-tRNA biosynthesis.” KEGG Kyoto Encyclopedia of Genes and Genomes, tRNA transfer ribonucleic acid. The colors of the nodes, varying from yellow to red, indicates the level of significance with red being highly significant; the size of the nodes depicts the impact of the pathway with larger nodes illustrating a higher impact. (DOCX 54 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Teitsma, X.M., Yang, W., Jacobs, J.W.G. et al. Baseline metabolic profiles of early rheumatoid arthritis patients achieving sustained drug-free remission after initiating treat-to-target tocilizumab, methotrexate, or the combination: insights from systems biology. Arthritis Res Ther 20, 230 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: