Next Article in Journal
Role of HMGB1 in the Interplay between NETosis and Thrombosis in Ischemic Stroke: A Review
Next Article in Special Issue
Probiotic Lactobacillus and Bifidobacterium Strains Counteract Adherent-Invasive Escherichia coli (AIEC) Virulence and Hamper IL-23/Th17 Axis in Ulcerative Colitis, but Not in Crohn’s Disease
Previous Article in Journal
The Emerging Role of Galectins and O-GlcNAc Homeostasis in Processes of Cellular Differentiation
Previous Article in Special Issue
Interplay between Cellular and Molecular Mechanisms Underlying Inflammatory Bowel Diseases Development—A Focus on Ulcerative Colitis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 9
Issue 8
10.3390/cells9081793
cells-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

Molecular Changes in the Non-Inflamed Terminal Ileum of Patients with Ulcerative Colitis

by
Ho-Su Lee
1,2,†,
Maaike Vancamelbeke
3,†,
Sare Verstockt
1,
Tom Wilms
1,
Bram Verstockt
3,4,
João Sabino
3,4,
Marc Ferrante
3,4,
Séverine Vermeire
3,4 and
Isabelle Cleynen
1,*
1
Laboratory of Complex Genetics, Department of Human Genetics, KU Leuven, 3000 Leuven, Belgium
2
Department of Biochemistry and Molecular Biology, University of Ulsan College of Medicine, Seoul 05505, Korea
3
Translational Research Center for Gastrointestinal Disorders (TARGID), Department Chronic Diseases, Metabolism & Ageing (CHROMETA), KU Leuven, 3000 Leuven, Belgium
4
Department of Gastroenterology and Hepatology, University Hospitals Leuven, KU Leuven, 3000 Leuven, Belgium
*
Author to whom correspondence should be addressed.
†
These authors have contributed equally to this work.
Cells 2020, 9(8), 1793; https://doi.org/10.3390/cells9081793
Submission received: 22 June 2020 / Revised: 22 July 2020 / Accepted: 24 July 2020 / Published: 28 July 2020
(This article belongs to the Special Issue The Molecular and Cellular Basis for Inflammatory Bowel Diseases (IBD) 2020)
Browse Figures
Versions Notes

Abstract

:
Ulcerative colitis is a chronic inflammatory disease confined to the colon. Although the etiopathogenesis remains unknown, small bowel dysfunctions like histological and permeability alterations have been described in ulcerative colitis. We evaluated the molecular gene signature in the non-inflamed terminal ileum of 36 ulcerative colitis patients (7 active, with Mayo endoscopic subscore ≥2, and 29 inactive) as compared to 15 non-inflammatory bowel disease controls. Differential gene expression analysis with DESeq2 showed distinct expression patterns depending on disease activity and maximal disease extent. We found 84 dysregulated genes in patients with active extensive colitis and 20 in inactive extensive colitis, compared to controls. There was an overlap of 5 genes: REG1B, REG1A, MUC4, GRAMD2, and CASP10. In patients with left-sided colitis, ileal gene expression levels were similar to controls. Based on gene co-expression analysis, ileal changes in active ulcerative colitis patients were related to immune functions. The ileal changes in the inactive ulcerative colitis subjects converged into the maintenance of the intestinal barrier through increased mitochondrial function and dampened immune functions. In conclusion, we identified molecular changes in the non-inflamed ileum of ulcerative colitis that are dependent on colonic inflammation.

1. Introduction

Inflammatory bowel diseases (IBD), comprising ulcerative colitis (UC) [1] and Crohn’s disease (CD), are chronic disorders of the gastrointestinal tract, characterized by relapsing gut inflammation [2]. UC is limited to the colon and is characterized by a continuous mucosal inflammation, whereas CD can affect all parts of the gastrointestinal tract and is determined by a non-continuous transmural inflammation [2]. From a pathophysiological point of view, there is currently no explanation for the spatial restriction of the inflammation in UC [3]. Small bowel dysfunction has, however, been reported in UC [3]. For example, histological alterations including structural alterations of villi and inflammatory cell infiltration [1,4,5,6,7], and increased intestinal permeability, have been detected in the small bowel of patients with UC [8,9]. There is also experimental evidence of small bowel dysfunction from animal models of colitis. During experimental colitis, alterations in small bowel permeability and brush border enzymatic activity are observed [7,10,11]. Moreover, in IL10-deficient mice which spontaneously develop colitis, colitis is mitigated by reducing small intestinal permeability upon treatment with the zonulin peptide inhibitor AT-1001 [12]. This again suggests a possible link between small bowel dysfunction and colonic inflammation. To date, however, the precise mechanisms underlying this dysfunction have not been elucidated.
The intestinal transcriptomes of IBD patients have been previously analysed to gain insights into the molecular mechanisms leading to disease development [13,14]. Likewise, genome-wide expression profiling of unaffected terminal ileum tissue from patients with UC can provide better insights into the molecular characteristics of small bowel dysfunction in these patients. To the best of our knowledge, only two studies report transcriptomic data of unaffected terminal ileum in UC, albeit on at least partially overlapping datasets, and with the primary goals not directed towards finding the molecular alterations in UC terminal ileum [15,16]. They found significant alterations in gene expression levels in the terminal ileum of UC patients, with Dual Oxidase Maturation Factor 2 (DUOXA2) and Dual Oxidase 2 (DUOX2) being among the most significant [15,16]. The DUOX2/DUOXA2 complex plays an important role in maintaining mucosal immune homeostasis [17,18,19] and in the host-microbial interaction in IBD [16]. They did not investigate whether these alterations in the non-inflamed small bowel are related to the disease activity and/or disease extent of UC.
We performed RNA sequencing on mucosal non-inflamed terminal ileum tissue samples from active and inactive UC patients, as well as non-IBD controls, to identify if and which genes and gene co-expression modules were dysregulated. The two primary aims of our study were (i) to investigate the molecular alterations in non-inflamed ileum tissue in UC, and (ii) to determine whether these alterations are dependent on colonic inflammation. The latter was assessed in terms of disease activity and disease extent.

2. Materials and Methods

2.1. Study Subjects and Samples

Patients with an established diagnosis of UC [20] were included. Terminal ileum mucosal biopsies were obtained during endoscopy from macroscopically unaffected sites in 37 UC patients. Colonic inflammation activity at endoscopy was determined in accordance with the Mayo endoscopic subscore, with active colitis defined by a value of 2 or 3 [21]. The maximal UC disease extent during follow-up and at the time of endoscopy was categorized in accordance with the Montreal classification and based on macroscopic findings [22]. Sixteen non-IBD controls were included for comparison. These underwent endoscopy for polyp screening, and had normal mucosa without evidence of gut inflammation during endoscopy. The ethics committee of the University Hospitals Leuven approved the study on 29 March 2012 (Approval number: B322201213950/S53684). All patients provided written informed consent.

2.2. RNA Isolation and Sequencing

Ileal biopsies were immediately placed in RNAlater (Ambion, Austin, TX, USA). The total RNA was extracted using the AllPrep DNA/RNA Mini kit (Qiagen, Hilden, Germany) in accordance with the manufacturer’s instructions. RNA integrity and quantity were assessed using a 2100 Bioanalyzer (Agilent, Waldbronn, Germany) and a Nanodrop ND-1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). Extracted RNA was stored at −80 °C until further processing. Single-end RNA sequencing was conducted using an Illumina HiSeq 4000 system (Illumina, San Diego, CA, USA) in accordance with the manufacturer’s instructions. Raw RNA sequencing data were aligned to the reference genome using Hisat2 version 2.1.0 [23], and absolute counts were generated using HTSeq [24]. Non-coding transcripts or low-expressed transcripts with an average count ≤10 were excluded from further analysis, leaving 14,025 transcripts. Principal component analysis (PCA) was used to identify outliers and to check for clustering of samples. We excluded two outlier samples (one UC patient and one control) based on PCA (Supplementary Figure S1A), leaving 36 UC patients (61% male; median age 52 years) and 15 non-IBD controls (40% male; median age 55 years) for further analysis (Table 1 and Supplementary Figure S1B). PCA was performed using the FactoMineR package [25]. Grouping phenotypes (case versus non-IBD control, with and without cases broken down in disease activity and/or disease extent) were added as supplementary variables in the analysis, in order to calculate their correlation with the obtained expression-based PCs.

