Investigation of microRNAs in mouse macrophage responses to lipopolysaccharide-stimulation by combining gene expression with microRNA-target information
© Chiu and Wu 2015
Published: 9 December 2015
Toll-like receptors, which stimulated by pathogen-associated molecular patterns such as lipopolysaccharides (LPS), induces the releasing of many kinds of proinflammatory cytokines to activate subsequent immune responses. Plenty of studies have also indicated the importance of TLR-signalling on the avoidance of excessive inflammation, tissue repairing and the return to homeostasis after infection and tissue injury. The significance of TLR-signalling attracts many attentions on the regulatory mechanisms since several years ago. However, as newly discovered regulators, how and how many different microRNAs (miRNAs) regulate TLR-signalling pathway are still unclear.
By integrating several microarray datasets and miRNA-target information datasets, we identified 431 miRNAs and 498 differentially expressed target genes in bone marrow-derived macrophages (BMDMs) with LPS-stimulation. Cooperative miRNA network were constructed by calculating targets overlap scores, and a sub-network finding algorithm was used to identify cooperative miRNA modules. Finally, 17 and 8 modules are identified in the cooperative miRNA networks composed of miRNAs up-regulate and down-regulate genes, respectively.
We used gene expression data of mouse macrophage stimulated by LPS and miRNA-target information to infer the regulatory mechanism of miRNAs on LPS-induced signalling pathway. Also, our results suggest that miRNAs can be important regulators of LPS-induced innate immune response in BMDMs.
The immune system is composed of innate and adaptive immune systems, in which the machinery behind innate immune system is very ingenious due to the conserved ability to recognize microbial components during evolution process. Toll-like receptors (TLRs) are a kind of specialized sensor responsible for recognizing molecular patterns of pathogens and bridging innate and adaptive immunity by activation of antigen-presenting cells such as dendritic cells and macrophages. Upon exposure of molecular patterns of pathogens, TLRs initiate a cascade of intracellular signalling transduction that leads to secretion of inflammatory mediators, including cytokines and chemokines, inducing expression changes of a broad spectrum of genes that regulate defense function against pathogens.
TLR4, the best characterized member of TLR family, can detect lipopolysaccharide (LPS) from Gram-negative bacteria. Once TLR4 is stimulated by LPS, two independent signalling pathways (TRIF-dependent and MyD88-dependent) will be initiated, and many import transcription factors (TNF-α, IFN-β, IL-1β, etc.) will be activated as well. Although several mechanisms which are responsible for the regulation of TLR4-induced signalling pathways have been already well-studied, including physical interactions, changes of conformation, phosphorylation, ubiquitylation, and etc [1–3], how an anti-inflammatory response is induced and the processes by which inflammation remain unclear . Among these regulators, miRNAs have drew enough notice as a newly identified regulators in last several years.
MicroRNAs (miRNAs), a class of small non-coding RNAs, has emerged as important regulators in the development of immune and inflammatory responses [4–8]. The importance of miRNAs in the modulation of normal and pathological immune function has been shown in multiple studies in which deregulation of miRNAs was demonstrated to accompany diseases associated with excessive or uncontrolled inflammation. As a critical part of regulatory networks in innate immunity, the dysregulation of some specific miRNAs would modulate the regulation of inflammatory in stimulated immune cells. According to the literature , many miRNAs are reported to be dys-regulated in activated immune cells, such as miR-146a.
Recently, many studies have attempted to develop methods to understand miRNA co-operativity. Most of these studies however did not actually considered coexpression profiles of mRNAs and miRNAs. Considering that most miRNAs exert their functions through interactions with other miRNAs, an understanding of a miRNA network context using both co-expression pattern and the sequence complementarity between miRNAs and mRNAs is essential to discover the cooperative regulation of miRNAs.
In this paper, we used two kinds of data sets to infer the regulatory mechanism of miRNAs on LPS-induced signalling pathway. mRNA expression data and miRNA-target information data are used to identify reliable regulatory relationship between miRNAs and mRNA, and then cooperative miRNAs can be identified using hypergeometric test. Besides, all differentially expressed target genes were used to perform coherent functional analysis, the possible function of miRNAs which actually regulate gene expression could be hence inferred. Using these cooperative miRNAs and their common target genes, a cooperative miRNA network could be constructed and the cooperative miRNA modules could be also identified.
Materials and methods
The work flow for identification of co-regulating miRNAs in mouse macrophages with stimulation of LPS. Microarray expression data and miRNA-target information data are used to identify reliable regulatory relationship between miRNAs and mRNA, and then cooperative microRNAs can be identified using hypergeometric test. Besides, all differentially expressed target genes were used to perform coherent functional analysis, the possible function of miRNAs which actually regulate gene expression could be hence inferred. Using these cooperative miRNAs and their common target genes, a cooperative miRNA network could be constructed and the cooperative miRNA modules could be also identified.
In this study, five expression profile datasets and four miRNA-target information datasets are used. Expression datasets, which were downloaded from NCBI GEO database , are from five experiments generated individually by Barish et al , El Chartouni et al , Ghigo et al , Mages et al  and Millien et al , which were downloaded from NCBI GEO database . Barish et al used Illumina MouseWG-6 v2.0 expression beadchip to compare the changes of gene expression between wildtype and Bcl6 KO bone marrow-derived macrophages of mouse in the absence or presence of LPS (GSE24792). El Chartouni et al used Agilent-014868 Whole Mouse Genome Microarray to study TLR4 induced transcriptional responses in interleukin 4 primed mouse bone marrow-derived macrophages (GSE21895). Ghigo et al used Agilent-014868 Whole Mouse Genome Microarray to compare the difference of gene expression of mouse bone marrow-derived macrophage response to T. whipplei and LPS (GSE20210). Magges et al used Affymetrix Mouse Genome 430 2.0 Array to study the development of a state of refractoriness to a second stimulation in mouse bone marrow-derived macrophages treated with LPS (GSE8621). Millien et al used Illumina MouseWG-6 v2.0 expression beadchip to study the effects of mouse bone marrow derived macrophages treated with PAO, IFN-gamma, or LPS (GSE48609). On the other hand, miRNA-target gene interaction datasets were downloaded from Mouse Genome Informatics (MGI) database . These datasets mainly come from four different data sources, in which miRTarBase  data are validated, but microt-cds , miRDB  and Pictar  data are all predicted. The list of mouse genetic markers which also downloaded from MGI database was used as the list of all genes and miRNA in mouse.
Expression data preprocessing and differentially expressed genes identification
In each expression dataset, only wild-type samples in the absence or presence of LPS were used in subsequent data analyses. The value distribution service of GEO2R in NCBI GEO database was used to determined if this data are normalized and cross-comparable, and the results revealed that only the datasets from Barish et al study and from Millien et al study are not needed to be normalized. Other three datasets were quantile-normalized using R package "affyPLM" . For every expression datasets, differential expression analyses were performed using R package "limma" . Benjamini-Hochberg method was used to control the false discovery rate (FDR), and genes with adjusted p-value less than 0.05 were identified as differentially expressed (DE) genes.
Preprocessing of miRNA-target gene interaction data
Due to a large part of miRNA-target gene interaction data consists of predicted data, there are many false positives in these datasets. Therefore, only the data from miRTarBase and the data appear at least in two of three predicted data sources were retained, and then we obtained a list with 822697 miRNA-target gene interaction data. We obtained 712 miRNAs from this list, and then remove the miRNAs which could be found in the list of mouse genetic markers. Finally, 703 miRNAs were retained.
Gene set enrichment analysis
We carried out the gene set enrichment analysis to identify the miRNAs which actually participate in the regulation of target genes. For each miRNA, a 2 × 2 contingency table was built (whether the genes are targets of this miRNA or not versus whether the genes are differentially expressed) to access the enrichment level using Fisher's exact test. After all miRNAs were tested, Benjamini-Hochberg procedure for multiple test correction was used to correct p-values. This miRNA would be considered to have a significant number of differentially expressed target genes if the adjusted p-value of Fisher's exact test is less than 0.05, and this miRNA could be thus considered a reliable regulator in mouse macrophages under LPS-stimulated condition.
Coherent functional analysis
Cytoscape v3.2.1 (http://cytoscape.org/) with ClueGO v2.1.6 plug-in  was used to perform GO and KEGG pathway analyses. ClueGO let us to obtain the distribution of the list of differentially expressed target genes across the GO terms and pathways. We used right-sided hypergeometric test to calculate enrichment level and used Benjamini-Hochberg adjustment to correct p-values in multiple test. A GO term or pathway would be considered to be enriched in differentially expressed target genes if its adjusted p-value is less than 0.05. We analyzed up-regulated target genes and down-regulated target genes using these settings in GO and KEGG pathway analyses, respectively.
Construct cooperative miRNA network
, where Targets(i) and Targets(j) represent the set of target genes of miRNA i and miRNA j, respectively. As a result, the list of cooperative miRNAs and the target overlap scores were used to construct CMN.
Extraction of cooperative miRNA modules (CMM)
In this paper, a sub-network detection algorithm, Molecular COmplex DEtection algorithm (MCODE) , was applied to find possible coherent modules in the global cooperative miRNA network. The MCODE, which is a clustering algorithm based on graph theorem, is specifically designed to find complexes by identifying densely connected subgraphs in networks. MCODE algorithm is composed of three stages: vertex weighting, complex prediction and an optional post-processing step. The core clustering coefficient is used as the weights of nodes. Once the weights are computed, the algorithm traverses the weighted graph using greedy algorithm to isolate densely connected regions. The post-processing step filters or adds nodes based on connectivity criteria. MCODE has a parameter that specifies the size of clusters returned. All parameters in MCODE are applied with default values. Sub-networks are filtered if they do not contain at least a 2-core (graph of minimum degree 2). This approach allows us to assign one miRNA to multiple clusters, considering a biological principle that miRNAs can have multiple functions.
Differentially expressed mRNA
GO and KEGG enrichment of differentially expressed genes
Construction and characterization of cooperative miRNA networks
Cooperative miRNA modules in which miRNAs up-regulate genes.
Mir106b, Mirlet7b, Mir98, Mir297c, Mir9-1, Mir9-2, Mir200b, Mir467d, Mir9-3, Mir669f, Mir467c, Mir466n, Mir297a-4, Mir669m-2, Mir669m-1, Mir302a, Mir467b, Mirlet7a-1, Mir200a, Mir141, Mir297a-3, Mir669l, Mirlet7c-2, Mirlet7f-1, Mir467a-10, Mir467f, Mir467a-9, Mir669d-2, Mir669d, Mir467e, Mir669b, Mir467a-4, Mir93, Mir467a-3, Mir467a-1, Mir467a-2, Mir297b, Mir302d, Mir17, Mir467a-8, Mir467a-7, Mir467a-6, Mir467a-5
Mir294, Mir295, Mir302c, Mir466c-2, Mir466c-1, Mir466p, Mir291b, Mir291a, Mir466, Mir466e, Mir467g, Mir466b-4, Mir466m, Mir292, Mir466b-3, Mir466b-7, Mir466b-6, Mir466b-5, Mir466b-8, Mir539, Mir493, Mir466o, Mir106a, Mir20bm Mir302b
Mir20a, Mir153, Mir202, Mir322, Mirlet7i, Mirlet7c-1, Mir466i, Mirlet7e, Mir7b, Mir22, Mir7-2, Mirlet7f-2, Mir1961, Mirlet7d, Mir488, Mir19b-1, Mir320, Mir669p-1, Mir196b, Mir669p-2, Mirlet7a-2
Mir503, Mir216a, Mir883a, Mir421, Mir883b, Mir551b, Mir5619, Mir30c-2, Mir3065, Mir3069, Mir410, Mir128-2, Mir5624, Mir103-1, Mir742, Mir344e, Mir344d-2, Mir133a-2, Mir3074-1, Mir181b-2, Mir133a-1, Mir499, Mir743, Mir497, Mir3082, Mir300, Mir344d-1, Mir344d-3, Mir124a-3, Mir30b, Mir124a-1, Mir124a-2, Mir495, Mir448, Mir101a, Mir188, Mir30c-1, Mir27b, Mir382, Mir301b, Mir374c, Mir323, Mir384
Mir494, Mir429, Mir15b, Mir128-1, Mir5616, Mir1907, Mir107, Mir30a, Mir33, Mir103-2, Mir140, Mir222, Mir181b-1, Mir3074-2, Mir7-1, Mir216b, Mir181d, Mir16-1, Mir19b-2, Mir468, Mir15a, Mir743b, Mir96, Mir181c, Mir1912, Mir214, Mir674, Mir1906-1, Mir582, Mir27a, Mir1906-2, Mir301
Mir465b-1, Mir152, Mir32, Mir367, Mir465, Mir218-1, Mir465c-2, Mir465c-1, Mir758, Mir129-1, Mir30d, Mir129-2, Mir3066, Mir706, Mir204, Mir669c, Mir19a, Mir186, Mir465b-2, Mir1192, Mir875, Mir92b, Mir16-2, Mir669h, Mir673
Mir694, Mir199a-2, Mir199a-1, Mir199b
Mir212, Mir351, Mir125b-1, Mir125b-2
Mir449a, Mir218-2, Mir449b, Mir143, Mir1930, Mir3059, Mir34b, Mir5104, Mir23b, Mir211, Mir23a, Mir672, Mir224
Mir466f-2, Mir466f-1, Mir466f-3
Mir29a, Mir29b-1, Mir29c
Mir1981, Mir340, Mir142
Mir1a-2, Mir148a, Mir148b
Mir3107, Mir486, Mir326
Mir544, Mir377, Mir1903
Mir692-3, Mir692-1, Mir692-2
Mir149, Mir3085, Mir5627
Cooperative miRNA modules in which miRNAs down-regulate genes.
Mir129-1, Mir19a, Mirlet7i, Mir466l, Mir1912, Mirlet7e, Mir19b-2, Mir1961, Mir98, Mir350, Mir466n, Mirlet7d, Mir129-2, Mir1192, Mirlet7f-1, Mirlet7c-1, Mirlet7a-1, Mirlet7a-2, Mir96, Mirlet7c-2, Mir669m-2, Mir669m-1, Mir291a, Mir7-1, Mir497, Mir335, Mirlet7f-2, Mir19b-1, Mir204, Mir875, Mir883a, Mirlet7g, Mir743b, Mir743
Mir351, Mir5615-2, Mir93, Mirlet7b, Mir876, Mir470, Mir17, Mir677, Mir1983, Mir3063, Mir9-2, Mir320, Mir495, Mir9-3, Mir203, Mir182, Mir672, Mir449a, Mir488, Mir26a-2, Mir9-1, Mir330, Mir669k, Mir20b, Mir106a, Mir128-1, Mir466d, Mir466i, Mir218-1, Mir539, Mir883b, Mir153, Mir20a, Mir292
Mir465c-1, Mir465c-2, Mir465, Mir465b-2, Mir465b-1
Mir362, Mir128-2, Mir147, Mir338, Mir412, Mir125a, Mir322, Mir1941, Mir200b, Mir125b-2, Mir3097, Mir34b, Mir370, Mir125b-1, Mir133a-1, Mir133a-2, Mir1298
Mir124a-3, Mir5624, Mir124a-1, Mir673, Mir124a-2
Mir667, Mir103-2, Mir3078, Mir107
Mir5627, Mir3090, Mir5623, Mir346, Mir3087, Mir3113
Mir24-2, Mir149, Mir24-1
LPS-induced signalling must be tightly regulated to avoid excessive inflammation and to allow for tissue repair and the return to homeostasis after infection and tissue injury. Though miRNAs are critical regulators of immune responses, whether they are involved in LPS-induced signalling pathway and how is their expression regulated in mouse macrophages are still unclear. We combine mRNA expression data and miRNA-target information to infer the regulatory mechanism of miRNAs on LPS-induced signalling pathway. Also, our results suggest that miRNAs can be important regulators of LPS-induced innate immune response in BMDMs.
This study was supported by National Cheng Kung University and Ministry of Science and Technology of Taiwan MOST-103-2221-E-006-174-MY2.
Publication charges of this articles was funded by National Cheng Kung University and Ministry of Science and Technology of Taiwan MOST-103-2221-E-006-174-MY2.
This article has been published as part of BMC Genomics Volume 16 Supplement 12, 2015: Joint 26th Genome Informatics Workshop and 14th International Conference on Bioinformatics: Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S12.
- He X, Jia H, Jing Z, Liu D: Recognition of pathogen-associated nucleic acids by endosomal nucleic acid-sensing toll-like receptors. Acta Biochim Biophys Sin (Shanghai). 2013, 45 (4): 241-258.View ArticleGoogle Scholar
- O'Neill LAJ: When signaling pathways collide: positive and negative regulation of toll-like receptor signal transduction. Immunity. 2008, 29 (1): 12-20.PubMedView ArticleGoogle Scholar
- Kondo T, Kawai T, Akira S: Dissecting negative regulation of toll-like receptor signaling. Trends Immunol. 2012, 33 (9): 449-458.PubMedView ArticleGoogle Scholar
- O'Neill LA, Sheedy FJ, McCoy CE: Micrornas: the fine-tuners of toll-like receptor signalling. Nat Rev Immunol. 2011, 11 (3): 163-175.PubMedView ArticleGoogle Scholar
- Alam MM, O'Neill LA: Micrornas and the resolution phase of inflammation in macrophages. Eur J Immunol. 2011, 41 (9): 2482-2485.PubMedView ArticleGoogle Scholar
- Contreras J, Rao DS: Micrornas in inflammation and immune responses. Leukemia. 2012, 26 (3): 404-413.PubMedView ArticleGoogle Scholar
- Dai R, Ahmed SA: Microrna, a new paradigm for understanding immunoregulation, inflammation, and autoimmune diseases. Transl Res. 2011, 157 (4): 163-179.PubMedView ArticleGoogle Scholar
- Xiao C, Rajewsky K: Microrna control in the immune system: basic principles. Cell. 2009, 136 (1): 26-36.PubMedView ArticleGoogle Scholar
- Liu G, Abraham E: Micrornas in immune response and macrophage polarization. Arterioscler Thromb Vasc Biol. 2013, 33 (2): 170-177. doi:10.1161/ATVBAHA.112.300068PubMedView ArticleGoogle Scholar
- Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, Robertson CL, Serova N, Davis S, Soboleva A: Ncbi geo: archive for functional genomics data sets-update. Nucleic Acids Res. 2013, 41 (Database issue): 991-995.View ArticleGoogle Scholar
- Barish GD, Yu RT, Karunasiri M, Ocampo CB, Dixon J, Benner C, Dent AL, Tangirala RK, Evans RM: Bcl-6 and nf-kappab cistromes mediate opposing regulation of the innate immune response. Genes Dev. 2010, 24 (24): 2760-2765.PubMedView ArticleGoogle Scholar
- El Chartouni C, Rehli M: Comprehensive analysis of tlr4-induced transcriptional responses in interleukin 4-primed mouse macrophages. Immunobiology. 2010, 215 (9-10): 780-787.PubMedView ArticleGoogle Scholar
- Ghigo E, Barry AO, Pretat L, Al Moussawi K, Desnues B, Capo C, Kornfeld H, Mege J-L: Il-16 promotes t. whipplei replication by inhibiting phagosome conversion and modulating macrophage activation. PLoS One. 2010, 5 (10): 13561.-View ArticleGoogle Scholar
- Mages J, Dietrich H, Lang R: A genome-wide analysis of lps tolerance in macrophages. Immunobiology. 2007, 212 (9-10): 723-737.PubMedView ArticleGoogle Scholar
- Millien VO, Lu W, Shaw J, Yuan X, Mak G, Roberts L, Song L-Z, Knight JM, Creighton CJ, Luong A, Kheradmand F, Corry DB: Cleavage of fibrinogen by proteinases elicits allergic responses through toll-like receptor 4. Science. 2013, 341 (6147): 792-796.PubMedView ArticleGoogle Scholar
- Eppig JT, Blake JA, Bult CJ, Kadin JA, Richardson JE, MGDG: The mouse genome database (mgd): facilitating mouse as a model for human biology and disease. Nucleic Acids Res. 2015, 43 (Database issue): 726-736.View ArticleGoogle Scholar
- Hsu S-D, Tseng Y-T, Shrestha S, Lin Y-L, Khaleel A, Chou C-H, Chu C-F, Huang H-Y, Lin C-M, Ho S-Y, Jian T-Y, Lin F-M, Chang T-H, Weng S-L, Liao K-W, Liu C-C, Huang H-D: mirtarbase update 2014: an information resource for experimentally validated mirna-target interactions. Nucleic Acids Res. 2014, 42 (Database issue): 78-85.View ArticleGoogle Scholar
- Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, Vergoulis T, Koziris N, Sellis T, Tsanakas P, Hatzigeorgiou AG: Diana-microt web server: elucidating microrna functions through target prediction. Nucleic Acids Res. 2009, 37 (Web Server issue): 273-276.View ArticleGoogle Scholar
- Wong N, Wang X: mirdb: an online resource for microrna target prediction and functional annotations. Nucleic Acids Res. 2015, 43 (Database issue): 146-152.View ArticleGoogle Scholar
- Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microrna target predictions. Nat Genet. 2005, 37 (5): 495-500.PubMedView ArticleGoogle Scholar
- Bolstad B, Collin F, Brettschneider J, Simpson K, Cope L, Irizarry R, Speed T: Quality assessment of affymetrix genechip data. Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Edited by: Gentleman R, Carey V, Huber W, Irizarry R, Dudoit S. 2005, Springer, New York, 3: 33-47.View ArticleGoogle Scholar
- Smyth G: Bioinformatics and computational biology solutions using r and bioconductor. Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Edited by: Gentleman R, Carey V, Huber W, Irizarry R, Dudoit S. 2005, Springer, New York, 23: 397-420.View ArticleGoogle Scholar
- Bindea G, Galon J, Mlecnik B: Cluepedia cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013, 29 (5): 661-663.PubMedView ArticleGoogle Scholar
- Na Y-J, Kim JH: Understanding cooperativity of micrornas via microrna association networks. BMC Genomics. 2013, 14 (Suppl 5): 17-View ArticleGoogle Scholar
- Bader GD, Hogue CWV: An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003, 4: 2-PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.