Open Access Highly Accessed Research article

Changes in colon gene expression associated with increased colon inflammation in interleukin-10 gene-deficient mice inoculated with Enterococcus species

Matthew PG Barnett1*, Warren C McNabb23, Adrian L Cookson1, Shuotun Zhu4, Marcus Davy5, Bianca Knoch16, Katia Nones7, Alison J Hodgkinson8 and Nicole C Roy1

Author Affiliations

1 Food, Metabolism & Microbiology Section, AgResearch Grasslands, Tennent Drive, Palmerston North 4442, New Zealand

2 Food & Textiles Group, AgResearch Grasslands, Tennent Drive, Palmerston North 4442, New Zealand

3 Riddet Institute, Massey University, Palmerston North, NZ

4 Faculty of Medical and Health Sciences, The University of Auckland, 85 Park Road, Grafton, Auckland 1023, New Zealand

5 Plant & Food Research, Ruakura Research Centre, East Street, Hamilton 3214, New Zealand

6 Institute of Food, Nutrition & Human Health, Massey University, Tennent Drive, Palmerston North 4442, New Zealand

7 Plant & Food Research, Batchelar Road, Palmerston North 4442, New Zealand

8 Dairy Science & Technology, AgResearch Ruakura, East Street, Hamilton 3214, New Zealand

For all author emails, please log on.

BMC Immunology 2010, 11:39  doi:10.1186/1471-2172-11-39

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2172/11/39


Received:26 January 2010
Accepted:15 July 2010
Published:15 July 2010

© 2010 Barnett et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Inappropriate responses to normal intestinal bacteria may be involved in the development of Inflammatory Bowel Diseases (IBD, e.g. Crohn's Disease (CD), Ulcerative Colitis (UC)) and variations in the host genome may mediate this process. IL-10 gene-deficient (Il10-/-) mice develop CD-like colitis mainly in the colon, in part due to inappropriate responses to normal intestinal bacteria including Enterococcus strains, and have therefore been used as an animal model of CD. Comprehensive characterization of changes in cecum gene expression levels associated with inflammation in the Il10-/- mouse model has recently been reported. Our aim was to characterize changes in colonic gene expression levels in Il10-/- and C57BL/6J (C57; control) mice resulting from oral bacterial inoculation with 12 Enterococcus faecalis and faecium (EF) strains isolated from calves or poultry, complex intestinal flora (CIF) collected from healthy control mice, or a mixture of the two (EF·CIF). We investigated two hypotheses: (1) that oral inoculation of Il10-/- mice would result in greater and more consistent intestinal inflammation than that observed in Il10-/- mice not receiving this inoculation, and (2) that this inflammation would be associated with changes in colon gene expression levels similar to those previously observed in human studies, and these mice would therefore be an appropriate model for human CD.

Results

At 12 weeks of age, total RNA extracted from intact colon was hybridized to Agilent 44 k mouse arrays. Differentially expressed genes were identified using linear models for microarray analysis (Bioconductor), and these genes were clustered using GeneSpring GX and Ingenuity Pathways Analysis software. Intestinal inflammation was increased in Il10-/- mice as a result of inoculation, with the strongest effect being in the EF and EF·CIF groups. Genes differentially expressed in Il10-/- mice as a result of EF or EF·CIF inoculation were associated with the following pathways: inflammatory disease (111 genes differentially expressed), immune response (209 genes), antigen presentation (11 genes, particularly major histocompatability complex Class II), fatty acid metabolism (30 genes) and detoxification (31 genes).

Conclusions

Our results suggest that colonic inflammation in Il10-/- mice inoculated with solutions containing Enterococcus strains is associated with gene expression changes similar to those of human IBD, specifically CD, and that with the EF·CIF inoculum in particular this is an appropriate model to investigate food-gene interactions relevant to human CD.

Background

The term 'Inflammatory Bowel Disease' (IBD) refers to a heterogeneous collection of conditions characterized by chronic inflammation of the gastrointestinal tract, and includes Crohn's Disease (CD) and Ulcerative Colitis (UC) [1]. While there is some overlap in disease pathology, CD and UC also have distinct pathologic features; CD can, for example, affect any part of the gastrointestinal tract, whereas UC is confined to the colon and rectum, often causing diarrhea. The inflammation seen in CD is typically discontinuous, segmental and involves all layers of the intestinal wall. In UC, inflammation tends to be continuous and superficial, only affecting the mucosal layer of the colonic wall [2].

The exact etiology and pathogenesis of IBD is still unclear, although there is strong epidemiological evidence for a genetic contribution to disease susceptibility. Several candidate genes for IBD susceptibility have been identified, including nucleotide-binding oligomerization domain containing 2 (NOD2) [3-5], tumour necrosis factor (TNF) [6], members of the toll-like receptor (TLR) family [7], IL-4 [8] and IL-18 [9], and a number of genes encoding transporter molecules, such as the ATP-binding cassette, sub-family B (MDR/TAP), member 1 (ABCB1) [10,11] and solute carrier family 22 (organic cation/ergothioneine transporter), member 4 (SLC22A4) genes [12,13].

The IL-10 gene deficient (Il10-/-) mouse has been used as a model of IBD [14-21]. These mice, when bred onto a C57BL/6J (C57) background, have been reported to develop CD-like colitis by 12 weeks of age when raised under conventional conditions [19], while female 129 Ola × C57Il10-/- mice have been shown to develop colitis from 20 weeks of age under specific pathogen free (SPF) conditions [21].

The precise mechanism that results in inflammation in Il10-/- mice is unclear, although, as is the case in human IBD, there is evidence of an inappropriate inflammatory response to normal intestinal flora [22].

Clinical isolates of Enterococcus faecalis have been shown to induce IBD-like symtoms in germ-free Il10-/- mice [14,23,24]. Enterococcus species are a common component of the intestinal flora of healthy humans and animals [25-27], comprising up to 1% of the adult microflora [28]. Enterococcus faecalis and Enterococcus faecium are the two species most commonly detected in the human bowel [29-31], and both are known to carry a variety of virulence factors (reviewed in [25]) which may play a role in the establishment of inflammation.

Based on these published studies, and on our own observations of only mild inflammation in 12 week old Il10-/- mice (C57 background) that were raised under conventional conditions (M. P. G. Barnett, "unpublished observations"), we decided to establish bacterially-inoculated Il10-/- mice as a model of IBD in order to test food-gene interactions associated with IBD. We tested two hypotheses: (1) that oral inoculation of Il10-/- mice with a mixture of pure Enterococcus isolates (both faecalis and faecium), alone or combined with conventional intestinal flora derived from healthy and conventionally raised C57 mice, would result in greater and more consistent intestinal inflammation than that observed in Il10-/- mice not receiving this inoculation, and (2) that this inflammation would be associated with changes in colon gene expression levels in key pathways similar to those previously observed in human studies, and these mice would therefore be an appropriate model for human CD [32,33].

We have previously published body weight, histology and preliminary gene expression data using this mouse model [34]. Here we describe in detail gene expression changes in colonic tissue in response to bacterial inoculation in Il10-/- on a C57 background, measured using high density oligonucleotide microarrays.

Results

Animal body weight

There was no difference in mean body weight between Il10-/- and C57 mice at the start of the experiment (Il10-/- 18.9 ± 1.0 g; C57 18.8 ± 0.8 g). As reported previously [34], Il10-/- mice in the conventional conditions (C) and EF·CIF groups gained less weight during the course of the trial, both when compared to similarly-inoculated C57 mice, and when compared to Il10-/- mice in the other three treatment groups.

Intestinal Histology

Total intestinal histology results have been reported elsewhere, in which colon was shown to be the intestinal section most susceptible to inflammation [34]. A more detailed analysis of the colon histology results showed that inflammatory cell infiltration was the most prominent feature of the observed inflammation in inoculated Il10-/- mice (accounting for approximately 80% of total colon HIS), and was significantly higher in these animals, when compared with C57 mice in the same group, or when compared with non-inoculated Il10-/- mice (i.e., SPF and C groups). Furthermore, Il10-/- mice inoculated with EF inocula (either EF or EF·CIF) showed significantly higher tissue destruction, both when compared with similarly-inoculated C57 mice, and when compared with Il10-/- under SPF and C conditions (Figure 1).

thumbnailFigure 1. Summary of the colon histological injury score (HIS) in Il10-/- and C57 mice in response to various treatments. Data represent mean ± SD for 5 mice fed the AIN-76A diet, and housed in specific pathogen free (SPF) or conventional (C) conditions, or housed in conventional conditions and orally inoculated with 12 Enterococcus strains (EF), with conventional intestinal flora from normal C57 mice (CIF), or a 1:1 combination (EF·CIF). Asterisks denote Il10-/- mice with values significantly different (*P < 0.05, ***P < 0.01) from C57 mice in the same treatment group, while the "†" symbol denotes Il10-/- mice with values significantly (P < 0.05) different from Il10-/- mice under SPF conditions for the same histological parameter. One Il10-/- mouse in the CIF treatment group died during the course of the trial. An autopsy was unable to unequivocally establish the cause of death.

Plasma Serum Amyloid A (SAA)

Overall, plasma SAA levels were higher in Il10-/- mice compared with C57 mice (Il10-/- 87 ± 100; C57 15 ± 42 μg/ml, P < 0.001), and this difference between mouse strains was also significant in each of the treatment groups (P < 0.05). The interaction between strain and group was not significant (P = 0.53), and no overall significant within-strain effect of treatment group was detected (P = 0.15). While the lack of a significant between-strain difference in mice inoculated with EF·CIF seems to be due to the presence of two outliers (one in the C57 group and one in Il10-/- group), there is no data from any other observations which would justify the removal of these values. These data are summarized in Table 1.

Table 1. Plasma serum amyloid A and cytokine data from Il10-/- and C57 mice

A small number of SAA values were negative (i.e. were at the limit of detection). Because these data were log transformed (which resulted in homogeneous variance), a constant was added to all values before transformation. This constant was 6.8, being the absolute value of the minimum recorded value (-6.3 μg/ml) plus 0.5.

Plasma Cytokines

For five of the cytokines measured (namely IL-10, IFNγ, TNFα, granulocyte monocyte colony-stimulating factor (GM-CSF), and IL-17), there were four or fewer non-zero values across the dataset, thus no statistical analysis was performed. Data from analysis of the remaining five plasma cytokines are summarized in Table 1.

Microarrays

We have previously described microarray data from these experiments comparing Il10-/- mice with C57 mice in each of the five treatment groups [34]. The total number of probes differentially expressed in any of the within-treatment (Il10-/- vs. C57) comparisons was 6,521, and this set was used for unsupervised hierarchical clustering. This analysis grouped the array slides into two main clusters: inoculated Il10-/- in one cluster, and C57 mice and non-inoculated Il10-/- mice in the other cluster (Figure 2).

thumbnailFigure 2. Hierarchical clustering of 20 microarray slides (representing RNA pooled from either two or three mice per slide) and 6,521 probe sets. A heat map of 6,521 transcripts and 2 dendograms that group gene probes (left) and microarray slides (top) together is shown. This represents all differentially expressed probes for the comparison Il10-/- vs. C57 mice for any of the five treatment groups (SPF, C, EF, CIF, EF·CIF). Each line is a probe, and each column is an array slide. Expression signal intensities are shown in red and blue, indicating high and low expression, respectively. Arrays for Il10-/- mice are represented in yellow and C57 mice in blue in the bar above the heat map.

We therefore decided to focus on genes differentially expressed when comparing inoculated Il10-/- mice (i.e. CIF, EF, and EF·CIF) with Il10-/- mice not receiving an inoculation, in order to better understand interactions or relationships between genes which may be leading to the increased inflammation observed in the inoculated mice (first hypothesis), and to compare these changes with those previously observed in human IBD (second hypothesis).

There were no significant differences in expression (q < 0.05) of any probes between Il10-/- mice in the C or SPF groups at 12 weeks of age, thus subsequent analysis of gene expression changes in inoculated Il10-/- mice was with those in the C group. Il10-/- in the CIF group had approximately 1,500 gene probes differentially expressed compared to those in group C, while Il10-/- mice inoculated with EF or EF·CIF had over 3,500 differentially expressed probes when compared with Il10-/- mice in group C (Table 2).

Table 2. Number of differentially expressed genes in colon tissue of Il10-/- mice compared to C57 mice

Comparison with human IBD

Thirty two genes are known to undergo changes in transcriptional activity in IBD [35]. In order to assess the relevance of the inoculation protocols used in our study to that of human disease, we compared changes in expression of these thirty-two genes across each of the treatment groups. These results, summarised in Table 3, show that the two inoculations containing EF strains were most effective in mimicking human IBD on the basis of gene expression changes, in particular with respect to changes in gene expression levels of cytokines and chemokines and their receptors.

Table 3. Differentially expressed genes in colon tissue of Il10-/- mice compared to C57 mice, and in Il10-/- mice orally inoculated with intestinal bacteria, in common with those observed in human studies

Ingenuity pathways analysis

Table 4 shows the top functions in each treatment group and the number of differentially expressed genes associated with each function. For all groups analyzed using IPA, the most significant biological functions were grouped into three categories: 1) Diseases and Disorders; 2) Molecular and Cellular Functions; and 3) Physiological System Development and Function (Table 4).

