Skip to main content

High-level integration of murine intestinal transcriptomics data highlights the importance of the complement system in mucosal homeostasis

Abstract

Background

The mammalian intestine is a complex biological system that exhibits functional plasticity in its response to diverse stimuli to maintain homeostasis. To improve our understanding of this plasticity, we performed a high-level data integration of 14 whole-genome transcriptomics datasets from samples of intestinal mouse mucosa. We used the tool Centrality based Pathway Analysis (CePa), along with information from the Reactome database.

Results

The results show an integrated response of the mouse intestinal mucosa to challenges with agents introduced orally that were expected to perturb homeostasis. We observed that a common set of pathways respond to different stimuli, of which the most reactive was the Regulation of Complement Cascade pathway. Altered expression of the Regulation of Complement Cascade pathway was verified in mouse organoids challenged with different stimuli in vitro.

Conclusions

Results of the integrated transcriptomics analysis and data driven experiment suggest an important role of epithelial production of complement and host complement defence factors in the maintenance of homeostasis.

Background

The mammalian gastrointestinal (GI) tract is crucial for the digestion and absorption of nutrients, energy metabolism, and homeostasis of the gut barrier and mucosal immunity. A number of specialized adaptations of the mammalian mucosal immune system have evolved to maintain a peaceful co-existence with the microbial symbionts while responding appropriately to prevent infection by enteric pathogens [1, 2]. Changes in external factors like the diet or intake of medication can influence microbiota ecology but also host metabolic processes [3,4,5,6,7,8,9,10,11,12,13,14,15]. The intestinal epithelium plays an important role in orchestrating innate defences [16, 17] and signalling to the numerous cells of the immune system located underneath the epithelial layer [18,19,20].

The gut functionalities described above are attributed to groups of genes organised into various functional pathways [21,22,23,24] responding to physiological changes. These pathways can be modulated by enteric infection, toxic compounds in food or produced by the microbiota, ionic and osmotic changes as well as substantial variations in nutrient availability. We hypothesized that several pathways are involved in maintaining homeostasis, the regulation of which depends on the type of perturbation. The diverse range of changing conditions encountered at the intestine would require a high-level of functional plasticity compared to other tissues. This theory is supported by the fact that a higher number of genes are specifically expressed in the gut mucosa than that in the heart, liver, kidney, and other organs that carry out a narrower range of functions [25,26,27]. Transcriptional responses of the intestinal mucosa to individual stimuli or perturbations have been documented widely in literature [15, 28, 29]. However, little is known about which key biological pathways provide functional plasticity in the intestinal mucosa. Interest in understanding this plasticity stems from the current trend to develop (dietary) interventions to optimise gut health and reduce the risk of disease.

Therefore, it is essential to investigate the functional plasticity of mucosal tissues at the functional genomic level in terms of pathways. Such an approach may aid to identify key sets of biosynthetic and signalling pathways involved in the mucosal responses, but also to identify the commonalities and differences in the expression of pathways responding to various environmental and physiological perturbations. To investigate this, we used publicly available gene expression datasets generated from mouse intestinal tissues exposed to orally administered challenges. From the results of the analysis on these datasets we identified the pathway “Regulation of Complement Cascade” that appears to play an important role in the functional plasticity of the intestinal epithelial response to different nutritional, microbial, and chemical challenges.

The complement system consists of several inactive pre-proteins produced in the liver that circulate in the blood which are crucial for efficient clearance of invading organisms. It is part of the innate immune system and activation of the complement cascade plays a key role in the opsonisation of micro-organisms to increase phagocytosis by macrophages and neutrophils at the sites of infection [30]. To avoid complement injury to autologous tissues, complement activation is controlled by a number of fluid-phase and cell surface proteins. Given the importance of this pathway in the GI system and in our high-level data integration we validated the results of our analysis with an in vitro experiment on mouse intestinal organoids. This experimental validation allowed us to, i) validate our data driven experimental design approach, ii) fulfil our aim of gaining more understanding of the functional plasticity of the GI tract.

Results

Classification of intestinal gene expression datasets

We identified 14 publicly available datasets meeting the search criteria. We classified them into three broad stimulation categories: Diet (7 experiments); Drug (3 experiments); and Immune Challenge (4 experiments) based on the type of intervention. Stimulations that were given as part of the feed of the animal were classified as ‘Diet’. Stimulation with a substance that is used as medication was classified into the ‘Drug’ category. Any substance that elicits a strong immune response was classified as an ‘Immune Challenge’. It could be argued that the DSS challenge belongs to the Drug category rather than Immune Challenge category as it can be used as a drug. However, since the effects of DSS likely result from a primary epithelial damage leading to the translocation of bacterial antigens, here we have chosen to classify it as an Immune Challenge.

These 14 experiments comprised of 37 experimental conditions (Additional file 1: Table S1), where a condition is defined as a unique combination of an inbred mouse strain, a specific intervention and an intestinal tissue sampled at a certain time point, as depicted in Fig. 1.

Fig. 1
figure1

Experimental datasets: The 37 conditions from 14 experiments (with 17 GEO accession numbers) used in this study are detailed in a timeline based on the age of the mice. Mice are selected to be part of an experiment based on weight, hence their age can vary within a range. The stars denote the start of the intervention, an empty star indicates the range of age when the intervention starts, when the age is not the same for all the animals in the group. Triangles denote the end of the interventions, an empty triangle indicates the start of the range of age of the animals. Challenges have been divided in three categories (colour coded): Diet, Drug and Immune Challenge. The names given for each dataset are abbreviated to show the challenge category in the first part of the name, the tissue sampled at the end (SI: small intestine; C: colon) and the text in the middle indicates the nature of challenge. Additional detailed explanations for the abbreviated condition names and the control conditions are given in Additional file 1: Table S1

Significant pathway results for all the datasets

In order to identify pathways specifically regulated by the challenges we used a modified version of the Centrality based Pathway Analysis (CePa) algorithm [31] and the ‘in-reach’ and ‘out-reach’ centrality options. Each analysis was performed on a single comparison, where a comparison is made between a stimulated condition versus the corresponding control in that experiment. We used the Reactome database [32] for the pathway information, which is arranged heirarchically from ‘root’ pathways (very broad), to more specific, ‘leaf’ pathways. We only work with the leaf pathways for the pathway analysis.

For all 37 conditions analysed, 710 pathways were significantly enriched, see Additional file 2: Table S2 for p-values of these pathways. The majority of these pathways (84%) were significant in both centrality measures whereas about 11% were significant only in the ‘in-reach’ centrality and 5% significant in only the ‘out-reach’ centrality. An overview of the pathway analysis results is given in Additional file 4: Figure S1. This figure shows that the responses are partially influenced by the tissue that was sampled since the red points (large intestine) mostly separate from the green points (small intestine). The number of signifcantly regulated pathways for each condition is highly variable, with the maximum being 377 from the experiment with Heme (Diet.Add.Heme.C) and the minimum being 37 in the experiment where 20% of the energy in the diet came from fat (Diet.Fat.20.P.SI). The average number of signifcantly regulated pathways is highest in the conditions belonging to the Diet category with 132 pathways, the second highest is Immune Challenge with 101 pathways followed by Drug with 87 pathways on average.

Comparison of the significant pathways in the three experimental challenge categories

