Structural Remodeling of the Human Colonic Mesenchyme in Inflammatory Bowel Disease

Summary Intestinal mesenchymal cells play essential roles in epithelial homeostasis, matrix remodeling, immunity, and inflammation. But the extent of heterogeneity within the colonic mesenchyme in these processes remains unknown. Using unbiased single-cell profiling of over 16,500 colonic mesenchymal cells, we reveal four subsets of fibroblasts expressing divergent transcriptional regulators and functional pathways, in addition to pericytes and myofibroblasts. We identified a niche population located in proximity to epithelial crypts expressing SOX6, F3 (CD142), and WNT genes essential for colonic epithelial stem cell function. In colitis, we observed dysregulation of this niche and emergence of an activated mesenchymal population. This subset expressed TNF superfamily member 14 (TNFSF14), fibroblastic reticular cell-associated genes, IL-33, and Lysyl oxidases. Further, it induced factors that impaired epithelial proliferation and maturation and contributed to oxidative stress and disease severity in vivo. Our work defines how the colonic mesenchyme remodels to fuel inflammation and barrier dysfunction in IBD.

[1]  Alex E. Lash,et al.  Gene Expression Omnibus: NCBI gene expression and hybridization array data repository , 2002, Nucleic Acids Res..

[2]  R. Shivdasani,et al.  Wnt Secretion from Epithelial Cells and Subepithelial Myofibroblasts Is Not Required in the Mouse Intestinal Stem Cell Niche In Vivo , 2014, Stem cell reports.

[3]  Måns Magnusson,et al.  MultiQC: summarize analysis results for multiple tools and samples in a single report , 2016, Bioinform..

[4]  Hans Clevers,et al.  Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and Barrett's epithelium. , 2011, Gastroenterology.

[5]  S. Koch.  Extrinsic control of Wnt signaling in the intestine. , 2017, Differentiation; research in biological diversity.

[6]  A. Regev,et al.  Spatial reconstruction of single-cell gene expression , 2015, Nature Biotechnology.

[7]  D. Sahoo,et al.  Identification and prospective isolation of a mesothelial precursor lineage giving rise to smooth muscle cells and fibroblasts for mammalian internal organs, and their vasculature , 2012, Nature Cell Biology.

[8]  I. Pinchuk,et al.  PD-1 ligand expression by human colonic myofibroblasts/fibroblasts regulates CD4+ T-cell activity. , 2008, Gastroenterology.

[9]  K. Kaestner,et al.  Foxl1-Expressing Mesenchymal Cells Constitute the Intestinal Stem Cell Niche , 2015, Cellular and molecular gastroenterology and hepatology.

[10]  S. Itzkovitz,et al.  Subepithelial telocytes are an important source of Wnts that supports intestinal crypts , 2018, Nature.

[11]  Cheng Li,et al.  Adjusting batch effects in microarray expression data using empirical Bayes methods. , 2007, Biostatistics.

[12]  Aaron T. L. Lun,et al.  Scater: pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in R , 2017, Bioinform..

[13]  R. Flavell,et al.  Fibroblasts and myofibroblasts of the intestinal lamina propria in physiology and disease. , 2016, Differentiation; research in biological diversity.

[14]  Wei Shi,et al.  featureCounts: an efficient general purpose program for assigning sequence reads to genomic features , 2013, Bioinform..

[15]  K. Csiszȧr,et al.  Lysyl oxidases: a novel multifunctional amine oxidase family. , 2001, Progress in nucleic acid research and molecular biology.

[16]  Renaud Gaujoux,et al.  A flexible R package for nonnegative matrix factorization , 2010, BMC Bioinformatics.

[17]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[18]  Sarah A Teichmann,et al.  Computational assignment of cell-cycle stage from single-cell transcriptome data. , 2015, Methods.

[19]  Luís Torgo,et al.  Data Mining with R: Learning with Case Studies , 2010 .

[20]  Laleh Haghverdi,et al.  Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors , 2018, Nature Biotechnology.

[21]  Steven L Salzberg,et al.  HISAT: a fast spliced aligner with low memory requirements , 2015, Nature Methods.

[22]  J. Marioni,et al.  Pooling across cells to normalize single-cell RNA sequencing data with many zero counts , 2016, Genome Biology.

[23]  Max Kuhn,et al.  Building Predictive Models in R Using the caret Package , 2008 .

[24]  Fabian J. Theis,et al.  destiny: diffusion maps for large-scale single-cell data in R , 2015, Bioinform..

[25]  W. Fibbe,et al.  Mesenchymal stromal cells: sensors and switchers of inflammation. , 2013, Cell stem cell.

[26]  N. Ashley,et al.  Myofibroblasts are distinguished from activated skin fibroblasts by the expression of AOC3 and other associated markers , 2016, Proceedings of the National Academy of Sciences.

[27]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[28]  Richard A. Flavell,et al.  The Stromal Intervention: Regulation of Immunity and Inflammation at the Epithelial-Mesenchymal Barrier , 2017, Cell.

[29]  Jack Satsangi,et al.  Regional variation in gene expression in the healthy colon is dysregulated in ulcerative colitis , 2008, Gut.

[30]  N. Barker.  Adult intestinal stem cells: critical drivers of epithelial homeostasis and regeneration , 2013, Nature Reviews Molecular Cell Biology.

[31]  Terry P Yamaguchi,et al.  Wnt5a Potentiates TGF-β Signaling to Promote Colonic Crypt Regeneration After Tissue Injury , 2012, Science.

[32]  Andy Liaw,et al.  Classification and Regression by randomForest , 2007 .

[33]  J. Burch,et al.  The serosal mesothelium is a major source of smooth muscle cells of the gut vasculature , 2005, Development.

[34]  H. Clevers,et al.  Generation of BAC Transgenic Epithelial Organoids , 2013, PloS one.

[35]  R. Russell,et al.  Intestinal label-retaining cells are secretory precursors expressing Lgr5 , 2013, Nature.

[36]  Xavier Robin,et al.  pROC: an open-source package for R and S+ to analyze and compare ROC curves , 2011, BMC Bioinformatics.

[37]  Ryan M. Anderson,et al.  Periostin potently promotes metastatic growth of colon cancer by augmenting cell survival via the Akt/PKB pathway. , 2004, Cancer cell.

[38]  K. Basler,et al.  GLI1-expressing mesenchymal cells form the essential Wnt-secreting niche for colon stem cells , 2018, Nature.

[39]  S. Dudoit,et al.  Normalization of RNA-seq data using factor analysis of control genes or samples , 2014, Nature Biotechnology.

[40]  Shuifang Zhu,et al.  Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads , 2014, BMC Bioinformatics.

[41]  Sarah L. Brown,et al.  Myd88-dependent positioning of Ptgs2-expressing stromal cells maintains colonic epithelial proliferation during injury. , 2007, The Journal of clinical investigation.

[42]  S. Ward,et al.  A novel population of subepithelial platelet-derived growth factor receptor α-positive cells in the mouse and human colon. , 2013, American journal of physiology. Gastrointestinal and liver physiology.

[43]  C. M. van der Loos.  Chromogens in Multiple Immunohistochemical Staining Used for Visual Assessment and Spectral Imaging: The Colorful Future , 2010 .