Table 4. Key functions associated with inflammation in Il10-/- mice identified using Ingenuity Pathways Analysis

Table 5 lists the Canonical Pathways showing significant differences between Il10-/- and C57 mice within each treatment group, calculated by IPA. Calculation was either according to ratio (the number of genes from the data set that map to the canonical pathway in question divided by the total number of genes that map to the same canonical pathway) or significance (Fischer's exact test was used to calculate a P-value determining the probability that the association between the genes in the dataset and the canonical pathway is explained by chance alone). A full list of the genes within the ten most significant Canonical Pathways is shown in Additional file 1 Table S1.

Table 5. Key canonical pathways associated with inflammation in Il10-/- mice

Additional file 1. Pathways with genes differentially expressed in Il10-/- mice in response to the 1) SPF, 2) C, 3) EF, 4) CIF and 5) EF·CIF treatments when compared with similarly-inoculated C57 mice. Fold-change and P-values for genes in the top 10 most significant canonical pathways as identified by Ingenuity Pathways Analysis (using Fischer's exact test, as shown in Table 5 within the main text). For those comparisons where there was a significant difference between Il10-/- and C57 mice, the text is shown in italics.

Format: DOC Size: 382KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

In order to examine those molecules directly interacting with IL-10 (and therefore gain insight into the effect of the lack of this protein in the Il10-/- mice), the Neighbourhood Explorer function of IPA was used. The Network Neighborhood of IL-10 consists of 405 molecules which have been identified as directly interacting with IL-10. The number of differentially expressed genes (Il10-/- vs. C57) within this network as a result of the treatments was as follows: SPF (7); C (8); EF (86); CIF (79); EF·CIF (82). Two genes of interest in this comparison were regenerating islet-derived 3 beta (Reg3b) (Table 3) and the polymeric immunoglobulin receptor (pIgR), as these were the only genes which showed significantly higher expression in Il10-/- mice in the SPF, C and CIF groups (compared to C57 mice), but for which there was no difference in the EF and EF·CIF groups. The expression of Reg3b mRNA was between 4- and 8-fold higher in the colon of Il10-/- mice in the SPF and C and CIF groups when compared to C57 mice in the same groups, but there was no differential expression of this molecule in the EF and EF·CIF groups (Table 3). In the case of the related Reg3g gene, while its expression was higher in Il10-/- mice in all experimental groups, the difference was greatest in the C and CIF groups (~9-fold), less in the SPF and EF·CIF groups (~5-fold) and lowest in the EF group (3-fold, Table 3). In the case of pIgR, expression was higher in the SPF (2.0-fold), C (2.2-fold) and CIF (1.7-fold) groups when comparing Il10-/- mice with C57 mice, but there was no difference in the EF or EF·CIF groups.

In order to better understand the effect of the EF·CIF inoculation in Il10-/- mice compared with the same inoculation in C57 mice, genes from the top five canonical pathways on the basis of IPA P-value (Oxidative Phosphorylation, Antigen Presentation Pathway, IL-10 Signaling, Interferon Signaling, LPS/IL-1 mediated inhibition of RXR function) were combined with those from the top five on the basis of calculated ratio (Antigen Presentation Pathway, Interferon Signaling, Oxidative Phosphorylation, IL-10 Signaling, Circadian Rhythm Signaling). A network was built by connecting these molecules according to their interactions, as determined by the IPA Knowledge base (Figure 3). A similar analysis was performed for EF·CIF-inoculated Il10-/- mice compared to Il10-/- mice which were not inoculated (Il10-/- (EF·CIF vs. C), Figure 4).

thumbnailFigure 3. Generation of a biological network of genes of the most significant Canonical Pathways for the EF·CIF (Il10-/- vs. C57) comparison. The network was generated by IPA using all molecules from significantly affected Canonical Pathways (Fatty Acid Metabolism; LPS/IL-1 mediated inhibition of RXR function; Tryptophan Metabolism; β-Alanine Metabolism; Valine, Leucine and Isoleucine Degradation; Antigen Presentation Pathway; Interferon Signaling; IL-10 Signaling; Fatty Acid Elongation in Mitochondria). Connections were applied based on known interactions between these genes within the Ingenuity Pathways Knowledge Base. Central genes and their direct interactions were identified and were supported by published information. *Genes that are detected 2 or more times on the array.

Genes or gene products are represented as nodes, and the biological relationship between two nodes is represented as a line (i.e. an edge). All edges are supported by at least 1 reference from the literature. Red and green colored nodes indicate genes with up- and down-regulated expression, respectively. The intensity of the colors specifies the degree of up- or down-regulation. Greater intensity represents a higher level of differential expression. Nodes and edges are displayed with various shapes and labels that present the functional class of genes and the nature of the relationships between the nodes, as shown in legend below the figure.

thumbnailFigure 4. Generation of a biological network of genes of the most significant Canonical Pathways for the Il10-/- (EF·CIF vs. C) comparison. The network was generated by IPA using all molecules from significantly affected Canonical Pathways (Oxidative Phosphorylation; Antigen Presentation Pathway; IL-10 Signaling; Interferon Signaling; LPS/IL-1 mediated inhibition of RXR function; Circadian Rhythm Signaling). Connections were applied based on known interactions between these genes within the Ingenuity Pathways Knowledge Base. Central genes and their direct interactions were identified and were supported by published information. *Genes that are detected 2 or more times on the array. Genes or gene products are represented as nodes, and the biological relationship between two nodes is represented as a line (i.e. an edge). All edges are supported by at least 1 reference from the literature. Red and green coloured nodes indicate genes with up- and down-regulated expression, respectively. The intensity of the colours specifies the degree of up- or down-regulation. Greater intensity represents a higher level of differential expression. Nodes and edges are displayed with various shapes and labels that present the functional class of genes and the nature of the relationships between the nodes, as shown in the legend to Figure 3.

Key molecules in the Il10-/- (EF·CIF vs. C) comparison were: IFNγ, TNF, STAT1, IL-1β and suppressor of cytokine signaling 1 (SOCS1) and SOCS3. These molecules also appear to be important when comparing EF·CIF-inoculated Il10-/- with C57 receiving the same inoculation (Figure3), as do peroxisome proliferator-activated receptor alpha (PPARα) and TLR4. Conversely, the expression levels of the pro-inflammatory IL-18 gene were reduced in the Il10-/- (EF·CIF vs. C) comparison, and in Figure 4, IL-18 is shown to interact directly with five other molecules in the network; in contrast, this molecule does not feature in the EF·CIF (Il10-/- vs. C57) comparison (Figure 3).

Quantitative real-time polymerase chain reaction (qRT-PCR)

For the three xenobiotic metabolism genes selected for qRT-PCR validation (Cyp2c40, Ces2 and Sult1a1) the reduced expression in inoculated Il10-/- mice observed in the microarray analysis was confirmed (both magnitude and direction) by qRT-PCR. In the case of the immune response genes Ifng and Ncf4, the direction of the gene expression change (an increase in both cases) was confirmed, but the magnitude of the change for Ifng was higher in the PCR analysis compared to the microarray results (Figure 5).

thumbnailFigure 5. Real time quantitative PCR validation of gene expression results from microarray analysis. Real time quantitative PCR was performed to determine the relative expression of Ifng (a), Sult1a1 (b), Cyp2c40 (c), Ces2 (d) and Ncf4 (e), five genes which showed differential expression (as determined using microarray analysis) in Il10-/- mice as a result of bacterial inoculation. Results for the differentially expressed genes were normalized against the geometric mean of CANX and HPRT. An asterisk indicates a significant difference in gene expression for the comparison of interest.

Discussion

Our study clearly shows that oral inoculation of 5 week old Il10-/- mice (C57 background) with a bacterial preparation (12 Enterococcus strains combined with complex intestinal flora from C57 mice raised under conventional conditions) results in consistent and increased colon inflammation at 12 weeks of age, compared with Il10-/- mice which were not inoculated (C group). These tissue changes were similar to those seen in human IBD, and were accompanied by significant changes in expression of genes involved in 'immune response', 'inflammatory diseases' and 'antigen presentation'.

Histological analysis

Several studies have reported the development of intestinal inflammation in Il10-/- mice when raised under either SPF or conventional conditions in the absence of bacterial inoculation [16,19,36]. It has been also reported that the induction of IBD-like colitis (and therefore relevance as a model) in the Il10-/- mouse is dependent on the background mouse strain (with C3H/HeJBir or BALB/C Il10-/- mouse being more susceptible, and C57BL/6 Il10-/- mice (the strain used in this study) being more resistant to the development of intestinal inflammation [1,19]) and the commensal bacteria to which the mice are exposed [23]. We observed no significant difference in intestinal inflammation between Il10-/- mice (C57 background) and C57 control mice when kept in either SPF or conventional conditions. In addition, comparison of colon gene expression levels between Il10-/- and C57 mice in these two treatment groups showed relatively few changes in the current study. The Il10-/- mice in these studies were brought in from The Jackson Laboratory (where they had been raised under SPF conditions), were free of all tested viruses, bacteria and mycoplasma and other opportunistic organisms, and were maintained under quarantine during the trial. It therefore seems likely that the failure of these mice to develop inflammation under conventional conditions is due at least in part to insufficient exposure to commensal intestinal bacteria. Our results clearly demonstrate that, in the absence of such bacterial exposure and in the conventional conditions used in our study, the Il10-/- mouse on a C57 background at 12 weeks of age is not a suitable model for human IBD.

Relevance of bacterial inoculation

All of the inoculation protocols we tested resulted in increased intestinal (particularly colon) inflammation in Il10-/- mice, both when compared with similarly-inoculated C57 mice and when compared with Il10-/- mice which received no inoculation (SPF or C groups). Although the CIF inoculation showed some evidence of colon inflammation, based on HIS the two inoculation preparations containing Enterococcus species (EF and EF·CIF) were more effective in the induction of inflammation. In particular, the EF·CIF inoculation resulted in a significant increase in overall HIS, with the greatest effect observed in the colon. The nature of the inflammation within the colon was transmural and discontinuous, affecting different layers of the intestinal wall, similar to human CD [2]. Our findings for the Il10-/- mice inoculated with Enterococcus species are in agreement with those of other studies where these bacterial species have been shown to induce inflammation in Il10-/- mice [14,23,37]. Our results support observations that inoculation of Il10-/- mice with Enterococcus species gives rise to characteristics of chronic inflammatory diseases, such as a dysregulated immune response acting in combination with inflammatory mechanisms that may lead to tissue damage.

This study used Enterococcus strains isolated from calves and poultry, and there were several reasons for this. First, these strains were readily available, and there is clear evidence in the literature that Enterococcus species induce inflammation in this mouse model [14,23,24]; establishing a relevant mouse model of IBD within our laboratory was the key goal of these studies. Furthermore, enterococci are ubiquitous as a gastrointestinal bacterium in warm-blooded animals, and using those strains isolated from calves ensured that we included enterococci from a mammalian host. Lastly, all strains were fecal isolates, and therefore derived from an appropriate environment.

We do not have specific information on the relationship of these isolates to those derived from mice. It is possible that some virulence factors in poultry may be different from calves, which may be different from mice. However, we have shown that the enterococci, either alone or in combination with CIF derived from control mice, trigger inflammation regardless of the fact that they are from a different host. Thus the original host may be immaterial, and there is some factor (or factors) common to Enterococcus strains which triggers inflammation in a genetically predisposed host. Other studies have used human oral Enterococcus isolates in Il10-/- mice [23,37], which also suggests that there is no host specificity.

Finally, while the ability of the bacteria within the various inocula to colonize, and to persist, within the murine gut was not established in this study, neither colonisation nor persistence is necessarily a pre-requisite if the appropriate inflammatory response is observed even from transiently dosed enterococci. Thus, although the long-term fate of the inoculated bacteria within the murine gut was not established, the consequences of the inoculation in terms of intestinal inflammation, which was the key outcome of the study, have been clearly established.

Gene expression profile in inoculated Il10-/- mice

Each of the inoculation procedures induced changes in the expression levels of a large number of genes, both when comparing Il10-/- mice with similarly inoculated C57 mice and when comparing the inoculated Il10-/- mice with those raised under conventional conditions.

It is apparent from the unsupervised hierarchical clustering analysis of all differentially expressed genes (Figure2) that the various groups of mice fall into two main clusters, with C57 mice and Il10-/- mice in the SPF and C groups (i.e. those which had not received any inoculation) in one, and inoculated Il10-/- mice in the other. The Il10-/- mice that did not receive any inoculation were clustered more closely with the C57 mice than with those Il10-/- mice which did receive a bacterial inoculation. This confirms the histological observation that, in the absence of bacterial inoculation, the Il10-/- mouse on a C57 background at 12 weeks of age is not a suitable model for human IBD in the conditions used for our study. Furthermore, bacterial inoculation led to many more changes in the expression of genes that would normally interact with IL-10.

Of the three inoculation protocols, CIF appears to be the least effective in triggering intestinal inflammation, showing the lowest number of relevant differentially expressed genes. This is in agreement with published studies showing Enterococcus species to be particularly effective in triggering inflammation in the Il10-/- mouse model.

Expression levels of mRNA encoding pro-inflammatory cytokines IFN-γ, TNF-α, and IL-1 in colon tissue were higher in inoculated Il10-/- mice compared with C57 mice, regardless of the inoculation, although mRNA expression of the pleiotropic pro-inflammatory cytokine IL-18 was lower in the colon of EF·CIF-inoculated mice. The pro-inflammatory IL-6 was detected in plasma samples at higher levels in both EF and EF·CIF-inoculated Il10-/- mice compared with similarly inoculated C57 mice and (in the case of EF) compared with non-inoculated Il10-/- mice, although there was no difference in the colon mRNA levels of IL-6 for these two inoculations as determined by microarray analysis. IL-6, a pleiotropic cytokine, has been shown to play a crucial role in the chronic inflammatory process in IBD [38]. This cytokine is produced by macrophages, lymphocytes and intestinal epithelial cells in response to intestinal microflora [39] and stimulates T-cell expansion via anti-apoptopic signaling [38]. T-cell accumulation increases IL-6 levels and thereby a vicious cycle is induced, leading to chronic inflammation. While there were also some differences in the concentrations of other plasma cytokines both within and across strains (e.g. IL-1α higher in EF-inoculated Il10-/- mice, IL-4 lower in EF- and EF·CIF-inoculated mice, both Il10-/- and C57), the levels of Th1/Th2 cytokines measured in plasma samples were generally low or undetectable. In future studies, inflammatory markers such as these Th1/Th2 cytokines could be measured at the site of inflammation to confirm observed gene expression changes, and further verify the relevance of this model of IBD. Our findings are in agreement with observations that the pathogenesis of IBD in humans is mediated in part through an imbalance of pro- and anti-inflammatory cytokines [40]. Our inoculated Il10-/- mice showed a cytokine profile (increased IFNγ, TNFα, IL-1, IL-6 gene expression levels in colon tissues) that is characteristic of the CD inflammatory response. CD is thought to be dominated by a T-helper (Th)-1 response [41], with increased production of IFNγ [42] and IL-2 [43]. TNF-α has been identified as a key pathogenic cytokine for immune-mediated inflammatory diseases such as CD [44].

As shown in Table 3, this is further supported by the fact that inoculation of Il10-/- mice with Enterococcus strains led to differential expression (compared with Il10-/- mice that did not receive an inoculation) of many of the thirty two genes (19 for EF·CIF, 18 for EF) previously identified in the literature as being relevant to IBD in human studies [35], including the ATP-binding cassette (ABC) transporter ABCB1A. Our study is also in agreement with other studies using Il10-/- mice, in which MHC class II mRNA expression levels in cecal samples was increased, and cytochrome-P450 expression was decreased [45]. We observed up-regulation of MHC class II molecule mRNA (antigen presentation pathway, Table 5) and down-regulation of mRNA expression of several cytochrome P450 genes (Cyp4b1, Cyp2c9, Cyp2c18), as well as a sulfotransferase (Sult1c1), a carboxylesterase (Ces1) and the ATP-binding cassette (ABC) transporters Abcb1b, Abcb4, Abcb6 and Abcb9 in the colon of EF·CIF-inoculated Il10-/- mice. Other xenobiotic metabolism genes such as aldo-keto reductases (Akr1c6), a flavin containing monooxygenase (Fmo5), and three of the genes selected for RT-PCR validation (Ces2, Sult1a1 and Cyp2c40) were also down-regulated. A decrease in the expression of detoxification enzymes and transporters may be an effect secondary to the occurrence of inflammation; both inflammation and infection are known to down-regulate the activity and expression levels of drug metabolizing enzymes and transporters [46].

Our results clearly show that expression of mRNA encoding the CD14 protein, a receptor expressed on the surface of macrophages, monocytes and neutrophils which is essential for LPS-dependent signal transduction via TLR4 [47], was up-regulated in Il10-/- mice as a result of bacterial inoculation, particularly associated with EF inoculation (both EF and EF·CIF treatments). This could be a mechanism by which inflammation is initiated in our mouse model of IBD. Quantitative trait loci and microarray analyses of intestinal tissues of two Il10-/- strains with divergent susceptibility to developing inflammation, C3H/HeJBir Il10-/- mice (most susceptible) and C57 Il10-/- mice (less susceptible), identified Cd14 as a candidate gene associated with this difference in susceptibility [1]. In addition, a polymorphism in the human CD14 promoter (resulting in increased expression of the CD14 receptor) was shown to be associated with both CD [48] and UC [49]. Both of these studies refer to genetic susceptibility to inflammation associated with CD14, while our study shows a change in the expression of Cd14 mRNA in inflamed mice, possibly caused by increased cell influx, thus they are not directly comparable. All of these studies do, however, highlight the potential importance of CD14 in the initiation of intestinal inflammation.

In the current study, Reg3b mRNA levels were higher in the colon of Il10-/- mice in the SPF and C (non-inflamed) and CIF (moderately inflamed) groups when compared to C57 mice in the same groups, but not in the EF and EF·CIF groups, which showed more severe signs of inflammation (Table 3); a similar pattern was observed for the Reg3g gene. Reg3b and Reg3g are expressed in mouse intestine [50], and encode murine orthologues of human pancreatitis-associated proteins (PAP) which may be involved in the innate immune response to bacterial colonisation of the intestinal tract [51] and can inhibit the inflammatory response by blocking NFκB activation [7].

The Reg3b and Reg3g genes are induced after development of intestinal inflammation in germ-free severe combined immunodeficiency (SCID) mice colonised with commensal bacteria [52], and bacterial-epithelial contact may drive Reg3g expression as a mechanism to limit microbial penetration and maintain mucosal surface integrity [53]. In our study, the fact that expression levels of both Reg3b and Reg3g were lower in Il10-/- mice inoculated with Enterococcus strains relative those in the C and CIF groups may reflect an inability of these mice to appropriately respond to this inoculation, and a failure to effectively deal with the additional enterococci introduced via inoculation [54]. Alternatively, lower overall expression levels of these two Reg3 genes in response to Enterococcus inoculation may result in insufficient suppression of the immune response, potentially a crucial early step in the initiation and development of intestinal inflammation in these mice. This is in agreement with studies in the Il2-/- mouse model, where increased expression levels of Reg3b and Reg3g genes might be associated with prevention of colitis triggered by colonization with commensal E. coli [51]. While the levels of Reg3b may be sufficient to explain the lack of phenotype in the SPF/C groups (and to a lesser extent the CIF group), the absence of IL10 (which has been shown to regulate Reg3b in rat cells in vitro [55]) in the Il10-/- mice may mean that insufficient Reg3b is expressed to suppress inflammation induced by enterococci. It is not clear from the literature what constitutes a 'suppressive' amount of Reg3b in a colitic phenotype, and further studies would be required to establish this, and to clarify functional roles of either Reg3b or Reg3g in Enterococcus-induced inflammation in this model.

The decreased expression of polymeric immunoglobulin receptor (Pigr) mRNA in Il10-/- mice in the EF and EF·CIF groups showed a similar pattern of expression to Reg3b in our study. The protein product of the Pigr gene is a glycoprotein expressed on the basolateral surface of secretory columnar and crypt epithelial cells, and is responsible for active transport of secretory antibodies such as IgA across the secretory epithelium that lines the mucosal surfaces. Studies with Pigr knockout (Pigr-/-) mice have produced strong evidence that innate secretory antibodies protect against invasion by pathogenic bacteria via immune exclusion; these mice are unable to bind and actively transport dimeric IgA to the mucosa and thus are more susceptible to infection with virulent strains such as Salmonella [56]. Furthermore, Pigr-/- mice exhibit profound immunopathological changes and clinical disease in response to induction of colitis with dextran sulfate sodium, suggesting that Pigr may play an important role in modulating inflammatory responses in the mucosa during active colitis [57]. The low Pigr gene expression in our Il10-/- mice receiving an inoculation with Enterococcus strains suggests a reduced capacity to prevent induction of a systemic immune response, a possible mechanism by which otherwise harmless commensal bacterial strains may have initiated the chronic inflammation observed in these two treatment groups.

In the current study, antigen presentation genes such as major histocompatibility complex (MHC) class II molecules (human leukocyte antigen (HLA) family members - H2-Ab1, H2-Aa, H2-Eb1 and Cd74) showed higher expression levels in colon in Il10-/- mice in all treatment groups (SPF, C, EF, CIF, EF·CIF) compared with C57 mice in the same treatment groups. The up-regulation of antigen presentation genes in our study is consistent with increased expression of HLA class II molecules that typically occurs in IBD in humans [58]. In addition, MHC class II antigen expression has been observed in the colon epithelial cells of Il10-/- mice as early as three weeks of age, with higher expression levels observed in mice at three to six months of age [19]. It has also been reported that IFNγ appears to play a role in development of intestinal colitis, but that neither IFNγ nor MHC class II expression are required for sustaining disease once it has become established [19]. In our study, in contrast to antigen presentation gene response, there was no difference in expression levels of colon Tnf or Ifng genes between Il10-/- and C57 mice in the C or SPF groups, but both of these genes were up-regulated in Il10-/- mice as a result of bacterial inoculation. This occurs when compared with similarly-inoculated C57 mice, and with Il10-/- mice in the C group. These results are further evidence that IFNγ does play a role in the development of intestinal colitis, whereas MHC class II expression alone is not sufficient to trigger an inflammatory response in this model.

Conclusions

Overall, our findings indicate that inoculation of Il10-/- mice with solutions containing intestinal bacteria increases colon inflammation, and that the use of Enterococcus strains in particular results in a more appropriate model of IBD compared with non-inoculated Il10-/- mice. High density oligonucleotide microarrays have identified gene expression changes in Il10-/- colonic tissue in response to bacterial inoculation that are consistent with the current knowledge of mechanisms responsible for human IBD. In addition, there is preliminary evidence for the inflammatory response in this model being initiated by a failure of the normal mechanisms which recognize commensal bacterial, through molecules such as pIgR and REG3A/REG3G. Our data suggest that in particular the EF·CIF inoculation, which results in exposure to a highly complex bacterial environment, gives an appropriate and relevant mouse model of human IBD in which a variety of food components could be tested to establish potential ameliorating effects, and to understand the mechanisms by which these effects may occur.

Methods

Animals and diet

This study was reviewed and approved by the AgResearch Ruakura Animal Ethics Committee in Hamilton, New Zealand according to the Animal Protection Act (1960) and Animal Protection Regulations (1987) and amendments. Twenty five male Il10-/- (C57 background, formal designation B6.129P2-Il10tm1Cgn/J) mice and twenty five male C57 control mice were received from The Jackson Laboratory (Bar Harbor, ME, USA) at approximately 5 weeks of age. Prior to this time, Il10-/- mice were raised under SPF conditions, and were therefore free of all tested viruses, bacteria and mycoplasma and other opportunistic organisms (including Helicobactor, Pasteurella, and Pseudomonas). Mice were maintained under quarantine throughout the trial, and were housed either in pairs or groups of three (5 mice per treatment) in shoebox-style cages containing untreated wood shavings and a plastic tube for environmental enrichment. The animal room was maintained at a temperature of ~22°C and humidity of ~50% with a 12-hour light/dark cycle. All mice had ad libitum access to water, which was refreshed twice a week. An AIN-76A powdered diet, prepared as previously described [59], was supplied twice a week, with sufficient provided to meet the daily intake of Il10-/- mice, as determined in a previous feeding trial (data not shown). The diet for all groups was sterilized by gamma irradiation (25 kGy, Schering-Plough, Wellington, New Zealand) to a level required for SPF conditions, to minimize the possibility of bacterial introduction to the SPF group of animals. All mice were weighed twice weekly and carefully monitored for disease symptoms (weight loss, soft faeces and inactivity).

Experimental design

Both Il10-/- and C57 mice were randomly divided into five treatment groups with five animals per group (Figure 6). One group of mice was housed in SPF conditions (isolator cages supplied with high efficiency particulate air (HEPA)- filtered air (Tecniplast SpA, Buguggiate, Italy)); a second group was maintained under conventional conditions (C), while the remaining groups were kept in conventional conditions and orally inoculated (200 μl) with solutions containing either the 12 strains of Enterococcus listed in Table 1 (EF; 1.2 × 108 colony forming units (CFU)), CIF derived from healthy age-matched C57BL/6 mice raised under conventional conditions, or a combination of the two (EF·CIF; 6.0 × 107 CFU from the EF inoculum), as described below. The CIF inoculation was included to better mimic the complete microbiota associated with the mouse gastrointestinal tract. Microbial ecology of the gastrointestinal tract is a complex interaction between microorganisms, the host and food components with potentially as few as 50% of microbes being able to be cultured in the laboratory, many requiring unknown or undefined conditions for growth [60]. It is likely that many of the uncultivable micro-organisms may have a role in gastrointestinal tract homeostasis in addition to eliciting inflammation.

thumbnailFigure 6. Overall study design. Both Il10-/- and C57 mice were fed AIN-76A diet from the time of arrival (34 ± 3 days of age), and were randomly divided into five treatment groups with five animals per group. One group of mice was housed in SPF conditions, a second group was maintained under conventional conditions (C), while the remaining groups were kept in conventional conditions and orally inoculated (200 μl) with solutions containing either the 12 strains of Enterococcus listed in Table 1 (EF), conventional intestinal flora (CIF) derived from healthy age-matched C57 mice raised under conventional conditions, or a combination of the two (EF·CIF). Inoculation was performed at approximately 5.5 weeks of age, and tissue sampling at 12 weeks of age.

Preparation of bacterial solutions for oral inoculation

Enterococcus strains were obtained from calf or poultry faeces using sterile Amies transfer swabs (Biolab Limited, Auckland, New Zealand) and inoculating directly onto Slanetz and Bartley agar (Oxoid, Hampshire, UK). Plates were incubated at 42°C for 48 hours with red or purple colonies subcultured and incubated using the same culture conditions. Putative enterococcus strains were stored as glycerol stocks at -85°C in brain heart infusion broth supplemented with 30% (v/v) glycerol. Speciation of the E. faecium and E. faecalis strains was performed using PCR primers specific for amplification of the D-ala D-ala ligase gene - ddlEF (5'-TAGAGACATTGAATATGCC-3') and ddlER (5'-AATCGCACCGGCTCAATC-3') - which were modified from a previous study [61].

Seventy-two hours prior to inoculation, colonies of the E. faecalis and E. faecium strains (Table 6) were sub-cultured from slope tubes onto fresh Slanetz & Bartley agar, and incubated at 42°C for 48 h. A single colony from each culture was subsequently transferred to 5 ml of Todd Hewitt Broth (Oxoid, Hampshire, UK), and incubated at 37°C for 24 h. Each Enterococcus culture was centrifuged to pellet the bacterial cells (3000 g, 10 min, 4°C), which were then re-suspended in 5 ml sterile PBS (pH 7.4) and pooled (EF inoculum).

Table 6. Strains of Enterococcus faecalis and E. faecium used in solutions for oral inoculation

The healthy, age-matched C57BL/6 mice from which flora were being collected were euthanized by CO2 asphyxiation and cervical dislocation and the gastrointestinal tract (from stomach to just below the caecum) removed. Digesta were collected from the intestine and caecum by gently washing with sterile PBS, pH 7.4, then suspended in a total of 30 ml PBS. After mixing by gentle inversion and a settling period of approximately 5 min, the suspension was collected (CIF inoculum). Inoculation solutions EF and CIF were mixed in a 1:1 (v/v) solution to obtain the EF·CIF inoculum.

Sample collection

At 12 weeks of age, samples were collected from all mice. To minimize time variation between the last food intake and sampling, mice were fasted overnight on the night before sampling. On the morning of sampling, food was returned for two hours, followed by a further two hour fast immediately prior to tissue sampling [62].

Mice from all treatments were euthanized using CO2 asphyxiation followed by cervical dislocation. Blood was sampled via cardiac puncture (0.5 to 1 ml), cells pelleted and the plasma snap-frozen and stored at -85°C for subsequent cytokine and SAA analysis. The intestine was quickly removed, cut open lengthwise and flushed with 0.9% sodium chloride to remove any traces of digesta, then carefully laid out on an ice-cold stainless steel tray. Sections of each intestinal region (duodenum, jejunum, ileum and colon) were frozen in liquid nitrogen before storage at -85°C, for gene profiling, and a sub-sample from each intestinal region was fixed in 10% phosphate buffered formalin and stored at room temperature until histological evaluation of inflammation.

Histology

Histological examination was performed as previously described [34,59]. Briefly, formalin-fixed samples were processed, sectioned, stained with haematoxylin and eosin and evaluated for inflammation under a light microscope, using a modification of a previously described scoring system [15,63]. A histological injury score (HIS) was assigned based on the presence of inflammatory lesions, tissue destruction and tissue repair.

Plasma SAA

Inflammation was also determined by analysis of SAA levels in plasma using a murine-specific Phase SAA ELISA kit, according to the manufacturer's protocol (Tridelta Development Ltd., Maynooth, County Kildare, Ireland). Pipetting of standards and samples was performed using an epMotion 5070 Liquid Handling Workstation (Eppendorf South Pacific Pty. Ltd., NSW, Australia). Briefly, standards or samples plus biotinylated monoclonal SAA antibody were incubated in microtitre plate wells pre-coated with capture monoclonal SAA antibody. In one step, SAA (in the standard or sample) was captured and labelled in a sandwich format. After washing to remove unbound material, wells were incubated with streptavidin-horseradish peroxidase prior to the addition of enzyme substrate (3,3',5,5'-tetramethylbenzidine, TMB). The reaction was stopped with the addition of 2 M sulphuric acid. Optical density in the wells was measured at 450 nm (630 nm as reference) using an automated plate reader (Versamax, Molecular Devices, CA, USA).

Plasma cytokines

Plasma samples were analyzed for the presence of Th1/Th2-related cytokines (IL-1α, IL-2, IL-5, IL-6, IL-10, IFNγ, TNFα, granulocyte monocyte colony-stimulating factor (GM-CSF), IL-4 and IL-17) using a FlowCytomix Multiplex (Th1/Th2 10plex) mouse kit (Bender MedSystems GmbH, Vienna, Austria) according to the manufacturer's instructions. Briefly, microbeads coated with antibodies to the specified cytokines were mixed with plasma samples. Biotinylated secondary antibodies were then added, and the mixture incubated at room temperature for 2 h with constant shaking (500 rpm). The amount of cytokine that bound to the antibodies was then detected using streptavidin conjugated to phycoerythrin, with end-point measurement in a flow cytometer (FACScan, Becton Dickinson, North Ryde, NSW, Australia). Levels of each cytokine were determined by comparison with a standard curve (concentration range for all cytokines 0-20000 pg/ml).

RNA isolation and synthesis of labeled cRNA

Total RNA from intact colon tissue was isolated by homogenizing the samples in TRIzol (Invitrogen, Auckland, New Zealand) according to the manufacturer's instructions. RNA was quantified with a Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA), and RNA integrity was assessed with an RNA 6000 Nano LabChip kit using the Agilent 2100 Bioanalyser (Agilent Technologies, Palo Alto, CA, USA). Extracted RNA was purified using RNeasy spin columns (QIAGEN, San Diego, CA, USA). Only total RNA with an OD 260/280 ratio > 2.0, a Bioanalyzer 28 s/18 s peak ratio ≥ 1.2 and an RNA integrity number ≥ 7.5 was used for microarray hybridization. An equimolar pool of total RNA extracts from colon tissues of two or three different mice per treatment was made, resulting in two pools per treatment. A reference design was used for microarray hybridization: colonic RNA extracts from all mice were pooled in an equimolar proportion and used as the reference sample.

The Low RNA Input Fluorescent Linear Amplification Kit (Agilent Technologies Inc., Palo Alto, CA, USA) was used to synthesize cDNA and fluorescent cRNA. Labelled cRNA was made on the same day for all pools, including the reference sample. cDNA was synthesized from 500 ng of purified total RNA from each pool according to the manufacturer's protocol. Cyanine 3-cytidine triphosphate (Cy3, PerkinElmer, Waltham, Massachusetts 02451, USA) was used to label sample groups, while the reference RNA was labeled with cyanine 5-CTP (Cy5, PerkinElmer, Waltham, MA 02451, USA). Hybridization was performed according to a reference design without dye swap.

Microarray hybridization and scanning

The in situ hybridization kit-plus (Agilent Technologies Inc., Palo Alto, CA, USA) was used to hybridize cRNA samples to Agilent Technologies Mouse G4121A - 44 k 60 mer oligo arrays. For each experimental group, two pools of cRNA were hybridized, thus in total twenty arrays were used. Cy3-labelled cRNA (0.75 μg) and Cy5-labelled cRNA (0.75 μg) were hybridised onto the microarray according to the manufacturer's protocol, as previously described [34]. Slides were scanned using a GenePix 4200A scanner (Molecular Devices Corporation, Sunnyvale, CA, USA) at a photomultiplier tube (PMT) setting of 450 V. Spot identification and quantification were performed using GenePix 6.0 software (Molecular Devices Corporation, Sunnyvale, CA, USA). The microarray data are available as accession GSE12223 in the Gene Expression Omnibus repository at the National Center for Biotechnology Information http://www.ncbi.nlm.nih.gov/geo/info/linking.html webcite.

Microarray data analysis

Statistical analysis was performed using linear models for microarray analysis (limma) within the Bioconductor framework [64]. Before analysis, poor quality spots were manually flagged and filtered out. Quality of the microarray data was assessed on diagnostic plots (boxplots and density plots) and spatial images generated from the raw (non-processed) data. All twenty arrays passed these strict criteria and were included in the analyses. Intensity ratio values for all microarray spots were normalized using a within-slide global Locally Weighted Scatterplot Smoothing procedure to remove the effect of systematic variation in the microarrays; no background correction was necessary due to homogeneous hybridization. The normalized data from the arrays of each treatment group were averaged. For each comparison, differentially expressed genes were identified using FDR control with a threshold of q < 0.05.

Further analysis of the differentially expressed genes (for example, clustering, self-organizing maps) was performed using Bioconductor and GeneSpring 7.3 (Agilent Technologies Inc., Palo Alto, CA, USA). To gain an overall idea of the pattern of gene expression changes across the various treatments, unsupervised hierarchical cluster analysis was performed on the group of all probes differentially expressed in any of the within-treatment Il10-/- vs. C57 comparisons using Bioconductor.

Network, pathways and functional analyses were generated using Ingenuity Pathways Analysis (IPA, Version 5.0 or Version 6.0, Ingenuity Systems, Redwood City, CA, USA; http://www.ingenuity.com webcite). For analysis using IPA, the full dataset (representing all IDs on the Agilent array) was uploaded for analysis, and the set of differentially expressed genes analyzed using the complete gene-list as a reference data set. For analyses where more than 800 differentially expressed genes were identified as network-eligible, a fold-change cut-off was applied to reduce this number to 800, as is recommended by IPA. The fold-change cut off for each list is reported in the results section. In the case of genes being replicated on the array, the median value of the replicates was used for IPA analysis.

Functional analysis in IPA identified biological functions and/or diseases that were most relevant to the data set. Genes from the dataset that met the FDR cutoff (q < 0.05) and were associated with biological functions and/or diseases in the IPA Knowledge Base were considered for the analysis. Fischer's exact test was used to calculate the probability that each biological function and/or disease assigned to that data set is due to chance alone.

Canonical pathways analysis identified those pathways from the IPA library that were most significant to the data set. Only those genes from the data set meeting the FDR cutoff (q < 0.05) and associated with a canonical pathway in the IPA Knowledge Base were considered for this analysis. The significance of the association between the data set and the canonical pathway was measured in two ways: 1) The ratio of the number of genes from the data set that map to the canonical pathway in question divided by the total number of genes that map to the same canonical pathway; 2) Fischer's exact test was used to calculate the probability that the association between the genes in the dataset and the canonical pathway is explained by chance alone.