There was a large overlap in the results between the three categories (Fig. 2). In addition, there were several pathways unique to each challenge category, the largest number of pathways were identified in the Diet category, which also had the largest number of experimental conditions in our study. There is a large number (212) of pathways that are shared among all challenge categories (Fig. 2). These 212 pathways belong to 24 of the 27 root pathways. The distribution of the 212 leaf pathways among the 24 root pathways of Reactome is shown in Table 1. The roots with the largest number of results are Metabolism, Disease, and Signalling pathways, but the proportion of common leaf pathways is similar to the proportion of all the leaf pathways in the roots in the Reactome database.

Fig. 2
figure2

Number of leaf pathways enriched in differentially expressed genes in the three challenge categories. The three circles are indicative of the significant leaf pathways in the datasets belonging to one of the three challenge categories Diet, Drug, or Immune Challenge. The 212 common pathways of all the three classes are indicated in the centre

Table 1 The 24 root pathways common to the three challenge categories. The names of the root pathways are given in the first column. The second column shows the number of leaf pathways in a particular root and this number is shown as a percentage of all the leaf pathways in the root pathway in the database

Regulation of pathways shared by the challenge categories

In order to investigate the most differently regulated pathways among the three conditions, the 212 common pathways were ranked based on a Difference Score. The Difference Score was calculated based on the mean node scores of the pathway nodes. The node scores of the pathway are a t-statistic which is a differential value between the experimental condition and the control. A mean of the node scores in the pathway under experimental conditions within one challenge category was calculated to end up with three scores for a pathway. The sum of the difference between the three mean values, the Difference Score was used to rank the list of common pathways. Table 2 shows the top 10 pathways ranked by this method and their corresponding Difference Scores.

Table 2 Top 10 pathways with the most difference in expression between the three perturbation classes. The second column shows the specific experimental conditions in the perturbation classes that were most different. The Difference Score was calculated using the t-values of the pathway nodes in the given conditions

The largest differences between the three challenge categories and gene expression of nodes were found in the innate immunity pathway, ‘Regulation of Complement Cascade’ (Fig. 3). All the genes in the Reactome pathway ‘Regulation of Complement Cascade’, with their differential expression values, are provided in Additional file 3: Table S3. The Regulation of Complement Cascade pathway was significantly regulated in 17 conditions in vivo (11 Diet, 3 of 7 Drug and 3 out of 5 Immune Challenge conditions). In the Diet conditions containing added fat, the complement factor genes mentioned above were also increased in expression (Fig. 3, Additional file 3: Table S3). In the Drug category, penicillin (Drug.Penicillin.SI), expression of the complement pathway genes highlighted in Fig. 3 were unaffected or mostly down-regulated whereas some were upregulated by the herbal drug conditions (Drug.HerbalDrug.SI, Drug.HerbalDrug.SPF.C). The biggest change was in the Immune Challenge category infection, where Salmonella upregulated the complement factors C4, C2, C3. These complement factors are required for activation of the complement cascade via the classical pathway and mannose lectin pathway (Fig. 3). C3 and factor B which are required for activation of the alternative pathway were also upregulated by some conditions in the Immune Challenge category (Fig. 3). Although C5 expression was only moderately upregulated by some of the Immune Challenge conditions, C6 was strongly upregulated. The other complement factors forming the membrane attack complex (MAC) on the surface of microorganisms were not strongly regulated under any of the conditions (Additional file 3: Table S3). The host protection factors (Fig. 3) CD55, CD46 and factor H, which are important for protection of host cell membranes when complement activation is triggered by microbes, displayed similar expression patterns as the complement factors (Fig. 3, Additional file 3: Table S3).

Fig. 3
figure3

The three complement pathways leading to enhanced phagocytosis of microorganisms. The binding of C3b to a receptor expressed on the surface of phagocytes and formation of the C5 convertase which generates chemotactic factors C5a and C3a, and the membrane attack complex for lysis of microbial membranes are depicted in the figure. Names of enzymatic products or complexes are shown. Common gene names are shown in boxes and are shaded in red when positively regulated in the datasets included in this study and green when negatively regulated. Each panel shows the complement pathway; however, it is superimposed with the pathway gene expression measurements for three different experimental conditions, one from each challenge category. Panel a represents the expression of the genes in the condition Diet.Fat.45.P.SI, panel b, Drug.Penicillin.SI and panel c, IC.STyph.4 day.SL1344.C. The legend shows the range of differential regulation. Complement-mediated defence mechanisms are shown in grey boxes. Modified from Microbes Online by Srijana Khanal [33]

Response of representative stimulants of the three challenge categories in mice intestinal organoids

To investigate the regulation of genes in the pathway ‘Regulation of Complement Cascade’ under similar conditions as tested in vivo, we performed experiments on adult stem cell derived ileal organoids from mouse. The organoids contained all the main epithelial cell lineages found in the tissue of origin, including comprehensive components of the complement cascade [34,35,36].

As defined stimuli, we used TNFα, an inflammatory cytokine-induced by infection or activation of inflammatory pathways, bacterial flagellin, an agonist of an innate immune receptor Toll Like Receptor 5 (TLR5) and a pharmacological agonist of PPARα a transcription factor and a major regulator of lipid metabolism.

The inflammatory cytokine TNFα which is induced by infection (e.g. with pathogenic Salmonella) induces expression of all complement related genes except C5, C8GH and CR2 (Fig. 4). In contrast, none of the genes were significantly altered in expression by incubation with flagellin, despite its ability to activate TLR5 signalling on HEK reporter cells expression TLR (data not shown). None of the genes related to the Regulation of Complement Cascade pathway were altered by the PPARα agonist. Instead, we observed the PPARα agonist significantly altered the expression of host receptors and CFI which are involved in protection of the host from complement activation.

Fig. 4
figure4

Expression of 10 chosen genes from the ‘Regulation of Complement Cascade’ pathways with significance calculated with ANOVA. Each graph contains information on different genes, the x-axis contains information on the treatment of the organoids and the y-axis has the fold change of the control genes. Data were analysed using Prism statistical software (v5.0, Graphpad, San Diego, US), measured for normality using the Kolmogorov-Smirnov test, and represented as Box and Whisker plots. A t-statistic test was performed on the RT-qPCR results of the 10 genes using the same methods as on the nodes in the gene set pathway analysis of CePa. All data were considered significantly different from the Blank (indicated in grey) when P < alpha (0.05) and indicated with * (P < 0.05 = *, P < 0.01 = **, P < 0.001 = ***)

Discussion

By integrating the results of experiments in which intestinal homeostasis was perturbed by completely different challenges, including probiotics, antibiotics, infectious agents, and major dietary components, we were able to investigate the plasticity of the GI tissue in terms of engaging various (biological) pathways. To the best of our knowledge, this type of study, has not yet been performed on this scale, focusing on one tissue and different types of challenges. By grouping the different challenge conditions in categories comprising drugs, dietary ingredients or potentially inflammatory agents, we grouped the responses of the mucosal tissue to facilitate broader comparisons. The results revealed pathways which are regulated in all categories. In addition, we observed large differences in the expression profiles of pathway genes between different exposure conditions, in some cases, irrespective of the challenge category. By focussing on the commonly regulated pathways, we show that the gut mucosa employs similar pathway systems. However, these pathway systems are used in different combinations and with different intra-pathway gene expression profiles, to respond to different exposures.

High-level data integration and pathway level analysis