2.3. Differential Gene Expression Analysis

Absolute counts were analysed for differential gene expression between cases and controls using the DESeq2 (version 1.20.0) package [26] in R (version 3.5.1), and adjusted for age and sex. Differentially expressed genes were defined by an absolute log2 fold-change ≥1 (i.e., fold-change ≥2) and a corrected p-value (Benjamini–Hochberg method [27]) below 0.05. We reported corrected p-values for differential gene expression analysis unless noted otherwise.

2.4. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

To validate RNA-seq data, gene expression for selected differentially expressed genes in ileal mucosal biopsies was studied trough a qRT-PCR analysis. One sample of active extensive UC was unavailable for qRT-PCR. cDNA was synthesised from 0.5 μg RNA using the RevertAid H Minus First Strand cDNA synthesis kit (Fermentas, St. Leon-Rot, Germany), following the manufacturer’s protocol. To facilitate the analysis of multiple targets, we preamplified the cDNA (1:10 dilution) with a multiplex cDNA preamplification using the TaqMan PreAmp Master Mix kit (Applied Biosystems, Foster City, CA), according to the manufacturer’s protocol. The selected genes included DUOXA2, DUOX2, REG1B, REG1A, MUC4, GRAMD2, and CASP10, with ACTB and GAPDH as endogenous reference genes. An overview of the validated genes and the corresponding TaqMan Assay ID are given in Supplementary Table S1 (Thermo Fisher Scientific, Massachusetts, MA). qPCR was performed in a final reaction volume of 10 µL on a QuantStudioTM 12K Flex Real-Time PCR system (Thermo Fisher Scientific, Massachusetts, MA). All samples were analysed in triplicate. mRNA levels were normalized to the geometric mean of the two housekeeping genes (ACTB and GAPDH) and quantified using the comparative (ΔΔ) Ct method.

2.5. Weighted Gene Co-Expression Network Analysis

To identify modules of co-expressed genes, we applied the Weighted Gene Correlation Network Analysis (WGCNA, version 1.66) [28] using R (version 3.5.1). We used the same design as applied for differential gene expression analysis for data input for WGCNA: briefly, we normalized the counts and corrected for age and sex using the DESeq2 package. We constructed co-expression modules by calculating a matrix of the Pearson correlations between all gene pairs. To identify modules of highly co-regulated genes, we used average linkage hierarchical clustering to group genes based on the topological overlap of their connectivity, followed by a dynamic cut-tree algorithm to dynamically separate clustering dendrogram branches into gene modules (using deepSplit = 2 and a minimal module size of 30). A soft-threshold power value of 4 was selected in this analysis, which corresponds to a 87.2% scale-free topology. Each module was assigned a unique colour identifier.
To identify which modules are most relevant in the context of colonic inflammation, each module was tested for correlation with phenotypes of interest (disease activity, disease extent…). For this analysis, each module is represented by its eigengene, i.e., the first principal component of the expression patterns of all genes within a given module, summarized into a single characteristic. We calculated the corrected p-value using the Benjamini–Hochberg method [27].

2.6. Pathway Analysis

Ingenuity Pathway Analysis (IPA, QIAGEN Inc., www.qiagen.com/ingenuity) and the Reactome pathway database (www.reactome.org, version 68) were used to identify upstream regulators and canonical pathways associated with the differentially expressed genes, and with genes within each WGCNA module. To further evaluate the biological relevance of the differentially expressed genes and the genes within the modules, we also performed a gene ontology (GO) enrichment analysis of biological processes using PANTHER tools version 15.0 (pantherdb.org/tools/) [29]. A corrected p-value < 0.05 was considered as significantly enriched.

2.7. Statistical Analysis

Continuous variables were expressed as median and interquartile range (IQR), and categorical variables as frequencies and percentages. Data were compared using the Mann–Whitney U test for continuous variables, and with Fisher’s exact or Chi2 test for categorical variables. A p-value < 0.05 was considered statistically significant. Statistical analyses were done using R (version 3.5.1).

3. Results

3.1. Baseline Characteristics

The baseline characteristics of the included patients and controls are listed in Table 1. Of the 36 UC patients, 7 (19.4%) had active colonic disease at the time of endoscopy (5 with an endoscopic Mayo score of 2, and 2 with an endoscopic Mayo score of 3), and 29 (80.6%) were inactive. The maximal extent of disease throughout follow-up in the active UC cases was left-sided colitis in 2 patients (28.6%) and extensive colitis in 5 patients (71.4%). At the time of endoscopy, among the active UC cases, 4 patients had left-sided colitis, and 3 had extensive colitis. The maximal disease extent among the inactive UC cases included 4 cases of proctitis, 11 left-sided colitis, and 14 extensive colitis. At time of endoscopy, 1 active UC patient (14.3%) and 13 inactive UC patients (44.8%) received biologics.

3.2. Gene Dysregulation in the Normal Terminal Ileum of UC Patients is Dependent on Colonic Inflammation

Based on PCA, which characterized the general gene expression differences among the 51 samples, we did not see any obvious clusters in the data (Supplementary Figure S1B). The first two principal components explained 35.6% of variance between the samples (26.3% for PC1, 9.3% for PC2). The highest correlation with PC1 was seen for the combination of disease activity and extent (r2 = 0.14, p = 1.21 × 10−1), while disease extent (irrespective of disease activity) was significantly correlated with PC2 (r2 = 0.13, p = 3.18 × 10−2). This suggested disease activity and/or disease extent contribute to the gene expression pattern in the terminal ileum in UC, albeit a small one.
We then separately compared gene expression in normal ileum of active and inactive UC patients with non-IBD controls (Figure 1, Table 2, Supplementary Figure S2A, and Supplementary Table S2). In active UC, we identified 18 differentially expressed genes, with DUOXA2 being the most significant (log2 fold-change = 4.9, p = 5.58 × 10−3). The cytokines IL22, STAT1, and IFN-γ were the top 3 predicted upstream regulators of some of these 18 genes: IL22 for MUC1 and NOS2; STAT1 for DUOX2, MUC1, and NOS2; and IFN-γ for DUOX2, DUOXA2, MUC1, and NOS2 (Supplementary Table S3). In inactive UC, we identified two dysregulated genes compared to non-IBD controls: REG1B (log2 fold-change = 2.8, p = 3.42 × 10−2) and CEBPD (log2 fold-change = 1.1, p = 2.12 × 10−2). REG1B was also one of the 18 dysregulated genes found in active UC, but with more pronounced dysregulation in active than in inactive UC (active UC: log2 fold-change = 4.2, p = 1.98 × 10−2; inactive UC: log2 fold-change = 2.8, p = 3.42 × 10−2) (Table 2, Supplementary Figure S2A).
We next studied if gene expression profiles in the terminal ileum depend on maximal disease extent. We therefore separately performed a differential gene expression analysis for extensive UC and left-sided colitis/proctitis, in active and inactive UC (Figure 1, Table 3, and Supplementary Table S2). We found 84 differentially expressed genes in active extensive UC compared to non-IBD controls, and 20 in inactive extensive UC, with an overlap of 5 genes (REG1B, REG1A, MUC4, GRAMD2, and CASP10; Figure 1, Table 3, Supplementary Figure S2B, and Supplementary Table S2). In line with the RNA sequencing results, qRT-PCR validation found significantly increased expression levels of each of these 5 genes in active extensive UC (Figure 2), as well as for DUOXA2 and DUOX2, which were also among the top dysregulated genes for active UC (Table 2). In inactive extensive UC, increased REG1B and MUC4 expression levels were validated using qRT-PCR (Figure 2). Although REG1A, GRAMD2, and CASP10 genes could not be confirmed with qRT-PCR, we observed borderline significance for REG1A and GRAMD2 (p-value = 0.09 and 0.13, respectively), and the observed p-values for all genes where concordant with their rank observed for RNA sequencing data (Table 3).
IL1A, NFkB (complex), and PPRC1 were the top 3 predicted upstream regulators of 7, 8, and 4 out of 84 genes, respectively (Supplementary Table S3). The top 3 predicted upstream regulators in inactive extensive colitis were miR-155-5p, SOCS1, and TREM1, regulating CXCL2 and SOCS1; CXCL2 and SOCS1; and CXCL2, HES4, and TIFA, respectively (Supplementary Table S3). Given the number of dysregulated genes, we only performed pathway analyses on the 84 differentially expressed genes in active extensive UC. Among these 84 genes, 6 genes were involved in “Endothelin-1 Signalling” (p = 1.51 × 10−4), and 5 in “Granulocyte Adhesion and Diapedesis” (p = 5.50 × 10−4) (Supplementary Table S4). Significantly enriched biological processes among the dysregulated 84 genes were the regulation of inflammatory response, cell wall disruption in another organism, and antimicrobial humoral immune response mediated by antimicrobial peptide, hereby implicating the role of the antimicrobial response in the ileum of UC patients (Supplementary Table S4).
Only LCT (Lactase) was dysregulated in active UC with left-sided colitis (log2 fold-change = −5.3, p = 2.25 × 10−2, Supplementary Figure S2C). It should however be noted that these were only two patients. No significant differences were seen between inactive UC patients with left-sided colitis/proctitis (n = 15) and non-IBD controls (Supplementary Figure S2C). The analysis of extensive UC at the time of endoscopy in active UC patients showed similar findings both for extensive UC and maximal disease extent.

