Pathway analysis of transcriptomic data shows immunometabolic effects of vitamin D
- Amadeo Muñoz Garcia1,2,
- Martina Kutmon1,3,
- Lars Eijssen1,
- Martin Hewison2,
- Chris T Evelo1,3 and
- Susan L Coort1⇑
- 1Department of Bioinformatics – BiGCaT, NUTRIM School of Nutrition and Metabolism in Translational Research, Maastricht University, Maastricht, The Netherlands
- 2Institute of Metabolism and Systems Research, The University of Birmingham, Birmingham, UK
- 3Maastricht Centre for System Biology (MaCSBio), Maastricht University, Maastricht, The Netherlands
- Correspondence should be addressed to S L Coort: susan.coort{at}maastrichtuniversity.nl
Abstract
Unbiased genomic screening analyses have highlighted novel immunomodulatory properties of the active form of vitamin D, 1,25-dihydroxyvitamin D (1,25(OH)2D). However, clearer interpretation of the resulting gene expression data is limited by cell model specificity. The aim of the current study was to provide a broader perspective on common gene regulatory pathways associated with innate immune responses to 1,25(OH)2D, through systematic re-interrogation of existing gene expression databases from multiple related monocyte models (the THP-1 monocytic cell line (THP-1), monocyte-derived dendritic cells (DCs) and monocytes). Vitamin D receptor (VDR) expression is common to multiple immune cell types, and thus, pathway analysis of gene expression using data from multiple related models provides an inclusive perspective on the immunomodulatory impact of vitamin D. A bioinformatic workflow incorporating pathway analysis using PathVisio and WikiPathways was utilized to compare each set of gene expression data based on pathway-level context. Using this strategy, pathways related to the TCA cycle, oxidative phosphorylation and ATP synthesis and metabolism were shown to be significantly regulated by 1,25(OH)2D in each of the repository models (Z-scores 3.52–8.22). Common regulation by 1,25(OH)2D was also observed for pathways associated with apoptosis and the regulation of apoptosis (Z-scores 2.49–3.81). In contrast to the primary culture DC and monocyte models, the THP-1 myelomonocytic cell line showed strong regulation of pathways associated with cell proliferation and DNA replication (Z-scores 6.1–12.6). In short, data presented here support a fundamental role for active 1,25(OH)2D as a pivotal regulator of immunometabolism.
- Received 16 November 2017
- Accepted 11 December 2017
- Made available online as an Accepted Preprint 12 December 2017
- © 2018 Society for Endocrinology