The most important criteria for dataset selection was the age at which the mice were sampled, because the mucosal immune system and intestinal microbiota of mammals is known to change dramatically around weaning [37]. This holds true for mice [38, 39], therefore we only included datasets where the mice were sampled two or more weeks after weaning. The inclusion criteria provided datasets that are comparable, but they still differed in many aspects such as the use of microarray platforms, in sampling of the tissue and in the control conditions. By using controls within an experiment, we expect to only eliminate differences caused by platforms and retain biological influences like the type of tissue and the perturbation. Therefore, we used a high-level data integration method that started with the identification of differently expressed pathways as detected within an individual experiment and/or experimental condition. Before the pathway analysis the discriminating factor between the datasets was the platform in which the data was measured (data not shown). After the high-level data integration, the main differences among the datasets is the sampled tissue and the type of perturbation, see Additional file 4: Figure S1.

We used the CePa algorithm that considers a pathway’s topology by using different network centrality measures. Based on the biological information behind the pathways, we decided to use two centrality calculations, ‘In-reach’ and ‘Out-reach’ to capture regulation of pathways down-stream (important for signalling pathways) and up-stream (for metabolic pathways) respectively.

Most of the results were significant in both the centrality calculations irrespective of the type of pathway. This apparent indifference to topology is also observed by Bayerlova et al. [40] in a different pathway database using a variety of algorithms. In the aforementioned study, among the algorithms that used pathway topology, the CePa GSA algorithm performed consistently well. But, as also discussed in Khatri et al. [41], in order to make the most use of pathway topology it is important to be able to better annotate the edges between pathway nodes.

Several pathways were regulated by all three challenge categories

The results of the integrated pathway analysis show a notable overlap in the pathway response between the three challenge categories. The Diet category has the highest number of leaf pathway results, and also the highest number of experimental conditions. Surprisingly, most of the results from Drug and Immune Challenges were also shared by the Diet class. The results demonstrate that there is a group of pathways that are commonly regulated by interventions which perturb homeostasis. These common pathways contribute to a major extent towards the capability of the intestinal mucosa to display a high-level of functional plasticity. Most of the transcriptomics data used in this study came from intestinal scrapings which is greatly enriched for different types of epithelial cells involved in innate immunity and cross-talk with the immune cells in the lamina propria [2]. The different functions of these epithelial cells contribute to the functional plasticity of the epithelium. The results of this study revealed another layer of plasticity which is based on the specific use of a common set of pathways. These common pathways are significantly differentially expressed compared to controls in at least one experimental condition in each category. Simplifying the regulatory output of a whole pathway is difficult due to issues like up-regulation of the expression of inhibitory molecules.

The pathway regulation of complement Cascade is regulated by multiple intestinal challenges

One of the 212 common pathways which responds differently to the challenge categories is the ‘Regulation of Complement Cascade’. This pathway showed the largest difference in node expression profiles between the three challenge categories as shown in Table 2. The local production of complement factors must be important in intestinal homeostasis as the pathway shows maximum difference between the challenge categories and is regulated in several experimental conditions. Although the complement cascade is mentioned in two of the experiments used in this study [42, 43], the effects of the experimental conditions on the complement cascade were not explored in detail.

Complement factors involved in complement activation by one of three pathways are produced in the liver and enter the circulation. Complement factors reach tissue sites of infection through acute inflammation which results in the exudation of fluid and plasma proteins and an emigration of leukocytes into the extravascular compartment. Our observation that key complement factors involved in the complement activation pathways (e.g. C2, C3, C4, factor B), complement control and host protection (CD55, CD46, CFI, CFH), are increased in expression by epithelial cells exposed to infectious challenge, suggests that local complement production may be needed in the intestine as an early defence mechanism against encounters with microorganisms due to a dysfunctional barrier or infection. This idea has been proposed before [44, 45] especially with regards to inflammatory conditions and here we have explored this hypothesis in our dataset. This hypothesis was supported by our finding that transcription of complement system genes was regulated in mouse organoids in response to agonists of different signalling pathways. Only complement factors involved in the early stages of pathway activation and opsonization were strongly upregulated under inflammatory challenge conditions. Activation of the complement pathways in the mucosal tissues would lead to early opsonization of microorganisms and production of chemokines such as C5a to attract immune cells. The observed upregulation of host factors involved in protection of autologous cell membranes from complement damage is also compatible with the hypothesis that complement activation occurs in the mucosal tissue. The upregulation of C2, C3, C4, and factor B but not C5 and the complement factors (C6-C9) required for generation of the membrane attack complex (MAC), which is an important effector protein, is consistent with other literature describing complement factor expression in colorectal carcinoma cell lines stimulated with various cytokines [46, 47]. C3 and C4 transcripts have been localised to intestinal crypts in biopsies from Crohn’s patients [48]. This is consistent with the notion that epithelial cells can be induced to express complement factors needed for opsonisation of invading bacteria, but not cell lysis.

In the Diet category of challenges high fat diets strongly upregulated transcription of complement pathway genes, which may be due to a low-grade inflammation and a hyperpermeable gut [49, 50]. Interestingly, depletion of the microbiota with penicillin reduced pathway expression, suggesting that the microbiota contribute to ‘tonic’ stimulation of the complement related pathways via stimulation of innate immunity.

To confirm that the intestinal epithelium could express complement factors in response to signalling pathways targeted by the dietary challenges, we stimulated small intestinal crypt derived mouse organoids with TNFα, flagellin, an agonist of PPARα, or culture medium as a control and compared the relative transcript abundance of selected complement pathway genes by reverse transcription polymerase chain reaction (RT-PCR). TNFα is known to be induced by invasive infection by enteric infection with Salmonella [51], which was used in the Immune Challenge category [52]. Furthermore, receptors for TNFα are present on intestinal cells and signal in response to TNFα [53]. TLR5 was chosen because flagellin was orally administered in one of the studies in immune challenge category and PPARα agonist was selected because it was administered in one of the studies in the diet category with the aim of altering lipid metabolism [42]. One of the characteristics of inflammation is compromised barrier function, leading to a cascade of events in the lamina propria, e.g. TNFα secretion by invading immune cells, or bacterial fragment translocation [2]. As some TLR receptors may only be active or signal via the basolateral membrane of enterocytes [2, 54] we stimulated intact organoids. Genes to be observed were selected based on biological significance and fluctuation in the datasets that were analysed and RT-qPCR was performed on those mice genes (written in italics to differentiate from human genes). Stimulation with the inflammatory cytokine TNFα increased transcription of complement factors C2, C3 and regulatory proteins Cfi, Cd46 and Cd55, whereas expression of C5, C8 and Cfh was not significantly changed and Cd59 was significantly down-regulated. The agonist of PPARα significantly up-regulated transcription of regulatory proteins Cd55, Cr2, Cd46, and Cfi and significantly down-regulated expression of Cd59. Surprisingly, flagellin (which was shown to activate TLR5 in an intestinal cancer cell line [55]) did not significantly alter expression of any of these genes. The reasons for this are unclear but may be due to aberrant expression or regulation of TLR in intestinal cancer cells which are known to display biological variations such as aneuploidy, chromosome rearrangements or mutations [56]. It has also been proposed that TLR5 signalling is tightly controlled in epithelial cells to avoid chronic inflammatory responses to bacterial MAMPs from the intestinal lumen, and that expression of this receptor is exclusively present on Paneth cells in the small intestine [54]. Having TLR5 exclusively expressed in Paneth cells could explain the low responsiveness to flagellin, since this cell type is abundant at low levels in organoids [57]. Interestingly, activation of the PPARα pathway also increased expression of protective factors Cd55, Cd46, Cfi and Cr2 which allows the complement system to play a role in B cell activation and maturation. This links PPARα to regulation of complement cascade in the gut and the effects of high fat diets on this pathway in vivo.

