Int J Biol Sci 2022; 18(5):1813-1828. doi:10.7150/ijbs.67112 This issue

Research Paper

Lipidomic Trajectories Characterize Delayed Mucosal Wound Healing in Quiescent Ulcerative Colitis and Identify Potential Novel Therapeutic Targets

Jacob Tveiten Bjerrum1 Corresponding address, Yulan Wang2 Corresponding address, Jingtao Zhang2, Lene Buhl Riis3, Ole Haagen Nielsen1, Jakob Benedict Seidelin1

1. Department of Gastroenterology, Medical Section, Herlev Hospital, University of Copenhagen, Denmark.
2. Singapore Phenome Centre, Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore.
3. Department of Pathology, Herlev Hospital, University of Copenhagen, Denmark.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Bjerrum JT, Wang Y, Zhang J, Riis LB, Nielsen OH, Seidelin JB. Lipidomic Trajectories Characterize Delayed Mucosal Wound Healing in Quiescent Ulcerative Colitis and Identify Potential Novel Therapeutic Targets. Int J Biol Sci 2022; 18(5):1813-1828. doi:10.7150/ijbs.67112. Available from https://www.ijbs.com/v18p1813.htm

File import instruction

Abstract

Graphic abstract

Improving the long-term prognosis of ulcerative colitis (UC) requires sustained deep mucosal colonic healing with histologic remission, making the study of colonic tissue regeneration essential. In experimental colitis models, lipid metabolites are recognized as pivotal components of this process. This study aimed to describe the kinetics of wound healing and lipid metabolites engaged in regeneration in the normal colonic mucosa and how they are affected in UC to reveal new therapeutic targets. Experimental colonic wounds were created endoscopically in quiescent UC (n=21) and controls (n=9), and the healing process was surveilled by serial endoscopies and cross-sectional wound biopsies post-wounding. Biopsies were analyzed by liquid chromatography coupled with mass spectrometry. Endoscopic wound scores were significantly higher in UC at day two (p=0.001) and seven (p<0.0001) post-wounding, demonstrating a prolonged wound healing process. The wound scores were correlated with lipid mediators crucial for normal regeneration and sustained UC-specific changes in key phospholipids and eicosanoids, i.e., lysophosphatidylcholine, phosphatidylcholine, lysophosphatidic acid, phosphatidylglycerol, phosphatidylinositol, prostaglandin D2, and prostaglandin E1, were observed. A prolonged wound healing process is identified in quiescent UC with altered disease specific lipidomic trajectories providing potential novel therapeutic avenues for stimulating mucosal regeneration as an add-on to the traditional immune suppression treatment.

Keywords: Eicosanoids, inflammatory bowel disease, metabolomics, metabonomics, mucosal healing, phospholipids

Introduction

Ulcerative colitis (UC) [1] is the most prevailing entity under the umbrella term of inflammatory bowel disease (IBD) with an increased incidence and prevalence worldwide, although a plateau seems to have been reached in the Western world [2]. The disease is characterized by inflammation and concomitant wounds in the colon. Affected individuals experience a substantial impact on the quality of life during phases of flares, characterized by bloody diarrhea and abdominal discomfort [3]. UC carries a cumulative risk of colectomy of 7.5% after 5 years [4], and certain subgroups (e.g., concomitant primary sclerosing cholangitis, early onset, extensive and long-standing activity, or a familial history of colorectal cancer (CRC)) have a slightly increased risk for development of CRC [5]. In recent years it has gradually become evident that the main risk factors for disease activity, colectomy, or CRC are sustained wound formation and inflammation, and that medication improving mucosal healing modifies the disease course favorably [6]. On the other hand, there are several lines of evidence showing that anti-inflammatory strategies might impair important pathways involved in mucosal tissue regeneration [7]. Accordingly, research on wound healing in the healthy colon and in patients with UC is critical in order to identify potential therapeutic targets that could modify the disease course through promotion of tissue regeneration.

Research within the lipidome - a subset of the metabolome - is receiving increased attention [8] as it represents lipid metabolites that are ubiquitous and with a surprising complex plethora of bioactive functions [9]. One such member of the lipidome, prostaglandin E2 (PGE2), a principal eicosanoid with multiple functions, has been shown to promote early wound repair in mice by inducing wound-associated epithelial (WAE) cells, forming an initial barrier over the wound [10]. Once a wound is closed by WAE cells, the level of PGE2 diminish, which is essential for the sequential wound healing process, as a continuously high level of PGE2 inhibits fibroblast migration, represses normal stem cell functions, and drives cells into a constant proinflammatory state as seen in vitro as well as in murine models of colitis [10-12]. After PGE2 levels subside during wound repair, the anti-inflammatory PGD2 and its metabolites rise, and the wound healing process enters the resolution phase [13]. Here resolvin E1, another eicosanoid and an endogenous lipid mediator derived from omega-3-fatty acids, inhibits the inflammatory process and promotes epithelial cell migration and proliferation, and ultimately wound repair in an experimental mouse model of mucosal wound healing [14]. Another key lipid is lysophosphatidic acid (LPA), which is an intermediate phospholipid released by activated platelets and fibroblasts upon injury. This metabolite promotes platelet aggregation and induces cellular tension and cell surface fibronectin assembly, which are key events in the wound healing process [15]. Thus, eicosanoids and phospholipids are considered pivotal components in the various stages of the wound healing process, and studies in both human [16] and experimental colitis models [14,17] have provided promising treatment results with some of these lipid metabolites. However, only very few of the lipids are well-described in human wound healing, and due to the lack of applicable human in vivo models, research on the lipidome and human colonic mucosal wound healing remains unexploited.