Because we have used an Il10-/- mouse model, the "Neighborhood Explorer" feature of IPA was used to investigate genes associated with Il10. The Neighborhood Explorer is a network which displays all molecules that directly interact with a molecule of interest (in this case mRNA encoded by the Il10 gene, which is non-functional in the Il10-/- mice), either regulating or being regulated by that molecule, or physically interacting with it. These interactions are based on information contained within the Ingenuity Pathways Knowledge Base.

Quantitative real-time polymerase chain reaction

Five genes showing differential expression in EF·CIF-inoculated Il10-/- mice (compared with similarly inoculated C57 mice) were selected for validation and their expression levels established using quantitative real-time polymerase chain reaction (qRT-PCR). cDNA was synthesized from the same total RNA samples used for the microarray analysis, using the Transcriptor First Strand cDNA Synthesis Kit (Roche Diagnostics, Mannheim, Germany). Reverse transcription was performed using 1.0 μg of total RNA and oligo-dT primers according to the manufacturer's instructions.

The PCR conditions were: 95°C for 5 min, 35 cycles at 95°C for 15 s, 63°C for 10 s and 72°C for 15 s. Melting curve analyses were performed by increasing the temperature (1°C/s) from 65°C to 95°C, with continuous fluorescence acquisition. With the exception of neutrophil cytosolic factor 4 (Ncf4) [59], primers for the selected genes (cytochrome P450, family 2, subfamily C, polypeptide 40 (Cyp2c40), carboxylesterase 2 (Ces2), sulfotransferase 1A (Sult1a1) and IFNγ (Ifng)) were designed using Primer 3.0 [65], with available public sequences. RefSeq IDs and primer sequences are as previously described [66]. Ifng and Ncf4 are genes associated with immune response, while Cyp2c40, Ces2 and Sult1a1 are xenobiotic metabolism genes. PCR conditions for all primers were optimized and amplicons were sequenced to confirm identity. Specificities of all PCR reactions were verified by melting curves analyses and agarose gel electrophoresis. Data were normalized against two reference genes (hypoxanthine guanine phosphoribosyl transferase 1 (Hprt1) and calnexin (Canx) [59]).