CD55 and CD46 have other functions, besides their role in protection of host cells from complement activation which may be relevant for intestinal homeostasis. CD55 binds to the neutrophil receptor CD97 expressed on neutrophils to promote neutrophil migration through the epithelium [58]. Binding of antibodies to CD46 on Caco-2 cells was shown to induce intracellular signalling and improved cell proliferation and wound healing [59].

Conclusions

In conclusion, high-level data integration of transcriptomics datasets from intestinal tissue from in vivo experiments was a valuable approach to identify common pathways associated with functional plasticity and intestinal homeostasis. The identified pathways are regulated in different combinations to generate different physiological responses, or genes within a pathway are differentially regulated contributing further to the plasticity. The “Regulation of Complement Cascade” pathway is one of many pathways regulated by multiple intestinal challenges suggesting it is an important mechanism in the periphery of the intestine, which might have poor access to circulating complement components from the blood. Epithelial expression of complement factors involved in opsonisation and chemotaxis of host phagocytes, but not formation of the MAC complex, indicates a primary function in opsonisation of microbes and chemotaxis of host immune cells. Coincident with the increased intestinal expression of complement factors is the expression of host factors involved in complement control and protection such as CD55 and CD46 which have secondary functions in innate immunity and wound healing.

Methods

Datasets

The R tool GEOmetadb [60] was used to search Gene Expression Omnibus (GEO) [57, 61] for publicly available datasets generated from intestinal samples of mice. We selected 14 transcriptomics experiments (17 GEO datasets) from the 450 available (as of 07-07-2015) which analysed any intestinal tissue section from an intervention in weaned mice. All selected experiments used single channel microarrays with at least 3 biological replicates and were published on GEO between 2006 and 2014. Three of the datasets were obtained using Illumina microarray platforms whereas the others were obtained using versions of the Affymetrix platform. Most experimental data were obtained from inbred C57BL6J mice, but there were two experiments that used BALB/c mice, one that used the 129S1/svlmj mouse strain and another that used IQI mice. Eleven experiments were performed on female mice and three experiments on male mice.

There were two experiments where the fat content of the diet was increased by reducing the carbohydrate portion of the diet [11, 62]. The former experiment provided three different levels of fat (20, 30 and 45% of total energy in the diet) and performed transcriptomics on three sections of the small intestine. The latter experiment provided 60% fat as total energy in the diet and measured gene expression in RNA isolated from two halves of the small intestine. One experiment tested five different fibres which were substituted for part of the corn starch in the diet and the transcriptomics data was generated from the colon [63]. The rest of the dietary interventions were additives or supplements to the mouse diet. In one experiment, dietary heme was added to a high fat diet and the response was measured in the colon [64]. In another experiment quercetin was added to the standard diet and the response was measured in two parts of the small intestine (jejunum and ileum), as well as in the colon [65]. Another experiment added a synthetic PPARα agonist to the diet and a microarray analysis was performed on the small intestine [42]. One experiment tested a probiotic on two strains of mice and two sections of the intestine (small and large intestine) [66]. There are three experiments where a drug was administered to the animals. In one of these experiments a herbal drug was tested on wild-type and specific pathogen free mice, and the response measured in the small intestine and colon [43]. The other 2 drug experiments involved administration of antibiotics. One study utilised a mix of several antibiotics to strongly deplete the abundance of gut bacteria, and investigated gene expression in the colon [67]. In the other experiment with antibiotics, a low dose of penicillin was administered daily to the animals from an early age and the diet was also changed before measuring the response in the small intestine [68]. There were 2 challenge experiments, one with Salmonella Typhimurium [52], where the colon was sampled andanother involved the parasite Giardia duodenalis [69] where the small intestine was sampled. The effect of flagellin from Salmonella enterica serovar Typhimurium [70] on the intestinal immune system and the severity of DSS-induced colitis was investigated [71] in the colon. More details on the experimental conditions and the control conditions can also be found in Additional file 1: Table S1.

Data pre-processing

Using GEOQuery [72], we downloaded normalised datasets of preselected experiments from GEO. In each experiment, the normalization was performed with one of the following methods, GCRMA, RMA, MAS5 or quantile normalization. Our high-level integration approach does not require uniform normalization, so in each case we preferred the method chosen by the authors of the original study. The probes were mapped to mouse Entrez identifiers using the annotation files from the platform that was used for microarray analysis. After a quality check using Principle Component Analysis plots, these mice gene identifiers were then mapped to their human homologs using the NCBI HomoloGene database [RRID:SCR_002924]. Code snippets for these steps can be found on the GitHub page https://github.com/nirupamaBenis/PathwayLevelDataIntegration.

Pathway database

The analysis of all the datasets was performed using pathways from the Reactome database [RRID:SCR_003485], a freely accessible and a manually curated database available in different formats. Pathways from Reactome were downloaded in the BioPAX [RRID:SCR_009881] (Biological Pathway Exchange) [73] format (version 51) from the official website. These pathways were then converted to a pathway catalogue object in R that can be used by the pathway analysis algorithm. This was accomplished by using the pathway2Graph function from the R package rBiopaxParser [RRID:SCR_002744] [74, 75].

All pathways in the Reactome database are arranged in a hierarchy, larger ‘root’ pathways consist of more and more specific pathways, ending in several ‘leaf’ pathways. This hierarchy is depicted in a simplified cartoon in the inset of Additional file 5: Figure S2. The main image in Additional file 5: Figure S2 shows a network of all the root pathways in Reactome version 51. These 27 root pathways contain 1639 pathways within their hierarchy, of which 950 are leaf pathways.

Pathway analysis

We used a modified version of the algorithm CePa (Centrality based Pathway analysis), which uses pathways as networks where the nodes in a pathway could be small molecules (compounds), macro-molecules (proteins or RNA) or complexes (more than one protein). The topological information of the pathway is used to assign weights to each node using centralities. The user can choose between one of four centrality measures, in-degree (number of edges that are directed towards the node), out-degree (number of edges that are directed outwards from the node), in-reach (longest path that brings information to the node) and out-reach (longest path that directs out of the node), with another option of giving equal weight to all nodes. This centrality information is used along with the expression data to give a list of significantly enriched pathways for given conditions vs their controls. There are two methods of using the expression data in the CePa package, Over-Representation Analysis (ORA) and Gene Set Analysis (GSA). ORA usually takes a list of differentially expressed genes which could be ranked with p-values or fold changes. GSA takes the entire matrix of expression values to find enriched pathways and this is the method we chose to use.

The gene expression data is mapped to the nodes of the pathways, when the node is a protein the expression value of the corresponding gene is used as such. When the node is a complex, the largest component from a Principle Component Analysis of the expression values of all the corresponding proteins is assigned as the node expression value. Subsequently, these expression values are inputted in a t-statistic to obtain a differential expression value for each node, which can be positive or negative based on the up- or down-regulation of that protein. This differential node value is multiplied with the centrality-based weight of the node to obtain a final node value. This calculation is performed for each of the nodes in the pathway and all these values are averaged to obtain a pathway level score. The pathway level score is then tested for significance by substituting random expression values in the same pathway calculations and comparing the obtained value with the original dataset and the value obtained with the randomized data. The fraction of the iterations on which a higher score is obtained with the randomized data is used to represent the p-value. This p-value calculation was modified from the original CePa function which randomizes the replicates of the tested conditions. As we work with a minimum of three samples per condition, we modified this calculation to be able to handle smaller sample sizes. The original algorithm randomized the data across samples in order to calculate the significance of a pathway score. We decided to randomize the expression values by genes so that there is a larger chance of the values being truly random and thus without a biological signal. Because the hierarchical nature of the pathway database implies dependence between the pathways, we decided not to perform a multiple testing correction.

We weighted the nodes with the in-reach and out-reach centrality calculations, because they assign higher weights to the nodes down-stream and up-stream of the pathway respectively. In this way, we can capture signalling pathways, where the effectors are more likely to be down-stream of the pathway. However, we did not rule out the metabolic pathways where the enzymes are generally up-stream in a pathway. The threshold of the p-values was set at 0.01 to compensate for the lack of multiple testing correction.

Intestinal organoid cultures

Three dimensional (3-D) crypt derived murine intestinal organoids were grown as described in literature [34, 76,77,78]. Briefly, a 2 cm duodenal section was opened longitudinally and washed in ice-cold phosphate-buffered saline solution (PBS). After scraping excess villi, the tissue was transferred to PBS containing 2.5 mM EDTA and incubated for 30 min. Following incubation, the sections were washed with PBS and remaining residue was passed on a 70 μm cell strainer, pelleted at 300 x g for 5 min, and suspended in matrigel basement membrane (Growth factor reduced, Corning) at a density of 50–100 crypts per 50 μl. After inversed polymerization at 37 °C for > 10 min, 600 μl basal culture medium (DMEM/F12) was added, enriched with mouse EGF, Hepes 1 M (Invitrogen), N-acetylcysteine (Sigma), B-27 (Thermo-Fisher), Noggin, and R-spondin. The culture was passaged 1:4 every 7 days by mechanical disruption and re-suspension in fresh Matrigel. All experiments were performed after 2 passages of the organoid cultures.

Stimulation of organoids and reverse transcriptase-quantitative PCR

The 3-D organoids were stimulated with TNFα (10 ng/ml), a PPARα agonist (WY14643 0.1% v/v), and flagellin (200 ng/ml) for 6 h before total RNA was extracted with the Qiagen mini-kit according to manufacturer’s instructions along with a 15 min DNAse step. Purity and integrity measurements were performed on a DS-11 spectrophotometer (DeNovix) and 1 μg total RNA was reverse transcribed into cDNA using a QScript kit (Quantabio). Quantification of gene expression (RT-qPCR) was performed using a Rotor-gene Q2 plex RT-cycler (Qiagen) on primers specified in Table 3 with the rotor-gene SYBR green PCR kit, also from Qiagen. These genes were selected based on their representative contribution to the pathway ‘Regulation of Complement Cascade’. Relative expression levels were calculated following methods described in [81] using individual amplification values, with 18S and β-Actin as endogenous control genes for normalization.

Table 3 Primers used for RT-qPCR. The information on the primers used to quantify 10 genes is given in each row along with the publication from which this sequence was obtained

Untreated 3-D organoids were used as control to obtain relative gene expression values of the 10 chosen complement pathway genes when stimulated by the three treatments. Genes encoding C3 and C5 were selected as they are key factors in the three complement activation pathways. Genes encoding CD46, CD55 and CD59 were chosen because they are involved in protection of host membranes when complement pathway is activated. The other five genes (C8GH, CFI, CFH, CR1, and C2) were chosen because their expression varied substantially under the 17 experimental conditions where the ‘Regulation of Complement Cascade’ pathway was significantly affected. Complement Factor H and Factor I are involved in the regulation of complement activation, C2 is a component required for activation of the classical and alternate pathways. The murine CR2 contains 25 exons; a common first exon is spliced to exon 2 and to exon 9 in transcripts encoding CR1 and CR2 which encode receptors binding complement complexes on host immune cells.

Availability of data and materials

All the datasets used in the analysis were retrieved from the online repository Gene Expression Omnibus and the identifiers for each dataset is provided in the manuscript as are the hyperlinks to publications based on the data. R code and detailed instructions to perform these analyses can be found at https://github.com/nirupamaBenis/PathwayLevelDataIntegration.

Abbreviations

BioPAX:

Biological Pathway eXchange

CePa:

CEntrality based Pathway Analysis

DMEM:

Dulbecco’s Modified Eagle Medium

DSS:

Dextran Sulfate Sodium

EDTA:

Ethylenediaminetetraacetic acid

EGF:

Epidermal Growth Factor

GCRMA:

Guanine Cytosine Robust Multi-array Analysis

GEO:

Gene Expression Omnibus

GI:

Gastrointestinal tract

GSA:

Gene Set Analysis

MAC:

Membrane Attack Complex

ORA:

Over-Representation Analysis

PBS:

Phosphate-Buffered Saline

RMA:

Robust Multi-array Averaging

RT-PCR:

Reverse Transcription Polymerase Chain Reaction