In this study, a novel human in vivo colonic wound healing assay was applied in which small experimental wounds were created with a biopsy forceps in the human colonic mucosa. The wound healing process was surveilled by serial endoscopies at day two and seven post-wounding during which a wound score was calculated and biopsies across the experimental injuries were obtained. Multiple initial injuries were made in order to sample only one wound healing biopsy from each injury site. Multiple sampling from the same injury site was thus not performed. The serial biopsies were subsequently analyzed by liquid chromatography coupled with mass spectrometry (LC-MS) producing a metabolic profile focusing on phospholipids and eicosanoids. Using this approach, the aims were: 1) to describe the kinetics of colonic mucosal wound healing in healthy individuals and patients with quiescent UC, 2) to characterize the lipid changes during wound healing, and thus 3) to identify the most likely new potential therapeutic targets that might be manipulated pharmacologically to enhance tissue regeneration among patients with UC.

Materials and methods

Population

Cohort 1, a test cohort for the initial histological evaluation of the wound healing process. Eight patients with UC in remission and eight healthy controls were included. The results were used to provide a rational timeframe for the observation points in the subsequent Cohort 2.

Cohort 2, the kinetics and lipidomics cohort. Sample size calculation is an inherent problem in omics analyses, however, based on our previous extensive experience with metabonomics we estimated that eight patients in each group were the absolute minimum to detect differences of relevance [18-20]. Thus, twenty-one patients with quiescent UC as well as nine healthy control subjects were included.

All patients were recruited from the outpatient clinic at the Dept. of Gastroenterology, Herlev Hospital, University of Copenhagen, Denmark, and only patients with an affirmed diagnosis of UC according to well-established criteria [1] were included. The disease had to be in complete endoscopically (endoscopic Mayo sub-score 0) [21], histological (Geboes score < 0.1) [22] and clinically remission (total Mayo score < 2) [21] at time of inclusion. On-going maintenance treatment with 5-aminosalicylic acid, thiopurine, and biologics was allowed but needed to be on a stable dosing at least 3 months prior to inclusion. The healthy controls were volunteer subjects without any known diseases and free of daily medication. No dietary restrictions were imposed.

All participants included were between 18 and 75 years old, without a history of cancer, and without any other medication affecting the immune system apart from the above-mentioned stable maintenance medication. Lactating or pregnant subjects as well as subjects with psychiatric or neurological disorders that would affect decision-making were excluded from participating in the study.

The project, including the in vivo wound assay on humans, was approved by the Scientific Ethics Committee of the Copenhagen Capital Region (Protocol No. H-18003763).

The human in vivo wound healing assay

A simple human in vivo colonic mucosal wound healing model (Figure 1A and B) and concomitant wound scores (Figure 1C) were applied as previously described [23]. In short, an experimental wound was generated in the rectosigmoid part of the colon by use of biopsy forceps (Radial Jawtm 4, 2.8 mm, Boston Scientific, Marlborough, MA, USA), which generated an elongated wound of approximately 7 × 2 mm and extended into the lamina muscularis mucosa only. To avoid unnecessary mechanical forces affecting the wound healing process and intersegmental variation, the wounds were placed in between mucosal folds at the rectosigmoid intersection. The wound and healing processes were documented by successive macroscopic high-definition imaging and recorded through a standard endoscope at the given time points, and all digitally recorded procedures were independently reviewed and scored by two gastroenterologists in a blinded fashion. Subsequent molecular characterization was made possible by obtaining biopsies across the experimental wounds at the pre-defined time points with the use of a biopsy forceps angled 90 degrees to the initial experimental wounds. Importantly, each experimental wound was sampled only once to avoid repeated sampling across a single wound. The model consequently allowed a description of the kinetics of wound healing in the human colonic mucosa on both the macroscopic and molecular levels.

Participants in Cohort 1 had a sigmoidoscopy performed at day zero, one, and two, whereas subjects in Cohort 2 had a sigmoidoscopy at day zero, two, and seven. Bowel preparation with enemas was performed prior to each sigmoidoscopy. At the index endoscopy, four biopsies were obtained; three biopsies were snap frozen and stored at -80 °C, and one biopsy fixed in formalin. At the two subsequent sigmoidoscopies, two biopsies were obtained exactly across the index biopsies (Figure 1B); the first biopsy was snap frozen and stored at -80 °C and the second fixed in formalin.

Histological assessment of inflammation

Inflammation was assessed using the Geboes score [22], which grades inflammation and structural changes associated with UC on a scale from 0.0 - 5.4. However, a modified Geboes score was applied leaving out grade of “ulceration”, as all of the samples had ulcers from the experimental wound. The modified Geboes score thus ranged from 0.0-4.4.

Lipidomics

Sample preparation

The tissue samples were weighed into their respective bead beating tubes. A total of 50 µL of phospholipid internal standard mixture, 5 µL of oxylipins internal standard, 10 µL of butylated hydroxytoluene (BHT) and 435 µL of methanol were added into the respective tubes [24]. Tissue samples were homogenized and incubated at 4 ºC for 1 hour. Samples were then centrifuged for 10 min at 10,000 g and 4 ºC. Five µL of the supernatant was diluted 500 times for phospholipid quantification in positive ion mode, and 350 µL of the supernatant was dried and reconstituted in 50 µL of methanol for quantification of phospholipid and eicosanoids in negative ion mode. Two µL of the reconstituted sample was pooled for phospholipid identification in negative ion mode, and a 100 times dilution of the pooled sample was used for phospholipid identification in positive ion mode.

Ultraperformance liquid chromatography-mass spectrometry (UPLC-MS) qualification of phospholipids and eicosanoids

An ACQUITY UPLC system coupled tandem quadrupole mass spectrometer (Xevo TQ-S, Waters, Milford, MA, USA) operated in electrospray ionization (ESI) mode and reverse phase UPLC column, ACQUITY UPLC BEH C18 column (2.1 mm × 100 mm, 1.7 µm, Waters), was employed for the separation and detection of these compounds. The data acquisition was performed in multiple reaction-monitoring (MRM) mode. For phospholipids, capillary voltage was set at 3.20 kV and 2.00 kV for the positive and negative ion mode runs, respectively. The phospholipids were separated on an isocratic flow of 0.4 mL/min for 6 min with A (methanol (98%), water (2%), ammonium acetate (5 mM), and formic acid (0.01%)) for positive detection, or 5 min with B (methanol (98%), water (2%), ammonium acetate (10 mM), and formic acid (0.01%)) for negative ion mode detection. The column temperature was maintained at 45 ºC, and the injection volume was 2 µL for both ion modes. The separation of eicosanoids were performed at a flow rate of 0.6 mL/min with an initial gradient of 25% B (isopropyl alcohol:acetonitrile (90:10)) and 75% A (water with formic acid (0.1%)) for 1 min, then a linear gradient to 95% B at 8 min, held until 8.5 min. The total injection volume was 5 µL, and the column temperature was maintained at 40 ºC. The capillary voltage was set at 2.5 kV, source temperature at 150 ºC, and a desolvation gas flow rate of 1000 l/hour at 600 ºC.