Threshold cycle (Ct) values were obtained in triplicates for each sample on the LightCycler 480 (Roche Diagnostics, Mannheim, Germany) using LightCycler 480 SYBR Green I Master (Roche Diagnostics, Mannheim, Germany) in 20 μl reactions, according to the manufacturer's protocol. Standard curves for all selected genes and reference genes were generated using serial dilutions of pooled cDNAs from all samples. LightCycler 480 Relative Quantification Software was used to calculate mRNA concentrations based on the appropriate standard curves and normalized ratios (target/reference).

Statistical analysis

Statistical analyses of live weight, HIS, qRT-PCR and plasma SAA and cytokines were performed using GenStat (VSN International, Hemel Hempstead, UK; 9th edition, 2006 or 10th edition, 2007). Differences in HIS both between and within mouse strains were analyzed using an ANOVA with pooled variance. As expected, there were a high number of zero histology scores for the C57 mice and so the data for the two strains were analyzed separately. Plasma concentrations of SAA and cytokines were analyzed by ANOVA using strain and treatment as factors. Log transformed values were used for statistical analysis of the histology and SAA data. Due to the skewed nature of the data, and the large proportion of zeros observed for some of the plasma cytokine measurements, a variety of different transformations (such as log, square-root, and non-parametric rank transformations) were required across these variables. Specific transformations used are described in the results section for the cytokine analyses.

