Next Article in Journal
Ammonium Fluoride as Suitable Additive for HILIC-Based LC-HRMS Metabolomics
Previous Article in Journal
Metabolic Fingerprint of Chronic Obstructive Lung Diseases: A New Diagnostic Perspective
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 9
Issue 12
10.3390/metabo9120291
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mucosal Metabolomic Profiling and Pathway Analysis Reveal the Metabolic Signature of Ulcerative Colitis

by
Joseph Diab
1,
Terkel Hansen
1,
Rasmus Goll
2,3,
Hans Stenlund
4,
Einar Jensen
1,
Thomas Moritz
4,5,
Jon Florholmen
2,3 and
Guro Forsdahl
1,*
1
Natural Products and Medicinal Chemistry Research Group, Department of Pharmacy, Faculty of Health Sciences, University of Tromsø The Arctic University of Norway, 9037 Tromsø, Norway
2
Research Group of Gastroenterology and Nutrition, Department of Clinical Medicine, Faculty of Health Sciences, University of Tromsø The Arctic University of Norway, 9037 Tromsø, Norway
3
Department of Medical Gastroenterology, University Hospital of North Norway, 9037 Tromsø, Norway
4
Swedish Metabolomics Center, Department of Molecular Biology, Umeå University, 90736 Umeå, Sweden
5
The Novo Nordisk Foundation Center for Basic Metabolic Research, 2200 Copenhagen, Denmark
*
Author to whom correspondence should be addressed.
Metabolites 2019, 9(12), 291; https://doi.org/10.3390/metabo9120291
Submission received: 15 October 2019 / Revised: 21 November 2019 / Accepted: 25 November 2019 / Published: 27 November 2019
Browse Figures
Versions Notes

Abstract

:
The onset of ulcerative colitis (UC) is characterized by a dysregulated mucosal immune response triggered by several genetic and environmental factors in the context of host–microbe interaction. This complexity makes UC ideal for metabolomic studies to unravel the disease pathobiology and to improve the patient stratification strategies. This study aims to explore the mucosal metabolomic profile in UC patients, and to define the UC metabolic signature. Treatment- naïve UC patients (n = 18), UC patients in deep remission (n = 10), and healthy volunteers (n = 14) were recruited. Mucosa biopsies were collected during colonoscopies. Metabolomic analysis was performed by combined gas chromatography coupled to time-of-flight mass spectrometry (GC-TOF-MS) and ultra-high performance liquid chromatography coupled with mass spectrometry (UHPLC-MS). In total, 177 metabolites from 50 metabolic pathways were identified. The most prominent metabolome changes among the study groups were in lysophosphatidylcholine, acyl carnitine, and amino acid profiles. Several pathways were found perturbed according to the integrated pathway analysis. These pathways ranged from amino acid metabolism (such as tryptophan metabolism) to fatty acid metabolism, namely linoleic and butyrate. These metabolic changes during UC reflect the homeostatic disturbance in the gut, and highlight the importance of system biology approaches to identify key drivers of pathogenesis which prerequisite personalized medicine.

Graphical Abstract

1. Introduction

Inflammatory bowel diseases (IBD) are chronic, relapsing inflammatory disorders in the gastrointestinal tract that affect around 0.3% of the population in Europe and North America with increasing worldwide incidence [1]. The two major forms of IBD, ulcerative colitis (UC) and Crohn’s disease (CD), are characterized by a dysregulated mucosal immune response triggered by several genetic and environmental factors in the context of host–microbe interaction [2]. The interaction between these components yield a network effect, defined as ‘IBD interactome’, which results in an overwhelming complexity [3]. This complexity cannot be solved by studying one component in isolation from the others. For instance, only 10% of IBD cases can be explained by genetic variance [4], while it remains unclear whether the alteration in the microbiota is primary or secondary to the chronic inflammation in IBD [5]. Recently, multiomics approaches were suggested as a tool to unravel the IBD interactome, and to improve the patient stratification strategies toward personalized medicine [3,6]. Metabolomics, defined as the comprehensive measurement of all metabolites (low-molecular-weight molecules) in a biological specimen, is perhaps the most closely linked to the phenotype [7]. It describes the pathophysiology of a disease at the molecular level, and provides predictive, prognostic and diagnostic markers of diverse disease states [8]. Previous published metabolomic data from IBD patients have given first hints on metabolic changes during the course of the disease. However, these data are either generated from serum [9,10,11,12,13,14], stool [12], or urine [15] samples. Moreover, these studies were restricted to pediatric patients [12,14], or treated patients [9,11,13,15].
This study aimed to explore the mucosal metabolomic profile in treatment-naïve UC patients compared to treated UC patients in deep remission and to healthy subjects. The high throughput metabolomic analysis was performed by a combined gas chromatography coupled to time-of-flight mass spectrometry (GC-TOF-MS) and ultra-high performance liquid chromatography coupled with mass spectrometry (UHPLC-MS) on biopsy samples. Our results maps the metabolic changes during IBD, and highlights the metabolic signatures of the IBD interactome network.

2. Results

2.1. Subjects Characteristics

Colon biopsies collected from newly diagnosed treatment-naïve UC patients (n = 18), UC patients in state of deep remission (n = 10), and healthy controls (n = 14) were included in this study. The study group characteristics are shown in Table 1. In addition, clinical data such as tumor necrosis factor (TNF) gene expression, levels of fecal calprotectin and C-reactive protein are provided. Furthermore, data on daily supplementation with omega-3 and previous treatment with antibiotics are included in Table 1.

2.2. Mucosal Metabolite Profiles in Treatment-Naïve UC Patients, UC Remission Patients and Controls

Mucosal metabolite profiles were compared to identify significant alteration in metabolite composition in treatment-naïve patients and UC deep remission patients compared with controls (Supplementary Table S1). The results are summarized as a Venn diagram in Figure 1. Among the 177 metabolites included in this study, the mucosal levels of 60 metabolites were altered in UC treatment-naïve patients compared with healthy controls. Among these metabolites, the mucosal levels of 38 metabolites were higher, and those of 22 metabolites were lower. Similarly, the mucosal levels of 21 metabolites were changed in UC remission patients compared with healthy controls. Accordingly, the mucosal levels of 10 metabolites were higher and those of 11 metabolites were lower. The most prominent changes among the study groups were in lysophospholipids, acyl carnitine, and amino acid profiles. In addition, 46 metabolites were changed in treatment-naïve UC patients compared to deep remission UC patients.

2.3. Discriminative Models for UC State