For internal standards, chemicals, and quantification analysis please see Supplementary material.

Data analysis

Univariate analysis

All continuous variables, i.e., wound score, Geboes score, biopsy weight, and lipids, were analyzed with nonparametric statistics using GraphPad Prism (v9.0.0, GraphPad, San Diego, CA, USA) and data presented as medians with interquartile ranges. The Mann-Whitney and Wilcoxon tests were applied when relevant with a significant level of p<0.05. A false discovery rate (FDR) below 5% was used for correction of multiple testing.

Multivariate data analysis

Multivariate statistical analyses were initiated with principal component analysis (PCA) to visualize the general structure of the mean centered data set and to identify any abnormalities or outliers (based on the principals of Hotelling´s T2) using SIMCA-P+ (v12.0, Umetrics, Umeå, Sweden). Class-belongings, i.e., control or UC day zero, two, or seven, and variables, were characterized by supervised models using the orthogonal partial least squares discriminant analysis (OPLS-DA) on data scaled to unit variance (UV). The models were validated with a seven-fold cross validation and permutation tests (200 permutations). The significance of the OPLS-DA models was additionally validated by the analysis of variance of the cross-validated residuals (CV-ANOVA), and a p-value < 0.05 indicated a valid model. The models were considered valid only if the permutation test and the CV-ANOVA test were satisfied at the same time. Variables holding strong differential power in the models were identified with the use of correlation coefficients. Similarly, correlation analyses and validations were performed between the colonic mucosal wound scores and phospholipids and eicosanoids, respectively. All data are deposited at MetaboLights.

 Figure 1 

The wound healing model. A and B) illustrative presentation of the model. C) The colonic mucosal composite wound healing score ranges from 1-5 with two variables: 1) a regeneration score ranging from 1-3 and 2) an inflammation score of 1 or 2 - a high wound score correlates with impaired wound healing, e.g. a, wound score 5 with no fibrin coverage (3 points), marked peripheral hyperemia (1 point) and edema (1 point), thus representing delayed wound healing, and b; wound score 2 with complete fibrin coverage (1 point) and edema (1 point) representing timely wound healing.

Int J Biol Sci Image

(View in new window)

Results

Clinical characteristics and samples

Clinical data for both cohorts are detailed in Table 1. All samples from day zero were obtained without any signs of histological inflammation; Geboes score < 0.1 [22], confirming the clinical and endoscopic Mayo scores [21], and thus, ensuring deep mucosal healing in patients with UC.

 Table 1 

All patients included with ulcerative colitis (UC) were in deep remission (endoscopic Mayo score 0 and histological normal colonic mucosa, Geboes score < 0.1), and were categorized in accordance with previous extent of the disease. E1: proctitis, E2: left-sided colitis, E3: extensive colitis. EIM: extraintestinal manifestation.

Table 1Cohort 1Cohort 2
CharacteristicsUlcerative colitis
n = 8
Healthy controls
n = 8
Ulcerative colitis
n = 21
Healthy controls
n = 9
Gender (male/female)4/44/49/124/5
Age, years (mean, range)52 (33-67)24 (22-26)48 (19-73)47 (23-56)
Age at diagnosis (≤25/> 25 years)0/8na8/13na
Years with disease (≤10/> 10 years)7/1na12/9na
Extension (E1, E2 or E3)2/6/0na5/10/6na
Smoking/former/non-smoking0/0/80/0/80/6/150/0/9
EIM (present/not present)0/8na3/18na
Previous glucocorticoid (independent/dependent)8/0na14/7na
Current daily medication, n
Systemic 5-aminosalicylic acid (1.6-3.2 g)
5na18na
Current topical 5-aminosalicylic acid (1 g)0na1na
Current azathioprine (100-150 mg)0na4na
Previous infliximab (5 mg/kg)0na2na
None of the above medications3819

On average, samples had a wet weight of 4.89 mg, however, significant differences were found between the average weight of index biopsies and the subsequent cross-sectional wound biopsies (Figure S1).

Cohort 1, Modified Geboes scores

No difference in the modified Geboes scores was observed within each group when day one and day two were compared (Figure 2A). However, when comparing the two groups at day one and two, a clear trend was found with higher modified Geboes scores (at day 2, 2.2 vs. 1.1; p = 0.05) among patients with UC indicating a more intense inflammatory response to wounding. The histological evaluation also revealed that cells resembling WAEs did not appear until 48h post-wounding in both UC and controls (Figure 2B). These observations indicate that the early inflammatory response to injury might be more rigorous and sustained in patients with UC as previously described [23], and that the regenerative healing process extends beyond day two. Thus, to ensure that the actual wound healing process was investigated, and not just the inflammatory response caused by the wounding, subjects in the subsequent cohort 2 had sigmoidoscopies performed at day zero, two, and seven.

Cohort 2, Wound scores

The wound scores at day two after the artificial wounding were significantly higher in UC than controls (3 vs. 2; p=0.002), primarily due to insufficient fibrin-like coverage and peripheral hyperemia (Figure 2B). The significant difference in wound scores between UC and controls persisted at day seven (p<0.0001). At day seven, all patients with UC had achieved complete white fibrin-like coverage, but peripheral hyperemia and edema were, however, on-going characteristics.

Lipidomics data

Univariate analysis