References

  1. 1.

    Hooper LV, Littman DR, Macpherson AJ. Interactions between the microbiota and the immune system. Science. 2012;336:1268–73.

    CAS  Article  Google Scholar 

  2. 2.

    Wells JM, Rossi O, Meijerink M, van Baarlen P. Epithelial crosstalk at the microbiota-mucosal interface. Proc Natl Acad Sci U S A. 2011;(Supplement 1):4607–14. https://doi.org/10.1073/pnas.1000092107.

    Article  Google Scholar 

  3. 3.

    Van den Abbeele P, Gérard P, Rabot S, Bruneau A, El Aidy S, Derrien M, et al. Arabinoxylans and inulin differentially modulate the mucosal and luminal gut microbiota and mucin-degradation in humanized rats. Environ Microbiol. 2011;13:2667–80. https://doi.org/10.1111/j.1462-2920.2011.02533.x.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Kelly CJ, Zheng L, Campbell EL, Saeedi B, Scholz CC, Bayless AJ, et al. Crosstalk between microbiota-derived short-chain fatty acids and intestinal epithelial HIF augments tissue barrier function. Cell Host Microbe. 2015;17:662–71. https://doi.org/10.1016/j.chom.2015.03.005.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Schokker D, Zhang J, Vastenhouw SA, Heilig HG, Smidt H, Rebel JM, et al. Long-lasting effects of early-life antibiotic treatment and routine animal handling on gut microbiota composition and immune system in pigs. PLoS One. 2015;10:e0116523.

    Article  Google Scholar 

  6. 6.

    Mayer EA, Savidge T, Shulman RJ. Brain-gut microbiome interactions and functional bowel disorders. Gastroenterology. 2014;146:1500–12. https://doi.org/10.1053/j.gastro.2014.02.037.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Collins SM, Bercik P. The relationship between intestinal microbiota and the central nervous system in Normal gastrointestinal function and disease. Gastroenterology. 2009;136:2003–14.

    Article  Google Scholar 

  8. 8.

    Burger-van Paassen N, Vincent A, Puiman PJ, van der Sluis M, Bouma J, Boehm G, et al. The regulation of intestinal mucin MUC2 expression by short-chain fatty acids: implications for epithelial protection. Biochem J. 2009;420:211–9. https://doi.org/10.1042/BJ20082222.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Atarashi K, Tanoue T, Oshima K, Suda W, Nagano Y, Nishikawa H, et al. Treg induction by a rationally selected mixture of clostridia strains from the human microbiota. Nature. 2013;500:232–6. https://doi.org/10.1038/nature12331.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Atarashi K, Tanoue T, Shima T, Imaoka A, Kuwahara T, Momose Y, et al. Induction of Colonic Regulatory T Cells by Indigenous Clostridium Species. Science. 2011;331:337–41. https://doi.org/10.1126/science.1198469.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    de Wit NJW, Boekschoten MV, Bachmair E-M, Hooiveld GJEJ, de Groot PJ, Rubio-Aliaga I, et al. Dose-dependent effects of dietary fat on development of obesity in relation to intestinal differential gene expression in C57BL/6J mice. PLoS One. 2011;6:e19145. https://doi.org/10.1371/2Fjournal.pone.0019145.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Serino M, Luche E, Gres S, Baylac A, Bergé M, Cenac C, et al. Metabolic adaptation to a high-fat diet is associated with a change in the gut microbiota. Gut. 2012;61:543–53. https://doi.org/10.1136/gutjnl-2011-301012.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Willing BP, Russell SL, Finlay BB. Shifting the balance: antibiotic effects on host-microbiota mutualism. Nat Rev Microbiol. 2011;9:233–43.

    CAS  Article  Google Scholar 

  14. 14.

    Niewold TA. The nonantibiotic anti-inflammatory effect of antimicrobial growth promoters, the real mode of action? A Hypothesis. Poult Sci. 2007;86:605–9 http://ps.oxfordjournals.org/content/86/4/605.abstract.

    CAS  Article  Google Scholar 

  15. 15.

    Schokker D, Zhang J, Zhang L, Vastenhouw SA, Heilig HGHJ, Smidt H, et al. Early-life environmental variation affects intestinal microbiota and immune development in new-born piglets. PLoS One. 2014;9:e100040. https://doi.org/10.1371/2Fjournal.pone.0100040.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Wells JM, Brummer RJ, Derrien M, MacDonald TT, Troost F, Cani PD, et al. Homeostasis of the gut barrier and potential biomarkers. Am J Physiol Gastrointest Liver Physiol. 2017;312:G171–93.

    Article  Google Scholar 

  17. 17.

    Peterson LW, Artis D. Intestinal epithelial cells: regulators of barrier function and immune homeostasis. Nat Rev Immunol. 2014;14:141–53. https://doi.org/10.1038/nri3608.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Wells JM, Loonen LMP, Karczewski JM. The role of innate signaling in the homeostasis of tolerance and immunity in the intestine. Int J Med Microbiol. 2010;300:41–8.

    CAS  Article  Google Scholar 

  19. 19.

    Akira S, Takeda K. Toll-like receptor signalling. Nat Rev Immunol. 2004;4:499–511. https://doi.org/10.1038/nri1391.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Sharma R, Young C, Neu J. Molecular modulation of intestinal epithelial barrier: contribution of microbiota. J Biomed Biotechnol. 2010;2010:305879.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Mootha VK, Lindgren CM, Eriksson K-F, Subramanian A, Sihag S, Lehar J, et al. PGC-1[alpha]-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73. https://doi.org/10.1038/ng1180.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9. https://doi.org/10.1038/75556.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    CAS  Article  Google Scholar 

  24. 24.

    Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2014;42(Database issue):D472–7.

    CAS  Article  Google Scholar 

  25. 25.

    Lindskog C, Linné J, Fagerberg L, Hallström BM, Sundberg CJ, Lindholm M, et al. The human cardiac and skeletal muscle proteomes defined by transcriptomics and antibody-based profiling. BMC Genomics. 2015;16:475. https://doi.org/10.1186/s12864-015-1686-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Tissue-based map of the human proteome. Science. 2015;347 http://science.sciencemag.org/content/347/6220/1260419.abstract.

    Article  Google Scholar 

  27. 27.

    Fagerberg L, Hallström BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, et al. Analysis of the human tissue-specific expression by genome-wide integration of Transcriptomics and antibody-based proteomics. Mol Cell Proteomics. 2014;13:397–406. https://doi.org/10.1074/mcp.M113.035600.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Kar SK, Jansman AJM, Benis N, Ramiro-Garcia J, Schokker D, Kruijt L. Dietary protein sources differentially affect microbiota, mTOR activity and transcription of mTOR signaling pathways in the small intestine. PLOS ONE. 2017;12:1–9. https://doi.org/10.1371/journal.pone.0188282.

    Article  Google Scholar 

  29. 29.

    Reikvam DH, Erofeev A, Sandvik A, Grcic V, Jahnsen FL, Gaustad P, et al. Depletion of murine intestinal microbiota: effects on gut mucosa and epithelial gene expression. PLoS One. 2011;6:e17996. https://doi.org/10.1371/252Fjournal.pone.0017996.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Ricklin D, Hajishengallis G, Yang K, Lambris JD. Complement: a key system for immune surveillance and homeostasis. Nat Immunol. 2010;11:785–97. https://doi.org/10.1038/ni.1923.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Gu Z, Wang J. CePa: an R package for finding significant pathways weighted by multiple network centralities. Bioinformatics. 2013;29:658–60.

    CAS  Article  Google Scholar 

  32. 32.

    Fabregat A, Sidiropoulos K, Garapati P, Gillespie M, Hausmann K, Haw R, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2016;44:D481–7.

    CAS  Article  Google Scholar 

  33. 33.

    Srijana Khanal. Complement pathways: types, Functions and Regulation – microbeonlin. 2017. https://microbeonline.com/complement-system-pathways-functions-regulation/. Accessed 1 May 2018.

    Google Scholar 

  34. 34.

    Sato T, Vries RG, Snippert HJ, van de Wetering M, Barker N, Stange DE, et al. Single Lgr5 stem cells build crypt–villus structures in vitro without a mesenchymal niche. Nature. 2009;459:262–5. https://doi.org/10.1038/nature07935.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Ootani A, Li X, Sangiorgi E, Ho QT, Ueno H, Toda S, et al. Sustained in vitro intestinal epithelial culture within a Wnt-dependent stem cell niche. Nat Med. 2009;15:701–6. https://doi.org/10.1038/nm.1951.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    van der Hee B, Madsen O, Smidt H, Wells JM. Congruence of location-specific transcriptional programs in intestinal organoids during long-term culture. bioRxiv. 2019:600940. https://doi.org/10.1101/600940.

  37. 37.

    Cera KR, Mahan DC, Cross RF, Reinhart GA, Whitmoyer RE. Effect of age, weaning and postweaning diet on small intestinal growth and jejunal morphology in young swine. J Anim Sci. 1988;66:574–84.

    CAS  Article  Google Scholar 

  38. 38.

    Gu X, Li D, She R. Effect of weaning on small intestinal structure and function in the piglet. Arch Anim Nutr. 2002;56:275–86.

    CAS  Google Scholar 

  39. 39.

    Kelly D, Smyth JA, McCracken KJ. Digestive development of the early-weaned pig. 2. Effect of level of food intake on digestive enzyme activity during the immediate post-weaning period. Br J Nutr. 1991;65:181–8.

    CAS  Article  Google Scholar 

  40. 40.

    Bayerlová M, Jung K, Kramer F, Klemm F, Bleckmann A, Beißbarth T. Comparative study on gene set and pathway topology-based enrichment methods. BMC Bioinformatics. 2015;16:334. https://doi.org/10.1186/s12859-015-0751-5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Khatri P, Sirota M, Butte AJ. Ten years of pathway analysis: current approaches and outstanding challenges. PLoS Comput Biol. 2012;8:e1002375. https://doi.org/10.1371/journal.pcbi.1002375.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Bünger M, van den Bosch HM, van der Meijde J, Kersten S, Hooiveld GJEJ, Müller M. Genome-wide analysis of PPARα activation in murine small intestine. Physiol Genomics. 2007;30:192–204 http://physiolgenomics.physiology.org/content/30/2/192.abstract.

    Article  Google Scholar 

  43. 43.

    Munakata K, Takashima K, Nishiyama M, Asano N, Mase A, Hioki K, et al. Microarray analysis on germfree mice elucidates the primary target of a traditional Japanese medicine juzentaihoto: acceleration of IFN-$α$ response via affecting the ISGF3-IRF7 signaling cascade. BMC Genomics. 2012;13:30. https://doi.org/10.1186/1471-2164-13-30.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Kopp ZA, Jain U, Van Limbergen J, Stadnyk AW. Do antimicrobial peptides and complement collaborate in the intestinal mucosa? Front Immunol. 2015;6:17. https://doi.org/10.3389/fimmu.2015.00017.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Sina C, Derer S. The intestinal complement system in inflammatory bowel disease: shaping intestinal barrier function. Semin Immunol. 2018;37:66–73. https://doi.org/10.1016/J.SMIM.2018.02.008.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Andoh A, Fujiyama Y, Bamba T, Hosoda S. Differential cytokine regulation of complement C3, C4, and factor B synthesis in human intestinal epithelial cell line, Caco-2. J Immunol. 1993;151:4239–47 http://www.ncbi.nlm.nih.gov/pubmed/8409399. Accessed 27 Oct 2017.

    CAS  PubMed  Google Scholar 

  47. 47.

    Bernet-Camard MF, Coconnier MH, Hudault S, Servin AL. Differential expression of complement proteins and regulatory decay accelerating factor in relation to differentiation of cultured human colon adenocarcinoma cell lines. Gut. 1996;38:248–53 http://www.ncbi.nlm.nih.gov/pubmed/8801206. Accessed 1 May 2018.

    CAS  Article  Google Scholar 

  48. 48.

    Laufer J, Oren R, Goldberg I, Horwitz A, Kopolovic J, Chowers Y, et al. Cellular localization of complement C3 and C4 transcripts in intestinal specimens from patients with Crohn’s disease. Clin Exp Immunol. 2000;120:30–7. https://doi.org/10.1046/J.1365-2249.2000.01168.X.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Bischoff SC, Barbara G, Buurman W, Ockhuizen T, Schulzke J-D, Serino M, et al. Intestinal permeability--a new target for disease prevention and therapy. BMC Gastroenterol. 2014;14:189.

    Article  Google Scholar 

  50. 50.

    Chassaing B, Gewirtz AT. Gut microbiota, low-grade inflammation, and metabolic syndrome. Toxicol Pathol. 2014;42:49–53.

    Article  Google Scholar 

  51. 51.

    Zhang Y-G, Wu S, Xia Y, Sun J. Salmonella -infected crypt-derived intestinal organoid culture system for host-bacterial interactions. Physiol Rep. 2014;2:e12147. https://doi.org/10.14814/phy2.12147.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Liu X, Lu R, Xia Y, Sun J. Global analysis of the eukaryotic pathways and networks regulated by Salmonella typhimurium in mouse intestinal infection in vivo. BMC Genomics. 2010;11:722. https://doi.org/10.1186/1471-2164-11-722.

    Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Rossi O, Karczewski J, Stolte EH, Brummer RJM, van Nieuwenhoven MA, Meijerink M, et al. Vectorial secretion of interleukin-8 mediates autocrine signalling in intestinal epithelial cells via apically located CXCR1. BMC Res Notes. 2013;6:431. https://doi.org/10.1186/1756-0500-6-431.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Price AE, Shamardani K, Lugo KA, Deguine J, Roberts AW, Lee BL, et al. A Map of Toll-like Receptor Expression in the Intestinal Epithelium Reveals Distinct Spatial, Cell Type-Specific, and Temporal Patterns. Immunity. 2018;49:560–575.e6. https://doi.org/10.1016/j.immuni.2018.07.016.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Gewirtz AT, Navas TA, Lyons S, Godowski PJ, Madara JL. Cutting edge: bacterial flagellin activates basolaterally expressed TLR5 to induce epithelial proinflammatory gene expression. J Immunol. 2001;167:1882–5. https://doi.org/10.4049/JIMMUNOL.167.4.1882.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Liu Y, Mi Y, Mueller T, Kreibich S, Williams EG, Van Drogen A, et al. Multi-omic measurements of heterogeneity in HeLa cells across laboratories. Nat Biotechnol. 2019;37:314–22. https://doi.org/10.1038/s41587-019-0037-y.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991–5. https://doi.org/10.1093/nar/gks1193.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Leemans JC, te Velde AA, Florquin S, Bennink RJ, de Bruin K, van Lier RAW, et al. The epidermal growth factor-seven transmembrane (EGF-TM7) receptor CD97 is required for neutrophil migration and host defense. J Immunol. 2004;172:1125–31 http://www.ncbi.nlm.nih.gov/pubmed/14707087. Accessed 1 May 2018.

    CAS  Article  Google Scholar 

  59. 59.

    Cardone J, Al-Shouli S, Kemper C. A novel role for CD46 in wound repair. Front Immunol. 2011;2:28. https://doi.org/10.3389/fimmu.2011.00028.

    Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Zhu Y, Davis S, Stephens R, Meltzer PS, Chen Y. GEOmetadb: powerful alternative search engine for the gene expression omnibus. Bioinforma. 2008;24:2798–800. https://doi.org/10.1093/bioinformatics/btn520.

    CAS  Article  Google Scholar 

  61. 61.

    Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10. https://doi.org/10.1093/nar/30.1.207.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Desmarchelier C, Dahlhoff C, Keller S, Sailer M, Jahreis G, Daniel H. C57Bl/6 N mice on a western diet display reduced intestinal and hepatic cholesterol levels despite a plasma hypercholesterolemia. BMC Genomics. 2012;13:84. https://doi.org/10.1186/1471-2164-13-84.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Lange K, Hugenholtz F, Jonathan MC, Schols HA, Kleerebezem M, Smidt H, et al. Comparison of the effects of five dietary fibers on mucosal transcriptional profiles, and luminal microbiota composition and SCFA concentrations in murine colon. Mol Nutr Food Res. 2015;59:1590–602.

    CAS  Article  Google Scholar 

  64. 64.

    IJssennagger N, Derrien M, van Doorn GM, Rijnierse A, van den Bogert B, Müller M, et al. Dietary Heme alters microbiota and mucosa of mouse Colon without functional changes in host-microbe Cross-talk. PLoS One. 2012;7:e49868. https://doi.org/10.1371/2Fjournal.pone.0049868.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Mutch DM, Crespy V, Clough J, Henderson CJ, Lariani S, Mansourian R, et al. Hepatic cytochrome P-450 reductase-null mice show reduced transcriptional response to quercetin and reveal physiological homeostasis between jejunum and liver. Am J Physiol Gastrointest Liver Physiol. 2006;291:G63–72 http://ajpgi.physiology.org/content/291/1/G63.abstract.

    CAS  Article  Google Scholar 

  66. 66.

    Mariman R, Tielen F, Koning F, Nagelkerken L. The probiotic mixture VSL#3 has differential effects on intestinal immune parameters in healthy female BALB/c and C57BL/6 mice. J Nutr. 2015;145:1354–61. https://doi.org/10.3945/jn.114.199729.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Reikvam DH, Derrien M, Islam R, Erofeev A, Grcic V, Sandvik A, et al. Epithelial-microbial crosstalk in polymeric Ig receptor deficient mice. Eur J Immunol. 2012;42:2959–70. https://doi.org/10.1002/eji.201242543.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Cox LM, Yamanishi S, Sohn J, Alekseyenko AV, Leung JM, Cho I, et al. Altering the intestinal microbiota during a critical developmental window has lasting metabolic consequences. Cell. 2015;158:705–21. https://doi.org/10.1016/j.cell.2014.05.052.

    CAS  Article  Google Scholar 

  69. 69.

    Tako EA, Hassimi MF, Li E, Singer SM. Transcriptomic analysis of the host response to Giardia duodenalis infection reveals redundant mechanisms for parasite control. MBio. 2013;4:e00660–13. https://doi.org/10.1128/mBio.00660-13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Carvalho FA, Nalbantoglu I, Aitken JD, Uchiyama R, Su Y, Doho GH, et al. Cytosolic flagellin receptor NLRC4 protects mice against mucosal and systemic challenges. Mucosal Immunol. 2012;5:288–98. https://doi.org/10.1038/mi.2012.8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Ahmad R, Chaturvedi R, Olivares-Villagomez D, Habib T, Asim M, Shivesh P, et al. Targeted colonic claudin-2 expression renders resistance to epithelial injury, induces immune suppression, and protects from colitis. Mucosal Immunol. 2014;7:1340–53. https://doi.org/10.1038/mi.2014.21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Davis S, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics. 2007;23:1846–7.

    Article  Google Scholar 

  73. 73.

    Demir E, Cary MP, Paley S, Fukuda K, Lemer C, Vastrik I, et al. The BioPAX community standard for pathway data sharing. Nat Biotechnol. 2010;28:935–42.

    CAS  Article  Google Scholar 

  74. 74.

    Benis N, Schokker D, Kramer F, Smits MA, Suarez-Diez M, Benis N, et al. Building pathway graphs from BioPAX data in R. F1000Research. 2016;5:2414. https://doi.org/10.12688/f1000research.9582.2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Kramer F, Bayerlová M, Klemm F, Bleckmann A, Beißbarth T. RBiopaxParser-an R package to parse, modify and visualize BioPAX data. Bioinformatics. 2013;29:520–2.

    CAS  Article  Google Scholar 

  76. 76.

    Sato T, Stange DE, Ferrante M, Vries RGJ, Van Es JH, Van den Brink S, et al. Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and Barrett’s epithelium. Gastroenterology. 2011;141:1762–72.

    CAS  Article  Google Scholar 

  77. 77.

    Dekkers JF, Wiegerinck CL, de Jonge HR, Bronsveld I, Janssens HM, de Winter-de Groot KM, et al. A functional CFTR assay using primary cystic fibrosis intestinal organoids. Nat Med. 2013;19:939–45. https://doi.org/10.1038/nm.3201.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    van der Hee B, Loonen LMP, Taverne N, Taverne-Thiele JJ, Smidt H, Wells JM. Optimized procedures for generating an enhanced, near physiological 2D culture system from porcine intestinal organoids. Stem Cell Res. 2018. https://doi.org/10.1016/j.scr.2018.02.013.

    CAS  Article  Google Scholar 

  79. 79.

    Li X, Ding F, Zhang X, Li B, Ding J. The expression profile of complement components in Podocytes. Int J Mol Sci. 2016;17:471. https://doi.org/10.3390/ijms17040471.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  80. 80.

    Nie F, Su D, Shi Y, Chen J, Wang H, Qin W, et al. A preliminary study on the role of the complement regulatory protein, cluster of differentiation 55, in mice with diabetic neuropathic pain. Mol Med Rep. 2015;11:2076–82. https://doi.org/10.3892/mmr.2014.2896.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8. https://doi.org/10.1038/nprot.2008.73.

    CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work has been financially supported by the Systems Biology Investment Programme of Wageningen University, KB-17-003.02-022. Funding from the European Union Seventh Framework Programme (FP7/2007–2013) project SysmedID under grant agreement 305564 is gratefully acknowledged. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data, and in writing the manuscript.

Author information

Affiliations

Authors

Contributions

NB performed the data acquisition and analysis and prepared the manuscript. JW and BH designed the validation experiment and contributed significantly to the content of the manuscript. MS helped with the direction and critical revision of the manuscript. SKK contributed to the design of the validation experiment and contributed to the biological interpretation and the manuscript. BH contributed to the design of the validation experiment and executed it besides providing input to the manuscript. VMdS and MSD guided the data analysis, helped with the interpretation of the results and contributed to the manuscript. DS helped with the data analysis, contributed to the biological interpretation of the results and contributed to the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Nirupama Benis.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Description of the 37 conditions. This table describes the abbreviated names of experimental conditions as used in the paper and the control conditions in each of the experiments.

Additional file 2: Table S2.

Significant pathways with p-values in each experimental condition both centralities. This table has the p-values of all the significant leaf pathways in the 37 experimental conditions in both the in-reach and out-reach centrality calculations. If the p-value of the pathway was above the threshold of 0.01 an empty space is shown in the table. The p-values of at least one of the centrality calculations have to be lower than the threshold.

Additional file 3: Table S3.

Differential expression of genes in “Regulation of Complement Cascade” pathway. This table contains the differential gene expression values of 24 genes in 17 experimental conditions. The differential values were obtained with a T test of the experimental condition vs the control in that experiment.

Additional file 4: Figure S1.

PCA of all the significant pathways over the experimental conditions. Green points represent p-values from the pathway analysis from the small intestine and the red ones from the large intestine. Circles represent experimental conditions from the Diet category, the triangles are from the Drug category and squares from the Immune Challenge category.

Additional file 5: Figure S2.

Network of Reactome root pathway. The nodes in this network represent the 27 root pathways as present in Reactome v51 and the edges indicate the ‘leaf’ pathways shared by connected root pathways. The thickness of the edges indicates the number of leaf pathways shared by the nodes. The nodes are labelled with the names of the root pathways and the number of enclosed leaf pathways is given between brackets. The inset shows a simplified example of root and leaf pathways, where the cartoon has one root pathway with three leaf pathways.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Benis, N., Wells, J.M., Smits, M.A. et al. High-level integration of murine intestinal transcriptomics data highlights the importance of the complement system in mucosal homeostasis. BMC Genomics 20, 1028 (2019). https://doi.org/10.1186/s12864-019-6390-x

Download citation

Keywords

  • Pathway analysis
  • Transcriptomics
  • Data integration
  • Intestine
  • Complement pathway
  • Homeostasis