Principle component analysis (PCA) was used as an unbiased multivariate analysis to have an overview of the variation within the data, to detect outliers, and to determine subgroups. The two main components explained 29% of the variability in the combined metabolomic data set (42 observations, 177 variables). Accordingly, the PCA t1/t2-scores plot (Figure 2A) revealed a distinct metabolomic profile in inflamed mucosa taken from treatment-naïve UC patients compared to noninflamed mucosa taken from UC remission patients and healthy controls. In addition, it was observed that the UC remission patients differed to a lesser extent from the healthy controls. Conversely, PCA did not show specific clustering patterns of the study subjects according to age, sex or activity score (Supplementary Figure S2).
A supervised orthogonal partial least squares projection to latent structures-discriminant analysis (OPLS-DA) model was built to identify the differential metabolites between active UC patients, remission UC patients and healthy controls. A significant OPLS-DA model (P-value from cross-validated analysis of variance (CV-ANOVA) was 6.15 × 10−7), with maximum separation between the study groups with good predictive ability Q2cum >0.5 was obtained (Figure 2B). Additionally, a permutation test (Supplementary Figure S2) indicated that the obtained OPLS-DA model was not influenced by overfitting. The metabolites were ranked according to the variables importance in projection (VIP) scores (Supplementary Table S2) to identify the most distinctive metabolites between the study groups at a VIP threshold > 1.5. Glutamic acid, asparagine, lysophosphatidylethanolamine LPE (O-18:0), hypoxanthine, lysophosphatidylcholine LPC (20:3), hydroxyl carnitine, and LPC (20:4) were identified as the most important metabolites in the model, and the mucosal levels of these metabolites among the study groups are represented in Figure 3.

2.4. Pathway Analysis

Integrated pathway analysis was performed to capture the metabolic pathways disruption during the active UC state, and to ease the biological interpretation. The annotated metabolites were mapped into 50 metabolic pathways (Supplementary Table S3) according to the Kyoto encyclopedia of genes and genomes KEEG database. As Figure 4 shows, several pathways were the most perturbed, ranging from amino acid metabolism (such as tryptophan metabolism, and alanine, aspartate and glutamate metabolism) to antioxidant defense pathway (glutathione pathway). Furthermore, the pathway analysis revealed a disruption in the long- and short-chain fatty acid (LCFA and SCFA) metabolism, namely linoleic metabolism and butyrate metabolism. The impact value of altered metabolic pathways, based on topology analysis, ranged from 0.01–0.66. A summary of significantly altered pathways is provided in Table 2. In addition, the complete result from the pathway analysis containing all 50 metabolic pathways is provided in Supplementary Table S3.

3. Discussion

This study provides a unique and detailed snapshot of the mucosal metabolite profile in clearly stratified UC patients (treatment-naïve, newly diagnosed, and deep remission patients). The reported 177 metabolites revealed a distinctive metabolic fingerprint in active UC patients compared with healthy controls. In addition, the metabolomic profiling coupled with pathway analysis provided a deeper understanding of the metabolome changes among UC patients with ongoing active inflammation. Several metabolic pathways were identified, including pathways related to amino acid metabolism, SCFA and LCFA metabolism, and glutathione metabolism.
To our knowledge, this is the first study of mucosal metabolomic profile in treatment-naïve and deep remission UC patients. In contrast, previous studies were restricted to bio-fluids. It is well established that tissues are under greater homeostatic regulation than plasma16. Thus, it provides highly consistent measurements among individuals [16], and better understanding of the molecular basis of diseases [17]. Moreover, previous studies included treated and untreated UC patients. In the current work, however, only treatment- naïve UC patients were represented in the active UC group. In addition, the state of remission was defined by strict criteria (endoscopy, histology, and normalized TNF gene expression). Notably, remission patients were excluded from the pathway analysis. This stratification of patients allows capturing key metabolic alterations that are exclusively associated with the UC onset. Furthermore, the combination of two analytical metabolomic platforms allowed analysing metabolites in different polarity and molecular weight ranges, and gaining a wider prospective of the metabolome [18].
According to the pathway analysis, the omega-6 linoleic acid (ω-6 LA 18:2) metabolism had the highest impact score in the pathway analysis. Ω-6 LA, is an essential fatty acid, which is metabolised to dihomo-γ-linolenic acid (ω-6 DGLA 20:3). The latter is converted by fatty acid desaturase 1 (FADS1) to ω-6 arachidonic acid (AA 20:4) [19]. DGLA and AA are esterified with glycerol in the phospholipids, such as LPC, in the cell membrane, and released by phospholipase A2 during inflammation [19]. The released AA and DGLA are metabolised to form bioactive pro- and anti-inflammatory mediators. In the current data, LA was found to be lower in active UC patients. In contrast, the mucosal levels of LPC (20:3) and LPC (20:4) were higher in treatment-naïve UC compared to healthy controls, and were considered among the top discriminant metabolites between the study groups. This finding supports evidence suggesting that the onset of IBD is characterized by an imbalance between pro- and anti-inflammatory mediators [20]. For instance, the mucosal levels of AA related pro-inflammatory metabolites were elevated in treatment-naïve UC patients [21]. In addition, variations in the FADS1 gene were found to be associated with higher susceptibility to IBD [22,23]. Therefore, it seems that the increased metabolism of LA to AA is a crucial step in the IBD pathology.
Another important finding is the alteration in the amino acid metabolism, namely the tryptophan (Trp) metabolism and the alanine, aspartate and glutamate metabolism. Recently, Trp emerged as the hub of host–microbiota crosstalk considering that Trp metabolism pathways leading to serotonin, kynurenine (Kyn), and indole derivatives are under the direct or indirect control of the microbiota [24]. It was shown that supplementation with Trp improves the clinical symptoms and reduces the pro-inflammatory cytokines production in experimental colitis [25]. Furthermore, indole derivatives act as ligands for the aryl hydrocarbon receptor (AHR) inducing local production of interleukin-22 (IL-22), which maintains intestinal homeostasis, promotes immune defense and tissue repair. In the current study, we report a decreased mucosal level of Trp and an increased level of Kyn. This is in alignment with previous studies, which have reported low serum level of Trp in UC patients [10,12]. Notably, a large cohort study consisting of 148 UC patients has concluded that a higher Trp metabolism rate is associated with UC activity [26].
Furthermore, the current data demonstrates several perturbation in amino acid metabolism during UC. For instance, the mucosal levels of glutamic acid and asparagine were low in healed mucosa, and were gradually elevated in UC remission patients and active UC patients. Accordingly, glutamic acid and asparagine were discriminative between treatment-naïve UC patients, UC patients in remission and healthy controls. Interestingly, in a previous study, high levels of amino acids were detected in stool samples from IBD patients, and were linked with the gut microbiota dysbiosis [27]. In addition, higher urinary level of asparagine and glutamic acid were reported [28]. Notably, previous study of mucosal amino acids profile in IBD patients demonstrated increased levels of several amino acid, such as aspartate, glutamine, and glutamic acid in active UC patients [29]. However, we cannot determine to which degree the reported changes in mucosal amino acid levels are caused by gut microbiota.
Altered butyrate metabolism is another evidence of the bacterial dysbiosis in UC. It is well documented that the alteration in butyrate and other short chain fatty acid (SCFA) production is a hallmark of active UC patients [30]. For instance, it was found that dysbiosis in IBD patients is characterized by a decrease in the number of SCFAs/butyrate-producing bacteria [31]. Another study has reported reduction of butyrate and propionate in stool samples of IBD patients [32]. Although the current data did not show significant changes in butyrate related metabolites in the mucosa, the decreased mucosal level of glutamine in UC patients might indicate that glutamine is being used as energy source instead of butyrate, as previously reported [33]. Interestingly, previous data have shown low abundance of proteins related to this specific utilization of butyrate in UC patients’ mucosa [34].
The variation in the acylcarnitine profile, demonstrated in the current data, could also indicate energy impairment. Acylcarnitine is a mediator that transfers catabolism products of fatty acids and amino acids into mitochondria for β-oxidation [35]. This is a key step in the process of energy production. Therefore, the accumulation of medium and long chain fatty acyl carnitine, according to the current data, provides further evidence of the mitochondrial dysfunction. However, it is unclear yet whether the mitochondrial dysfunction in IBD is caused by a dysbiosis or if it is induced by the pro-inflammatory cytokines, such as TNF [36].
Although the inclusion criteria for remission patients was mucosal healing and immunological remission [37], the present work reveals a distinct metabolome in UC deep remission patients with respect to healthy controls and active UC patients. This comes in alignment with previously published data which reports a distinct mucosal lipid composition fingerprint in UC deep remission patients compared with healthy controls and treatment-naïve UC patients [38]. Consequently, from a clinical point of view, these findings supports the emerging importance of ‘Omics’ analysis in improving the current scoring system, monitoring the disease progression and improving the treatment strategies [39].
The relatively small sample size in the current study preclude subgroup analysis according to the severity of the disease. Hence, the reported results are exploratory and need to be validated by a larger cohort, which include inflamed and non-inflamed mucosa from UC patients. In addition, to further get insight in the mechanistic behind the alteration in the metabolic pathways, gene expression and/or protein data, preferably from the same patients, should be studied. Combining such multi omics data might also underline metabolite changes caused by the gut microbiota. Furthermore, we suggest the absolute quantification and identification of metabolites involved in the pathways of interest, especially tryptophan and butyrate pathways using targeted analysis. This is especially of interest for future evaluation of clinical validity, where absolute quantitative levels is a necessity. Suggestively, future studies also need to explore the relationship between metabolic changes, microbiota dysbiosis, and the activity of IBD. This approach will provide key insight into the disease outcome and response to treatment.