A probability value of less than 0.05 was considered significant while a probability value greater than 0.05 but lower than 0.10 was considered a trend. It must be noted that one Il10-/- mouse from the CIF inoculation group died during the experiment (52 days of age) for unknown reasons. This animal was not scored for histological signs of inflammation, and no data from this animal (e.g. live weight) were included in any of the analyses. Unless otherwise stated, data are presented in the text as mean ± SD.

Authors' contributions

MPGB participated in the study design, carried out the animal trials and Pathways Analysis, and drafted the manuscript. WCM participated in the design of the study and discussion of the results. ALC supplied bacterial cultures, participated in this aspect of the study design, and generated relevant discussion. SZ performed all histological evaluations. MD performed all Bioconductor analyses and generated heat maps, BK assisted with interpretation and discussion around microarrays, KN ran RT-PCR reactions, and AJH carried out cytokine analyses. NCR supervised the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors wish to thank Bruce Sinclair, Yvonne Dommels, Emma Bermingham, Janine Cooney and Shalome Bassett for critical evaluation of the manuscript, Kim Oden for assistance with the measurement of plasma SAA levels, Kate Broadley for assistance with the microarray hybridization, Ric Broadhurst for advice and assistance with animal care and tissue sampling, Harold Henderson and Zaneta Park-Ng for statistical help and advice, and Pieter Demmers for assistance with microarray data submission. M. Barnett was funded by Postdoctoral fellowship AGRX0504 from the New Zealand Foundation for Research, Science and Technology (FRST) during the course of these studies. The work described in this paper is part of Nutrigenomics New Zealand http://www.nutrigenomics.org.nz, a collaboration between AgResearch Limited, Plant & Food Research and The University of Auckland, which is largely funded by FRST.

References

  1. de Buhr MF, Mahler M, Geffers R, Hansen W, Westendorf AM, Lauber J, Buer J, Schlegelberger B, Hedrich HJ, Bleich A: Cd14, Gbp1, and Pla2g2a: three major candidate genes for experimental IBD identified by combining QTL and microarray analyses.

    Physiol Genomics 2006, 25(3):426-434. PubMed Abstract | Publisher Full Text OpenURL

  2. Yantiss RK, Odze RD: Diagnostic difficulties in inflammatory bowel disease pathology.

    Histopathology 2006, 48(2):116-132. PubMed Abstract | Publisher Full Text OpenURL

  3. Hampe J, Cuthbert A, Croucher PJ, Mirza MM, Mascheretti S, Fisher S, Frenzel H, King K, Hasselmeyer A, MacPherson AJ, et al.: Association between insertion mutation in NOD2 gene and Crohn's disease in German and British populations.

    Lancet 2001, 357(9272):1925-1928. PubMed Abstract | Publisher Full Text OpenURL

  4. Ogura Y, Bonen DK, Inohara N, Nicolae DL, Chen FF, Ramos R, Britton H, Moran T, Karaliuskas R, Duerr RH, et al.: A frameshift mutation in NOD2 associated with susceptibility to Crohn's disease.

    Nature 2001, 411(6837):603-606. PubMed Abstract | Publisher Full Text OpenURL

  5. Hugot JP, Chamaillard M, Zouali H, Lesage S, Cezard JP, Belaiche J, Almer S, Tysk C, O'Morain CA, Gassull M, et al.: Association of NOD2 leucine-rich repeat variants with susceptibility to Crohn's disease.

    Nature 2001, 411(6837):599-603. PubMed Abstract | Publisher Full Text OpenURL

  6. van Heel DA, Udalova IA, De Silva AP, McGovern DP, Kinouchi Y, Hull J, Lench NJ, Cardon LR, Carey AH, Jewell DP, et al.: Inflammatory bowel disease is associated with a TNF polymorphism that affects an interaction between the OCT1 and NF(-kappa)B transcription factors.

    Hum Mol Genet 2002, 11(11):1281-1289. PubMed Abstract | Publisher Full Text OpenURL

  7. Pierik M, Joossens S, Van Steen K, Van Schuerbeek N, Vlietinck R, Rutgeerts P, Vermeire S: Toll-like receptor-1, -2, and -6 polymorphisms influence disease extension in inflammatory bowel diseases.

    Inflamm Bowel Dis 2006, 12(1):1-8. PubMed Abstract | Publisher Full Text OpenURL

  8. Klein W, Tromm A, Griga T, Fricke H, Folwaczny C, Hocke M, Eitner K, Marx M, Duerig N, Epplen JT: Interleukin-4 and interleukin-4 receptor gene polymorphisms in inflammatory bowel diseases.

    Genes Immun 2001, 2(5):287-289. PubMed Abstract | Publisher Full Text OpenURL

  9. Glas J, Torok HP, Tonenchi L, Kapser J, Schiemann U, Muller-Myhsok B, Folwaczny M, Folwaczny C: Association of polymorphisms in the interleukin-18 gene in patients with Crohn's disease depending on the CARD15/NOD2 genotype.

    Inflamm Bowel Dis 2005, 11(12):1031-1037. PubMed Abstract | Publisher Full Text OpenURL

  10. Potocnik U, Ferkolj I, Glavac D, Dean M: Polymorphisms in multidrug resistance 1 (MDR1) gene are associated with refractory Crohn disease and ulcerative colitis.

    Genes Immun 2004, 5(7):530-539. PubMed Abstract | Publisher Full Text OpenURL

  11. Ardizzone S, Maconi G, Bianchi V, Russo A, Colombo E, Cassinotti A, Penati C, Tenchini ML, Bianchi Porro G: Multidrug resistance 1 gene polymorphism and susceptibility to inflammatory bowel disease.

    Inflamm Bowel Dis 2007, 13(5):516-523. PubMed Abstract | Publisher Full Text OpenURL

  12. Leung E, Hong J, Fraser AG, Merriman TR, Vishnu P, Krissansen GW: Polymorphisms in the organic cation transporter genes SLC22A4 and SLC22A5 and Crohn's disease in a New Zealand Caucasian cohort.

    Immunol Cell Biol 2006, 84(2):233-236. PubMed Abstract | Publisher Full Text OpenURL

  13. Silverberg MS, Duerr RH, Brant SR, Bromfield G, Datta LW, Jani N, Kane SV, Rotter JI, Philip Schumm L, Hillary Steinhart A, et al.: Refined genomic localization and ethnic differences observed for the IBD5 association with Crohn's disease.

    Eur J Hum Genet 2007, 15(3):328-335. PubMed Abstract | Publisher Full Text OpenURL

  14. Balish E, Warner T: Enterococcus faecalis induces inflammatory bowel disease in interleukin-10 knockout mice.

    Am J Pathol 2002, 160(6):2253-2257. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Kennedy RJ, Hoper M, Deodhar K, Erwin PJ, Kirk SJ, Gardiner KR: Interleukin 10-deficient colitis: new similarities to human inflammatory bowel disease.

    Br J Surg 2000, 87(10):1346-1351. PubMed Abstract | Publisher Full Text OpenURL

  16. Kuhn R, Lohler J, Rennick D, Rajewsky K, Muller W: Interleukin-10-deficient mice develop chronic enterocolitis.

    Cell 1993, 75(2):263-274. PubMed Abstract | Publisher Full Text OpenURL

  17. Madsen KL, Malfair D, Gray D, Doyle JS, Jewell LD, Fedorak RN: Interleukin-10 gene-deficient mice develop a primary intestinal permeability defect in response to enteric microflora.

    Inflamm Bowel Dis 1999, 5(4):262-270. PubMed Abstract | Publisher Full Text OpenURL

  18. Aharoni R, Kayhan B, Brenner O, Domev H, Labunskay G, Arnon R: Immunomodulatory therapeutic effect of glatiramer acetate on several murine models of inflammatory bowel disease.

    J Pharmacol Exp Ther 2006, 318(1):68-78. PubMed Abstract | Publisher Full Text OpenURL

  19. Berg DJ, Davidson N, Kuhn R, Muller W, Menon S, Holland G, Thompson-Snipes L, Leach MW, Rennick D: Enterocolitis and colon cancer in interleukin-10-deficient mice are associated with aberrant cytokine production and CD4(+) TH1-like responses.

    J Clin Invest 1996, 98(4):1010-1020. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Cohen SL, Moore AM, Ward WE: Flaxseed oil and inflammation-associated bone abnormalities in interleukin-10 knockout mice.

    J Nutr Biochem 2005, 16(6):368-374. PubMed Abstract | Publisher Full Text OpenURL

  21. McCarthy J, O'Mahony L, O'Callaghan L, Sheil B, Vaughan EE, Fitzsimons N, Fitzgibbon J, O'Sullivan GC, Kiely B, Collins JK, et al.: Double blind, placebo controlled trial of two probiotic strains in interleukin 10 knockout mice and mechanistic link with cytokine balance.

    Gut 2003, 52(7):975-980. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Sydora BC, Tavernini MM, Wessler A, Jewell LD, Fedorak RN: Lack of interleukin-10 leads to intestinal inflammation, independent of the time at which luminal microbial colonization occurs.

    Inflamm Bowel Dis 2003, 9(2):87-97. PubMed Abstract | Publisher Full Text OpenURL

  23. Kim SC, Tonkonogy SL, Albright CA, Tsang J, Balish EJ, Braun J, Huycke MM, Sartor RB: Variable phenotypes of enterocolitis in interleukin 10-deficient mice monoassociated with two different commensal bacteria.

    Gastroenterology 2005, 128(4):891-906. PubMed Abstract | Publisher Full Text OpenURL

  24. Kim SC, Tonkonogy SL, Karrasch T, Jobin C, Sartor RB: Dual-association of gnotobiotic IL-10-/- mice with 2 nonpathogenic commensal bacteria induces aggressive pancolitis.

    Inflamm Bowel Dis 2007, 13(12):1457-1466. PubMed Abstract | Publisher Full Text OpenURL

  25. Jett BD, Huycke MM, Gilmore MS: Virulence of enterococci.

    Clin Microbiol Rev 1994, 7(4):462-478. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Sellon RK, Tonkonogy S, Schultz M, Dieleman LA, Grenther W, Balish E, Rennick DM, Sartor RB: Resident enteric bacteria are necessary for development of spontaneous colitis and immune system activation in interleukin-10-deficient mice.

    Infect Immun 1998, 66(11):5224-5231. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Eaton TJ, Gasson MJ: Molecular screening of Enterococcus virulence determinants and potential for genetic exchange between food and medical isolates.

    Appl Environ Microbiol 2001, 67(4):1628-1635. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Sghir A, Gramet G, Suau A, Rochet V, Pochart P, Dore J: Quantification of bacterial groups within human fecal flora by oligonucleotide probe hybridization.

    Appl Environ Microbiol 2000, 66(5):2263-2266. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Finegold SM, Sutter VL, Mathisen GE: Normal indigenous intestinal flora. In Human intestinal microflora in health and disease. Edited by Hentges DJ. New York: Academic Press; 1983:3-31. OpenURL

  30. Tannock GW, Cook G: Enterococci as members of the intestinal microflora of humans. In The enterococci: pathogenesis, molecular biology, and antibiotic resistance. Edited by Gilmore MS. Washington, DC.: ASM Press; 2002:101-132. OpenURL

  31. Noble CJ: Carriage of group D streptococci in the human bowel.

    J Clin Pathol 1978, 31(12):1182-1186. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Parronchi P, Romagnani P, Annunziato F, Sampognaro S, Becchio A, Giannarini L, Maggi E, Pupilli C, Tonelli F, Romagnani S: Type 1 T-helper cell predominance and interleukin-12 expression in the gut of patients with Crohn's disease.

    Am J Pathol 1997, 150(3):823-832. PubMed Abstract | PubMed Central Full Text OpenURL

  33. Monteleone G, Monteleone I, Fina D, Vavassori P, Del Vecchio Blanco G, Caruso R, Tersigni R, Alessandroni L, Biancone L, Naccari GC, et al.: Interleukin-21 enhances T-helper cell type I signaling and interferon-gamma production in Crohn's disease.

    Gastroenterology 2005, 128(3):687-694. PubMed Abstract | Publisher Full Text OpenURL

  34. Roy NC, Barnett MPG, Knoch B, Dommels YEM, McNabb WC: Nutrigenomics applied to an animal model of Inflammatory Bowel Diseases: transcriptomic analysis of the effects of eicosapentaenoic acid- and arachidonic acid-enriched diets.

    Mutat Res 2007, 622(1-2):103-116. PubMed Abstract | Publisher Full Text OpenURL

  35. Te Velde AA, De Kort F, Sterrenburg E, Pronk I, Ten Kate FJW, Hommes DW, Van Deventer SJH: Comparative analysis of colonic gene expression of three experimental colitis models mimicking inflammatory bowel disease.

    Inflammatory Bowel Diseases 2007, 13(3):325-330. PubMed Abstract | Publisher Full Text OpenURL

  36. Rennick DM, Fort MM: Lessons from genetically engineered animal models. XII. IL-10-deficient (IL-10(-/-) mice and intestinal inflammation.

    Am J Physiol Gastrointest Liver Physiol 2000, 278(6):G829-833. PubMed Abstract | Publisher Full Text OpenURL

  37. Ruiz PA, Shkoda A, Kim SC, Sartor RB, Haller D: IL-10 gene-deficient mice lack TGF-beta/Smad signaling and fail to inhibit proinflammatory gene expression in intestinal epithelial cells after the colonization with colitogenic Enterococcus faecalis.

    J Immunol 2005, 174(5):2990-2999. PubMed Abstract | Publisher Full Text OpenURL

  38. Mudter J, Neurath MF: Il-6 signaling in inflammatory bowel disease: pathophysiological role and clinical relevance.

    Inflamm Bowel Dis 2007, 13(8):1016-1023. PubMed Abstract | Publisher Full Text OpenURL

  39. Guo W, Magnotti LJ, Ding J, Huang Q, Xu D, Deitch EA: Influence of gut microflora on mesenteric lymph cytokine production in rats with hemorrhagic shock.

    J Trauma 2002, 52(6):1178-1185.

    disciussion 1185

    PubMed Abstract | Publisher Full Text OpenURL

  40. Sanchez-Munoz F, Dominguez-Lopez A, Yamamoto-Furusho JK: Role of cytokines in inflammatory bowel disease.

    World J Gastroenterol 2008, 14(27):4280-4288. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Podolsky DK: Inflammatory bowel disease.

    N Engl J Med 2002, 347(6):417-429. PubMed Abstract | Publisher Full Text OpenURL

  42. Bouma G, Strober W: The immunological and genetic basis of inflammatory bowel disease.

    Nat Rev Immunol 2003, 3(7):521-533. PubMed Abstract | Publisher Full Text OpenURL

  43. Silva MA, Menezes J, Wizman S, Gendron R, Oligny L, Seidman EG: Cytokine tissue levels as markers of disease activity in pediatric Crohn disease.

    Pediatr Res 2003, 54(4):456-461. PubMed Abstract | Publisher Full Text OpenURL

  44. Tincani A, Andreoli L, Bazzani C, Bosiso D, Sozzani S: Inflammatory molecules: a target for treatment of systemic autoimmune diseases.

    Autoimmun Rev 2007, 7(1):1-7. PubMed Abstract | Publisher Full Text OpenURL

  45. Hansen JJ, Holt L, Sartor RB: Gene expression patterns in experimental colitis in IL-10-deficient mice.

    Inflamm Bowel Dis 2009, 15(6):890-899. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Aitken AE, Richardson TA, Morgan ET: Regulation of drug-metabolizing enzymes and transporters in inflammation.

    Annu Rev Pharmacol Toxicol 2006, 46:123-149. PubMed Abstract | Publisher Full Text OpenURL

  47. Wright SD, Ramos RA, Tobias PS, Ulevitch RJ, Mathison JC: CD14, a receptor for complexes of lipopolysaccharide (LPS) and LPS binding protein.

    Science 1990, 249(4975):1431-1433. PubMed Abstract | Publisher Full Text OpenURL

  48. Klein W, Tromm A, Griga T, Fricke H, Folwaczny C, Hocke M, Eitner K, Marx M, Duerig N, Epplen JT: A polymorphism in the CD14 gene is associated with Crohn disease.

    Scand J Gastroenterol 2002, 37(2):189-191. PubMed Abstract | Publisher Full Text OpenURL

  49. Obana N, Takahashi S, Kinouchi Y, Negoro K, Takagi S, Hiwatashi N, Shimosegawa T: Ulcerative colitis is associated with a promoter polymorphism of lipopolysaccharide receptor gene, CD14.

    Scand J Gastroenterol 2002, 37(6):699-704. PubMed Abstract | Publisher Full Text OpenURL

  50. Narushima Y, Unno M, Nakagawara K, Mori M, Miyashita H, Suzuki Y, Noguchi N, Takasawa S, Kumagai T, Yonekura H, et al.: Structure, chromosomal localization and expression of mouse genes encoding type III Reg, RegIII alpha, RegIII beta, RegIII gamma.

    Gene 1997, 185(2):159-168. PubMed Abstract | Publisher Full Text OpenURL

  51. Bohn E, Bechtold O, Zahir N, Frick JS, Reimann J, Jilge B, Autenrieth IB: Host gene expression in the colon of gnotobiotic interleukin-2-deficient mice colonized with commensal colitogenic or noncolitogenic bacterial strains: common patterns and bacteria strain specific signatures.

    Inflamm Bowel Dis 2006, 12(9):853-862. PubMed Abstract | Publisher Full Text OpenURL

  52. Keilbaugh SA, Shin ME, Banchereau RF, McVay LD, Boyko N, Artis D, Cebra JJ, Wu GD: Activation of RegIIIbeta/gamma and interferon gamma expression in the intestinal tract of SCID mice: an innate response to bacterial colonisation of the gut.

    Gut 2005, 54(5):623-629. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Cash HL, Whitham CV, Behrendt CL, Hooper LV: Symbiotic bacteria direct expression of an intestinal bactericidal lectin.

    Science 2006, 313(5790):1126-1130. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Brandl K, Plitas G, Mihu CN, Ubeda C, Jia T, Fleisher M, Schnabl B, DeMatteo RP, Pamer EG: Vancomycin-resistant enterococci exploit antibiotic-induced innate immune deficits.

    Nature 2008, 455(7214):804-807. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Folch-Puy E, Granell S, Dagorn JC, Iovanna JL, Closa D: Pancreatitis-associated protein I suppresses NF-kappa B activation through a JAK/STAT-mediated mechanism in epithelial cells.

    J Immunol 2006, 176(6):3774-3779. PubMed Abstract | Publisher Full Text OpenURL

  56. Wijburg OL, Uren TK, Simpfendorfer K, Johansen FE, Brandtzaeg P, Strugnell RA: Innate secretory antibodies protect against natural Salmonella typhimurium infection.

    J Exp Med 2006, 203(1):21-26. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Murthy AK, Dubose CN, Banas JA, Coalson JJ, Arulanandam BP: Contribution of polymeric immunoglobulin receptor to regulation of intestinal inflammation in dextran sulfate sodium-induced colitis.

    J Gastroenterol Hepatol 2006, 21(9):1372-1380. PubMed Abstract | Publisher Full Text OpenURL

  58. Selby WS, Poulter LW, Hobbs S, Jewell DP, Janossy G: Heterogeneity of HLA-DR-positive histiocytes in human intestinal lamina propria: a combined histochemical and immunohistological analysis.

    J Clin Pathol 1983, 36(4):379-384. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Dommels YEM, Butts C, Zhu ST, Davy M, Martell S, Hedderly D, Barnett MPG, Broadley K, Roy NC: Characterization of inflammation in mdr1a-/- mice and identification of related gene expression changes in colonic epithelial cells.

    Genes Nutr 2007, 2(2):209-223. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Tannock GW: Molecular assessment of intestinal microflora.

    Am J Clin Nutr 2001, 73(2 Suppl):410S-414S. PubMed Abstract | Publisher Full Text OpenURL

  61. Dutka-Malen S, Evers S, Courvalin P: Detection of glycopeptide resistance genotypes and identification to the species level of clinically relevant enterococci by PCR.

    J Clin Microbiol 1995, 33(1):24-27. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Park EI, Paisley EA, Mangian HJ, Swartz DA, Wu MX, O'Morchoe PJ, Behr SR, Visek WJ, Kaput J: Lipid level and type alter stearoyl CoA desaturase mRNA abundance differently in mice with distinct susceptibilities to diet-influenced diseases.

    J Nutr 1997, 127(4):566-573. PubMed Abstract | Publisher Full Text OpenURL

  63. Kruschewski M, Foitzik T, Perez-Canto A, Hubotter A, Buhr HJ: Changes of colonic mucosal microcirculation and histology in two colitis models: an experimental study using intravital microscopy and a new histological scoring system.

    Dig Dis Sci 2001, 46(11):2336-2343. PubMed Abstract | Publisher Full Text OpenURL

  64. Smyth GK: Limma: linear models for microarray data. In Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. New York: Springer; 2005:397-420. Publisher Full Text OpenURL

  65. Ramakers C, Ruijter JM, Deprez RH, Moorman AF: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data.

    Neurosci Lett 2003, 339(1):62-66. PubMed Abstract | Publisher Full Text OpenURL

  66. Nones K, Dommels YE, Martell S, Butts C, McNabb WC, Park ZA, Zhu S, Hedderley D, Barnett MP, Roy NC: The effects of dietary curcumin and rutin on colonic inflammation and gene expression in multidrug resistance gene-deficient (mdr1a-/-) mice, a model of inflammatory bowel diseases.

    Br J Nutr 2009, 101(2):169-181. PubMed Abstract | Publisher Full Text OpenURL