3.3. Co-Expression Modules Are Associated with Colonic Inflammation

Besides looking at the differential gene expression of individual genes, we assessed whether we could find modules of genes with a similar ileal expression pattern. The WGCNA analysis found 39 co-expression modules (Supplementary Figure S3) that ranged in size from 60 to 1194 genes. None of these modules were significantly correlated with the studied traits after correction for multiple testing (corrected p-value < 0.05). Six modules showed a nominal significant (uncorrected p-value < 0.05) correlation with disease: three with inactive UC and three with active UC (Supplementary Figure S3). Of the three modules, one had a significant correlation with inactive UC, while two were positively correlated (colour code: purple and grey60) and included genes involved in mitochondrial translation and actin cytoskeleton, respectively; and one was negatively correlated (colour code: violet) and included genes involved in signal regulatory protein family interactions and lymphotoxin β receptor signalling (Supplementary Table S4 and Supplementary Table S5). These modules were found to correlate irrespective of maximal disease extent during follow-up (left-sided colitis/proctitis or extensive colitis). Of the 20 differentially expressed genes in inactive extensive UC, 6 were assigned to the purple module (Supplementary Table S5).
The three modules most correlated with active UC (colour code: saddlebrown, steelblue, and brown) were all positively correlated and mainly involved in immune functions (e.g., IFN and cytokine signalling, NFkB signalling, and antigen presentation) (Table 4, Supplementary Table S4). The most significantly enriched GO biological processes were the defence response to virus, protein deubiquitination, and regulation of biological quality, respectively. The frequently enriched GO pathways were the microbial response and immune reactions (Supplementary Table S4). Importantly, unlike in patients with inactive disease, this correlation seemed to be driven by maximal disease extent, as a significant correlation was only observed for extensive disease and not with left-sided disease extent during follow-up (Table 4, Supplementary Figure S3). The small sample size of these subgroups should however preclude from drawing firm conclusions. Similar results were obtained when we applied the disease extent at the time of endoscopy for patients with active disease (Supplementary Figure S3). Among 84 dysregulated genes in active extensive UC, 29 genes were assigned to the saddlebrown (n = 18) and brown (n = 11) modules (Table 4).

4. Discussion

We found molecular alterations in the non-inflamed ileum of UC patients, which were dependent on the presence of colonic inflammation and disease extent, and thus point to a cross-talk between colon and ileum (see Figure 3 for an overview of the main findings). While patients with inactive colitis showed a similar gene expression profile in the ileum as non-IBD controls, more alterations in ileal gene expression were seen in patients with active colitis at the time of endoscopy. Patients who had previously experienced extensive disease showed a distinct gene expression pattern at the ileal mucosa, regardless of disease activity in the colon.
When interpreting the results, we should keep in mind the potential limitations of our study. First, the current analysis is limited by the small sample size of our cohort, particularly for active left-sided disease, where we only have two samples. Further studies with a larger sample size are required to verify the present observations. Second, we acknowledge the possibility of confounders such as smoking or different medication influencing our results, although in our study we did not observe an influence of the use of biologics on PCA clusters. We also conducted a differential gene expression analysis including the use of biologics as a covariate and found that the obtained results were largely the same as the current results (data not shown). Third, it would have been interesting to have included the paired colonic mucosal tissue samples from the included patients, as these provide insight into the relationship between molecular alterations at the terminal ileum and colonic inflammation. These samples, however, were not available. Fourth, although the results of our study are supported by previous multiple lines of evidence that the identified genes do have biological relevance for IBD, we lack the functional analysis to specifically address the biological consequences of the identified genes. More mechanistic studies are needed to further demonstrate the precise mechanism by which the dysregulated genes influence the development of UC. There are also some strengths that we want to highlight. To the best of our knowledge, we are the first to specifically report on the molecular changes in the unaffected terminal ileum in UC patients and the relationship of these changes with the presence and extent of colonic inflammation, based on high-throughput mucosal transcriptional profiling. Besides, in addition to performing simple and well-established differential gene expression analyses, we used an unsupervised bioinformatics approach (WGCNA) to analyse gene co-expression networks.
We demonstrated that both the presence of colonic inflammation and the disease extent were associated with transcriptional alterations of the terminal ileum. Comparative analyses of gene expression levels in ileal biopsies of active UC patients versus non-IBD controls identified the significantly different expression of 18 genes. By contrast, patients with inactive colitis had a similar ileal gene expression signature to non-IBD controls. Importantly though, in both active and inactive colitis patients, the strongest dysregulations were seen in patients with extensive disease as maximal disease extent during follow-up. In the group of patients with active and extensive disease, 84 genes were significantly dysregulated in the terminal ileum, and with the strongest fold changes. Interestingly, the patients who previously had extensive colitis, though without colonic inflammation at the time of endoscopy, presented a distinct expression pattern from patients who previously had left-sided colitis. This suggests that the non-inflamed terminal ileum of UC patients does not restore to the non-IBD control expression pattern after the macroscopic resolution of extensive inflammation. Rather, these patients still show transcriptional alterations in their terminal ileum, although they are largely different from the alterations seen in active extensive disease patients, with an overlap of only 5 of the 20 genes. In particular, CEBPD showed a significant up-regulation in inactive extensive UC patients compared with controls. CEBPD is known to be induced in many inflammation-related diseases [30] and plays a possible anti-inflammatory role [31,32,33]. In a Cebpd-knockout mouse model, significant increases in oxidative stress and mitochondrial dysfunction were induced by radiation [34,35]. However, the role of these alterations in the context of chronic colonic inflammation in UC remains unclear.
The most significantly differentially expressed gene in the terminal ileum of UC patients with active colitis was DUOXA2. This gene is known to be required for the maturation and function of DUOX2, which was also differentially expressed in active UC. This increased DUOXA2 and DUOX2 expression confirmed earlier findings where both genes were found to be increased in the ileum of both UC and CD [15]. The DUOX2/DUOXA2 complex plays a role in maintaining mucosal immune homeostasis by the hydrogen peroxide-based epithelial defence [17,18,19] and in the interaction between host and microbiome [16]. In addition, the potential role of DUOX2 in maintaining the epithelial barrier in the small bowel has been described in a Duoxa-knockout mouse model [36], and in patients with CD [37]. Another interesting pair of genes from our analyses is REG1A/REG1B. We observed a universal upregulation of REG1A/REG1B across subtypes of UC, including both active extensive colitis and inactive extensive colitis. REG1 proteins are involved in cell differentiation and proliferation within the gastrointestinal tract [38,39], and have been implicated in intestinal inflammation. For example, they showed upregulation during intestinal inflammation in infectious colitis [40]; were found to be up-regulated in metaplastic Paneth cells in both un-inflamed and inflamed colonic mucosa from IBD patients, with increased expression during inflammation [41,42]; and REG1A/REG1B was also found to be widely expressed in a newly identified T-cell subtype isolated from CD patients using a single-cell analysis of intestinal T cells [43]. Finally, genes MUC1 and MUC4 encode membrane-bound mucins from epithelial cells which aid in the intestinal mucosal barrier function [44,45,46]. Previously, our group found a possible critical role of epithelial barrier genes, including MUC1 and MUC4, in IBD [47]. Particularly, MUC1 and MUC4 showed increased expression levels compared with controls not only in active CD and UC, but also in uninflamed CD ileum and IBD colon. Furthermore, previous genome-wide association studies [48] have identified the MUC1 locus as a susceptibility locus for IBD, with additional evidence from expression quantitative trait loci data (GTEx Analysis Release V8, dbGaP Accession phs000424.v8.p2) and experimental data [49]. Taken together, for these three pairs of genes (DUOXA2/DUOX2, REG1A/REG1B, and MUC1/MUC4), there is biologically relevant evidence in the context of IBD, and we here speculate that they also play a role in the non-inflamed ileum of UC patients, indicating the need for further studies.
An important remaining question is if the observed molecular alterations are specific to the location (ileum versus colon) or the disease (UC versus CD). While we unfortunately did not have colonic tissue available from the included patients, the identified dysregulated genes largely overlapped with those described for colonic tissue in UC [50,51,52,53], or ileum in CD [15,52,54] (Supplementary Table S6). Of the top 20 dysregulated genes in active extensive colitis, 14 showed similar findings in inflamed ileal tissue of CD [15,52,54], as well as in colonic tissue from patients with UC (13 genes) [50,51,52,53]. In addition, previous studies have reported that MUC1, DUOX2, and DMBT1 are altered in the unaffected ileum of CD patients [55,56]. Genes dysregulated in the terminal ileum of inactive extensive colitis were also previously found to be differentially expressed in colonic tissue from patients with UC (20 of 20 genes) [50,51,52,53], and to a lesser extent in inflamed ileal tissue of CD (11 of 20 genes) [15,52,54] (Supplementary Table S6). Taken together, this suggests that some of the molecular alterations we observed in the normal ileum of active UC patients may be a shared feature of IBD, independent of the phenotype or tissue location, and thus do not represent an explanation for the spatial restriction of the inflammation in UC. It should also be noted that these molecular alterations might be associated with general intestinal inflammation itself, and not necessarily specific to UC or IBD. Additional research including the comparison with, for example, infectious colitis is needed to explore the role of these alterations in the specific context of IBD. Nevertheless, given the aberrant molecular properties of the non-inflamed ileum in UC, it is clear there is a cross-talk between the colon and ileum in UC, although we cannot currently infer a causal and/or temporal relationship between these changes and colonic inflammation.
Even in the absence of large gene expression differences, gene co-expression network analysis can identify modules of co-expressed genes [57]. In inactive UC, we found that the correlated modules seemed to converge into intestinal barrier maintenance and mitochondrial function gene expression, and dampened immune functions. Together with multiple lines of evidence of alterations in mitochondrial function in patients with UC [51,58,59,60], this finding suggests a potential role for mitochondrial function in maintaining intestinal homeostasis. Remarkably, in general, mitochondrial function is described as decreasing in the colon of UC patients, while we found an increased function in the non-inflamed ileum of inactive UC patients. In active UC, the identified modules converged into increased immune functions. These associations can either precede (and possibly promote) colonic inflammation, be a consequence of colonic inflammation (i.e., a secondary phenomenon), or both. In conclusion, this transcriptomic analysis has identified significant molecular alterations in the non-inflamed ileum in UC. The observed changes are dependent on UC disease activity and extent (as illustrated in Figure 3), suggesting a cross-talk between the ileum and colon. With the current study, however, we cannot infer the causal and/or temporal relationship. Future longitudinal prospective studies with paired ileal and colonic samples including non-IBD controls with and without gut inflammation are therefore needed.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/9/8/1793/s1, Figure S1: Principal component analysis plot, Figure S2. Gene expression analysis, Figure S3: Weighted gene co-expression network analysis, Table S1: Selected genes for the qRT-PCR analysis, Table S2: Significantly differentially expressed genes, Table S3: The predicted upstream regulators of differentially expressed genes, Table S4: The highly ranked pathways and gene ontology biological process, Table S5: Significantly correlated co-expression modules in inactive ulcerative colitis, Table S6: Evidence from previous studies for identified dysregulated genes in Table 3.

Author Contributions

H.-S.L., M.V.: statistical and data analysis, writing; M.V.: experimental procedures; S.V. (Sare Verstockt), B.V.: support, statistical and data analysis, and experimental procedures; T.W.: experimental procedures; J.S., M.F., S.V. (Séverine Vermeire): patient and sample inclusion, study design; I.C.: study design, writing. All authors have read and agreed to the published version of the manuscript.

Funding

H.-S. Lee was supported by a National Research Foundation of Korea (NRF) MRC grant funded by the Korean government (MSIT) (2018R1A5A2020732). B. Verstockt is a doctoral fellow, and M. Ferrante and S. Vermeire are Senior Clinical Investigators of the Research Foundation Flanders (FWO), Belgium.

Acknowledgments

We would like to thank the patients who contributed to this study. The authors would also like to thank Vera Ballet, Eline Vandeput, Nooshin Ardeshir Davani, Helene Blevi, Tamara Coopmans, Sophie Organe, and Willem-Jan Wollants for an excellent job in maintaining the Biobank database and sample processing; and Vanessa Brys, Jens Van Bouwel, Wim Meert, Alvaro Cortes Calabuig, Céline Helsmoortel, and Wouter Bossuyt (Genomics Core Facility, UZ/KU Leuven, Belgium) for the technical and analytical assistance with the RNA sequencing (library preparation, data generation, raw data processing).

Conflicts of Interest

B. Verstockt received financial support for research from Pfizer; lecture fees from Abbvie, Ferring Pharmaceuticals, Janssen, R-biopharm, and Takeda; consultancy fees from Janssen and Sandoz. J. Sabino acted as a speaker for Abbvie and Nestlé Health Sciences. M. Ferrante received financial support for research from Amge, Biogen, Pfizer, Takeda, and Janssen; lecture fees from Ferring, Boehringer-Ingelheim, Chiesi, Merck Sharpe & Dohme, Tillotts, Janssen Biologics, Abbvie, Takeda, Mitsubishi Tanabe, and Zeria; consultancy fees from Abbvie, Boehringer-Ingelheim, Ferring, Janssen Biologics, Merck Sharpe & Dohme, and Takeda. S. Vermeire received financial support for research from MSD, Abbvie, Janssen, and Pfizer; lecture fees from Abbott, Abbvie, Merck Sharpe & Dohme, Ferring Pharmaceuticals, and UCB Pharma; consultancy fees from Pfizer, Ferring Pharmaceuticals, Shire Pharmaceuticals Group, Merck Sharpe & Dohme, and AstraZeneca Pharmaceuticals. H.-S.L., M.V., Sa.V, and I.C. declare no conflicts of interest.

References

  1. Ishii, H.; Isomoto, H.; Shikuwa, S.; Hayashi, T.; Inoue, N.; Yamaguchi, N.; Ohnita, K.; Nanashima, A.; Ito, M.; Nakao, K.; et al. Peyer’s Patches in the Terminal Ileum in Ulcerative Colitis: Magnifying Endoscopic Findings. J. Clin. Biochem. Nutr. 2010, 46, 111–118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Xavier, R.J.; Podolsky, D.K. Unravelling the pathogenesis of inflammatory bowel disease. Nature 2007, 448, 427–434. [Google Scholar] [CrossRef]
  3. Mourad, F.H.; Barada, K.A.; Saade, N.E. Impairment of Small Intestinal Function in Ulcerative Colitis: Role of Enteric Innervation. J. Crohns Colitis 2017, 11, 369–377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Salem, S.N.; Truelove, S.C. Small-intestinal and gastric abnormalities in ulcerative colitis. Br. Med. J. 1965, 1, 827–831. [Google Scholar] [CrossRef] [Green Version]
  5. Jankey, N.; Price, L.A. Small intestinal histochemical and histological changes in ulcerative colitis. Gut 1969, 10, 267–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ferguson, R.; Allan, R.N.; Cooke, W.T. A study of the cellular infiltrate of the proximal jejunal mucosa in ulcerative colitis and Crohn’s disease. Gut 1975, 16, 205–208. [Google Scholar] [PubMed] [Green Version]
  7. Sedghi, S.; Barreau, F.; Morilla, I.; Montcuquet, N.; Cazals-Hatem, D.; Pedruzzi, E.; Rannou, E.; Treton, X.; Hugot, J.P.; Ogier-Denis, E.; et al. Increased Proliferation of the Ileal Epithelium as a Remote Effect of Ulcerative Colitis. Inflamm. Bowel Dis. 2016, 22, 2369–2381. [Google Scholar] [CrossRef]
  8. Buning, C.; Geissler, N.; Prager, M.; Sturm, A.; Baumgart, D.C.; Buttner, J.; Buhner, S.; Haas, V.; Lochs, H. Increased small intestinal permeability in ulcerative colitis: Rather genetic than environmental and a risk factor for extensive disease? Inflamm. Bowel Dis. 2012, 18, 1932–1939. [Google Scholar] [CrossRef]
  9. Chang, J.; Leong, R.W.; Wasinger, V.C.; Ip, M.; Yang, M.; Phan, T.G. Impaired Intestinal Permeability Contributes to Ongoing Bowel Symptoms in Patients With Inflammatory Bowel Disease and Mucosal Healing. Gastroenterology 2017, 153, 723–731. [Google Scholar] [CrossRef]
  10. Amit-Romach, E.; Reifen, R.; Uni, Z. Mucosal function in rat jejunum and ileum is altered by induction of colitis. Int. J. Mol. Med. 2006, 18, 721–727. [Google Scholar] [CrossRef] [Green Version]
  11. Mazzon, E.; Cuzzocrea, S. Absence of functional peroxisome proliferator-activated receptor-alpha enhanced ileum permeability during experimental colitis. Shock 2007, 28, 192–201. [Google Scholar] [CrossRef] [PubMed]
  12. Arrieta, M.C.; Madsen, K.; Doyle, J.; Meddings, J. Reducing small intestinal permeability attenuates colitis in the IL10 gene-deficient mouse. Gut 2009, 58, 41–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Peters, L.A.; Perrigoue, J.; Mortha, A.; Iuga, A.; Song, W.M.; Neiman, E.M.; Llewellyn, S.R.; Di Narzo, A.; Kidd, B.A.; Telesco, S.E.; et al. A functional genomics predictive network model identifies regulators of inflammatory bowel disease. Nat. Genet. 2017, 49, 1437–1449. [Google Scholar] [CrossRef] [PubMed]
  14. Kugathasan, S.; Denson, L.A.; Walters, T.D.; Kim, M.O.; Marigorta, U.M.; Schirmer, M.; Mondal, K.; Liu, C.; Griffiths, A.; Noe, J.D.; et al. Prediction of complicated disease course for children newly diagnosed with Crohn’s disease: A multicentre inception cohort study. Lancet 2017, 389, 1710–1718. [Google Scholar] [CrossRef] [Green Version]
  15. Haberman, Y.; Tickle, T.L.; Dexheimer, P.J.; Kim, M.O.; Tang, D.; Karns, R.; Baldassano, R.N.; Noe, J.D.; Rosh, J.; Markowitz, J.; et al. Pediatric Crohn disease patients exhibit specific ileal transcriptome and microbiome signature. J. Clin. Investig. 2014, 124, 3617–3633. [Google Scholar] [CrossRef] [Green Version]
  16. Lloyd-Price, J.; Arze, C.; Ananthakrishnan, A.N.; Schirmer, M.; Avila-Pacheco, J.; Poon, T.W.; Andrews, E.; Ajami, N.J.; Bonham, K.S.; Brislawn, C.J.; et al. Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases. Nature 2019, 569, 655–662. [Google Scholar] [CrossRef]
  17. Geiszt, M.; Witta, J.; Baffi, J.; Lekstrom, K.; Leto, T.L. Dual oxidases represent novel hydrogen peroxide sources supporting mucosal surface host defense. FASEB J. 2003, 17, 1502–1504. [Google Scholar] [CrossRef]
  18. Ha, E.M.; Oh, C.T.; Bae, Y.S.; Lee, W.J. A direct role for dual oxidase in Drosophila gut immunity. Science 2005, 310, 847–850. [Google Scholar] [CrossRef]
  19. Rigoni, A.; Poulsom, R.; Jeffery, R.; Mehta, S.; Lewis, A.; Yau, C.; Giannoulatou, E.; Feakins, R.; Lindsay, J.O.; Colombo, M.P.; et al. Separation of Dual Oxidase 2 and Lactoperoxidase Expression in Intestinal Crypts and Species Differences May Limit Hydrogen Peroxide Scavenging During Mucosal Healing in Mice and Humans. Inflamm. Bowel Dis. 2017, 24, 136–148. [Google Scholar] [CrossRef] [Green Version]
  20. 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]
  21. Schroeder, K.W.; Tremaine, W.J.; Ilstrup, D.M. Coated oral 5-aminosalicylic acid therapy for mildly to moderately active ulcerative colitis. A randomized study. N. Engl. J. Med. 1987, 317, 1625–1629. [Google Scholar] [CrossRef] [PubMed]
  22. Silverberg, M.S.; Satsangi, J.; Ahmad, T.; Arnott, I.D.; Bernstein, C.N.; Brant, S.R.; Caprilli, R.; Colombel, J.F.; Gasche, C.; Geboes, K.; et al. Toward an integrated clinical, molecular and serological classification of inflammatory bowel disease: Report of a Working Party of the 2005 Montreal World Congress of Gastroenterology. Can. J. Gastroenterol. 2005, 19, 5A–36A. [Google Scholar] [CrossRef] [PubMed]
  23. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  25. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef] [Green Version]
  26. Love, M.; Anders, S.; Huber, W. Differential analysis of count data–the deseq2 package. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  27. 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]
  28. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [Green Version]
  29. Mi, H.; Muruganujan, A.; Ebert, D.; Huang, X.; Thomas, P.D. PANTHER version 14: More genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019, 47, D419–D426. [Google Scholar] [CrossRef]
  30. Balamurugan, K.; Sterneck, E. The many faces of C/EBPdelta and their relevance for inflammation and cancer. Int. J. Biol. Sci. 2013, 9, 917–933. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, Y.W.; Tseng, H.P.; Chen, L.C.; Chen, B.K.; Chang, W.C. Functional cooperation of simian virus 40 promoter factor 1 and CCAAT/enhancer-binding protein beta and delta in lipopolysaccharide-induced gene activation of IL-10 in mouse macrophages. J. Immunol. 2003, 171, 821–828. [Google Scholar] [CrossRef] [Green Version]
  32. Moore, F.; Santin, I.; Nogueira, T.C.; Gurzov, E.N.; Marselli, L.; Marchetti, P.; Eizirik, D.L. The transcription factor C/EBP delta has anti-apoptotic and anti-inflammatory roles in pancreatic beta cells. PLoS ONE 2012, 7, e31062. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Rustenhoven, J.; Scotter, E.L.; Jansson, D.; Kho, D.T.; Oldfield, R.L.; Bergin, P.S.; Mee, E.W.; Faull, R.L.; Curtis, M.A.; Graham, S.E.; et al. An anti-inflammatory role for C/EBPdelta in human brain pericytes. Sci. Rep. 2015, 5, 12132. [Google Scholar] [CrossRef] [Green Version]
  34. Pawar, S.A.; Shao, L.; Chang, J.; Wang, W.; Pathak, R.; Zhu, X.; Wang, J.; Hendrickson, H.; Boerma, M.; Sterneck, E.; et al. C/EBPδ Deficiency Sensitizes Mice to Ionizing Radiation-Induced Hematopoietic and Intestinal Injury. PLoS ONE 2014, 9, e94967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Banerjee, S.; Aykin-Burns, N.; Krager, K.J.; Shah, S.K.; Melnyk, S.B.; Hauer-Jensen, M.; Pawar, S.A. Loss of C/EBPdelta enhances IR-induced cell death by promoting oxidative stress and mitochondrial dysfunction. Free Radic. Biol. Med. 2016, 99, 296–307. [Google Scholar] [CrossRef] [Green Version]
  36. Grasberger, H.; Gao, J.; Nagao-Kitamoto, H.; Kitamoto, S.; Zhang, M.; Kamada, N.; Eaton, K.A.; El-Zaatari, M.; Shreiner, A.B.; Merchant, J.L.; et al. Increased Expression of DUOX2 Is an Epithelial Response to Mucosal Dysbiosis Required for Immune Homeostasis in Mouse Intestine. Gastroenterology 2015, 149, 1849–1859. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Levine, A.P.; Pontikos, N.; Schiff, E.R.; Jostins, L.; Speed, D.; Lovat, L.B.; Barrett, J.C.; Grasberger, H.; Plagnol, V.; Segal, A.W. Genetic Complexity of Crohn’s Disease in Two Large Ashkenazi Jewish Families. Gastroenterology 2016, 151, 698–709. [Google Scholar] [CrossRef] [Green Version]
  38. Miyaoka, Y.; Kadowaki, Y.; Ishihara, S.; Ose, T.; Fukuhara, H.; Kazumori, H.; Takasawa, S.; Okamoto, H.; Chiba, T.; Kinoshita, Y. Transgenic overexpression of Reg protein caused gastric cell proliferation and differentiation along parietal cell and chief cell lineages. Oncogene 2004, 23, 3572–3579. [Google Scholar] [CrossRef] [Green Version]
  39. Ose, T.; Kadowaki, Y.; Fukuhara, H.; Kazumori, H.; Ishihara, S.; Udagawa, J.; Otani, H.; Takasawa, S.; Okamoto, H.; Kinoshita, Y. Reg I-knockout mice reveal its role in regulation of cell growth that is required in generation and maintenance of the villous structure of small intestine. Oncogene 2007, 26, 349–359. [Google Scholar] [CrossRef] [Green Version]
  40. Peterson, K.M.; Guo, X.; Elkahloun, A.G.; Mondal, D.; Bardhan, P.K.; Sugawara, A.; Duggal, P.; Haque, R.; Petri, W.A., Jr. The expression of REG 1A and REG 1B is increased during acute amebic colitis. Parasitol. Int. 2011, 60, 296–300. [Google Scholar] [CrossRef] [Green Version]
  41. Van Beelen Granlund, A.; Ostvik, A.E.; Brenna, O.; Torp, S.H.; Gustafsson, B.I.; Sandvik, A.K. REG gene expression in inflamed and healthy colon mucosa explored by in situ hybridisation. Cell Tissue Res. 2013, 352, 639–646. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Granlund, A.; Beisvag, V.; Torp, S.H.; Flatberg, A.; Kleveland, P.M.; Ostvik, A.E.; Waldum, H.L.; Sandvik, A.K. Activation of REG family proteins in colitis. Scand. J. Gastroenterol. 2011, 46, 1316–1323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Uniken Venema, W.T.; Voskuil, M.D.; Vila, A.V.; van der Vries, G.; Jansen, B.H.; Jabri, B.; Faber, K.N.; Dijkstra, G.; Xavier, R.J.; Wijmenga, C.; et al. Single-Cell RNA Sequencing of Blood and Ileal T Cells From Patients With Crohn’s Disease Reveals Tissue-Specific Characteristics and Drug Targets. Gastroenterology 2019, 156, 812–815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hanson, R.L.; Hollingsworth, M.A. Functional Consequences of Differential O-glycosylation of MUC1, MUC4, and MUC16 (Downstream Effects on Signaling). Biomolecules 2016, 6, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Peterson, L.W.; Artis, D. Intestinal epithelial cells: Regulators of barrier function and immune homeostasis. Nat. Rev. Immunol. 2014, 14, 141–153. [Google Scholar] [CrossRef]
  46. Andrianifahanana, M.; Moniaux, N.; Batra, S.K. Regulation of mucin expression: Mechanistic aspects and implications for cancer and inflammatory diseases. Biochim. Biophys. Acta 2006, 1765, 189–222. [Google Scholar] [CrossRef]
  47. Vancamelbeke, M.; Vanuytsel, T.; Farré, R.; Verstockt, S.; Ferrante, M.; Van Assche, G.; Rutgeerts, P.; Schuit, F.; Vermeire, S.; Arijs, I.; et al. Genetic and Transcriptomic Bases of Intestinal Epithelial Barrier Dysfunction in Inflammatory Bowel Disease. Inflamm. Bowel Dis. 2017, 23, 1718–1729. [Google Scholar] [CrossRef]
  48. De Lange, K.M.; Moutsianas, L.; Lee, J.C.; Lamb, C.A.; Luo, Y.; Kennedy, N.A.; Jostins, L.; Rice, D.L.; Gutierrez-Achury, J.; Ji, S.G.; et al. Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat. Genet. 2017, 49, 256–261. [Google Scholar] [CrossRef] [Green Version]
  49. Breugelmans, T.; Van Spaendonk, H.; De Man, J.G.; De Schepper, H.U.; Jauregui-Amezaga, A.; Macken, E.; Lindén, S.K.; Pintelon, I.; Timmermans, J.P.; De Winter, B.Y.; et al. In depth study of transmembrane mucins in association with intestinal barrier dysfunction during the course of T cell transfer and DSS-induced colitis. J. Crohns Colitis 2020. [Google Scholar] [CrossRef]
  50. Taman, H.; Fenton, C.G.; Hensel, I.V.; Anderssen, E.; Florholmen, J.; Paulssen, R.H. Transcriptomic Landscape of Treatment-Naive Ulcerative Colitis. J. Crohns Colitis 2018, 12, 327–336. [Google Scholar] [CrossRef]
  51. Haberman, Y.; Karns, R.; Dexheimer, P.J.; Schirmer, M.; Somekh, J.; Jurickova, I.; Braun, T.; Novak, E.; Bauman, L.; Collins, M.H.; et al. Ulcerative colitis mucosal transcriptomes reveal mitochondriopathy and personalized mechanisms underlying disease severity and treatment response. Nat. Commun. 2019, 10, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Arijs, I.; De Hertogh, G.; Lemaire, K.; Quintens, R.; Van Lommel, L.; Van Steen, K.; Leemans, P.; Cleynen, I.; Van Assche, G.; Vermeire, S.; et al. Mucosal gene expression of antimicrobial peptides in inflammatory bowel disease before and after first infliximab treatment. PLoS ONE 2009, 4, e7984. [Google Scholar] [CrossRef] [PubMed]
  53. Vanhove, W.; Peeters, P.M.; Staelens, D.; Schraenen, A.; Van der Goten, J.; Cleynen, I.; De Schepper, S.; Van Lommel, L.; Reynaert, N.L.; Schuit, F.; et al. Strong Upregulation of AIM2 and IFI16 Inflammasomes in the Mucosa of Patients with Active Inflammatory Bowel Disease. Inflamm. Bowel Dis. 2015, 21, 2673–2682. [Google Scholar] [CrossRef]
  54. Verstockt, S.; De Hertogh, G.; Van der Goten, J.; Verstockt, B.; Vancamelbeke, M.; Machiels, K.; Van Lommel, L.; Schuit, F.; Van Assche, G.; Rutgeerts, P.; et al. Gene and Mirna Regulatory Networks During Different Stages of Crohn’s Disease. J. Crohns Colitis 2019, 13, 916–930. [Google Scholar] [CrossRef] [PubMed]
  55. Hamm, C.M.; Reimers, M.A.; McCullough, C.K.; Gorbe, E.B.; Lu, J.; Gu, C.C.; Li, E.; Dieckgraefe, B.K.; Gong, Q.; Stappenbeck, T.S.; et al. NOD2 status and human ileal gene expression. Inflamm. Bowel Dis. 2010, 16, 1649–1657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Zhang, T.; Song, B.; Zhu, W.; Xu, X.; Gong, Q.Q.; Morando, C.; Dassopoulos, T.; Newberry, R.D.; Hunt, S.R.; Li, E. An ileal Crohn’s disease gene signature based on whole human genome expression profiles of disease unaffected ileal mucosal biopsies. PLoS ONE 2012, 7, e37139. [Google Scholar] [CrossRef] [Green Version]
  57. Zhang, B.; Horvath, S. A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol. 2005, 4, Article17. [Google Scholar] [CrossRef]
  58. Santhanam, S.; Rajamanickam, S.; Motamarry, A.; Ramakrishna, B.S.; Amirtharaj, J.G.; Ramachandran, A.; Pulimood, A.; Venkatraman, A. Mitochondrial electron transport chain complex dysfunction in the colonic mucosa in ulcerative colitis. Inflamm. Bowel Dis. 2012, 18, 2158–2168. [Google Scholar] [CrossRef] [PubMed]
  59. Mottawea, W.; Chiang, C.K.; Muhlbauer, M.; Starr, A.E.; Butcher, J.; Abujamel, T.; Deeke, S.A.; Brandel, A.; Zhou, H.; Shokralla, S.; et al. Altered intestinal microbiota-host mitochondria crosstalk in new onset Crohn’s disease. Nat. Commun. 2016, 7, 13419. [Google Scholar] [CrossRef]
  60. Rath, E.; Moschetta, A.; Haller, D. Mitochondrial function—Gatekeeper of intestinal epithelial cell homeostasis. Nat. Rev. Gastroenterol. Hepatol. 2018, 15, 497–516. [Google Scholar] [CrossRef]
Cells 09 01793 g001 550
Figure 1. Overview of significantly dysregulated genes in the different comparisons. Venn diagrams of the differentially expressed genes in the comparative analysis of (A) active or inactive ulcerative colitis (UC) patients and non- inflammatory bowel disease (IBD) controls; and of (B) Extensive UC or left-sided colitis/proctitis and non-IBD controls, in active and inactive UC separately. The number of significantly differentially expressed genes in each comparative analysis is shown in circles. The number of overlapping genes between comparisons is also indicated.
Figure 1. Overview of significantly dysregulated genes in the different comparisons. Venn diagrams of the differentially expressed genes in the comparative analysis of (A) active or inactive ulcerative colitis (UC) patients and non- inflammatory bowel disease (IBD) controls; and of (B) Extensive UC or left-sided colitis/proctitis and non-IBD controls, in active and inactive UC separately. The number of significantly differentially expressed genes in each comparative analysis is shown in circles. The number of overlapping genes between comparisons is also indicated.
Cells 09 01793 g001
Cells 09 01793 g002 550
Figure 2. Validation of expression levels of seven selected mRNAs in ileal mucosal biopsies. Boxplots of the relative expression level of (A) DUOXA2, (B) DUOX2, (C) REG1B, (D) REG1A, (E) MUC4, (F) GRAMD2, and (G) CASP10 in controls (n = 15), active extensive UC (n = 4), and inactive extensive UC (n = 14) patients, as assessed by qRT-PCR (box, 25–75%; whisker, min-max). The expression levels are normalized to the geometric mean of ACTB and GAPDH. Data were compared using the one-tailed Mann–Whitney U test (solid line with down angle, * p < 0.05; ** p < 0.01; *** p < 0.001).
Figure 2. Validation of expression levels of seven selected mRNAs in ileal mucosal biopsies. Boxplots of the relative expression level of (A) DUOXA2, (B) DUOX2, (C) REG1B, (D) REG1A, (E) MUC4, (F) GRAMD2, and (G) CASP10 in controls (n = 15), active extensive UC (n = 4), and inactive extensive UC (n = 14) patients, as assessed by qRT-PCR (box, 25–75%; whisker, min-max). The expression levels are normalized to the geometric mean of ACTB and GAPDH. Data were compared using the one-tailed Mann–Whitney U test (solid line with down angle, * p < 0.05; ** p < 0.01; *** p < 0.001).
Cells 09 01793 g002
Cells 09 01793 g003 550
Figure 3. Schematic summary of the current study. Ileal expression from ulcerative colitis (UC) patients is dependent on UC disease activity and extent. The strongest dysregulations were seen in patients with extensive disease as maximal disease extent during follow-up. Co-expression module analysis found that ileal changes in active extensive UC are mainly related to the immune function (colour code for identified co-expression modules: saddlebrown, steelblue, and brown). Ileal changes in inactive UC, on the other hand, seem to be functioning to maintain the intestinal barrier with increased mitochondrial functions (colour code for identified co-expression modules: purple and grey60) and dampened immune functions (colour code for identified co-expression module: violet).
Figure 3. Schematic summary of the current study. Ileal expression from ulcerative colitis (UC) patients is dependent on UC disease activity and extent. The strongest dysregulations were seen in patients with extensive disease as maximal disease extent during follow-up. Co-expression module analysis found that ileal changes in active extensive UC are mainly related to the immune function (colour code for identified co-expression modules: saddlebrown, steelblue, and brown). Ileal changes in inactive UC, on the other hand, seem to be functioning to maintain the intestinal barrier with increased mitochondrial functions (colour code for identified co-expression modules: purple and grey60) and dampened immune functions (colour code for identified co-expression module: violet).
Cells 09 01793 g003
Table
Table 1. Baseline characteristics.
Table 1. Baseline characteristics.
Active UC (n = 7)Inactive UC (n = 29)Control (n = 15)
Male (%)5 (71.4)17 (58.6)6 (40.0)
Median age at endoscopy (years, IQR)52 (36–68)52 (33–60)55 (46–60)
Median disease duration (years, IQR)12 (6–13)10 (8–14)
Maximal disease extent (%)
 Proctitis04 (13.8)
 Left-sided colitis2 (28.6)11 (37.9)
 Extensive colitis5 (71.4)14 (48.3)
Disease extent at endoscopy (%)
 Left-sided colitis4 (57.1)N/A
 Extensive colitis3 (42.9)N/A
Concomitant medication (%)
 Corticosteroids1 (14.3)0
 Immunomodulators05 (17.2)
 Biologics1 (14.3)13 (44.8)
UC, ulcerative colitis; IQR, interquartile range; N/A, not applicable.
Table
Table 2. Significantly differentially expressed genes in the terminal ileum of UC patients.
Table 2. Significantly differentially expressed genes in the terminal ileum of UC patients.
UC (n = 36)
vs Control (n = 15)		Active UC (n = 7)
vs Control (n = 15)		Inactive UC (n = 29)
vs Control (n = 15)
Genelog2FCpcorrectedlog2FCpcorrectedlog2FCpcorrected
DUOXA22.920.064.865.58 × 10−32.380.23
REG1B3.134.85 × 10−34.150.022.750.03
PDZK1IP11.120.081.870.020.940.23
DUOX22.390.094.020.021.390.63
CASR−0.990.15−1.870.02−0.850.30
CBR30.700.131.150.030.580.31
DMBT11.480.112.400.031.140.34
MUC11.100.121.840.030.890.32
MTRNR2L10.381.00−5.940.030.840.95
GPR1101.960.183.370.041.590.41
REG1A1.450.152.140.041.220.23
NOS21.370.232.540.040.990.61
CFB0.810.271.520.040.590.61
GPR37L10.950.301.770.040.700.63
SPINK10.610.401.300.040.390.79
TSPO21.130.181.820.040.880.44
CLDN50.920.071.310.050.830.21
SSTR2−0.830.12−1.360.05−0.740.23
CEBPD1.070.020.730.621.140.02
Log2 fold-change and pcorrected for the genes that were significantly dysregulated. Significantly dysregulated genes are indicated in bold. Results in ascending order of p-values of the comparison between the patients with active colitis and controls. UC, ulcerative colitis; vs, versus; FC, fold-change.
Table
Table 3. Significantly dysregulated genes in the terminal ileum of ulcerative colitis patients with extensive colitis (as maximal disease extent).
Table 3. Significantly dysregulated genes in the terminal ileum of ulcerative colitis patients with extensive colitis (as maximal disease extent).
Active Extensive UC (n = 5)
vs Control (n = 15)		Inactive Extensive UC (n = 14)
vs Control (n = 15)
Genelog2FCpcorrectedlog2FCpcorrected
Top 20 dysregulated genes in active extensive UC
SPINK11.732.60 × 10−40.310.77
DUOXA25.462.81 × 10−42.790.07
PDZK1IP12.312.81 × 10−40.810.37
GPR1104.228.57 × 10−42.300.08
NOS23.008.57 × 10−41.580.10
GPR37L12.078.57 × 10−41.110.11
SLC37A11.028.57 × 10−40.390.35
TFEC−1.208.57 × 10−4−0.350.58
CCL281.921.47 × 10−30.510.67
CFB1.871.85 × 10−30.760.32
DMBT12.753.10 × 10−31.640.05
MUC12.113.44 × 10−31.200.08
PLA2G161.633.44 × 10−30.640.41
FUT31.143.44 × 10−30.490.30
NFKBIZ1.463.83 × 10−30.710.19
TIMD4−6.704.01 × 10−3−0.580.92
DUOX24.584.49 × 10−31.440.64
REG1A *2.444.49 × 10−31.650.03
TSPO22.134.49 × 10−31.180.13
ATP10B1.494.49 × 10−30.920.06
Dysregulated 20 genes in inactive extensive UC (REG1A is listed in the top 20 above)
MESDC10.450.371.082.89 × 10−5
CEBPD0.640.401.466.22 × 10−4
CARD141.550.182.119.90 × 10−3
REG1B *4.255.31 × 10−33.269.90 × 10−3
C2CD4B1.050.402.069.90 × 10−3
NCOA70.850.251.309.90 × 10−3
CD550.380.641.230.02
C2CD4A1.720.182.190.02
FAM83A1.630.232.230.02
HES40.830.471.710.02
GRAMD2 *1.970.021.530.02
HIC10.200.831.110.02
CXCL22.170.072.080.03
SOCS10.520.671.660.03
MUC4 *2.130.021.650.03
TIFA0.810.441.500.04
C4BPB1.030.131.090.05
CASP10 *1.350.031.030.05
EVA1B0.240.821.170.05
Log2 fold-change and pcorrected are given. Significantly dysregulated genes are indicated in bold. * Overlapping genes between active extensive UC vs control and inactive extensive UC vs control. FC, fold-change; UC, ulcerative colitis; vs, versus.
Table
Table 4. Significantly correlated co-expression modules in active ulcerative colitis.
Table 4. Significantly correlated co-expression modules in active ulcerative colitis.
Saddlebrown
(115 genes)		Steelblue
(113 genes)		Brown
(781 genes)
Active UC vs COr = 0.47r = 0.49r = 0.44
p = 0.03p = 0.02p = 0.04
pcorrected = 0.22pcorrected = 0.19pcorrected = 0.28
Extensive vs COr = 0.80r = 0.77r = 0.67
p = 2.78 × 10−05p = 6.32 × 10−05P = 1.17 × 10−3
pcorrected = 7.35 × 10−3pcorrected = 8.37 × 10−3pcorrected = 0.05
Lt-sided vs COr = −0.42r = −0.15r = −0.34
p = 0.09p = 0.58p = 0.18
pcorrected = 0.37pcorrected = 0.73pcorrected = 0.52
Top enriched pathway in ReactomeCytokine Signalling in Immune systemMetabolism of proteinsER-Phagosome pathway
Top enriched pathway in IPACell Death and SurvivalProtein Ubiquitination PathwayPhagosome Maturation
Top enriched GO biological processDefence response to virusProtein deubiquitinationRegulation of biological quality
Upstream Regulator *IFNγ, IFNα, STAT1EGF, EGFR, PDCD6NFE2L2, TCR, HSPA5
Overlapping dysregulated genes18/84 genes -11/84 genes
Correlation strengths r of each module were calculated for each ulcerative colitis subgroup with uncorrected/corrected p-values. The filled colours represent the strength of the association. A positive correlation (marked as red) means an upregulation in disease as compared to controls, while a negative correlation (marked as blue) refers to a downregulation in disease. The most highly scoring canonical pathways (according to p-value) according to Reactome, IPA, and GO enrichment analysis are given. * Top 3 ranked upstream regulators by IPA. UC, ulcerative colitis; CO, controls; Lt-sided, left-sided colitis as maximal disease extent; vs, versus.