4. Materials and Methods

4.1. Patients and Biopsy Collection

Mucosal biopsies were collected from newly diagnosed treatment-naïve UC patients (n = 18) and UC patients in deep remission (n = 10). The UC diagnosis was made upon clinical, endoscopic and histological criteria established by the European Crohn and Colitis Organization (ECCO) guidelines [40]. The degree of inflammation was endoscopic evaluated by the scoring system of ulcerative colitis disease activity index (UCDAI); UCDAI score of 3–5 is defined as mild, 6–8 as moderate, and 9–12 as severe UC [41]. TNF-α mRNA expression levels were measured by real-time PCR in mucosal biopsies to evaluate the UC activity [42]. The state of deep remission was achieved after treatment with anti-TNF-α monoclonal antibody biologics. Deep remission was defined as endoscopic healed mucosa by ECCO 2017 consensus (Mayo score = 0) [43] and, additionally, normalized mucosal TNF-α level [44]. Subjects performing endoscopy for colonic cancer screening, with normal findings (no ulcer, no redness) and normal colonic histological examination, served as healthy controls (n = 14).
All biopsies were acquired from the rectum or sigmoid colon (Table 1). In active UC patients, biopsies were obtained from the most inflamed mucosa. The dry weight of the biopsies ranged from 2–8 mg. All biopsies were dry-frozen immediately at −80 °C, and kept at this temperature until further analysis. The Regional Committee of Medical Ethics of North Norway and the Norwegian Social Science Data Services approved the study and the storage of biological material under the number (REK NORD 2012/1349).
In addition, all enrolled subjects have signed an informed consent form, and the study was conducted in accordance with the Declaration of Helsinki.

4.2. Chemicals and Reagents

Detailed information of chemicals used for GC-MS and UHPLC-MS analysis is provided in the supplementary data section.

4.3. Sample Preparation

Metabolite extraction was carried out as previously described [45]. Briefly, each biopsy was transferred to an Eppendorf tube and kept on ice. Then, the extraction solution (methanol:water (8:1)) with all internal standards was added to the biopsy in a solid-to-solvent ratio of 1:15 (w/v). The final concentration of UPLC-MS standards and GC-MS standards was 0.625 ng/mL and 5 ng/μL respectively. Two tungsten beads were added to each tube, and the samples were shaken at 30 Hz for 3 min in a MM301 Vibration Mill (Retsch GmbH & Company KG). The beads were removed, and the samples were further centrifuged at 14,000 rpm and 4 °C for 3 min. Finally, the supernatant was transfer to a micro vials for UHPLC-MS and GC-MS analysis, 200 µL was used for UHPLC-MS analysis and 150 µL for GC-MS analysis. Samples were dried using a vacuum concentrator (MIVac, SP, Warminster, PA, USA). Quality Control (QC) samples were prepared by pooling 10 μL from each extract. Extracts were stored at −80 °C until analysis.

4.4. UHPLC-MS Analysis

On the day of analysis, samples were reconstituted in 20 µL of methanol:water (1:1) solution. The UHPLC-MS analysis was performed with an Infinity 1290 Agilent (Agilent Technologies, Santa Clara, CA, USA) ultra-high performance liquid chromatograph coupled with tandem mass spectrometry (UHPLC-MS-MS) as previously described [46]. Briefly, 1 µL of each extract was injected into the UHPLC system equipped with an Acquity column (HSS T3, 2.1 × 50 mm, 1.8 µm C18) in combination with a 2.1 mm × 5 mm, 1.7 µm VanGuard charged-surface hybrid (CSH) precolumn (Waters Corporation, Milford, MA, USA), held at 60 °C. Mobile phases used were MilliQ water with 0.1% formic acid (A) and 75:25 acetonitrile: 2-propanol with 0.1% formic acid (B). The following gradient was used: 10% B for 2 min, then B was increased to 99% in 5 min and held at 99% for 2 min. Subsequently, B was decreased to 0.1% in 0.3 min and the flow-rate was increased to 0.8 mL min−1 for 0.5 min. These conditions were held for 0.9 min, after which the flow-rate was reduced to 0.5 mL min−1 for 0.1 min before the next injection. Samples were randomly injected. The first parallel of extracts was analyzed in positive mode. Then, the instrument was switched to negative mode and the second parallel of extracts was injected. Blank samples with only methanol:water (1:1) solution were run prior and after each samples set. MS parameters were kept identical between the modes, with exception of the capillary voltage. The exact masses of metabolites were detected with an Agilent 6550 Q-TOF mass spectrometer equipped with an iFunnel jet stream electrospray ion source (Agilent Technologies, Santa Clara, CA, USA). The flow gas temperature was set at 150 °C, the drying gas flow at 16 L min−1 and the nebulizer pressure at 35 psi. The sheath gas temperature was set at 350 °C and the sheath gas flow was 11 L min−1. The capillary voltage was set at 4000 V for the positive mode and 4500 V for the negative mode. The m/z range was 70–1700, and data were collected in centroid mode with an acquisition rate of 4 scans/s. The QC-samples were a part of the quality control of the analysis were run in the beginning of the sample set. Auto MS/MS acquisition was used when running QC samples to generate MS/MS data.

4.5. GC-MS Analysis

Prior to injection, derivatization was performed as previously described [47]. Briefly, 30 µL of a methoxyamine solution in pyridine (15 μg μL−1) was added to the dry extract, and then shaken for 15 min on a shaking table. Derivatization was carried out at 70 °C for 1 h followed by room temperature for 16 h. Afterwards, the samples were trimethylsilylated (TMS) with 30 μL methyl-N-(trimethylsilyl) trifluoroacetamide MSTFA at room temperature for 1 h. Finally, 30 μL of heptane (including 15 ng methylstearate μL−1) was added and the vials were vortexed before 1 µL was injected splitless by a CTC Combi Pal autosampler (CTC Analytics AG, Switzerland) into an Agilent 6890 GC equipped with a fused silica capillary column (10 m × 0.18 mm I.D.) with a chemically bonded 0.18 µm DB 5-MS stationary phase (J&W Scientific, Folsom, CA, USA). Samples were randomly injected. Blank samples with only heptane were run prior and after the samples set. The injector temperature was 270 °C and the purge flow-rate was 20 mL min−1. The column temperature was set to 70 °C for 2 min, then increased to 320 °C by a rate of 40 °C min−1, and held there for 2 min using a gas flow rate of 1 mL min−1. The GC was coupled to the ion source of a Pegasus III TOF-MS (Leco Corp., St Joseph, MI, USA). The transfer line and MS instrument settings were as follows: Transfer lines and ion source temperature were set to 300 and 350 respectively. The mass detecting range was set to 50 to 800 m/z. An alkane series (C10-C40) was run together with all samples.

4.6. Metabolites Identification and Data Processing

Targeted feature extraction of the acquired UHPLC–MS data was performed using the Profinder™ software package, version B.08.00 (Agilent Technologies Inc., Santa Clara, CA, USA). In-house libraries with exact masses and experimental retention times were used for identification. The libraries contained metabolites from the following chemical classes: acylcarnitines, amino acids, carbohydrates, fatty acids, bile acids, nucleotides, small peptides, and lysophospholipids, namely lysophosphatidylcholine (LPC) and lysophosphatidylethanolamine (LPE). The allowed ion species for metabolites identification were +H, +Na, +K, and +NH4 in positive ionization mode, and –H, +HCOO in negative ionization mode. The mass tolerance was 10 ppm and the retention time tolerance 0.1 min. Only one charge for each metabolite was allowed. The extracted peaks were aligned and matched between samples, and then each compound was manually checked for mass and retention time agreement with the library. A two-step filtering approach was used for peak quality control: First, peaks with bad characteristics (e.g., overloaded, sample noise, non-Gaussian) were excluded from the analysis. Second, only peaks present in at least 75% of at least one study group were included.
Raw GC–MS data files were exported in NetCDF format to a MATLAB 8.3 (R2014a) (Mathworks, Natick, MA) based in-house script for baseline correction, chromatogram alignment, and peak deconvolution. Metabolite annotation was performed based on the retention index (RI) values and MS spectra from the in-house mass spectra library established by the Swedish Metabolomics Centre (Umeå, Sweden). The total number of annotated metabolites by UHPLC-MS and GC-MS was 128 and 66 respectively. Seventeen metabolites were detected with both methodologies, the signal detected with the UHPLC method was included in the statistical analyses. The UHPLC–MS metabolites were normalized by the total peak areas, whereas GC–MS metabolites were normalized by internal standards as described before [48]. A combined data set containing 177 metabolites was submitted to statistical analysis.

4.7. Statistical Analysis

Statistical analysis was carried out using RStudio: Integrated Development Environment (version 1.0.143). Undetectable Metabolites, which represented 0.6% of the total reported metabolites, were assigned a value corresponding to half of the minimum positive value in the original data. Shapiro–Wilk test of normality was applied, and the data was not found normally distributed. Kruskal–Wallis one way analysis of variance test was performed to compare the mean concentration of metabolites between treatment-naïve UC, remission UC, and control groups. Acquired p-values were adjusted using Benjamini and Hochberg FDR method [49]. Dunn’s test [50] was applied as a post-hoc test, and significant p-value cut-off was corrected to 0.017 by Bonferroni multiple comparison method [51]. Multivariate analysis was carried out using SIMCA software (version 14.0.0.135559; Sartorius AB, Umea, Sweden). The metabolites were auto-scaled and mean-centered in order to adjust the importance of high and low abundance metabolites to an equal level [52]. Unsupervised PCA was first performed to assess the unicity of the metabolome for each of the study groups. Then, supervised OPLS-DA [53] was employed and metabolites were classified according to corresponding regression coefficients to identify the most important metabolites in discriminating between the study groups. The parameters of the OPLS-DA model were described by R2Xcum, R2Ycum and Q2cum, whereas, R2Xcum is the cumulative modeled variation in X, R2Ycum is the amount of variation in X correlated to Y (response matrix) and Q2cum is the cumulative predicted ability of the model [54]. The validity and degree of overfitting of the OPLS-DA model was assessed by conducting analysis of variance testing of cross-validated predictive residuals (CV-ANOVA), and permutation analyses.
Pathway analysis was performed using MetaboAnalyst 4.0, a web tool for metabolomics data analysis (http://www.metaboanalyst.ca/) [55]. First, all 177 metabolites were annotated according to ‘Human Metabolome Database’ (HMDB) [56] and linked to a metabolic pathway according to KEGG database [57]. Secondly, powerful pathway enrichment analysis coupled with pathway topology analysis was carried out to identify the altered metabolic pathways in active UC compared with healthy state. The enrichment analysis was based on a global test [58] while, the node/metabolite importance was measured by relative betweenness centrality [59]. Obtained P-values from the enrichment analysis were adjusted by Holm method [60]. Adjusted P-values lower than 0.05 were considered significant.

5. Conclusions

The present report provides an in-depth description of the mucosal metabolome in UC via a high-throughput metabolomic analysis of colon biopsies taken from UC treatment-naïve patients, UC patients in state of deep remission, and healthy subjects. The study of mucosal metabolites revealed the main metabolic signatures in active UC, and reflects the homeostatic disturbance in the gut. The reported metabolites were identified by searching the human-only metabolites database, and only human metabolic pathways were included in the pathway analysis. However, the gut microbiota seems to be heavily involved in altering several metabolic pathways in the colon mucosa. This highlights the importance of integrating IBD-ome compartments by system biology approaches to identify key drivers of pathogenesis that require personalized treatment.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/9/12/291/s1, Figure S1: Multivariate analysis of the mucosal metabolomic profiles according to the subjects’ sex, age, and activity score., Figure S2: Permutation test of the OPLS-DA model, Table S1: Kruskal Wallis analysis comparing the mucosal metabolomic profile among the study groups, Table S2: The variables importance in projection (VIP) scores of the 177 metabolites included in this study according to the OPLS-DA model, Table S3: Pathway analysis of the 177 metabolites included in this study in active UC compared with healthy controls.

Author Contributions

J.D. and G.F. prepared the manuscript. J.D., G.F., E.J., T.M. and J.F. planned the experiments. J.D. performed the lipidomics analysis and data processing. J.D., T.H. and H.S. performed the data analysis. J.F. and R.G. recruited patients, provided colon mucosa biopsies, reported the baseline characteristics, and provided clinical interpretation. All authors revised and approved the final version of the manuscript.

Funding

This work belongs to the Advanced Study of Inflammatory Bowel disease (ASIB) study supported by the Northern Norway Regional Health Authority [SFP-1134-13], and the University of Tromsø, The Arctic University of Norway. The publication fees for this article have been funded by a grant from the publication fund of UiT—The Arctic University of Norway.

Acknowledgments

We thank Renate Meyer for administrating the patient samples. Swedish Metabolomics Centre is acknowledged for help with the metabolomics analysis.

Conflicts of Interest

All authors certify that they have no financial interests or connections, direct or indirect, or other situations that might raise the question of bias in the work reported. This includes pertinent commercial or other sources of funding for the individual author(s) or for the associated department(s) or organization(s), personal relationships, or direct academic competition.

References

  1. Ng, S.C.; Shi, H.Y.; Hamidi, N.; Underwood, F.E.; Tang, W.; Benchimol, E.I.; Panaccione, R.; Ghosh, S.; Wu, J.C.Y.; Chan, F.K.L.; et al. Worldwide incidence and prevalence of inflammatory bowel disease in the 21st century: A systematic review of population-based studies. Lancet 2017, 390, 2769–2778. [Google Scholar] [CrossRef]
  2. Neurath, M.F. Targeting immune cell circuits and trafficking in inflammatory bowel disease. Nat. Immunol. 2019, 20, 970–979. [Google Scholar] [CrossRef] [PubMed]
  3. De Souza, H.S.P.; Fiocchi, C.; Iliopoulos, D. The IBD interactome: An integrated view of aetiology, pathogenesis and therapy. Nat. Rev. Gastroenterol. Hepatol. 2017, 14, 739–749. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, J.Z.; van Sommeren, S.; Huang, H.; Ng, S.C.; Alberts, R.; Takahashi, A.; Ripke, S.; Lee, J.C.; Jostins, L.; Shah, T.; et al. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nat. Genet. 2015, 47, 979–986. [Google Scholar] [CrossRef] [PubMed]
  5. Manichanh, C.; Borruel, N.; Casellas, F.; Guarner, F. The gut microbiota in IBD. Nat. Rev. Gastroenterol. Amp Hepatol. 2012, 9, 599. [Google Scholar] [CrossRef] [PubMed]
  6. Orazio, P.; Tommaso, M.; Stefano, C.; Panza, A.; Latiano, T.; Corritore, G.; Andriulli, A.; Latiano, A. Inflammatory bowel disease meets systems biology: A multi-omics challenge and frontier. OMICS J. Integr. Biol. 2016, 20, 692–698. [Google Scholar]
  7. Trivedi, D.K.; Hollywood, K.A.; Goodacre, R. Metabolomics for the masses: The future of metabolomics in a personalized world. New Horiz. Transl. Med. 2017, 3, 294–305. [Google Scholar] [CrossRef]
  8. Beger, R.D.; Dunn, W.; Schmidt, M.A.; Gross, S.S.; Kirwan, J.A.; Cascante, M.; Brennan, L.; Wishart, D.S.; Oresic, M.; Hankemeier, T.; et al. Metabolomics enables precision medicine: “A White Paper, Community Perspective”. Metab. Off. J. Metab. Soc. 2016, 12, 149. [Google Scholar] [CrossRef]
  9. Scoville, E.A.; Allaman, M.M.; Brown, C.T.; Motley, A.K.; Horst, S.N.; Williams, C.S.; Koyama, T.; Zhao, Z.; Adams, D.W.; Beaulieu, D.B.; et al. Alterations in lipid, amino acid, and energy metabolism distinguish crohn’s disease from ulcerative colitis and control subjects by serum metabolomic profiling. Metabolomics 2018, 14, 17. [Google Scholar] [CrossRef]
  10. Kohashi, M.; Nishiumi, S.; Ooi, M.; Yoshie, T.; Matsubara, A.; Suzuki, M.; Hoshi, N.; Kamikozuru, K.; Yokoyama, Y.; Fukunaga, K.; et al. A novel gas chromatography mass spectrometry-based serum diagnostic and assessment approach to ulcerative colitis. J. Crohns Colitis 2014, 8, 1010–1021. [Google Scholar] [CrossRef]
  11. Lai, Y.; Xue, J.; Liu, C.W.; Gao, B.; Chi, L.; Tu, P.; Lu, K.; Ru, H. Serum metabolomics identifies altered bioenergetics, signaling cascades in parallel with exposome markers in Crohn’s disease. Molecules 2019, 24, 449. [Google Scholar] [CrossRef] [PubMed]
  12. Kolho, K.L.; Pessia, A.; Jaakkola, T.; de Vos, W.M.; Velagapudi, V. Faecal and serum metabolomics in paediatric inflammatory bowel disease. J. Crohns Colitis 2017, 11, 321–334. [Google Scholar] [CrossRef] [PubMed]
  13. Murgia, A.; Hinz, C.; Liggi, S.; Denes, J.; Hall, Z.; West, J.; Santoru, M.L.; Piras, C.; Manis, C.; Usai, P.; et al. Italian cohort of patients affected by inflammatory bowel disease is characterised by variation in glycerophospholipid, free fatty acids and amino acid levels. Metabolomics 2018, 14, 140. [Google Scholar] [CrossRef] [PubMed]
  14. Daniluk, U.; Daniluk, J.; Kucharski, R.; Kowalczyk, T.; Pietrowska, K.; Samczuk, P.; Filimoniuk, A.; Kretowski, A.; Lebensztejn, D.; Ciborowski, M. Untargeted metabolomics and inflammatory markers profiling in children with Crohn’s disease and ulcerative colitis—A preliminary study. Inflamm. Bowel Dis. 2019, 25, 1120–1128. [Google Scholar] [CrossRef]
  15. Stephens, N.S.; Siffledeen, J.; Su, X.; Murdoch, T.B.; Fedorak, R.N.; Slupsky, C.M. Urinary NMR metabolomic profiles discriminate inflammatory bowel disease from healthy. J. Crohns Colitis 2013, 7, e42–e48. [Google Scholar] [CrossRef]
  16. Viant, M.R.; Ludwig, C.; Günther, U. 1D and 2D NMR spectroscopy: From metabolic fingerprinting to profiling. RSC Biomol. Sci. 2007, 44–70. [Google Scholar] [CrossRef]
  17. Gowda, G.A.N.; Zhang, S.; Gu, H.; Asiago, V.; Shanaiah, N.; Raftery, D. Metabolomics-based methods for early disease diagnostics. Expert Rev. Mol. Diagn. 2008, 8, 617–633. [Google Scholar] [CrossRef]
  18. Shulaev, V. Metabolomics technology and bioinformatics. Brief. Bioinform. 2006, 7, 128–139. [Google Scholar] [CrossRef]
  19. Hanna, V.S.; Hafez, E.A.A. Synopsis of arachidonic acid metabolism: A review. J. Adv. Res. 2018, 11, 23–32. [Google Scholar] [CrossRef]
  20. Das, U.N. Inflammatory bowel disease as a disorder of an imbalance between pro- and anti-inflammatory molecules and deficiency of resolution bioactive lipids. Lipids Health Dis. 2016, 15, 11–18. [Google Scholar] [CrossRef]
  21. Diab, J.; Al-Mahdi, R.; Gouveia-Figueira, S.; Hansen, T.; Jensen, E.; Goll, R.; Moritz, T.; Florholmen, J.; Forsdahl, G. A quantitative analysis of colonic mucosal oxylipins and endocannabinoids in treatment-naïve and deep remission ulcerative colitis patients and the potential link with cytokine gene expression. Inflamm. Bowel Dis. 2018, 25, 490–497. [Google Scholar] [CrossRef] [PubMed]
  22. Costea, I.; Mack, D.R.; Lemaitre, R.N.; Israel, D.; Marcil, V.; Ahmad, A.; Amre, D.K. Interactions between the dietary polyunsaturated fatty acid ratio and genetic factors determine susceptibility to pediatric Crohn’s disease. Gastroenterology 2014, 146, 929–931.e3. [Google Scholar] [CrossRef] [PubMed]
  23. Mesbah-Uddin, M.; Elango, R.; Banaganapalli, B.; Shaik, N.A.; Al-Abbasi, F.A. In-silico analysis of inflammatory bowel disease (IBD) GWAS loci to novel connections. PLoS ONE 2015, 10, e0119420. [Google Scholar] [CrossRef] [PubMed]
  24. Agus, A.; Planchais, J.; Sokol, H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe 2018, 23, 716–724. [Google Scholar] [CrossRef]
  25. Kim, C.J.; Kovacs-Nolan, J.A.; Yang, C.; Archbold, T.; Fan, M.Z.; Mine, Y. l-Tryptophan exhibits therapeutic function in a porcine model of dextran sodium sulfate (DSS)-induced colitis. J. Nutr. Biochem. 2010, 21, 468–475. [Google Scholar] [CrossRef]
  26. Nikolaus, S.; Schulte, B.; Al-Massad, N.; Thieme, F.; Schulte, D.M.; Bethge, J.; Rehman, A.; Tran, F.; Aden, K.; Häsler, R.; et al. Increased tryptophan metabolism is associated with activity of inflammatory bowel diseases. Gastroenterology 2017, 153, 1504–1516.e2. [Google Scholar] [CrossRef]
  27. Ni, J.; Shen, T.-C.D.; Chen, E.Z.; Bittinger, K.; Bailey, A.; Roggiani, M.; Sirota-Madi, A.; Friedman, E.S.; Chau, L.; Lin, A.; et al. A role for bacterial urease in gut dysbiosis and Crohn’s disease. Sci. Transl. Med. 2017, 9, eaah6888. [Google Scholar] [CrossRef]
  28. Martin, F.P.; Su, M.M.; Xie, G.X.; Guiraud, S.P.; Kussmann, M.; Godin, J.P.; Jia, W.; Nydegger, A. Urinary metabolic insights into host-gut microbial interactions in healthy and IBD children. World J. Gastroenterol. 2017, 23, 3643–3654. [Google Scholar] [CrossRef]
  29. Bjerrum, J.T.; Nielsen, O.H.; Hao, F.; Tang, H.; Nicholson, J.K.; Wang, Y.; Olsen, J. Metabonomics in ulcerative colitis: Diagnostics, biomarker identification, and insight into the pathophysiology. J. Proteome Res. 2010, 9, 954–962. [Google Scholar] [CrossRef]
  30. Parada Venegas, D.; De la Fuente, M.K.; Landskron, G.; González, M.J.; Quera, R.; Dijkstra, G.; Harmsen, H.J.M.; Faber, K.N.; Hermoso, M.A. Short Chain Fatty Acids (SCFAs)-mediated gut epithelial and immune regulation and its relevance for inflammatory bowel diseases. Front. Immunol. 2019, 10, 277. [Google Scholar] [CrossRef]
  31. Machiels, K.; Joossens, M.; Sabino, J.; De Preter, V.; Arijs, I.; Eeckhaut, V.; Ballet, V.; Claes, K.; Van Immerseel, F.; Verbeke, K.; et al. A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. Gut 2014, 63, 1275–1283. [Google Scholar] [CrossRef] [PubMed]
  32. Huda-Faujan, N.; Abdulamir, A.S.; Fatimah, A.B.; Anas, O.M.; Shuhaimi, M.; Yazid, A.M.; Loong, Y.Y. The impact of the level of the intestinal short chain Fatty acids in inflammatory bowel disease patients versus healthy subjects. Open Biochem. J. 2010, 4, 53–58. [Google Scholar] [CrossRef] [PubMed]
  33. Roediger, W.E.W. The colonic epithellium in ulcerative colitis: An energy-deficiancy disease? Lancet 1980, 316, 712–715. [Google Scholar] [CrossRef]
  34. Schniers, A.; Goll, R.; Pasing, Y.; Sørbye, S.W.; Florholmen, J.; Hansen, T. Ulcerative colitis: Functional analysis of the in-depth proteome. Clin. Proteom. 2019, 16, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Eaton, S.; Bartlett, K.; Pourfarzam, M. Mammalian mitochondrial beta-oxidation. Biochem. J. 1996, 320 Pt 2, 345–357. [Google Scholar] [CrossRef]
  36. Jackson, D.N.; Theiss, A.L. Gut bacteria signaling to mitochondria in intestinal inflammation and cancer. Gut Microbes 2019, 1–20. [Google Scholar] [CrossRef] [Green Version]
  37. Florholmen, J. Mucosal healing in the era of biologic agents in treatment of inflammatory bowel disease. Scand. J. Gastroenterol. 2015, 50, 43–52. [Google Scholar] [CrossRef]
  38. Diab, J.; Hansen, T.; Goll, R.; Stenlund, H.; Ahnlund, M.; Jensen, E.; Moritz, T.; Florholmen, J.; Forsdahl, G. Lipidomics in ulcerative colitis reveal alteration in mucosal lipid composition associated with the disease state. Inflamm. Bowel Dis. 2019, 25, 1780–1787. [Google Scholar] [CrossRef]
  39. Polytarchou, C.; Koukos, G.; Iliopoulos, D. Systems biology in inflammatory bowel diseases: Ready for prime time. Curr. Opin. Gastroenterol. 2014, 30, 339–346. [Google Scholar] [CrossRef] [Green Version]
  40. Stange, E.F.; Travis, S.P.; Vermeire, S.; Reinisch, W.; Geboes, K.; Barakauskiene, A.; Feakins, R.; Fléjou, J.F.; Herfarth, H.; Hommes, D.W.; et al. European evidence-based Consensus on the diagnosis and management of ulcerative colitis: Definitions and diagnosis. J. Crohns Colitis 2008, 2, 1–23. [Google Scholar] [CrossRef] [Green Version]
  41. Marteau, P.; Probert, C.S.; Lindgren, S.; Gassul, M.; Tan, T.G.; Dignass, A.; Befrits, R.; Midhagen, G.; Rademaker, J.; Foldager, M. Combined oral and enema treatment with Pentasa (mesalazine) is superior to oral therapy alone in patients with extensive mild/moderate active ulcerative colitis: A randomised, double blind, placebo controlled study. Gut 2005, 54, 960–965. [Google Scholar] [CrossRef] [PubMed]
  42. Olsen, T.; Goll, R.; Cui, G.; Husebekk, A.; Vonen, B.; Birketvedt, G.S.; Florholmen, J. Tissue levels of tumor necrosis factor-alpha correlates with grade of inflammation in untreated ulcerative colitis. Scand. J. Gastroenterol. 2007, 42, 1312–1320. [Google Scholar] [CrossRef] [PubMed]
  43. Magro, F.; Gionchetti, P.; Eliakim, R.; Ardizzone, S.; Armuzzi, A.; Barreiro-de Acosta, M.; Burisch, J.; Gecse, K.B.; Hart, A.L.; Hindryckx, P.; et al. Third European evidence-based consensus on diagnosis and management of ulcerative colitis. Part 1: Definitions, diagnosis, extra-intestinal manifestations, pregnancy, cancer surveillance, surgery, and ileo-anal pouch disorders. J. Crohns Colitis 2017, 11, 649–670. [Google Scholar] [CrossRef] [PubMed]
  44. Johnsen, K.M.; Goll, R.; Hansen, V.; Olsen, T.; Rismo, R.; Heitmann, R.; Gundersen, M.D.; Kvamme, J.M.; Paulssen, E.J.; Kileng, H.; et al. Repeated intensified infliximab induction—Results from an 11-year prospective study of ulcerative colitis using a novel treatment algorithm. Eur. J. Gastroenterol. Hepatol. 2017, 29, 98–104. [Google Scholar] [CrossRef] [PubMed]
  45. Trygg, J.; Gullberg, J.; Johansson, A.I.; Jonsson, P.; Antti, H.; Marklund, S.L.; Moritz, T. Extraction and GC/MS analysis of the human blood plasma metabolome. Anal. Chem. 2005, 77, 8086–8094. [Google Scholar]
  46. Lindahl, A.; Forshed, J.; Nordstrom, A. Overlap in serum metabolic profiles between non-related diseases: Implications for LC-MS metabolomics biomarker discovery. Biochem. Biophys. Res. Commun. 2016, 478, 1472–1477. [Google Scholar] [CrossRef] [Green Version]
  47. Karimpour, M.; Surowiec, I.; Wu, J.; Gouveia-Figueira, S.; Pinto, R.; Trygg, J.; Zivkovic, A.M.; Nording, M.L. Postprandial metabolomics: A pilot mass spectrometry and NMR study of the human plasma metabolome in response to a challenge meal. Anal. Chim. Acta 2016, 908, 121–131. [Google Scholar] [CrossRef]
  48. Redestig, H.; Fukushima, A.; Stenlund, H.; Moritz, T.; Arita, M.; Saito, K.; Kusano, M. Compensation for systematic cross-contribution improves normalization of mass spectrometry based metabolomics data. Anal. Chem. 2009, 81, 7974–7980. [Google Scholar] [CrossRef]
  49. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B (Methodol.) 1995, 57, 289–300. [Google Scholar] [CrossRef]
  50. Dunn, O.J. Multiple comparisons among means. J. Am. Stat. Assoc. 1961, 56, 52–64. [Google Scholar] [CrossRef]
  51. Bonferroni, C.E. Teoria statistica delle classi e calcolo delle probabilita. Pubblicazioni del R Istituto Superiore di Scienze Economiche e Commerciali di Firenze 1936, 8, 3–62. [Google Scholar]
  52. Van den Berg, R.A.; Hoefsloot, H.C.; Westerhuis, J.A.; Smilde, A.K.; van der Werf, M.J. Centering, scaling, and transformations: Improving the biological information content of metabolomics data. BMC Genom. 2006, 7, 142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wiklund, S.; Johansson, E.; Sjöström, L.; Mellerowicz, E.J.; Edlund, U.; Shockcor, J.P.; Gottfries, J.; Moritz, T.; Trygg, J. Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using OPLS class models. Anal. Chem. 2008, 80, 115–122. [Google Scholar] [CrossRef] [PubMed]
  54. Ni, Y.; Su, M.; Lin, J.; Wang, X.; Qiu, Y.; Zhao, A.; Chen, T.; Jia, W. Metabolic profiling reveals disorder of amino acid metabolism in four brain regions from a rat model of chronic unpredictable mild stress. FEBS Lett. 2008, 582, 2627–2636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2017, 46, D608–D617. [Google Scholar] [CrossRef]
  57. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, 353–361. [Google Scholar] [CrossRef] [Green Version]
  58. Goeman, J.J.; van de Geer, S.A.; de Kort, F.; van Houwelingen, H.C. A global test for groups of genes: Testing association with a clinical outcome. Bioinformatics 2004, 20, 93–99. [Google Scholar] [CrossRef] [Green Version]
  59. Aittokallio, T.; Schwikowski, B. Graph-based methods for analysing networks in cell biology. Brief. Bioinform. 2006, 7, 243–255. [Google Scholar] [CrossRef]
  60. Holm, S. A simple sequentially rejective multiple test procedure. Scand. J. Stat. 1979, 6, 65–70. [Google Scholar]
Metabolites 09 00291 g001 550
Figure 1. Venn diagram summarizing the comparison between the mucosal levels of metabolites in treatment-naïve ulcerative colitis (UC) patients and UC remission patients with healthy controls. Significantly altered metabolites were determined by Kruskal–Wallis test (False discovery rate (FDR) corrected P < 0.05), followed by the Dunn post-hoc test (Bonferroni adjusted P < 0.017). In total, the levels of 60 and 21 metabolites were changed in treatment-naïve UC and deep-remission UC, respectively, compared with healthy controls. The number of up/down regulated metabolites is indicated next to up/down green arrows. For simplicity, only the full names of significantly altered metabolites at a cut-off twofold change are presented. The red up/down arrows correspond to the direction of change (up/down regulation).
Figure 1. Venn diagram summarizing the comparison between the mucosal levels of metabolites in treatment-naïve ulcerative colitis (UC) patients and UC remission patients with healthy controls. Significantly altered metabolites were determined by Kruskal–Wallis test (False discovery rate (FDR) corrected P < 0.05), followed by the Dunn post-hoc test (Bonferroni adjusted P < 0.017). In total, the levels of 60 and 21 metabolites were changed in treatment-naïve UC and deep-remission UC, respectively, compared with healthy controls. The number of up/down regulated metabolites is indicated next to up/down green arrows. For simplicity, only the full names of significantly altered metabolites at a cut-off twofold change are presented. The red up/down arrows correspond to the direction of change (up/down regulation).
Metabolites 09 00291 g001
Metabolites 09 00291 g002 550
Figure 2. Multivariate analysis of the mucosal metabolomic profiles. Each subject was labeled according to the corresponding study group. (A) Principle component analysis (PCA) t1/t2-scores plots. The variation explained by PC1 and PC2 were 17.3% and 11.7%, respectively. t1 is the first component, which explains the largest variation, t2 is independent of t1 and explains second largest variation. (B) The t1/t2-score plot of the orthogonal partial least squares projection to latent structures-discriminant analysis (OPLS-DA) model (two predictive components and one orthogonal component) built from the mucosal metabolites profile of UC treatment-naïve patients, UC remission patients and healthy controls. t1 and t2 show the direction of class separation. The performance parameters R2Xcum, R2Ycum and Q2cum were 0.33, 0.77 and 0.53, respectively.
Figure 2. Multivariate analysis of the mucosal metabolomic profiles. Each subject was labeled according to the corresponding study group. (A) Principle component analysis (PCA) t1/t2-scores plots. The variation explained by PC1 and PC2 were 17.3% and 11.7%, respectively. t1 is the first component, which explains the largest variation, t2 is independent of t1 and explains second largest variation. (B) The t1/t2-score plot of the orthogonal partial least squares projection to latent structures-discriminant analysis (OPLS-DA) model (two predictive components and one orthogonal component) built from the mucosal metabolites profile of UC treatment-naïve patients, UC remission patients and healthy controls. t1 and t2 show the direction of class separation. The performance parameters R2Xcum, R2Ycum and Q2cum were 0.33, 0.77 and 0.53, respectively.
Metabolites 09 00291 g002
Metabolites 09 00291 g003 550
Figure 3. Jitter box plots of the mucosal level of the most discriminant metabolites (Variables importance of projection VIP score > 1.5 in OPLS-DA) between treatment-naïve UC, remission UC, and healthy controls. The levels of the metabolites were autoscaled for visualization. *P-value ≤ 0.017 versus healthy control was obtained by a Dunn post-hoc test.
Figure 3. Jitter box plots of the mucosal level of the most discriminant metabolites (Variables importance of projection VIP score > 1.5 in OPLS-DA) between treatment-naïve UC, remission UC, and healthy controls. The levels of the metabolites were autoscaled for visualization. *P-value ≤ 0.017 versus healthy control was obtained by a Dunn post-hoc test.
Metabolites 09 00291 g003
Metabolites 09 00291 g004 550
Figure 4. Pathway analysis, combining pathway enrichment and pathway topology analysis, of annotated metabolites in UC treatment-naïve patients and healthy controls. The x-axis marks the pathway impact and the y-axis represents the pathway enrichment. Each node marks a pathway, with larger sizes and darker colours represent higher pathway impact values and higher pathway enrichment.
Figure 4. Pathway analysis, combining pathway enrichment and pathway topology analysis, of annotated metabolites in UC treatment-naïve patients and healthy controls. The x-axis marks the pathway impact and the y-axis represents the pathway enrichment. Each node marks a pathway, with larger sizes and darker colours represent higher pathway impact values and higher pathway enrichment.
Metabolites 09 00291 g004
Table
Table 1. Description of study group characteristics.
Table 1. Description of study group characteristics.
CharacteristicsActive UC (Debut)UC RemissionHealthy Controls
Number of Subjects181014
Age, years (mean, range, P-value *)40 (20–68) 0.0948 (31–77) 0.1855 (26–83)
Gender (Female/Male)6/124/64/10
UCDAI Score (Mild, Moderate, Severe)12/2/4
Biopsy sampling site (Rectum/sigmoid)3/155/54/10
TNF-α, copies/μg of total RNA (mean, range, P-value *)18,122 (4600–31,700) 0.014675 (800–7300) 0.115478 (1800–11,300)
Fecal calprotectin, μg/g (mean, range, P-value *)828 (25–1970) < 0.0153 (25–150) 0.1546 (25–180)
C-Reactive protein, mg/L (mean, range, P-value *)16.5 (5–92) 0.085.6 (5–11) 0.315.2 (5–11)
Smoking/non-smoking1/171/93/11
Omega-3 daily supplementation (Yes/No)6/123/77/7
Antibiotic in the last 6 months prior to the biopsy (Yes/No)3/150/102/12
* computed P-value from the comparison of two means versus healthy controls group.
Table
Table 2. Altered metabolic pathways according to pathway analysis.
Table 2. Altered metabolic pathways according to pathway analysis.
KEEG PathwayNumb. MetabolitesMatched Metabolites from the Metabolomics DataAdjusted P-value **Impact ***
Linoleic Acid Metabolism15Linoleic acid *<0.0010.66
Alanine, Aspartate and Glutamate Metabolism24N-Acetyl-L-aspartic acid *; L-Asparagine *; L-Glutamine *; L-Glutamic acid *; Gamma-Aminobutyric acid; Fumaric acid; Succinic acid0.0140.53
Tryptophan Metabolism79L-Tryptophan *; 5-Hydroxyindoleacetic acid *; L-Kynurenine *; Picolinic acid; Quinolinic acid*<0.0010.15
Butyrate Metabolism40Gamma-Aminobutyric acid; L-Glutamic acid *; Fumaric acid0.0060.05
Glutathione Metabolism38L-Glutamic acid *; Cysteinylglycine; Pyroglutamic acid *; Ornithine *<0.0010.01
* Altered metabolites (P-value ≤ 0.017 versus healthy control obtained by Dunn post-hoc test). ** P-values were calculated from the enrichment analysis then adjusted by Holm method. *** Impact is the pathway impact score calculated from pathway topology analysis.

Share and Cite

MDPI and ACS Style

Diab, J.; Hansen, T.; Goll, R.; Stenlund, H.; Jensen, E.; Moritz, T.; Florholmen, J.; Forsdahl, G. Mucosal Metabolomic Profiling and Pathway Analysis Reveal the Metabolic Signature of Ulcerative Colitis. Metabolites 2019, 9, 291. https://doi.org/10.3390/metabo9120291

AMA Style

Diab J, Hansen T, Goll R, Stenlund H, Jensen E, Moritz T, Florholmen J, Forsdahl G. Mucosal Metabolomic Profiling and Pathway Analysis Reveal the Metabolic Signature of Ulcerative Colitis. Metabolites. 2019; 9(12):291. https://doi.org/10.3390/metabo9120291

Chicago/Turabian Style

Diab, Joseph, Terkel Hansen, Rasmus Goll, Hans Stenlund, Einar Jensen, Thomas Moritz, Jon Florholmen, and Guro Forsdahl. 2019. "Mucosal Metabolomic Profiling and Pathway Analysis Reveal the Metabolic Signature of Ulcerative Colitis" Metabolites 9, no. 12: 291. https://doi.org/10.3390/metabo9120291

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop