Skip to main content

Metabonomics uncovers a reversible proatherogenic lipid profile during infliximab therapy of inflammatory bowel disease



One-third of inflammatory bowel disease (IBD) patients show no response to infliximab (IFX) induction therapy, and approximately half of patients responding become unresponsive over time. Thus, identification of potential treatment response biomarkers are of great clinical significance. This study employs spectroscopy-based metabolic profiling of serum from patients with IBD treated with IFX and healthy subjects (1) to substantiate the use of spectroscopy as a semi-invasive diagnostic tool, (2) to identify potential biomarkers of treatment response and (3) to characterise the metabolic changes during management of patients with tumour necrosis factor-α inhibitors.


Successive serum samples collected during IFX induction treatment (weeks 0, 2, 6 and 14) from 87 IBD patients and 37 controls were analysed by 1H nuclear magnetic resonance (NMR) spectroscopy. Data were analysed with principal components analysis and orthogonal projection to latent structures discriminant analysis using SIMCA-P+ v12 and MATLAB.


Metabolic profiles were significantly different between active ulcerative colitis and controls, active Crohn’s disease and controls, and quiescent Crohn’s disease and controls. Metabolites holding differential power belonged primarily to lipids and phospholipids with proatherogenic characteristics and metabolites in the pyruvate metabolism, suggestive of an intense inflammation-driven energy demand. IBD patients not responding to IFX were identified as a potentially distinct group based on their metabolic profile, although no applicable response biomarkers could be singled out in the current setting.


1H NMR spectroscopy of serum samples is a powerful semi-invasive diagnostic tool in flaring IBD. With its use, we provide unique insights into the metabolic changes taking place during induction treatment with IFX. Of distinct clinical relevance is the identification of a reversible proatherogenic lipid profile in IBD patients with active disease, which partially explains the increased risk of cardiovascular disease associated with IBD.

Peer Review reports


Ulcerative colitis (UC) [1] and Crohn’s disease (CD) [2] are the two most prevailing entities of inflammatory bowel disease (IBD), characterised by chronic or recurrent episodes of intestinal inflammation. Both UC and CD present as multifactorial diseases believed to occur in genetically predisposed individuals due to an abnormal immunologic response to environmental and microbial components [3]. UC and CD share genetic and phenotypic characteristics; however, they are distinct diseases requiring precise differentiation to achieve optimal treatment regimens [4]. The diagnosis of IBD and the differentiation of UC and CD are based on a time-consuming and costly multidisciplinary approach (i.e. clinical history, endoscopy, radiology, histology, microbiology and haematology) that is hampered by the fact that approximately 10% of patients are left with the diagnosis of IBD unclassified [5]. Once a patient is diagnosed, selecting an individualised treatment regimen becomes the next challenging step, which is important because it involves costly biologic agents such as tumour necrosis factor-α (TNF-α) inhibitors, i.e. infliximab (IFX) [6, 7].

IFX is effective in inducing and maintaining remission of UC and luminal or fistulising CD. However, one-third of IBD patients show no response at all to IFX induction (i.e. primary non-responders) and up to half of patients responding to IFX become unresponsive over time (i.e. loss of response and thus withdrawal from treatment), whereas the remaining patients achieve long-term remission with continued IFX therapy [8,9,10,11]. Consequently, it is of great clinical and socioeconomic importance to identify biomarkers of efficacy of IFX before therapy is initiated. During ineffective IFX therapy, irremediable disease progression may take place, subsequently leading to severe outcomes. Further, the total cost of such unsuccessful treatments is an enormous burden on public health expenses.

In this respect, surprisingly few studies have focused on identifying predictive response biomarkers [12,13,14]. However, with the advent and feasibility of proton nuclear magnetic resonance (1H NMR) spectroscopy-based metabolic profiling [15, 16] of serum in IBD [17,18,19,20,21], metabonomics presents as an obvious tool for differential diagnoses and identification of predictive response biomarkers.

This study is a longitudinal cohort study of 1H NMR spectroscopy-based metabolic profiling of serum from IBD patients treated with IFX (1) to identify potential diagnostic biomarkers that hold differential power with respect to UC and CD patients and control subjects, (2) to provide insight into the disordered metabolism during active and quiescent IBD, and (3) to identify metabolic changes during treatment with IFX in order to explore markers of favourable outcomes.


Patient population

From May 2009 onward, blood samples were obtained from IBD patients prior to IFX infusion as part of a local drug monitoring program at the Department of Gastroenterology, Medical Section, Herlev Hospital, Denmark. All eligible UC (n = 49) and CD (n = 38) patients were naive to biologics, all had their diagnosis verified in accordance with well-established criteria [22, 23] and, at the time of IFX infusion, they were all graded in accordance with the Mayo score [24] (a score of 0–1: inactive UC, 2–4: mild UC, 5–8: moderate UC, and 9–12: severe UC), Harvey–Bradshaw (HB) score [25] (a score of 0–4: inactive CD, 5–8: mild CD, 9–16: moderate CD, and > 16: severe CD), and/or the perianal disease activity index [26] (PDAI; a score from 0 to 20, with a higher score indicating more severe disease). The healthy volunteer subjects (n = 37) were recruited locally among the staff at Herlev Hospital and were free of any daily medications (see Tables 1 and 2 for clinical details). Exclusion criteria were age above 80 or below 18 years; clinical evidence of any infections; recent (within 14 days) use of antibiotics, probiotics or high-dose prednisolone (≥60 mg/day); pregnancy; severe mental illness; and special food regimens such as a diet of fermentable oligosaccharides, disaccharides, monosaccharides and polyols, a diabetic diet or a gluten-free diet.

Table 1 Clinical details
Table 2 Clinical details

Classification of response to IFX

The outcome of IFX treatment was determined in accordance with previous studies [27, 28]:

  1. 1.

    Remission (Rem) was defined as a favourable clinical response to IFX induction (Mayo score < 2, HB score < 5, PDAI score < 5), followed by a sustained clinical remission at the initiation of maintenance therapy, i.e. at week 14.

  2. 2.

    Response (Res) was defined as a beneficial clinical response to IFX induction (reduced Mayo score, HB score and/or PDAI score) but without complete clinical remission (Mayo score ≥ 2, HB score ≥ 5, PDAI score ≥ 5) at initiation of maintenance therapy, i.e. at week 14 and despite subsequent dose optimisation.

  3. 3.

    Non-response (NRes) was defined as no clinical response to IFX induction therapy at weeks 2, 6 or 14.

Serum sample collection and preparation

Blood samples were collected during the induction treatment with IFX and obtained as trough levels with sampling 30 min prior to an IFX infusion. Hence, samples were available from time-point 0 (before first infusion of IFX), 2 weeks after the initial dose (before the second infusion), 6 weeks after the initial dose (before the third infusion) and 14 weeks after the initial dose (before continuing maintenance therapy, i.e. the fourth infusion). Patients with severe disease and no initial response to IFX treatment often never received their third and fourth infusions of IFX. One sample was available from each control subject. In total, 359 samples were available for analysis from non-fasting patients and controls. Within 3 h of sampling, the serum was collected after centrifugation (2500 × g for 5 min at ambient temperature) and stored at −80 °C until analysis.

The serum was thawed and prepared for 1H NMR spectroscopy by mixing 180 μL of serum with 60 μL of 400 mM phosphate buffer in an Eppendorf tube to achieve a final concentration of 100 mM. A total of 180 μL from each of these mixed serum samples was transferred into 3-mm sample tubes, briefly centrifuged with a hand centrifuge to remove air bubbles, and subsequently placed in the high-throughput SampleJet robotic system used for the 1H NMR experiment.

1H NMR spectroscopy

All 1H NMR experiments were performed at 288.1 K using a Bruker 600-MHz spectrometer operating at 599.35 MHz for protons and equipped with an inverse detection 5-mm cryogenic probe (BioSpin, Bruker, Rheinstetten, Germany). For all samples, the 1H NMR spectra were acquired using a Carr–Purcell–Meiboom–Gill (CPMG) pulse sequence with an echo time of 160 ms. The 90-degree pulse length was approximately 8 μs, as calculated automatically for each sample. A total of 128 scans were collected into 64,000 data points with a spectral width of 14 ppm (8370 Hz; steady-state scans = 8; acquisition time = 3.9 s). Metabolites were assigned according to previously published data [29,30,31] and online databases for metabonomics research such as Human Metabolome Database and Madison Metabolomics Consortium Database.

1H NMR spectroscopy procession

The free-induction decays for one-dimensional data were zero-filled to 64,000 data points and multiplied by an exponential function with a line-broadening factor of 1 Hz prior to Fourier transformation. All one-dimensional 1H NMR spectra were manually corrected for phase using TOPSPIN 3.1 (BioSpin, Bruker) and corrected for baseline distortions, and referenced to the left signal peak of α-glucose at 5.236 ppm using Mnova 8.1 (Mestrelab Research, Santiago de Compostela, Spain). The 1H NMR spectral region 0.5–9.0 ppm was binned with a width of 0.004 ppm (2 Hz) using the Mnova software. Water regions (4.23–5.18 ppm) were removed to avoid imperfect water saturation. The regions 1.43–1.50, 6.98–7.10, and 7.70–7.84 ppm, were aligned to avoid the effect of chemical-shift drift. The 1H NMR spectra were normalised to the total sum of the spectral integrals to compensate for sample concentration differences.

Multivariate data analysis

The applied approach to multivariate data analysis has previously been described in detail by Bjerrum et al. [32].

To visualise the general structure of each 1H NMR dataset and to identify any abnormalities or outliers (based on the principles of Hotelling’s T 2) within the data, principal components analyses were completed on mean-centred data (SIMCA-P+ software v12.0, Umetrics, Umeå, Sweden). Projection to latent structures discriminant analysis (PLS-DA) and orthogonal projection to latent structures discriminant analysis (O-PLS-DA) [33] were subsequently applied to the 1H NMR spectral data scaled to unit variance to discover metabolic differences. O-PLS is an extension of the partial least squares regression method, including an integrated orthogonal signal-correction filter [34]. As described by Cloarec et al. [35], analysis of the model is facilitated by a back-scaled transformation of the loadings, with incorporated colour-coded correlation coefficients (MATLAB v7.1, MathWorks, Natick, MA, USA) of the metabolites accountable for the differentiation. Briefly, each back-scaled loading is plotted as a function of the respective chemical shift with a colour code that signifies the weights of the discriminatory variables. A hot colour (e.g. red) signifies metabolites that are markedly different between groups, whereas a cool colour (e.g. blue) corresponds to no differences between groups.

To minimise the risk of overfitting and to validate the PLS models, a sevenfold cross-validation procedure was performed, involving iterative construction of models by repeatedly leaving out one-seventh of the samples and predicting them back into the model. This resulted in a cross-validation parameter, Q 2, indicating the predictability of the model in relation to its statistical validity. Furthermore, a cross-validation procedure, a permutation test, was completed for each model, wherein 200 models were constructed with the use of randomised classification of the samples, and Q 2 values were generated from these models and compared with the Q 2 values of the real model. If the maximum value of \( {\mathsf{Q}}_{\mathsf{max}}^2 \) from the permutation test was smaller than the Q 2 of the real model, the model was regarded as a predictable model. Similarly, R 2 was used to evaluate possibly overfitted models. Finally, with the use of an analysis of variance of the cross-validated residuals (CV-ANOVA) [36] test, a significance test was implemented to verify model validity. Only if the permutation test and the CV-ANOVA test were fulfilled at the same time were the models considered valid. The actual predictive value of each model was estimated with the use of receiver operating characteristics.

Correlation analyses between metabolites were performed with the use of the respective resonance peaks from the NMR spectra.


Classification of clinical outcome

Patients with UC and CD were primarily classified into three major groups stated above –Rem, Res and NRes. These classifications were based on clinical scores (Mayo or HB score and PDAI score) and concomitant biochemical parameters (Figs. 1 and 2). Patients with CD who went into remission had significant decreases in HB and/or PDAI scores and C-reactive protein (CRP), and an increase in haemoglobin. CD responders also had significant decreases in HB score and CRP, but no significant changes were seen in PDAI score or haemoglobin. CD non-responders experienced no significant changes at all (Fig. 1). Similarly, UC patients who went into remission had a significant decrease in Mayo score and an increase in haemoglobin and a trend toward decreased CRP. UC responders had a significant decline in Mayo score but without significant changes in haemoglobin or CRP. UC non-responders had no changes at all (Fig. 2).

Fig. 1
figure 1

Markers of disease activity in patients with Crohn’s disease (CD). Changes in Harvey–Bradshaw (HB) scores, perianal disease activity index (PDAI) and biochemical markers (C-reactive protein and haemoglobin levels) of disease activity during treatment of CD patients with infliximab (IFX) measured at the time of IFX initiation (0) and at the fourth infusion (14). Rem remission, Res responder, NRes non-responder. * P < 0.05; *** P < 0.001; **** P < 0.0001

Fig. 2
figure 2

Markers of disease activity in patients with ulcerative colitis (UC). Changes in Mayo scores and biochemical markers (C-reactive protein and haemoglobin levels) of disease activity during treatment of UC patients with infliximab (IFX) measured at the time of IFX initiation (0) and at the fourth infusion (14). Rem remission, Res responder, NRes non-responder. ** P < 0.01; *** P < 0.001; **** P < 0.0001

1H NMR spectra of serum

Typical 1H NMR spectra of serum obtained from patients with CD and UC as well as control subjects are shown in Additional file 1: Figure S1. A range of endogenous metabolites observed in the spectra is similar to that of previously reported metabolites [17,18,19,20,21] and consisted of different amino acids (e.g. valine, leucine, isoleucine, alanine, arginine, lysine, glutamine, glutamate, histidine, glycine, tyrosine and phenylalanine), short-chain fatty acids (e.g. acetate and isobutyrate), energy metabolism-related molecules (e.g. lactate, citrate and creatine), and membrane metabolites, including glycerophosphocholine (GPC) and choline. Notably, a range of lipids and lipoproteins (e.g. high-density lipoproteins (HDLs), low-density lipoproteins (LDLs) and very low-density lipoproteins (VLDL)) was also detected in the NMR spectra.

A number of unknown metabolites (U1, U2 and U3) was identified, but further characterisation was impossible in spite of correlation analysis with STOCSY (statistical total correlation spectroscopy) and database search.

Multivariate data analysis of 1H NMR spectra

Unsupervised analysis and outlier detection

Principal components analyses was initially applied to the NMR spectra acquired from the serum using data scaled to unit variance. Based on the principals of Hotelling’s T 2 (95% confidence limit), 16 samples were identified as outliers and subsequently excluded. The outliers were reviewed for any demographic abnormalities, but none could be identified. Thus, the remaining 343 samples were accessible for subsequent supervised multivariate data analysis (Additional file 2: Table S1).

Supervised analysis of phenotypes and identification of metabolites

To classify the samples, PLS-DA models comparing control subjects with either UC or CD patients, UC versus CD, treatment response types (i.e. Rem, Res and Nres), and different disease states (i.e. weeks 0, 2, 6 and 14) during the 14 weeks of treatment with IFX were generated, with CPMG NMR data as an X matrix and class information as the Y variables [35]. To uncover the metabolic changes holding differential power, the O-PLS-DA strategy was subsequently applied to each model. The O-PLS-DA models were constructed with one PLS component and one orthogonal component. The models were subsequently validated with the use of a permutation test (for the PLS-DA models) and a CV-ANOVA test (for the O-PLS-DA models).

As seen in Table 3, no significant models could be created for CD(0) versus UC(0), indicating that no differences in the metabolic profiles of active CD and UC could be identified. Both CD(0) and UC(0), however, were easily distinguished from control subjects, which was also true for CD Rem at all time points during IFX treatment (Table 3). In contrast, UC Rem became inseparable from control subjects during the 14 weeks of IFX treatment as patients went into remission (Table 3). The corresponding score plots and their back-scaled loading plots of the significant models are shown in Fig. 3, where a clear separation between each class is seen in the score plots. The corresponding back-scaled loading plots reflect the class differences of the NMR spectra and indicate relatively increased or decreased intensities of metabolites. These back-scaled loading plots were converted into a comprehensive list of markedly up- and down-regulated metabolites in Table 4. The actual prediction performance estimates are presented in Table 3 as area under the curve (AUC) and correlate with the results of the OPLS-DA models; all significant models produce AUC values above 90. The AUCs are presented in Additional file 3: Figure S2.

Table 3 Validation of PLS-DA and O-PLS-DA models
Fig. 3
figure 3

O-PLS-DA score plots and loading plots. O-PLS-DA score plots and corresponding coefficient-coded loading plots obtained from metabolic profiles of 1H NMR spectra of the serum samples. The score plots display the first PLS component and one orthogonal component for each model. A two-way separation of the samples is demonstrated in all plots. The corresponding back-scaled loading plots reflect the class differences in the NMR spectra. Upright peaks indicate a relatively increased intensity of metabolites, and downright peaks a decreased intensity of metabolites. The colours shown on the plot are associated with the significance of metabolites in separating the samples; red indicating significance at a level of P < 0.05. (0), before first infusion of infliximab; (2), before second infusion; (6), before third infusion; (14), before fourth infusion. Glc glucose, Gln glutamine, GPC glycerophosphocholine, Gly glycine, HDL high-density lipoproteins, His histidine, Lac lactate, NAG N-acetyl glycoprotein, Phe phenylalanine, U1, U2, U3 unknown metabolite, Val valine, VLDL very-low density lipoproteins, CD Crohn’s disease, UC ulcerative colitis, Rem remission

Table 4 Significant up and down-regulated metabolites

Similar assessments were made between Res patients and control subjects, between NRes patients and control subjects, and between the different response types (i.e. Rem, Res and NRes patients; Additional file 4: Table S2). CD Res patients and UC Res patients could be distinguished from control subjects throughout the IFX induction period, whereas CD NRes patients and UC NRes patients could not. Furthermore, the comparisons between the different response types resulted in insignificant models.

Next, the potential metabolic changes that take place during the 14 weeks of IFX treatment were investigated by comparing the consecutive time points within each response type (i.e. Rem, Res and NRes in UC and CD patients separately). Surprisingly, no significant models could be identified (Additional file 5: Table S3; Additional file 6: Table S4).

Finally, a correlation analysis was performed between phenylalanine and tyrosine as phenylalanine is related to the tyrosine metabolism and was found up-regulated in IBD. This resulted in a significant positive correlation (r = 0.447, P < 0.001, Additional file 7: Figure S3).

Supervised analysis of subphenotypes

A final analysis was performed on the different subphenotypes where statistically possible; these were sex, age at diagnosis (≤30, 31–40, ≥ 41 years), years with diagnosis (≤10 years/> 10 years), Mayo and HB scores (mild, moderate and severe), fistulising disease (no/yes), stenosis (no/yes), extent of inflammation (left-sided colitis/pancolitis), surgery (no/yes), smoking status (no/yes), extraintestinal manifestations (present/not present) and glucocorticoids (dependent/independent/responder). This was completed by creating PLS-DA models for each patient category, with NMR data as an X matrix and phenotype information as the Y variables. None of the models, however, turned out to be predictive (Additional file 5: Table S3 and Additional file 6: Table S4).


This 1H NMR spectroscopy-based metabonomics study of serum samples demonstrates the differential power and metabolic differences present between IBD patients and control subjects and the changes taking place during 14 weeks of IFX treatment in different treatment response types (i.e. Rem, Res or primary NRes).

Differentiation of CD and UC patients, control subjects and IFX response types

The study confirms previously published studies [17,18,19,20,21] and displays significant differential diagnostic power in models comparing active UC patients versus control subjects, active CD patients versus control subjects, and remission CD patients versus control subjects but no significant discrimination for active UC versus active CD patients, active UC versus remission UC patients, active CD versus remission CD patients, and remission UC patients versus control subjects (Table 3 and Additional file 5: Table S3 and Additional file 6: Table S4). No previous studies have actually compared the active and quiescent disease stages of CD and UC patients, although Dawiskiba et al. [17] did compare active IBD and remission IBD patients and found significant metabolic profiles that differentiated these two congregated cohorts. They also compared remission IBD patients and control subjects, but without any significant results. This, however, contrasts with the findings of Williams et al. [20], who were able to differentiate between remission CD patients and control subjects, as in the present study, and remission UC patients and control subjects, in contrast to the present study. The methodology and statistical procedures in the study by Williams et al. [20] and those used herein are comparable, and no obvious explanation has been identified for the discrepancy in outcome.

The sustained metabolic changes at 14 weeks of IFX treatment in the remission CD patients versus control subjects might be an indication of the more profound inflammation compared with the relatively superficial inflammation seen in UC patients; in Table 4, the number of significant up- or down-regulated metabolites decreases as treatment proceeds from the initial IFX infusion (CD Rem(0) vs. control) to the fourth infusion (CD Rem(14) vs. control). In contrast, patients with UC have fewer significant metabolic changes at the time of treatment initiation (UC Rem(0) vs. control), and these changes are already absent after 6 weeks of treatment. This circumstance might also explain why the metabolic profiles from the successive serum samples (i.e. weeks 0, 2, 6 and 14) within each response type (i.e. Rem, Res and NRes) did not show any significant discrimination (Additional file 5: Table S3 and Additional file 6: Table S4); the metabolic changes seen from week 0 through week 14 in patients with CD entering remission are simply too few (Table 4). Future studies need to collect serum samples within an expanded time frame, at least in CD patients. In contrast, patients with UC entering remission have very few significantly deregulated metabolites at the time of treatment initiation (Table 4, UC Rem(0) vs. control), and the subsequent metabolic changes may consequently be insufficient to create significant discrimination between successive serum samples.

As could be expected, the metabolic profiles of CD Res and UC Res patients stayed significantly different from that of control subjects (Additional file 4: Table S2) throughout the induction period. Surprisingly, this was not the case for CD NRes and UC NRes patients. The obvious explanation is the relatively low number of samples in these groups (n = 9 and 10, respectively), but it might also be explained by an entirely different metabolic profile in the serum of patients who are primary non-responders to IFX. However, comparing the different response types (i.e. Rem, Res and NRes) did not produce any significant discrimination. Thus, in order to answer this question, lager cohorts are needed in future studies. Furthermore, the lack of discriminant power might also be explained by the poor correlation in general between clinical disease activity indices and the actual burden of disease, although we did incorporate biochemical parameters to substantiate the clinical scores.

Metabolic profiles of IFX-treated IBD patients

An abnormal lipid metabolism was identified in both CD and UC patients compared with control subjects. LDL data could not be separated from the broad lipoprotein NMR peak, but the identified low levels of HDL and high levels of VLDL seem to be the consequence of malabsorption, increased transit times and the inflammatory environment (e.g. TNF-α, interleukin-1 and -6). TNF-α in particular has been shown to increase the release of free fatty acids from adipocytes, stimulate the production of triglycerides in the liver, and inhibit lipoprotein lipase and thereby the hydrolysis of triglycerides in VLDL, leading to an increase in VLDL and a decrease in HDL cholesterols [37]. This proatherogenic lipid profile correlates with the increased incidence of cardiovascular morbidity seen in IBD patients, particularly in women and adolescents [38]. As seen in this study (Table 4) and others [39, 40], treatment with IFX ameliorates this proatherogenic lipid profile in both UC and CD patients. Thus, women and young adults with chronically active disease or frequent flares might be considered candidates for concomitant statin treatment if other atherogenic risk factors are present. However, the biobank from which the samples originated did not contain patient data on quantitative lipid measurements or other potential atherogenic risk factors except for sex, age and smoking.

This study also demonstrated decreased levels of the membrane metabolite GPC (Table 4). GPC is a derivative of the membrane metabolite choline and one of the major forms of choline storage in the cytosol. The decreased levels of GPC may be a consequence of the compromised integrity of the intestinal mucosal membrane [15] with an inflammation-driven increased rate of apoptosis and cellular turnover [41]. Interestingly, sustained choline deprivation stimulates cell survival through nuclear factor-kappa B activation, leading to sustained inflammation and oncogenesis [42]. Furthermore, low levels of GPCs have been correlated with cardiovascular disease [43], underlining the proatherogenic nature of the inflammation.

In contrast to other studies [17, 19,20,21], low levels of glucose and lactate were found in both UC and CD patients compared with control subjects, and these changes persisted in CD patients achieving remission after 14 weeks of IFX treatment (Table 4). The low level of glucose could be explained by impaired intestinal absorption due to the inflammation, which is more profound in CD than UC. However, the concomitant low level of lactate and high level of pyruvate indicate the need for scavengers of reactive oxygen species [44] and an inflammation-driven high energy demand with high glycolytic activity and low lactate dehydrogenase activity and hence accumulation of pyruvate but no production of lactate [45].

Phenylalanine was persistently up-regulated and positively correlated with tyrosine in IBD patients even during induction of remission, especially in patients with CD (Table 4). This is in contrast to the other observed amino acids (i.e. valine, glutamine, glycine and histidine; Table 4). However, the pattern of increased aromatic amino acids (e.g. phenylalanine and tyrosine) and a decrease in branched chain amino acids (e.g. valine) is well known during catabolic conditions like sepsis and liver failure [46, 47]. The inflammation-driven process causes muscle proteolysis with a subsequent decrease in branched chain amino acids due to muscle energy consumption and an increase in aromatic amino acids due to hepatic insufficiency with impairment of the hepatic enzyme phenylalanine-4-hydroxylase [48].

Formate is an essential intermediary metabolite in virtually all living organisms, but in this study, formate was found in low levels in IBD patients (Table 4). This potentially reflects the high demand of formate locally in the inflamed intestines, where formate is tightly related to folate metabolism [49]. A previous study [15] has indicated an increased phospholipid metabolism (e.g. GPC and choline) and glutathione synthesis (e.g. glutamine and glycine) in inflamed colonic tissue from patients with UC, and both of these processes are linked via the methionine-homocysteine cycle, where formate and folate are essential and central elements [49].

N-Acetyl glycoprotein is a well-known acute-phase protein, which is in accordance with the high levels seen in this study (Table 4), but due to a long half-life, it is not applicable as a useful clinical tool [50]. Finally, three unidentified metabolites were consistently down-regulated in both active and inactive patients with IBD (Table 4).


1H NMR spectroscopy-based metabonomics on successive serum samples from patients with IBD has proven to be a potentially powerful semi-invasive diagnostic tool in flaring UC and CD, providing unique insights into the metabolic changes taking place during induction treatment with IFX. Thus, the identification of a reversible proatherogenic lipid profile in IBD patients with active disease might be of clinical importance in terms of statin treatment of IBD patients with atherogenic risk factors and continuous inflammation or frequent flares. Unfortunately, no discriminant biomarkers could be associated with the differentiation of CD and UC, nor could actual response biomarkers be identified, although IBD patients not responding to IFX were characterised as a potentially distinct group. However, with the increasing number of biological treatment options (e.g. TNF-α inhibitors and integrin antagonists) [51] and novel small molecules (e.g. Janus kinase inhibitors) [52], it will become imperative to develop predictive response biomarkers in order establish relational individualised medical treatment regimes.



area under the curve


Crohn’s disease




C-reactive protein


analysis of variance of the cross-validated residuals






high-density lipoproteins


inflammatory bowel disease




low-density lipoproteins


nuclear magnetic resonance




orthogonal projection to latent structures discriminant analysis


perianal disease activity index


projection to latent structures discriminant analysis






tumour necrosis factor-α


ulcerative colitis


very low-density lipoproteins


  1. Ungaro R, Mehandru S, Allen PB, Peyrin-Biroulet L, Colombel JF. Ulcerative colitis. Lancet. 2017;389:1756–70.

    Article  PubMed  Google Scholar 

  2. Baumgart DC, Sandborn WJ. Crohn's disease. Lancet. 2012;380:1590–605.

    Article  PubMed  Google Scholar 

  3. Ananthakrishnan AN. Epidemiology and risk factors for IBD. Nat Rev Gastroenterol Hepatol. 2015;12:205–17.

    Article  PubMed  Google Scholar 

  4. Cleynen I, Boucher G, Jostins L, Schumm LP, Zeissig S, Ahmad T, Andersen V, Andrews JM, Annese V, Brand S, et al. Inherited determinants of Crohn's disease and ulcerative colitis phenotypes: a genetic association study. Lancet. 2016;387:156–67.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Geboes K, Colombel J-F, Greenstein A, Jewell DP, Sandborn WJ, Vatn MH, Warren B, Riddell RH, Pathology Task Force of the International Organization of Inflammatory Bowel Diseases. Indeterminate colitis: a review of the concept--what's in a name? Inflamm Bowel Dis. 2008;14:850–57.

    Article  PubMed  Google Scholar 

  6. Vande Casteele N, Ferrante M, Van Assche G, Ballet V, Compernolle G, Van Steen K, Simoens S, Rutgeerts P, Gils A, Vermeire S. Trough concentrations of infliximab guide dosing for patients with inflammatory bowel disease. Gastroenterology. 2015;148:1320–29.

    Article  CAS  PubMed  Google Scholar 

  7. Steenholdt C, Brynskov J, Thomsen O, Munck LK, Fallingborg J, Christensen LA, Pedersen G, Kjeldsen J, Jacobsen BA, Oxholm AS, et al. Individualised therapy is more cost-effective than dose intensification in patients with Crohn's disease who lose response to anti-TNF treatment: A randomised, controlled trial. Gut. 2014;63:919–27.

    Article  PubMed  Google Scholar 

  8. Ford AC, Sandborn WJ, Khan KJ, Hanauer SB, Talley NJ, Moayyedi P. Efficacy of biological therapies in inflammatory bowel disease: systematic review and meta-analysis. Am J Gastroenterol. 2011;106:644–59.

    Article  CAS  PubMed  Google Scholar 

  9. Gisbert JP, Panés J. Loss of response and requirement of infliximab dose intensification in Crohn's disease: a review. Am J Gastroenterol. 2009;104:760–67.

    Article  CAS  PubMed  Google Scholar 

  10. Hanauer SB, Feagan BG, Lichtenstein GR, Mayer LF, Schreiber S, Colombel JF, Rachmilewitz D, Wolf DC, Olson A, Bao W, et al. Maintenance infliximab for Crohn's disease: the ACCENT I randomised trial. Lancet. 2002;359:1541–49.

    Article  CAS  PubMed  Google Scholar 

  11. Nielsen OH, Ainsworth MA. Tumor necrosis factor inhibitors for inflammatory bowel disease. N Engl J Med. 2013;369:754–62.

    Article  CAS  PubMed  Google Scholar 

  12. Leal RF, Planell N, Kajekar R, Lozano JJ, Ordás I, Dotti I, Esteller M, Masamunt MC, Parmar H, Ricart E, et al. Identification of inflammatory mediators in patients with Crohn's disease unresponsive to anti-tnfα therapy. Gut. 2015;64:233–42.

    Article  CAS  PubMed  Google Scholar 

  13. Arijs I, Li K, Toedter G, Quintens R, Van Lommel L, Van Steen K, Leemans P, De Hertogh G, Lemaire K, Ferrante M, et al. Mucosal gene signatures to predict response to infliximab in patients with ulcerative colitis. Gut. 2009;58:1612–19.

    Article  CAS  PubMed  Google Scholar 

  14. Arijs I, Quintens R, Van Lommel L, Van Steen K, De Hertogh G, Lemaire K, Schraenen A, Perrier C, Van Assche G, Vermeire S, et al. Predictive value of epithelial gene expression profiles for response to infliximab in Crohn's disease. Inflamm Bowel Dis. 2010;16:2090–98.

    Article  PubMed  Google Scholar 

  15. Bjerrum JT, Nielsen OH, Hao F, Tang H, Nicholson JK, Wang Y, Olsen J. Metabonomics in ulcerative colitis: Diagnostics, biomarker identification, and insight into the pathophysiology. J Proteome Res. 2010;9:954–62.

    Article  CAS  PubMed  Google Scholar 

  16. Bjerrum JT. Metabonomics: analytical techniques and associated chemometrics at a glance. Methods Mol Biol. 2015;1277:1–14.

    Article  CAS  PubMed  Google Scholar 

  17. Dawiskiba T, Deja S, Mulak A, Ząbek A, Jawień E, Pawełka D, Banasik M, Mastalerz-Migas A, Balcerzak W, Kaliszewski K, et al. Serum and urine metabolomic fingerprinting in diagnostics of inflammatory bowel diseases. World J Gastroenterol. 2014;20:163–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Fathi F, Oskouie AA, Tafazzoli M, Naderi N, Sohrabzedeh K, Fathi S, Norouzinia M, Nejad Rostami M. Metabonomics based NMR in Crohn's disease applying PLS-DA. Gastroenterol Hepatol Bed Bench. 2013;6:S82–86.

    PubMed  PubMed Central  Google Scholar 

  19. Zhang Y, Lin L, Xu Y, Lin Y, Jin Y, Zheng C. (1)H nmr-based spectroscopy detects metabolic alterations in serum of patients with early-stage ulcerative colitis. Biochem Biophys Res Commun. 2013;433:547–51.

    Article  CAS  PubMed  Google Scholar 

  20. Williams HRT, Willsmore JD, Cox IJ, Walker DG, Cobbold JFL, Taylor-Robinson SD, Orchard TR. Serum metabolic profiling in inflammatory bowel disease. Dig Dis Sci. 2012;57:2157–65.

    Article  CAS  PubMed  Google Scholar 

  21. Schicho R, Shaykhutdinov R, Ngo J, Nazyrova A, Schneider C, Panaccione R, Kaplan GG, Vogel HJ, Storr M. Quantitative metabolomic profiling of serum, plasma, and urine by (1)H NMR spectroscopy discriminates between patients with inflammatory bowel disease and healthy individuals. J Proteome Res. 2012;11:3344–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Van Assche G, Dignass A, Panes J, Beaugerie L, Karagiannis J, Allez M, Ochsenkühn T, Orchard T, Rogler G, Louis E, et al. The second European evidence-based consensus on the diagnosis and management of Crohn's disease: definitions and diagnosis. J Crohns Colitis. 2010;4:7–27.

    Article  PubMed  Google Scholar 

  23. Dignass A, Eliakim R, Magro F, Maaser C, Chowers Y, Geboes K, Mantzaris G, Reinisch W, Colombel J-F, Vermeire S, et al. Second European evidence-based consensus on the diagnosis and management of ulcerative colitis part 1: Definitions and diagnosis. J Crohns Colitis. 2012;6:965–90.

    Article  PubMed  Google Scholar 

  24. Schroeder KW, Tremaine WJ, Ilstrup DM. Coated oral 5-aminosalicylic acid therapy for mildly to moderately active ulcerative colitis. A randomized study. N Engl J Med. 1987;317:1625–29.

    Article  CAS  PubMed  Google Scholar 

  25. Harvey RF, Bradshaw JM. A simple index of Crohn's-disease activity. Lancet. 1980;1:514.

    Article  CAS  PubMed  Google Scholar 

  26. Irvine EJ. Usual therapy improves perianal Crohn's disease as measured by a new disease activity index. McMaster IBD study group. J Clin Gastroenterol. 1995;20:27–32.

    Article  CAS  PubMed  Google Scholar 

  27. Afif W, Loftus EV, Faubion WA, Kane SV, Bruining DH, Hanson KA, Sandborn WJ. Clinical utility of measuring infliximab and human anti-chimeric antibody concentrations in patients with inflammatory bowel disease. Am J Gastroenterol. 2010;105:1133–39.

    Article  CAS  PubMed  Google Scholar 

  28. Steenholdt C, Bendtzen K, Brynskov J, Thomsen O, Ainsworth MA. Cut-off levels and diagnostic accuracy of infliximab trough levels and anti-infliximab antibodies in Crohn's disease. Scand J Gastroenterol. 2011;46:310–18.

    Article  CAS  PubMed  Google Scholar 

  29. Dong F, Zhang L, Hao F, Tang H, Wang Y. Systemic responses of mice to dextran sulfate sodium-induced acute ulcerative colitis using (1)H NMR spectroscopy. J Proteome Res. 2013;12:2958–66.

    Article  CAS  PubMed  Google Scholar 

  30. Zhang X, Wang Y, Hao F, Zhou X, Han X, Tang H, Ji L. Human serum metabonomic analysis reveals progression axes for glucose intolerance and insulin resistance statuses. J Proteome Res. 2009;8:5188–95.

    Article  CAS  PubMed  Google Scholar 

  31. Nicholson JK, Foxall PJ, Spraul M, Farrant RD, Lindon JC. 750 mhz 1H and 1H-13C NMR spectroscopy of human blood plasma. Anal Chem. 1995;67:793–811.

    Article  CAS  PubMed  Google Scholar 

  32. Bjerrum JT, Wang Y, Hao F, Coskun M, Ludwig C, Günther U, Nielsen OH. Metabonomics of human fecal extracts characterize ulcerative colitis, Crohn's disease and healthy individuals. Metabolomics. 2015;11:122–33.

    Article  CAS  PubMed  Google Scholar 

  33. Bylesjö M, Rantalainen M, Cloarec O, Nicholson JK, Holmes E, Trygg J. OPLS discriminant analysis: combining the strengths of PLS-DA and SIMCA classification. J Chemom. 2006;20:341–51.

    Article  Google Scholar 

  34. Trygg J, Wold S. Orthogonal projections to latent structures (O-PLS). J Chemom. 2002;16:119–28.

    Article  CAS  Google Scholar 

  35. Cloarec O, Dumas ME, Trygg J, Craig A, Barton RH, Lindon JC, Nicholson JK, Holmes E. Evaluation of the orthogonal projection on latent structure model limitations caused by chemical shift variability and improved visualization of biomarker changes in 1H NMR spectroscopic metabonomic studies. Anal Chem. 2005;77:517–26.

    Article  CAS  PubMed  Google Scholar 

  36. Eriksson L, Trygg J, Wold S. CV-ANOVA for significance testing of PLS and OPLS models. J Chemom. 2008;22:594–600.

    Article  CAS  Google Scholar 

  37. Popa C, Netea MG, van Riel PLCM, van der Meer JWM, Stalenhoef AFH. The role of tnf-alpha in chronic inflammatory conditions, intermediary metabolism, and cardiovascular risk. J Lipid Res. 2007;48:751–62.

    Article  CAS  PubMed  Google Scholar 

  38. Singh S, Kullo IJ, Pardi DS, Loftus EV. Epidemiology, risk factors and management of cardiovascular diseases in IBD. Nat Rev Gastroenterol Hepatol. 2014;12:26–35.

    Article  PubMed  Google Scholar 

  39. Koutroubakis IE, Oustamanolakis P, Malliaraki N, Karmiris K, Chalkiadakis I, Ganotakis E, Karkavitsas N, Kouroumalis EA. Effects of tumor necrosis factor alpha inhibition with infliximab on lipid levels and insulin resistance in patients with inflammatory bowel disease. Eur J Gastroenterol Hepatol. 2009;21:283–88.

    Article  CAS  PubMed  Google Scholar 

  40. Parmentier-Decrucq E, Duhamel A, Ernst O, Fermont C, Louvet A, Vernier-Massouille G, Cortot A, Colombel J-F, Desreumaux P, Peyrin-Biroulet L. Effects of infliximab therapy on abdominal fat and metabolic profile in patients with Crohn's disease. Inflamm Bowel Dis. 2009;15:1476–84.

    Article  PubMed  Google Scholar 

  41. Pedersen J, LaCasse EC, Seidelin JB, Coskun M, Nielsen OH. Inhibitors of apoptosis (iaps) regulate intestinal immunity and inflammatory bowel disease (IBD) inflammation. Trends Mol Med. 2014;20:652–65.

    Article  CAS  PubMed  Google Scholar 

  42. Holmes-McNary MQ, Baldwin AS, Zeisel SH. Opposing regulation of choline deficiency-induced apoptosis by p53 and nuclear factor kappab. J Biol Chem. 2001;276:41197–1204.

    Article  CAS  PubMed  Google Scholar 

  43. Syme C, Czajkowski S, Shin J, Abrahamowicz M, Leonard G, Perron M, Richer L, Veillette S, Gaudet D, Strug L, et al. Glycerophosphocholine metabolites and cardiovascular disease risk factors in adolescents: a cohort study. Circulation. 2016;134:1629–36.

    Article  CAS  PubMed  Google Scholar 

  44. Liu J, Litt L, Segal MR, Kelly MJS, Pelton JG, Kim M. Metabolomics of oxidative stress in recent studies of endogenous and exogenously administered intermediate metabolites. Int J Mol Sci. 2011;12:6469–501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Voet D, Voet JG, Pratt CW. Principles of biochemistry. Hoboken, NJ, Chichester: John Wiley; 2012.

    Google Scholar 

  46. Morgan MY, Milsom JP, Sherlock S. Plasma ratio of valine, leucine and isoleucine to phenylalanine and tyrosine in liver disease. Gut. 1978;19:1068–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Freund HR, Ryan JA, Fischer JE. Amino acid derangements in patients with sepsis: Treatment with branched chain amino acid rich infusions. Ann Surg. 1978;188:423–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Scholl-Bürgi S, Schroecksnadel S, Jenny M, Karall D, Fuchs D. Chronic immune stimulation may cause moderate impairment of phenylalanine 4-hydroxylase. Pteridines. 2011;22:120–25.

    Article  Google Scholar 

  49. Lamarre SG, Morrow G, Macmillan L, Brosnan ME, Brosnan JT. Formate: an essential metabolite, a biomarker, or more? Clin Chem Lab Med. 2013;51:571–78.

    Article  CAS  PubMed  Google Scholar 

  50. Vermeire S, Van Assche G, Rutgeerts P. Laboratory markers in IBD: useful, magic, or unnecessary toys? Gut. 2006;55:426–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Danese S, Vuitton L, Peyrin-Biroulet L. Biologic agents for IBD: practical insights. Nat Rev Gastroenterol Hepatol. 2015;12:537–45.

    Article  CAS  PubMed  Google Scholar 

  52. Coskun M, Vermeire S, Nielsen OH. Novel targeted therapies for inflammatory bowel disease. Trends Pharmacol Sci. 2017;38:127–42.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


We acknowledge NMR access to spectrometers at HWB-NMR (Birmingham), supported by the European Union, FP7-INFRA-2010-1.1.8-261863 Bio-NMR. Michelle A. C. Reed was also supported in part by a COSMOS EU grant (FP7- INFRA-2012-1-312941).

Availability of data and materials

All data are deposited on MetaboLights.

Author information

Authors and Affiliations



Study concept and design: JTB, CS, MA and OHN; acquisition of data: JTB, MAR, KA and UG; data analysis and interpretation: JTB, MAR, KA, UG, FH and YW; technical and material support: MAR, KA and UG; manuscript drafting: JTB, CS and OHN; critical revision of the manuscript: MA, MAR, KA, UG, FH and YW. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jacob Tveiten Bjerrum.

Ethics declarations

Ethics approval and consent to participate

All patients participated in this study in confirmation with the principles outlined in the Declaration of Helsinki and with the approval of the Scientific Ethics Committee of the Copenhagen Capital Region (H-D-2009-055).

Consent for publication

Not applicable.

Competing interests

During the last 2 years, C. Steenholdt has served as a speaker for Abbvie, Pfizer and MSD, and as a consultant for Pfizer and Takeda Pharmaceuticals. The remaining authors have no competing interests to disclose.

Publisher's Note

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

Additional files

Additional file 1: Figure S1.

1H NMR spectra of serum and metabolite assignment. (DOCX 116 kb)

Additional file 2: Table S1.

Number of subjects and serum samples. (DOCX 18 kb)

Additional file 3: Figure S2.

Area under the receiver operating characteristic (ROC) curve. (DOCX 211 kb)

Additional file 4: Table S2.

Validation of PLS-DA and O-PLS-DA models. (DOCX 24 kb)

Additional file 5: Tables S3.

Validation of PLS-DA and O-PLS-DA models of UC phenotypes and treatment response. (DOCX 23 kb)

Additional file 6: Tables S4.

Validation of PLS-DA and O-PLS-DA models of CD phenotypes and treatment response. (DOCX 24 kb)

Additional file 7: Figure S3.

Correlation analysis between phenylalanine and tyrosine. (DOCX 49 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

Bjerrum, J.T., Steenholdt, C., Ainsworth, M. et al. Metabonomics uncovers a reversible proatherogenic lipid profile during infliximab therapy of inflammatory bowel disease. BMC Med 15, 184 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: