Transcriptome Profiling of Human Ulcerative Colitis Mucosa Reveals Altered Expression of Pathways Enriched in Genetic Susceptibility Loci

Human colonic mucosa altered by inflammation due to ulcerative colitis (UC) displays a drastically altered pattern of gene expression compared with healthy tissue. We aimed to understand the underlying molecular pathways influencing these differences by analyzing three publically-available, independently-generated microarray datasets of gene expression from endoscopic biopsies of the colon. Gene set enrichment analysis (GSEA) revealed that all three datasets share 87 gene sets upregulated in UC lesions and 8 gene sets downregulated (false discovery rate <0.05). The upregulated pathways were dominated by gene sets involved in immune function and signaling, as well as the control of mitosis. We applied pathway analysis to genotype data derived from genome-wide association studies (GWAS) of UC, consisting of 5,584 cases and 11,587 controls assembled from eight European-ancestry cohorts. The upregulated pathways derived from the gene expression data showed a highly significant overlap with pathways derived from the genotype data (33 of 56 gene sets, hypergeometric P = 1.49×10^–19). This study supports the hypothesis that heritable variation in gene expression as measured by GWAS signals can influence key pathways in the development of disease, and that comparison of genetic susceptibility loci with gene expression signatures can differentiate key drivers of inflammation from secondary effects on gene expression of the inflammatory process.     

Expression of heat shock protein 32 (hemoxygenase-1) in the normal and inflamed human stomach and colon: an immunohistochemical study

Heat shock protein 32 (Hsp32, hemoxygenase-1) is induced by reactive oxygen metabolites (ROM) and degrades heme leading to the formation of antioxidant bilirubin. Increased mucosal generation of ROM occurs in gastritis and inflammatory bowel disease. We aimed to assess mucosal expression of Hsp32 in normal stomach and colon and to test the hypothesis that disease-related differential expression occurs in inflamed tissue. Gastric body and antral mucosal biopsies were obtained from 33 patients comprising Helicobacter pylori-negative normal controls (n = 8), H pylori-negative gastritis patients (n = 11), and H pylori-positive gastritis patients (n = 14). Forty-seven archival colonic mucosal biopsies selected comprised normal histology (n = 10), active ulcerative colitis (UC) (n = 9), inactive UC (n = 8), active Crohn's disease (CD) (n = 8), inactive CD (n = 6), and other colitides (n = 6). Hsp32 expression in formalin-fixed sections was assessed by avidin-biotin peroxidase immunohistochemistry using a polyclonal rabbit anti-Hsp32 as the primary antibody. Immunohistochemical staining identified Hsp32 in all groups. Diffuse cytoplasmic staining was seen in gastric and colonic epithelial and lamina proprial inflammatory cells. Staining scores for Hsp32 were higher in antral H pylori-positive (P = 0.002) and H pylori-negative (P = 0.02) gastritis than in controls and in body H pylori-positive gastritis than in the other 2 groups (P < 0.01). Expression of Hsp32 was increased in active UC compared with inactive disease (P = 0.03) and normal controls (P = 0.02). In conclusion, Hsp32 is expressed constitutively in normal gastric and colonic mucosa, and differential expression occurs in these tissues when they are inflamed. Upregulation of Hsp32 may be an adaptive response to protect mucosa from oxidative injury in patients with gastritis and inflammatory bowel disease.     

Interleukin 35 (IL-35) and IL-37: Intestinal and peripheral expression by T and B regulatory cells in patients with Inflammatory Bowel Disease

The aim of the study was to characterize and to quantify peripheral and tissue. IL-35— and IL-37—producing cells in Inflammatory Bowel Disease (IBD) patients. We studied a total of 38 active UC, 31 inactive UC, 17 active CD, and 13 inactive CD and 50 non-inflamed tissues as control group. Gene expression was measured by real time polymerase chain reaction (RT-PCR) and protein expression was evaluated in tissue by immunohistochemistry and in peripheral blood mononuclear cells by flow cytometry. Higher levels of IL-35 was produced by intestinal regulatory B cells and circulating regulatory CD4⁺ and CD8⁺ T cells in active vs. inactive disease or healthy donors (P < 0.05). The IL-37 was conspicuously synthesized by circulating B cells, active natural killer cells and monocytes. These results suggest that down-regulation of inflammation in active IBD patients might be based on the increased expression of IL-35 and IL-37.     

Characterization of carboxypeptidase A6, an extracellular matrix peptidase

Carboxypeptidase A6 (CPA6) is a member of the M14 metallocarboxypeptidase family that is highly expressed in the adult mouse olfactory bulb and broadly expressed in embryonic brain and other tissues. A disruption in the human CPA6 gene is linked to Duane syndrome, a defect in the abducens nerve/lateral rectus muscle connection. In this study the cellular distribution, processing, and substrate specificity of human CPA6 were investigated. The 50-kDa pro-CPA6 is routed through the constitutive secretory pathway, processed by furin or a furin-like enzyme into the 37-kDa active form, and secreted into the extracellular matrix. CPA6 cleaves the C-terminal residue from a range of substrates, including small synthetic substrates, larger peptides, and proteins. CPA6 has a preference for large hydrophobic C-terminal amino acids as well as histidine. Peptides with a penultimate glycine or proline are very poorly cleaved. Several neuropeptides were found to be processed by CPA6, including Met- and Leu-enkephalin, angiotensin I, and neurotensin. Whereas CPA6 converts enkephalin and neurotensin into forms known to be inactive toward their receptors, CPA6 converts inactive angiotensin I into the biologically active angiotensin II. Taken together, these data suggest a role for CPA6 in the regulation of neuropeptides in the extracellular environment within the olfactory bulb and other parts of the brain.     

Dysregulation of Human β-Defensin-2 Protein in Inflammatory Bowel Disease

Background

Human β-defensin-2 (HBD2) is an antimicrobial peptide implicated in the pathogenesis of inflammatory bowel disease (IBD). Low copy number and concomitant low mRNA expression of the HBD2 gene have been implicated in susceptibility to colonic Crohn''s Disease (CD). We investigated the colonic distribution of HBD2 mRNA expression, and the contributions of genetic and environmental factors on HBD2 protein production.

Methodology/Principal Findings

We examined HBD2 mRNA expression at three colonic locations by microarray analysis of biopsies from 151 patients (53 CD, 67 ulcerative colitis [UC], 31 controls). We investigated environmental and genetic influences on HBD2 protein production using ex vivo cultured sigmoid colon biopsies from 69 patients (22 CD, 26 UC, 21 controls) stimulated with lipopolysaccharide (LPS) and/or nicotine for 24 hours. HBD2 and cytokines were measured in culture supernatants. Using DNA samples from these patients, regions in the HBD2 gene promoter were sequenced for NF-κB binding-sites and HBD2 gene copy number was determined. HBD2 mRNA expression was highest in inflamed (vs. uninflamed p = 0.0122) ascending colon in CD and in inflamed (vs. uninflamed p<0.0001) sigmoid colon in UC. HBD2 protein production was increased in inflamed UC biopsies (p = 0.0078). There was no difference in HBD2 protein production from unstimulated biopsies of CD, UC and controls. LPS-induced HBD2 production was significantly increased in CD (p = 0.0375) but not UC (p = 0.2017); this LPS-induced response was augmented by nicotine in UC (p = 0.0308) but not CD (p = 0.6872). Nicotine alone did not affect HBD2 production. HBD2 production correlated with IL8 production in UC (p<0.001) and with IL10 in CD (p<0.05). Variations in the HBD2 promoter and HBD2 gene copy number did not affect HBD2 production.

Significance/Conclusions

Colonic HBD2 was dysregulated at mRNA and protein level in IBD. Inflammatory status and stimulus but not germline variations influenced these changes.     

The role of adipose tissue-associated macrophages and T lymphocytes in the pathogenesis of inflammatory bowel disease

Inflammatory bowel diseases (IBDs) are chronic inflammatory disorders of the gastrointestinal tract that affect more than 3 million people worldwide, but the pathological etiology is still unknown. The overall purpose of our investigations was to elucidate the possibility of pathological causes of IBD, and therefore, we determined the difference of inflammatory cytokine profiles in adipose tissue macrophages (ATMs) and T lymphocytes (ATTs) obtained near active lesions of IBD; investigated whether the alteration in ATM activation induces genes involved in collagen formation; and evaluated the effects of fatty acid oxidation inhibitors on factors involved in inflammation and collagen production by ATMs in IBD. Adipose tissues (ATs) were collected near active lesions and also at the margin of resected segments of the bowel from IBD patients with ulcerative colitis (UC) and CD (n = 14/group). Normal appearing ATs from control subjects (n = 14) who had colon resection for adenocarcinoma were collected as far away from the cancer lesion as possible to rule out possible changes. Compared with inactive disease lesions, ATMs and ATTs from active lesions released more IL-6, IL-4 and IL-13. Treatments of cytokine IL-4 and/or IL-13 to ATMs reduced iNOS expression but increased Arg-I expression which were exacerbated when treated with T cell- and adipocyte-conditioned medium. However, fatty acid oxidation inhibitors prevented the effects of cytokines IL-4 and/or IL-13 on iNOS and Arg-I expressions. This study was the first to show the effect of IL-4 and IL-13 on collagen formation, through iNOS and Arg-I expressions, that was exacerbated in a condition that mimics in vivo condition of active lesions. Moreover, our study was the first to provide potential benefits of fatty acid oxidation inhibitors to ATMs on preventing collagen formation; thus, providing therapeutic implications for individuals with intestinal fibrosis and stricture lesions, although future study should be guaranteed to elucidate the underlying mechanisms.     

Computation of significance scores of unweighted Gene Set Enrichment Analyses

Background  

Gene Set Enrichment Analysis (GSEA) is a computational method for the statistical evaluation of sorted lists of genes or proteins. Originally GSEA was developed for interpreting microarray gene expression data, but it can be applied to any sorted list of genes. Given the gene list and an arbitrary biological category, GSEA evaluates whether the genes of the considered category are randomly distributed or accumulated on top or bottom of the list. Usually, significance scores (p-values) of GSEA are computed by nonparametric permutation tests, a time consuming procedure that yields only estimates of the p-values.     

Integrated miRNA and mRNA Expression Profiling in Inflamed Colon of Patients with Ulcerative Colitis

Background

Ulcerative colitis (UC) is associated with differential colonic expression of genes involved in immune response (e.g. IL8) and barrier integrity (e.g. cadherins). MicroRNAs (miRNAs) are regulators of gene expression and are involved in various immune-related diseases. In this study, we investigated (1) if miRNA expression in UC mucosa is altered and (2) if any of these changes correlate with mucosal mRNA expression. Integration of mRNA and miRNA expression profiling may allow the identification of functional links between dysregulated miRNAs and their target mRNA.

Methodology

Colonic mucosal biopsies were obtained from 17 UC (10 active and 7 inactive) patients and 10 normal controls. Total RNA was used to analyze miRNA and mRNA expression via Affymetrix miRNA 2.0 and Affymetrix Human Gene 1.0ST arrays, respectively. Both miRNA and gene expression profiles were integrated by correlation analysis to identify dysregulated miRNAs with their corresponding predicted target mRNA. Microarray data were validated with qRT-PCR. Regulation of IL8 and CDH11 expression by hsa-miR-200c-3p was determined by luciferase reporter assays.

Results

When comparing active UC patients vs. controls, 51 miRNAs and 1543 gene probe sets gave significantly different signals. In contrast, in inactive UC vs. controls, no significant miRNA expression differences were found while 155 gene probe sets had significantly different signals. We then identified potential target genes of the significantly dysregulated miRNAs and genes in active UC vs. controls and found a highly significant inverse correlation between hsa-miR-200c-3p and IL8, an inflammatory marker, and between hsa-miR-200c-3p and CDH11, a gene related to intestinal epithelial barrier function. We could demonstrate that hsa-miR-200c-3p directly regulates IL8 and CDH11 expression.

Conclusion

Differential expression of immune- and barrier-related genes in inflamed UC mucosa may be influenced by altered expression of miRNAs. Integrated analysis of miRNA and mRNA expression profiles revealed hsa-miR-200c-3p for use of miRNA mimics as therapeutics.     

Comparative proteomic studies on the pathogenesis of human ulcerative colitis

Ulcerative colitis (UC) is a chronic inflammatory disorder primarily affecting the colon mucosa. Its etiology and pathogenesis remain unclear. We used 2-DE and MS to identify differentially expressed proteins among the UC active, UC inactive, nonspecific colitis, and normal colon mucosa. Thirteen down-regulated and six up-regulated proteins were identified. Of the down-expressed proteins, eight (heat-shock protein 90 (HSPA9B), heat-shock protein 60 (HSPD1), H+-transporting two-sector ATPase (ATP5B), prohibitin (PHB), mitochondrial malate dehydrogenase (MDH2), voltage-dependent anion-selective channel protein 1 (VDAC1), thioredoxin peroxidase (PRDX1), and thiol-specific antioxidant (PRDX2)) were mitochondrial proteins, three (ATP5B, MDH2, triosephosphate isomerase) were involved in energy generation, three (PRDX1, PRDX2, SELENBP1) were cellular antioxidants, and six (HSPD1, HSPA9B, PRDX1, PRDX2, PHB, VDAC1) were stress-response proteins. Transmission electron microscopy revealed pathological alterations of mitochondrial ultrastructures even before the global colonocyte changes in the UC colon mucosa. PHB, an essential mitochondrial component protein, was down-expressed in the disease active as well as inactive colon mucosa from the patients of UC, indicative of an early event of mitochondrial changes during UC development. In contrast, aberrant activation of NFAT and ectopic expression of potential immunogenic proteins (tumor rejection antigen 1 and poliovirus receptor related protein 1) were found in the UC-diseased colon mucosa. Our findings suggest the implications of colonocyte mitochondrial dysfunction and perturbed mucosa immune regulation in the pathogenesis of UC and provide potential targets for the development of a new therapy.     

PAGE: Parametric Analysis of Gene Set Enrichment

Background  

Gene set enrichment analysis (GSEA) is a microarray data analysis method that uses predefined gene sets and ranks of genes to identify significant biological changes in microarray data sets. GSEA is especially useful when gene expression changes in a given microarray data set is minimal or moderate.     

Combined analysis of RNA‐sequence and microarray data reveals effective metabolism‐based prognostic signature for neuroblastoma

The relationship between metabolism reprogramming and neuroblastoma (NB) is largely unknown. In this study, one RNA‐sequence data set (n = 153) was used as discovery cohort and two microarray data sets (n = 498 and n = 223) were used as validation cohorts. Differentially expressed metabolic genes were identified by comparing stage 4s and stage 4 NBs. Twelve metabolic genes were selected by LASSO regression analysis and integrated into the prognostic signature. The metabolic gene signature successfully stratifies NB patients into two risk groups and performs well in predicting survival of NB patients. The prognostic value of the metabolic gene signature is also independent with other clinical risk factors. Nine metabolism‐related long non‐coding RNAs (lncRNAs) were also identified and integrated into the metabolism‐related lncRNA signature. The lncRNA signature also performs well in predicting survival of NB patients. These results suggest that the metabolic signatures have the potential to be used for risk stratification of NB. Gene set enrichment analysis (GSEA) reveals that multiple metabolic processes (including oxidative phosphorylation and tricarboxylic acid cycle, both of which are emerging targets for cancer therapy) are enriched in the high‐risk NB group, and no metabolic process is enriched in the low‐risk NB group. This result indicates that metabolism reprogramming is associated with the progression of NB and targeting certain metabolic pathways might be a promising therapy for NB.     

Development and validation of a five-lncRNA signature with prognostic value in colon cancer

Dysregulation of long noncoding RNAs (lncRNAs) has been found in a large number of human cancers, including colon cancer. Therefore, the implementation of potential lncRNAs biomarkers with prognostic prediction value are very much essential. GSE39582 data set was downloaded from database of Gene Expression Omnibus. Re-annotation analysis of lncRNA expression profiles was performed by NetAffx annotation files. Univariate and multivariate Cox proportional analyses helped select prognostic lncRNAs. Algorithm of random survival forest-variable hunting (RSF-VH) together with stepwise multivariate Cox proportional analysis were performed to establish lncRNA signature. The log-rank test was carried out to analyze and compare the Kaplan-Meier survival curves of patients’ overall survival (OS). Receiver operating characteristic (ROC) analysis was used for comparing the survival prediction regarding its specificity and sensitivity based on lncRNA risk score, followed by calculating the values of area under the curve (AUC). The single-sample GSEA (ssGSEA) analysis was used to describe biological functions associated with this signature. Finally, to determine the robustness of this model, we used the validation sets including GSE17536 and The Cancer Genome Atlas data set. After re-annotation analysis of lncRNAs, a total of 14 lncRNA probes were obtained by univariate and multivariate Cox proportional analysis. Then, the RSF-VH algorithm and stepwise multivariate Cox analysis helped to build a five-lncRNA prognostic signature for colon cancer. The patients in group with high risk showed an obviously shorter survival time compared with patients in group with low risk with AUC of 0.75. In addition, the five-lncRNA signature can be used to independently predict the survival of patients with colon cancer. The ssGSEA analysis revealed that pathways such as extracellular matrix-receptor interaction was activated with an increase in risk score. These findings determined the strong power of prognostic prediction value of this five-lncRNA signature for colon cancer.     

Bilateral sympathetic stellate ganglionectomy attenuates myocardial remodelling and fibrosis in a rat model of chronic volume overload

Reducing sympathetic neurohormone expression is a key therapeutic option in attenuating cardiac remodelling. Present study tested the feasibility of attenuating cardiac remodelling through reducing sympathetic neurohormone level by partial cardiac sympathetic denervation in a rat model of chronic volume overload. Male Sprague‐Dawley rats were randomized into sham group (S, n = 7), aortocaval fistula group (AV, n = 7), and aortocaval fistula with bilateral sympathetic stellate ganglionectomy group (AD, n = 8). After 12 weeks, myocardial protein expression of sympathetic neurohormones, including tyrosine hydroxylase, neuropeptide Y, growth associated protein 43, and protein gene product 9.5, were significantly up‐regulated in AV group compared to S group, and down‐regulated in AD group. Cardiac remodelling was aggravated in AV group compared to S group and attenuated in AD group. The myocardial deposition of extracellular matrix, including collagen I and III, was enhanced in AV group, which was reduced in AD group. Myocardial angiotensin II and aldosterone expressions were significantly up‐regulated in AV group and down‐regulated in AD group. Our results show that bilateral sympathetic stellate ganglionectomy could attenuate cardiac remodelling and fibrosis by down‐regulating sympathetic neurohormones expression in this rat model of chronic volume overload.     

The Effect of pH and Ion Channel Modulators on Human Placental Arteries

Chorionic plate arteries (CPA) are located at the maternofetal interface where they are able to respond to local metabolic changes. Unlike many other types of vasculature, the placenta lacks nervous control and requires autoregulation for controlling blood flow. The placental circulation, which is of low-resistance, may become hypoxic easily leading to fetal acidosis and fetal distress however the role of the ion channels in these circumstances is not well-understood. Active potassium channel conductances that are subject to local physicochemical modulation may serve as pathways through which such signals are transduced. The aim of this study was to investigate the modulation of CPA by pH and the channels implicated in these responses using wire myography. CPA were isolated from healthy placentae and pre-contracted with U46619 before testing the effects of extracellular pH using 1 M lactic acid over the pH range 7.4 - 6.4 in the presence of a variety of ion channel modulators. A change from pH 7.4 to 7.2 produced a 29±3% (n = 9) relaxation of CPA which increased to 61±4% at the lowest pH of 6.4. In vessels isolated from placentae of women with pre-eclampsia (n = 6), pH responses were attenuated. L-methionine increased the relaxation to 67±7% (n = 6; p<0.001) at pH 6.4. Similarly the TASK 1/3 blocker zinc chloride (1 mM) gave a maximum relaxation of 72±5% (n = 8; p<0.01) which compared with the relaxation produced by the TREK-1 opener riluzole (75±5%; n = 6). Several other modulators induced no significant changes in vascular responses. Our study confirmed expression of several ion channel subtypes in CPA with our results indicating that extracellular pH within the physiological range has an important role in controlling vasodilatation in the human term placenta.