A total of 98 phospholipid structures were detected and identified (Table S1). Comparing UC and controls, a significant lower amount of (16:0/18:2) phosphatidylglycerol (PG) was revealed at day seven. Further, a striking pattern was observed in the heatmap (Figure 3A) with significant differences within all examined classes of phospholipids at day two and seven as compared with day zero for UC. In contrast, only few and transient changes were observed in controls with no significant differences between day zero and seven. This pattern became especially evident when looking at the total amount of the different phospholipids measured (Figure 4) in which significant (p<0.05) changes were observed at day two within all phospholipid classes in UC. These changes were sustained at day seven, although only at a significant level for phosphatidic acid (PA) and LPA. The control group, on the other hand, only had a significant and transient reduction in phosphatidylcholine (PC).

A total of 33 eicosanoids were detected and identified (Table S2). Prominent and significant (p<0.05) reductions were observed within each group at day two and seven when compared to day zero, of which some were specific for UC, i.e., PGE1, 13,14-2OH-15keto-PGD2, PGJ2, 5(S)-hydroxyeicosapentanenoic acid (5(S)-HEPE), 13-hydroxyoctadecadienoic acid (13-HODE), and 14,15-dihydroxyeicosatetraenoic acid (14,15-di-HETE), whereas other lipid metabolites were specific for controls, i.e., 5(S)-HETE, 8-HETE, and 5(6)-epoxyeicosatrienoic acids (5(6)-EET). Nonetheless, no differences were noted when comparing UC and controls at the three time points (Figure 3B).

 Figure 2 

Wounds and histology. A) Modified Geboes and wound scores, median and interquartile range. The majority of subjects had a wound score of 3 (UC) or 2 (Control) at day 7 explaining the lack of range. B) Illustrative example of the histological appearance of the injury with acute inflammation. Black arrows, the edges of epithelial cells; green arrows, cuboidal regenerative epithelial cells; red arrows, flattened elongated wound associated epithelial (WAEs) cells. WAEs first appear after 48 hours.

Int J Biol Sci Image

(View in new window)

Multivariate analysis

PCA score plots from phospholipid and eicosanoid profiles containing all 90 samples identified five outliers. All five participants, from which these outliers originated, were reviewed in respect to demographics, clinical, and paraclinical parameters, including lipidomic data, but no obvious explanation could be identified, and the samples were consequently included in subsequent analyses.

OPLS-DA models based on phospholipids were valid at all three time points during wound healing in both UC and controls, and corresponding significant changes were identified in several phospholipids with both increased and decreased levels (Figure 5A-F, H). Especially UC had a high number of significantly (p<0.05) affected phospholipids, and the difference between UC and controls reached significance at day seven (Figure 5G), where levels of various species of phosphatidylserine (PS), phosphatidylinositol (PI), PG, PC, PA as well as LPA were significantly lower in patients with UC. The metabolic trajectory of this process is illustrated in Figure 5I, where patients with UC and controls had a similar starting point, but subsequently the paths diverged towards day two, and the return (i.e., day seven during wound healing) was notably different.

OPLS-DA models for the eicosanoids were valid when day two and seven were compared to day zero for both UC and controls (Figure 6A-E). Again, UC had a higher number of significantly (p<0.05) affected eicosanoids as compared to controls, but no significant differences were found between UC and controls. Thus, significant temporal changes in the levels of eicosanoids were identified during wounding and wound healing, but the metabolic trajectory did not differ significantly between UC and controls.

 Figure 3 

Heat maps of phospholipids and eicosanoids. Columns represent ulcerative colitis (UC) and healthy controls at day 0, 2, and 7, respectively. Rows represent each measured A) phospholipid or B) eicosanoid. Red and blue indicates relative increased or decreased concentrations.

Int J Biol Sci Image

(View in new window)

 Figure 4 

Measured total concentration of phospholipids (pmol/mg) during wound healing in ulcerative colitis and controls. Median and interquartile range. * Significant with a False Discovery Rate p<0.05. Lysophosphatidic acid (LPA), lysophosphatidylcholine (LPC), phosphatidic acid (PA), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), phosphatidylinositol (PI), phosphatidylserine (PS), sphingomyelin (SM), ulcerative colitis (UC).

Int J Biol Sci Image

(View in new window)

 Figure 5 

OPLS-DA models based on phospholipids. A-G: Phospholipids with correlation coefficients above 0.41 (UC) and 0.60 (controls) indicating significant differential power in the corresponding valid OPLS-DA models. Red, increased levels. Blue, decreased levels. H: Q2, cross-validation parameter indicating the predictability of the model. R2X, the fraction of the variation of the variables explained by the model. P-values (significant below 0.05) are based on the CV-ANOVA. I: Metabolic trajectory generated from the average of PCA scores of controls (orange color) and UC (blue) during wound healing. Lysophosphatidic acid (LPA), lysophosphatidylcholine (LPC), phosphatidic acid (PA), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), phosphatidylinositol (PI), phosphatidylserine.

Int J Biol Sci Image

(View in new window)

Int J Biol Sci Image

(View in new window)

 Figure 6 

OPLS-DA models based on eicosanoids. A-D: Eicosanoids with correlation coefficients above 0.41 (UC) and 0.60 (controls) indicating significant differential power in the corresponding valid OPLS-DA models. All measured eicosanoids were found at lower levels at day two and seven. E: Q2, cross-validation parameter indicating the predictability of the model. R2X, the fraction of the variation of the variables explained by the model. P-values (significant below 0.05) are based on the CV-ANOVA. Eicosatetraenoic acid (ETE), epoxy-docosapentaenoic acid (EpDPE), epoxyeicosatrienoic acids (EET), epoxy-octadecenoic acid (EpOME), hydroxyeicosapentanenoic acid (HEPE), hydroxyeicosatetraenoic acids (HETE), hydroxy-octadecadienoic acid (HODE), hydroperoxy-eicosapentanenoic acid (HpEPE), octadecadienoic acid (ODE), prostaglandin (PG), thromboxane (TX), ulcerative colitis (UC).

Int J Biol Sci Image

(View in new window)

 Figure 7 

Correlation between the wound score and phospholipids. Correlation analysis between wound scores (X-matrix, w*c[1]) and phospholipids (Y-matrix, w*c[2]) based on an OPLS-DA model and corresponding loadings. Q2, cross-validation parameter indicating the predictability of the model in relation to its statistical validity. R2, the fraction of the variation of the variables explained by the model. Lysophosphatidic acid (LPA), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), sphingomyelin (SM).

Int J Biol Sci Image

(View in new window)

Finally, a correlation analysis based on an OPLS-DA model and corresponding loadings was performed between wound scores (X-matrix) and phospholipids (Y-matrix) (Figure 7). As seen in this figure, wound scores were highly correlated with the phospholipid profile - especially with changes in sphingomyelin (SM) and phosphatidylethanolamine (PE) species - which were positively correlated with the wound scores, whereas levels of PG and LPA were negatively correlated. A similar model was built between wound scores and eicosanoids, but this model turned out to be invalid, i.e., failed cross-validation and the permutation test, and hence no correlation could be identified between the wound score and eicosanoids (data not shown).

Discussion

The human in vivo wound healing assay demonstrated that patients with even quiescent UC have an increased inflammatory response to intestinal wounding. This is in alignment with recently published data where inflammatory markers in an in vivo injury model showed rapid engagement of innate cytokines like interleukin (IL)- 1a/b and adaptive IL-17A along with recruitment of neutrophils and innate lymphoid cells in quiescent UC [23], and with a correlation between the histological grade of inflammatory infiltrate and wound scores. This UC-specific pattern of disrupted intestinal wounding and wound healing correlated in the current study with distinctive changes in the mucosal lipidome. Thus, both univariate, multivariate, and correlation analyses characterized the control group with only few and temporal changes contrasting the more profound and sustained chances in UC.

Previous studies have uncovered significant decreased amounts of LPC and PC in the colonic mucus of patients with UC [25-27], and human lipidomic studies on serum/plasma [28-30] as well as on colonic biopsies [31] have identified extensive alterations in the lipid composition with changed LPC and PC species. The low levels of LPC, and especially PC as seen in UC, might simply be due to an insufficient re-establishment of the mucus layer post-wounding, but might also be directly causative as different species and isomers of both LPC and PC act as pro- or anti-inflammatory stimuli [32]. A cell culture-based study suggests that PC inhibits the inflammatory response caused by tumor necrosis factor (TNF)-α through changes in the properties if the cell membrane, thus interfering with the signal transduction [33]. Human functional studies are, however, scarce, but a phase III placebo-controlled trial of a modified release of PC has shown to significantly alleviate disease activity in mesalazine-refractory UC [16], Whether this effect is a consequence of the anti-inflammatory properties; improvement of the mucus barrier, or both is, however, still not clarified.

A similar therapeutic effect has been achieved with local LPA in experimental colitis models by modulating intestinal epithelial migration, proliferation, and survival [17,34]. Oral application of LPA mitigates colonic wound closure and inhibits TNF-mediated intestinal barrier defects in mice [35,36]. No functional or interventional human studies are currently available, but in this study, LPA species are found at significant lower levels in UC at day seven (Figure 5G, Table S1), due to an aberrant inability to increase these LPA species as a response to wounding. LPA is produced through different pathways [37]; one involves the hydrolysis of the membrane derived PA, whereas another requires the elimination of the choline moiety from LPC and ultimately PC. Thus, the relatively low level of LPA in UC may be the consequence of transient or sustained suppressed levels of PA and PC, as demonstrated in the present study (Figure 5G), and might be part of the pathophysiology of compromised wound healing in UC.

Significantly suppressed levels of PS species are seen in UC at day seven compared to controls (Figure 5G). Along with its metabolically related membrane phospholipid PE, PS act as an essential metabolite in a balanced initial phase of wound healing, i.e., coagulation [38] and inflammation [39], a balance that seems to be off in UC wound healing as documented in the present study as well as in biopsy material from patients with flaring UC [31].

The majority of PG species remain unaffected by wounding in controls (Figure 5D-F, Table S1), which contrasts the findings in UC, where the PG level was significantly reduced compared to control (Figure 5G). Upon wounding, PG acts as a regulator of differentiation [40], and PG liposome treatment in murine wound models accelerates wound healing [41]. A somewhat similar pattern is observed with PI where most affected PI species are significantly decreased after wounding in UC (Figure 5G). Systemic or local treatment with PI inhibits T-cell proliferation and differentiation and ameliorates intestinal inflammation in a 2,4,6-trinitrobenzenesulfonic acid (TNBS)-induced colitis model [42,43]. Thus, these pronounced differences in PG and PI between controls and UC might explain the observed delay in wound healing observed in the present study. No human studies currently exist, but future therapies with supplementation of PG and PI could accelerate wound healing in UC by promotion of differentiation and anti-inflammatory effects.

Sphingomyelin, which is the most abundant sphingolipid and an important structural component of the cell membrane, was found in reduced levels, which correlated with high wound scores (Figure 7). Unfortunately, SM does not seem to offer an opportunity for therapeutic intervention, as it also has proinflammatory capabilities [44].

The majority of eicosanoids measured were significantly decreased at day two and/or seven in both UC and controls (Figure 3B and Figure 6A-D), but no differences could be found between day two and seven in UC or controls, or between UC and controls at any given time point (Figure 6E). The remaining eicosanoids were either unaffected or unmeasurable at several time points due to very low concentrations (Table S2).

Only few eicosanoids are well-described in relation to mucosal wound healing, and only so in experimental animal models. As previously described, PGE2 promotes early wound repair by inducing WAE cells in mice [10,11], whereas a subsequent increase in PGD2 and its metabolites function as a part of the anti-inflammatory response [45]. This does, however, not seem to apply to humans, as we find PGE2 and PGD2, including its metabolites, i.e., PGJ2, and 15-deoxy-Δ12,14-PGD2, to be continuously reduced during regeneration, especially in UC (Figure 6A, B). These findings are novel and detrimental to what is known from experimental animal models, which stress the importance and need for human experiments, and obviously merits further analyses and functional tests in future human studies. However, it could be hypothesized that this might represent a point of therapeutic intervention, as animal models suggest that PGD2 metabolites ameliorates inflammation [46]. Similarly, the sustained reduction of PGE1, 13-HODE, and 13-oxo-ODE in UC, which are all anti-inflammatory eicosanoids and improve the microcirculation [47], may be substituted, as animal [48] and human [49] studies have shown increased wound healing times with the infusion of PGE1.

The disease specific changes found in lipid abundance could be driven by differences in mucosal cell type composition during wounding and regeneration or by changes in the mucosa-near microbiome and cellular lipid metabolism. Accumulating evidence suggests a microbe-host-lipid co-metabolism [50], which seems especially evident for the eicosanoids [51]. Our previous data [23] indicates that the mucosal niche microbiota are changed during colonic mucosal injury and closely resembles the changes seen in active UC in terms of reduced α-diversity, bacterial composition, and bacterial load. These changes are similar in both patients with UC and healthy controls, which might explain why the majority of eicosanoids are reduced in the current study, and equally so in both UC and controls. However, this needs to be elucidated further via functional assays, just as the proposed potential therapeutic targets need to be tested and validated in interventional studies. Moreover, adjunctive comprehensive analyses of the wound niche microbiome, the transcriptome and proteome along with further lipidomic analyses could add significantly insight to the assessment of the wound healing model. Especially the role of specialized pro-resolving lipid mediators derived from polyunsaturated fatty acids like docosahexaenoic acid and eicosapentaenoic acid, which recently have been linked to protective actions in intestinal inflammation and to promote wound healing [14], might be of significant interest in future lipidomic studies.

It is important to note that the applied wound healing assay does not mimic the actual inflammation during regular flares of UC, but it is considered a tool to investigate the initial response to mucosal injury and subsequent mucosal healing process. Finally, no age-match was made between patients with UC and healthy control, nor were any restrictions imposed in terms of medication except that it had to be stable for at least 3 months. The latter needs to be kept in mind when interpreting the results, as the majority of patients with UC received maintenance treatment with 5-aminosalicylic acid, which may activate peroxisome proliferator-activated receptors with multiple effects on lipid metabolism, inflammatory processes, and cell proliferation [52].

Conclusions

In summary, 1) this study describes for the first time the kinetics of disrupted metabolic pathways in the wound healing process in the colon mucosa among patients with quiescent UC compared to healthy individuals. Thus, the analyses identify a persistent inflammatory response and delayed epithelial restitution during wound healing in UC with underlying UC-specific perturbations in lipid mediator engagement. Moreover, 2) the lipidomic characterization presented identifies significant temporal changes of several phospholipids and eicosanoids with significant differences to what is found in comparable animal models. Ultimately, 3) these results challenge our traditional treatment approach of IBD and may guide future medical treatment regimens. Therapeutic interventions ameliorating the low levels of especially LPC, PC, LPA PG, PI, PGD2, and PGE1 in inflamed colonic tissue may represent a novel potential therapeutic avenue for recovering flaring UC and subsequent epithelial healing, which, however, requires future functional and interventional human studies. In this way lipidomics may provide useful novel insights into specific disease phenotypes, potentially leading towards improved clinical practice, which will enable the development of personalized medicine in a wide range of diseases, including UC.

Abbreviations

BHT: butylated hydroxytoluene; CRC: colorectal cancer; CV-ANOVA: analysis of variance of the cross-validated residuals; EET: epoxyeicosatrienoic acids; EpDPE: epoxy-docosapentaenoic acid; EpOME: epoxy-octadecenoic acid; ESI: electrospray ionization; ETE: eicosatetraenoic acid; FDR: false discovery rate; HEPE: hydroxyeicosapentanenoic acid; HETE: hydroxyeicosatetraenoic acids; HODE: hydroxy-octadecadienoic acid; HpEPE: hydroperoxy-eicosapentanenoic acid; IBD: inflammatory bowel disease; IQR: interquartile range; LPA: lysophosphatidic acid; LPC: lysophosphatidylcholine; LPE: lysophosphatidylethanolamine; LC-MS: liquid chromatography mass spectrometry; LT: leukotriene; MRM: multiple reaction-monitoring; ODE: octadecadienoic acid; OPLS-DA: orthogonal projection to latent structure-discriminant analysis; PA: phosphatidic acid; PC: phosphatidylcholine; PCA: principal component analysis; PE: phosphatidylethanolamine; PG: phosphatidylglycerol; PGA/D/E/J: prostaglandin A/D/E/J; PI: phosphatidylinositol; PS: phosphatidylserine; SM: sphingomyelin; TNBS: trinitrobenzenesulfonic acid; TNF: tumor necrosis factor; TX: thromboxane; UC: ulcerative colitis; UV: unit variance; UPLC-MS: ultraperformance liquid chromatography-mass spectrometry; WAE: wound associated epithelial.

Supplementary Material

Supplementary figures and tables, information.

Attachment

Acknowledgements

Funding

This study received financial support from Aase og Ejnar Danielsens Fond, Civilingeniør Frode V. Nygaard og Hustrus Fond, Aage og Johanne Louis-Hansens Fond, Colitis-Crohn Foreningen, and Frimodt-Heineke Fonden. The funders were not involved in the study in any way.

Author Contribution

Study design J.T.B., J.B.S. Patient recruitment and sample processing J.T.B., J.B.S. Experimental work Y.W., J.Z., L.B.R, O.H.N. Data analysis, Y.W., J.Z., J.T.B. First draft of the manuscript J.T.B. All authors were involved in the subsequent critical review, editing, revision, and approval of the final manuscript.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Kobayashi T, Siegmund B, Le Berre C, Wei SC, Ferrante M, Shen B. et al. Ulcerative colitis. Nat Rev Dis Prim. 2020;6:74

2. Windsor JW, Kaplan GG. Evolving epidemiology of IBD. Curr Gastroenterol Rep. 2019;21:40

3. Le Berre C, Ananthakrishnan AN, Danese S, Singh S, Peyrin-Biroulet L. Ulcerative colitis and Crohn's disease have similar burden and goals for treatment. Clin Gastroenterol Hepatol. 2020;18:14-23

4. Rungoe C, Langholz E, Andersson M, Basit S, Nielsen NM, Wohlfahrt J. et al. Changes in medical treatment and surgery rates in inflammatory bowel disease: a nationwide cohort study 1979-2011. Gut. 2014;63:1607-16

5. Olén O, Erichsen R, Sachs MC, Pedersen L, Halfvarson J, Askling J. et al. Colorectal cancer in ulcerative colitis: a Scandinavian population-based cohort study. Lancet. 2020;395:123-31

6. Dal Buono A, Roda G, Argollo M, Peyrin-Biroulet L, Danese S. Histological healing: should it be considered as a new outcome for ulcerative colitis?. Expert Opin Biol Ther. 2020;20:407-12

7. Brazil JC, Quiros M, Nusrat A, Parkos CA. Innate immune cell-epithelial crosstalk during wound repair. J Clin Invest. 2019;129:2983-93

8. Bjerrum JT, Wang YL, Seidelin JB, Nielsen OH. IBD metabonomics predicts phenotype, disease course, and treatment response. EBioMedicine. 2021;71:103551

9. O'Donnell VB, Ekroos K, Liebisch G, Wakelam M. Lipidomics: current state of the art in a fast moving field. Wiley Interdiscip Rev Syst Biol Med. 2020;12:e1466

10. Miyoshi H, VanDussen KL, Malvin NP, Ryu SH, Wang Y, Sonnek NM. et al. Prostaglandin E2 promotes intestinal repair through an adaptive cellular response of the epithelium. EMBO J. 2017;36:5-24

11. Manieri NA, Drylewicz MR, Miyoshi H, Stappenbeck TS. Igf2bp1 is required for full induction of Ptgs2 mRNA in colonic mesenchymal stem cells in mice. Gastroenterology. 2012;143:110-21

12. Li Y, Soendergaard C, Bergenheim FH, Aronoff DM, Milne G, Riis LB. et al. COX-2-PGE2 signaling impairs intestinal epithelial regeneration and associates with TNF inhibitor responsiveness in ulcerative colitis. EBioMedicine. 2018;36:497-507

13. Diab J, Al-Mahdi R, Gouveia-Figueira S, Hansen T, Jensen E, Goll R. et al. 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. 2019;25:490-7

14. Quiros M, Feier D, Birkl D, Agarwal R, Zhou DW, García AJ. et al. Resolvin E1 is a pro-repair molecule that promotes intestinal epithelial wound healing. Proc Natl Acad Sci U S A. 2020;117:9477-82

15. Olorundare OE, Peyruchaud O, Albrecht RM, Mosher DF. Assembly of a fibronectin matrix by adherent platelets stimulated by lysophosphatidic acid and other agonists. Blood. 2001;98:117-24

16. Karner M, Kocjan A, Stein J, Schreiber S, von Boyen G, Uebel P. et al. First multicenter study of modified release phosphatidylcholine “LT-02” in ulcerative colitis: a randomized, placebo-controlled trial in mesalazine-refractory courses. Am J Gastroenterol. 2014;109:1041-51

17. Sturm A, Sudermann T, Schulte KM, Goebell H, Dignass AU. Modulation of intestinal epithelial wound healing in vitro and in vivo by lysophosphatidic acid. Gastroenterology. 1999;117:368-77

18. Bjerrum JT, Steenholdt C, Ainsworth M, Nielsen OH, Reed MAC, Atkins K. et al. Metabonomics uncovers a reversible proatherogenic lipid profile during infliximab therapy of inflammatory bowel disease. BMC Med. 2017;15:184

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

20. Bjerrum JT, Nielsen OH, Hao F, Tang H, Nicholson JK, Wang Y. et al. Metabonomics in ulcerative colitis: diagnostics, biomarker identification, and insight into the pathophysiology. J Proteome Res. 2010;9:954-62

21. Schroeder KW, Tremaine WJ, Ilstrup DM. Coated oral 5-aminosalicylic acid therapy for mildly to moderately active ulcerative colitis. N Engl J Med. 1987;317:1625-9

22. Geboes K, Riddell R, Ost A, Jensfelt B, Persson T, Löfberg R. A reproducible grading scale for histological assessment of inflammation in ulcerative colitis. Gut. 2000;47:404-9

23. Seidelin JB, Bahl MI, Licht TR, Mead BE, Karp JM, Johansen JV. et al. Acute experimental barrier injury triggers ulcerative colitis-specific innate hyperresponsiveness and ulcerative colitis-type microbiome changes in humans. Cell Mol Gastroenterol Hepatol. 2021;12:1281-96

24. Wolfer AM, Gaudin M, Taylor-Robinson SD, Holmes E, Nicholson JK. Development and validation of a high-throughput ultrahigh-performance liquid chromatography-mass spectrometry approach for screening of oxylipins and their precursors. Anal Chem. 2015;87:11721-31

25. Braun A, Treede I, Gotthardt D, Tietje A, Zahn A, Ruhwald R. et al. Alterations of phospholipid concentration and species composition of the intestinal mucus barrier in ulcerative colitis: a clue to pathogenesis. Inflamm Bowel Dis. 2009;15:1705-20

26. Ehehalt R, Wagenblast J, Erben G, Lehmann WD, Hinz U, Merle U. et al. Phosphatidylcholine and lysophosphatidylcholine in intestinal mucus of ulcerative colitis patients. A quantitative approach by nanoElectrospray-tandem mass spectrometry. Scand J Gastroenterol. 2004;39:737-42

27. Ding H, Dupont AW, Singhal S, Scott LD, Guha S, Younes M. et al. In vivo analysis of mucosal lipids reveals histological disease activity in ulcerative colitis using endoscope-coupled Raman spectroscopy. Biomed Opt Express. 2017;8:3426

28. Murgia A, Hinz C, Liggi S, Denes J, Hall Z, West J. 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

29. Daniluk U, Daniluk J, Kucharski R, Kowalczyk T, Pietrowska K, Samczuk P. et al. Untargeted metabolomics and inflammatory markers profiling in children with Crohn's disease and ulcerative colitis—a preliminary study. Inflamm Bowel Dis. 2019;25:1120-8

30. Guan S, Jia B, Chao K, Zhu X, Tang J, Li M. et al. UPLC-QTOF-MS-based plasma lipidomic profiling reveals biomarkers for inflammatory bowel disease diagnosis. J Proteome Res. 2020;19:600-9

31. Diab J, Hansen T, Goll R, Stenlund H, Ahnlund M, Jensen E. et al. Lipidomics in ulcerative colitis reveal alteration in mucosal lipid composition associated with the disease state. Inflamm Bowel Dis. 2019;25:1780-7

32. Liu P, Zhu W, Chen C, Yan B, Zhu L, Chen X. et al. The mechanisms of lysophosphatidylcholine in the development of diseases. Life Sci. 2020;247:117443

33. Treede I, Braun A, Sparla R, Kühnel M, Giese T, Turner JR. et al. Anti-inflammatory effects of phosphatidylcholine. J Biol Chem. 2007;282:27155-64

34. Sturm A, Zeeh J, Sudermann T, Rath H, Gerken G, Dignass AU. Lisofylline and lysophospholipids ameliorate experimental colitis in rats. Digestion. 2002;66:23-9

35. Lee S-J, Leoni G, Neumann P-A, Chun J, Nusrat A, Yun CC. Distinct phospholipase C-β isozymes mediate lysophosphatidic acid receptor 1 effects on intestinal epithelial homeostasis and wound closure. Mol Cell Biol. 2013;33:2016-28

36. Lin S, Han Y, Jenkin K, Lee S-J, Sasaki M, Klapproth J-M. et al. Lysophosphatidic acid receptor 1 is important for intestinal epithelial barrier function and susceptibility to colitis. Am J Pathol. 2018;188:353-66

37. Yun CC, Kumar A. Diverse roles of LPA signaling in the intestinal epithelium. Exp Cell Res. 2015;333:201-7

38. Bevers EM, Williamson PL. Getting to the outer leaflet: physiology of phosphatidylserine exposure at the plasma membrane. Physiol Rev. 2016;96:605-45

39. Elliott MR, Ravichandran KS. Clearance of apoptotic cells: implications in health and disease. J Cell Biol. 2010;189:1059-70

40. Qin H, Zheng X, Zhong X, Shetty AK, Elias PM, Bollag WB. Aquaporin-3 in keratinocytes and skin: its role and interaction with phospholipase D2. Arch Biochem Biophys. 2011;508:138-43

41. Arun SN, Xie D, Howard AC, Zhong Q, Zhong X, McNeil PL. et al. Cell wounding activates phospholipase D in primary mouse keratinocytes. J Lipid Res. 2013;54:581-91

42. van Dieren JM, Simons-Oosterhuis Y, Raatgeep HCR, Lindenbergh-Kortleve DJ, Lambers MEH, van der Woude CJ. et al. Anti-inflammatory actions of phosphatidylinositol. Eur J Immunol. 2011;41:1047-57

43. Fabia R, Ar'Rajab A, Willén R, Andersson R, Ahrén B, Larsson K. et al. Effects of phosphatidylcholine and phosphatidylinositol on acetic-acid-induced colitis in the rat. Digestion. 1992;53:35-44

44. Fischbeck A, Leucht K, Frey-Wagner I, Bentz S, Pesch T, Kellermeier S. et al. Sphingomyelin induces cathepsin D-mediated apoptosis in intestinal epithelial cells and increases inflammation in DSS colitis. Gut. 2011;60:55-65

45. Ajuebor MN, Singh A, Wallace JL. Cyclooxygenase-2-derived prostaglandin D2 is an early anti-inflammatory signal in experimental colitis. Am J Physiol Liver Physiol. 2000;279:G238-44

46. Li J, Guo C, Wu J. 15-deoxy-∆- 12,14 -prostaglandin J2 (15d-PGJ2), an endogenous ligand of PPAR-γ: function and mechanism. PPAR Res. 2019;2019:1-10

47. Braune S, Küpper J-H, Jung F. Effect of prostanoids on human platelet function: an overview. Int J Mol Sci. 2020;21:9020

48. Ruiz-Luque V, Parra-Membrives P, Escudero-Severín C, Aguilar-Luque J. Effect of alprostadil on colorectal anastomoses under relative ischemia. J Surg Res. 2019;236:230-7

49. De Caridi G, Massara M, Stilo F, Spinelli F, Grande R, Butrico L. et al. Effectiveness of prostaglandin E1 in patients with mixed arterial and venous ulcers of the lower limbs. Int Wound J. 2016;13:625-9

50. Lamichhane S, Sen P, Alves MA, Ribeiro HC, Raunioniemi P, Hyötyläinen T. et al. Linking gut microbiome and lipid metabolism: moving beyond associations. Metabolites. 2021;11:55

51. Wallace JL. Eicosanoids in the gastrointestinal tract. Br J Pharmacol. 2019;176:1000-8

52. Caioni G, Viscido A, D'Angelo M, Panella G, Castelli V, Merola C. et al. Inflammatory bowel disease: new insights into the interplay between environmental factors and PPARγ. Int J Mol Sci. 2021;22:985

Author contact

Corresponding address Corresponding authors: J.T. Bjerrum, Department of Gastroenterology 54R1M, Herlev Hospital, Denmark, Borgmester Ib Juuls Vej 1, DK-2730 Herlev, Denmark. Mail: bjerrum.jacobcom or Y. Wang, Singapore Phenome Centre, Lee Kong Chian School of Medicine, Nanyang Technological University, 59 Nanyang Drive, Level 5, Singapore 636921. Mail: yulan.wangedu.sg.


Received 2021-9-17
Accepted 2022-1-10
Published 2022-2-14