Abstract

Inflammatory bowel disease (IBD) comprises ulcerative colitis (UC) and Crohn disease (CD). Similar symptoms, but different treatment procedures for both diseases require precise diagnosis. MicroRNAs (miRNAs) are major posttranscriptional players that regulate the expression of genes during the inflammation and thus could be appropriate biomarkers for differentiation between UC and CD. For this purpose, we analyzed the expression of miR-21-3p, miR-31-3p, miR-125b-1-3p, miR-146a-3p, miR-155-5p, and E-cadherin (CDH1) genes associated with IBD, in 67 tissue samples: 28 inflamed mucosa samples ( UC, CD), 28 adjacent normal colonic mucosa ( UC, CD), and 11 normal mucosa from healthy patients using reverse transcription real-time RT-PCR. We found all analyzed miRNAs were significantly overexpressed in UC tissue as compared to adjacent normal tissue of patients with UC, as well as to normal mucosa from healthy controls. Four miRNAs (except miR-125b-1-3p) were significantly upregulated in CD lesions as compared to adjacent normal tissue of patients with CD, and four miRNAs, except miR-146a-3p, were significantly higher in CD samples compared to normal mucosa from healthy individuals. In the CD group, we found an inverse correlation between miR-155-5p or miR-146a-3p expressions and CDH1expression in inflamed mucosa. This type of correlation was also detected for miR-213p in adjacent normal tissue and CDH1 in inflamed mucosa, as well as between miR-155-5p and CDH1 in adjacent normal tissue. Elevated miRNA expression is characteristic for IBD-mediated inflammation process and inversely correlated with CDH1 gene expression, which suggest involvement of epithelial to mesenchymal transition (EMT) in IBD development.

1. Introduction

IBD is a group of chronic disorders that occur in the form of ulcerative colitis (UC) and Crohn disease (CD) [1]. Inflammatory changes in CD involve all the layers of the intestinal wall and can develop in any part of the digestive tract, predominantly in the final section of the small intestine and in the large intestine [2]. UC is characterized by continuous inflammatory changes within the mucous membrane and affects only the large intestine [3]. More than 2 million Europeans and 1.5 million North Americans are suffering with IBD [4]. The incidence of IBD is increasing in newly industrialized nations, which is associated with nutrition consumption of processed food, refined sugars, dairy, and lower intake of plant-based fibers [5]. However, the etiology of IBD still remains unknown. Genome-wide association studies (GWAS) support a role of heredity in the pathogenesis of IBD that accounts 23% and 16% in CD and UC, respectively [6]. In individuals with genetic risk, abnormal interactions between the host immune system and gut flora, as well as dysregulation of cellular responses such as autophagy and endoplasmic reticulum stress, induce an immune response in the gut resulting in intestinal inflammation [7]. Other possible reasons for IBD development include influence of other genetic, and especially epigenetic, changes including microRNA-mediated regulation of gene expression [8]. MicroRNAs (miRNAs) are an abundant class of endogenous, small, single-stranded, noncoding RNAs composed of 18-25 nucleotides formed from double stranded precursors. These small molecular players repress protein translation through binding to complementary 3 untranslated region (UTR), 5 UTR, or coding sequence of targeted mRNAs [9]. miRNAs play a major role in a wide range of developmental process including cell proliferation, cell cycle, cell differentiation, metabolism, apoptosis, developmental timing, neuronal cell fate, neuronal gene expression, brain morphogenesis, muscle differentiation, and stem cell division [10]. Many miRNAs display organ-specific expression patterns suggesting cell-type-specific functions [11]. Abnormal expression of miRNAs can lead to pathological processes involving the immune system such as carcinogenesis and autoimmunity [12]. In addition, miRNAs are stable, easily isolated from different types of tissues, as well as from blood, which allow to consider them as an ideal biological marker [13]. The aim of our study was to determine and compare the expression of five microRNAs (miR-155-5p, miR-21-3p, miR-31-3p, miR-146a-3p, and miR-125b-1-3p) in tissue specimens of patients with active CD and UC. These miRNAs were selected for analysis, because of their involvement in the regulation of gene expression responsible for the development and progression of IBD [1416]. We evaluated the usefulness of these miRNAs as potential biomarkers for inflammatory process during IBD and possibly to differentiate between UC and CD. For that purpose, we applied miRWalk algorithm analysis to assess the possible miRNA-target interactions (http://mirwalk.umm.uni-heidelberg.de) with genes involved in IBD pathology. We determined the expression of E-cadherin (CDH1), a major transmembrane glycoprotein important for the maintenance for intestinal epithelial homeostasis [17]. E-cadherin is essential for proper morphogenesis of intestine, because it is a core component of adherent junctions, the major cell-cell adhesion structure, in which disruption may lead to inflammation and cancer [18]. Polymorphisms of CDH1 result in abnormal protein localization and may be involved in pathogenesis of IBD which disturb intestinal homeostasis [19]. We aimed to correlate expression of CDH1 gene with miRNA expression patterns in order to evaluate their potential role in IBD development. Since there is still lack of standard methods that help to distinguish between UC and CD, our present study could contribute to the improvement of the accuracy of the diagnosis together with different endoscopic techniques.

2. Materials and Methods

2.1. Patients

The selection of patients from the CD and UC groups follows the endoscopic evaluation of the severity of inflammatory lesions according to the Mayo endoscopic score for UC and the Simple Endoscopic Score for CD (SES-CD). The research was conducted in accordance with the Declaration of Helsinki (1964), and the local ethical committee approved the study design (No. KE-0254/2013). The biological material was collected during colonoscopies performed in the Endoscopic Unit of the Department of Gastroenterology, Medical University of Lublin, Poland.

The inclusion criteria of patients enrolled in our study are as follows: diagnosed and histopathologically confirmed UC and CD, age 18-65, and active inflammatory changes assessed >2 points on the Mayo scale for UC and >3 points on the SES-CD for CD. The exclusion criteria are as follows: gastrointestinal infections caused by Clostridium difficile, Salmonella, Shigella, Campylobacter, and Yersinia; ischemic changes, inflammation following X-ray therapy; and oral administration of antibiotics and probiotics.

The patients were divided into following groups: the group of 12 patients diagnosed with CD, and the group of 16 patients with UC, and two biopsy specimens: from the inflammation-altered tissue and the normal tissue collected from each patient. The control group consisted of 11 patients without any inflammatory lesions found on colonoscopies performed within the screening program. In this group, the biopsy specimens of the normal tissue were collected. Patients’ characteristic is depicted in Table 1.

2.2. Sample Handling, Isolation of RNA, and Synthesis of cDNA

Immediately after the collection of tissue biopsies, the specimens were placed in RNAlater™ (Sigma Aldrich, St. Louis, MO) and stored at -80°C for further RNA isolation.

Total RNA including miRNAs were isolated using miRCURY RNA Isolation Kit Tissues (Exiqon, Vedbaek, Denmark) according to the manufacturer’s instructions. The quantity and purity of RNA were analyzed spectrophotometrically, and 100 ng of RNA was used for cDNA synthesis, which was synthesized using the Universal cDNA Synthesis Kit II (Exiqon). Spike-in UniSp6 was added to each sample, which allowed to monitor the quality of reaction with reverse transcriptase. Three negative controls of this reaction were prepared using the following reaction mixtures—without reverse transcriptase (no enzyme sample), without RNA template (no template), and using MS2 bacteriophage RNA (Roche, Cat. No. 10165948001) as the template (Mock).

One microgram of total RNA isolated from tissue samples, as mentioned above, was reverse transcribed using the Transcriptor Universal cDNA Master (Roche Diagnostics, Mannheim, Germany) and subsequently used for the evaluation of the expression of CDH1. Negative control samples were prepared with the omission of an enzyme and RNA template.

2.3. MicroRNA Expression Analysis

Each cDNA sample was diluted with nuclease-free water and mixed with SYBR Green Master Mix with LNA™ primers for the defined miRNAs (Exiqon) according to manufacturer’s instructions.

In all tissue samples, the expression of UniSp6 was assessed as a marker of efficiency of cDNA synthesis. Each sample positively verified was qualified for further analysis. The relative quantification (RQ) of five selected miRNAs: miR-21-3p, miR-31-3p, miR-125b-1-3p, miR-146a-3p, and miR-155-5p, was determined in triplicate using Light Cycler 480® II Instrument and Light Cycler 480® Software 1.5 (Roche Diagnostics). From the potential two reference miRNAs: miR-103a-3p and miR-423-3p, the latter was chosen as a reference gene due to a constant level of expression without significant differences between analyzed groups ().

2.4. CDH1 Expression Analysis

The amplification of each cDNA sample from tissue was performed using LightCycler®480 SYBR Green I Master and sets of primers for CDH1 (forward: CAGGCTCAAGCTATCCTTGC, reverse: AGTCATGCGTAGTGGTGCAT) according to the manufacturer’s protocol. GAPDH was used as the reference gene for data normalization (primer forward: CTCTGCTCCTCTGTTCGAC, reverse: GCCCAATACGACCAAATCC).

Relative quantification (RQ) was evaluated in triplicate using Light Cycler 480® II Instrument and Light Cycler 480® Software 1.5 (Roche Diagnostics).

2.5. Statistical Analysis

The normality of distribution of miRNA expression in the studied biological samples was assessed using histograms and the Kołmogorov-Smirnov and Shapiro-Wilk tests. Since the distribution was not normal, the intergroup differences in miRNA expression were analyzed using nonparametric tests (Kruskal-Wallis ANOVA by ranks, Mann–Whitney and Wilcoxon rank-sum tests). Correlations between variables were analyzed using Spearman’s test. The differences in expression were considered significant at . Statistical calculations were performed using Statistica, version 9.0 software (StatSoft, Inc., Tulsa, OK 74104, USA).

3. Results

3.1. Comparative Analysis of miRNA Expression Profiles in CD, UC, and Healthy Individuals

All miRNAs were expressed in examined tissue samples. In biopsy specimens obtained from patients with CD, miR-21-3p, miR-31-3p, miR-146a-3p, and miR-155-5p were overexpressed (except miR-125b-1-3p) in TD tissues in comparison to TAC samples (Table 2). Our findings indicate that expression of miR-21-3p, miR-31-3p, miR-125b-1-3p, miR-146a-3p, and miR155-5p was significantly higher in biopsy specimens from inflammation-altered tissue in comparison to the adjacent tissue of patients with UC (Table 3). When analyzing expression of selected miRNA in inflammation-altered tissue of CD or UC patients with biopsy specimens of healthy patients, statistically significant differences were found (Table 4), whereas miRNA analysis did not reveal statistically significant differences in analyzed tissues between the CD and UC patient groups (both, between inflamed mucosa, and adjacent tissues, data not shown). Given the fact that particular miRNAs maintain the intestinal epithelium intact through inhibiting EMT (epithelial to mesenchymal transition) and promoting proliferation of intestinal cells [20, 21] we aimed to test if the expression of the analyzed miRNAs could correlate with the expression of E-cadherin coding gene (CDH1), which is a hallmark of EMT [22]. We noted negative correlations between miR-155-5p and CDH1 in both tissues—adjacent normal tissue () and inflamed mucosa () from patients with CD (Figures 1(a) and 1(b)). Similarly, negative correlation was also found between miR-146a-3p and CDH1 in inflamed mucosa only () in the same patient’s group (Figure 2). A negative correlation was noted between miR-21-3p () in adjacent normal tissue and expression of CDH1 in the inflamed mucosa of patients with CD () (Figure 3).

Taken together, in inflamed mucosa from CD patients, three miRNAs (miR-21-3p, miR-146a3p, and miR-155-5p) were highly expressed, whereas CDH1 expression was downregulated. Such correlations were not detected in UC patients. In these patients’ group (UC), we only found that the duration of disease was significantly and negatively correlated with the expression level of miR146a-3p (Spearman’s correlation test, , , Figure 4). We did not reveal any other correlations between the miRNA expression and patient’s clinical status.

4. Discussion

Despite years of research, the exact cause of nonspecific inflammatory bowel diseases (IBD) has not been fully explained, and the diagnosis and differentiation between CD and UC remains often difficult. Understanding the molecular mechanisms that modulate the expression of inflammation-related genes is crucial for the detection of new biomarkers of IBD. miRNAs have a huge potential as biological markers, because there are key molecular players in regulating a multitude of cellular processes, including IBD development [23, 24]. Specific miRNA signatures have been observed in IBD, associated with canonical signaling pathways regulating autophagy, inflammation, fibrosis, or EMT [25]. One of the analyzed miRNAs in the present paper is miR-155 which was shown to inhibit the suppressor of cytokine signaling (SOCS1)—a negative regulator of the Janus kinase/signal transducers and activators of transcription (JAK/STAT signaling); thus, elevating the expression of miR-155 could enhance inflammation in the intestine of IBD patients by immune system perturbation [25, 26]. In our study, we noticed an increased expression of miR-155-5p in both inflamed tissue and, most important, also adjacent tissue of patients with CD, suggesting that elevated miR-155-5p expression contributes to CD development. We further show that tissues near to CD-related lesions displayed increased miR-21-3p expression, and given the fact that involvement of elevated expression of miR-21 resulted in loss of tight junction proteins and increased barrier permeability [27], the role of these miRNA in pathogenesis of IBD tissues seems to be evident [28]. The molecular background of IBD development was linked to the EMT process with specific miRNA signatures affecting signaling, transcription factors, and mesenchymal or epithelial marker expression [25]. We determined the expression of one of these markers—E-cadherin (CDH1 gene)—in relation to the miRNA expression pattern in tissues from CD and UC patients.

In our study, we confirmed the lower expression of CDH1 only in patients with CD, but not in patients with UC, which suggest the involvement of this protein in CD pathogenesis. We noticed that downregulated CDH1 was negatively correlated with miR-155-5p expression in both the inflamed mucosa and adjacent mucosa of such patients, which can assume that this expression pattern may differentiate CD from UC. miR-155 is primarily expressed in the thymus and spleen and minimally detected in normal physiologic conditions, but its overexpression has been observed in IBD and colorectal cancer (CRC) [29].

Based on our results and the fact that if miR-155 expression is induced by inflammatory cytokines and toll-like receptor ligands that leads to dysregulations in tissue development and immune response, we can conclude the elevated level of miR-155-5p in the adjacent mucosa of patients with CD can reflect the initiation of inflammatory process in the tissue microenvironment. Other overexpressed miRNAs, which were negatively correlated with the expression of CDH1, were miR-146a-3p and miR-21-3p, found in inflamed tissue specimens in patients with CD in our studies. Until now, the influence of selected miRNAs (miR-9, miR-25, miR-92a) for E-cadherin expression has been linked to cancer cell migration and invasion [30] or development of metastasis only [31, 32]. Thereby, our results point out a potential role of the analyzed miRNAs—miR-21-3p, miR-146a-3p, and miR-155-5p—in EMT-mediated CD development or progression. However, further investigation and functional studies are required to confirm this hypothesis.

From the clinical point of view, we noticed that the expression of miR-146a-3p was negatively correlated with the duration of disease in UC patients, suggesting that miR-146a-3p could be used for monitoring the disease progress, rather than the biomarker which differentiates UC from CD. miR-146a is a nitric oxide-responsive miRNA, which amplifies inflammatory responses, including the NFκB pathway [33]. On the other hand, it was reported that the NFκB transcription factor after LPS-stimulation upregulates miR-125b-1 expression in cultured human epithelial cells [34]. In our studies, we compared expression of miR-125b-1-3p in adjacent tissue and inflamed tissue in CD and UC, and we noticed significantly higher expression of this miRNA in UC, but not in CD. Our results are in agreement with the findings of Valmiki et al. [16, 35] who discovered the significantly higher relative expression of miR-125b in UC biopsies compared to non-IBD controls. The authors showed that miR-125b targets and negatively regulates TRAF6 and A20, the two molecules of NFκB signaling pathway that modulate several genes involved in inflammation. It was suggested that suppression of TRAF6 and A20 genes mediated by miR-125b could be one of the factors contributing to the chronic inflammation in UC [35].

Our study has some limitations. Firstly, the study is based on bioinformatic prediction of potential interactions between selected miRNA and expression of E-cadherin; thus, results have to be confirmed by functional tests. Secondly, the number of patients is limited due to several inclusion and exclusion criteria; thereby, our results need to be validated on a larger population. However, even on analyzing the small cohort of patients, we emphasized the confirmation of correlations between miRNA and CDH1 expression, which prompt us to continue this study in the future.

To sum up, paying special attention to the relationship between miRNAs and E-cadherin as a core component of colonic epithelium and signaling epicenter that regulates cell behavior [18] will enable developing a new therapeutic strategy to stop or slow down the progression of IBD.

5. Conclusions

Our study clearly confirmed that miR-155-5p, miR-146a-3p, miR-125b-1-3p, and miR-21-3p are associated with IBD, and their correlations with CDH1 expression deserve special attention suggesting the involvement of epithelial to mesenchymal transition (EMT) in CD development, but not in UC.

Data Availability

The data that support the findings of this study are available from the corresponding author [MG], upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

The study was designed by KC, MG, and TD. KC and AS were involved in planning and supervised the project. MG wrote the manuscript with support from KC and AS. The draft was corrected by AS and KC. The samples and clinical data were collected by TD and JK. RNA was isolated by MG. The RT-qPCR was performed by WJ. MC. performed the statistical analysis. The acquisition of the financial support was done by KC. Results were interpreted by MG, WJ, and MC. All authors read and approved the final version of the manuscript.

Acknowledgments

The study was supported by Medical University of Lublin DS440 grant. We would like to thank Agnieszka Styczyńska for the editorial assistance and proofreading.