Share and Cite

MDPI and ACS Style

Lee, H.-S.; Vancamelbeke, M.; Verstockt, S.; Wilms, T.; Verstockt, B.; Sabino, J.; Ferrante, M.; Vermeire, S.; Cleynen, I. Molecular Changes in the Non-Inflamed Terminal Ileum of Patients with Ulcerative Colitis. Cells 2020, 9, 1793. https://doi.org/10.3390/cells9081793

AMA Style

Lee H-S, Vancamelbeke M, Verstockt S, Wilms T, Verstockt B, Sabino J, Ferrante M, Vermeire S, Cleynen I. Molecular Changes in the Non-Inflamed Terminal Ileum of Patients with Ulcerative Colitis. Cells. 2020; 9(8):1793. https://doi.org/10.3390/cells9081793

Chicago/Turabian Style

Lee, Ho-Su, Maaike Vancamelbeke, Sare Verstockt, Tom Wilms, Bram Verstockt, João Sabino, Marc Ferrante, Séverine Vermeire, and Isabelle Cleynen. 2020. "Molecular Changes in the Non-Inflamed Terminal Ileum of Patients with Ulcerative Colitis" Cells 9, no. 8: 1793. https://doi.org/10.3390/cells9081793

APA Style

Lee, H. -S., Vancamelbeke, M., Verstockt, S., Wilms, T., Verstockt, B., Sabino, J., Ferrante, M., Vermeire, S., & Cleynen, I. (2020). Molecular Changes in the Non-Inflamed Terminal Ileum of Patients with Ulcerative Colitis. Cells, 9(8), 1793. https://doi.org/10.3390/cells9081793

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