BACKGROUND
the so-called "hygiene hypothesis" suggests that reduced exposure of children to microbes is associated with increased prevalence of common allergies in developed countries  <cit> .
at birth, the gastrointestinal tract is sterile and the neonatal immune response is characterized by a polarized t helper  <dig>  cytokine profile  <cit> .
during gut colonization by commensal microorganisms, the gut immune system is constantly challenged by a myriad of bacterial and food antigens.
gut colonization apparently plays a major role in driving the initial th2-skewed immune response toward a more finely balanced th1/th <dig> response, by boosting counterregulatory th <dig> immune responses  <cit> .
numerous studies using animal models have suggested the possible involvement of immunoregulatory lymphocytes, e.g., regulatory t cells  and/or interleukin- <dig>  producing b cells and cytokines  in intestinal homeostasis, which are driven by the intestinal bacterial burden  <cit>
many epidemiological studies suggest that there is an inverse relationship between infections in early childhood and the subsequent development of allergic diseases  <cit> .
therapy with broad-spectrum antibiotics is frequently performed in pediatric practice and children receiving this therapy within their first year of life are particularly prone to develop allergic diseases later in life  <cit> .
studies in animal models, such as gf animals, have also suggested that microbiota play a critical role in normal development of the immune system  <cit> .
oyama et al.
<cit>  reported that antibiotic use during infancy in mice promotes a shift in the th1/th <dig> balance toward th2-dominant immunity.
further, they demonstrated that gf mice do not develop oral tolerance, which was restored by microbial reconstitution in neonatal  mice but not in older mice  <cit> .
impairment of immune tolerance has been shown to augment disease in various models of allergy and/or autoimmunity, including diabetes onset in nod mice  <cit> , collagen-induced arthritis  <cit> , and experimental colitis  <cit> .
taken together, these findings provide new perspectives on the pathogenesis and recurrence of these diseases  <cit>
in the present study, to elucidate the impact of microbes on the immune system during the neonatal period, we performed microarray analysis of li, si, liv, and spl of mice with or without enteric microbiota, and of gf mice reconstituted with microbiota at different ages.
in order to extract useful information from the massive amount of gene expression data obtained by microarray assay, we employed a gene set enrichment analysis approach in the present study.
this type of analysis uses predetermined aggregations of genes  rather than individual genes to assess for coordinated expression in the samples.
single gene analysis may miss important effects on signaling because cellular processes often affect sets of genes acting in concert, with moderate effects on the strength of expression.
subtle but orchestrated changes of internally-related genes have often been found to be more important than a dramatic increase/decrease of a single gene.
therefore, it may be important to evaluate the statistical significance of changes in a gene aggregate or gene set, rather than in an individual gene.
since subramanian et al..  <cit>  initially proposed “gene set enrichment analysis ”, a number of algorithms optimized for this type of analysis have been developed  <cit> , and in the present study, we employed “metagene profiler ” developed by gupta et al.
<cit> .
unlike most existing methods, the main characteristic of mgp is that it evaluates statistical data for a set of genes independently from data of other gene sets.
because of the advantage of the analysis, it is logical to compare the results of tests for the same set of genes observed under different conditions, such as case–control experiments with multiple cases and time-course experiments, because the statistical evidence is evaluated using the same standard.
mgp was therefore suitable for the present study, which required multiple cross comparisons of the same samples  <cit> .
this study is part of a research project using gf and spf mice of the iqi strain, which has been established as an inbred strain from icr mice at the central institute for experimental animals .
we previously reported some of the results of microarray analyses conducted for the project.
these studies indicated that 1) activation of the ifn-α system in li differs significantly between gf and spf mice  <cit> , and 2) there is a striking commonality in transcriptome profiles between gf li and si, while the profiles of spf li and si share almost nothing in common  <cit> .
the present study, despite the use of different cohorts of animals, different versions of genechip arrays, and completely different strategies and algorithms for bioinformatics analysis from previous studies, has given support to our earlier findings.
furthermore, we demonstrated that the lack of neonatal encounter with commensal microorganisms may result in profound alteration of certain signaling pathways including tlrs, rac  <dig>  and type- <dig> ifn, which cannot be restored by later exposure to microbiota.
these findings may provide important insights into the molecular events underlying the interaction between neonatal immune systems and commensal microbiota.
RESULTS
generation of gnotobiotic mice and targets of microarray analysis
in the present study we prepared  <dig> groups of mice with different status of enteric flora: gf, spf, 0wexgf, and 5wexgf .
bacterial reconstitution was executed not by oral administration of fecal suspension but by cohabitation of pregnant gf mice with female spf mice of the same age, because 1) contact with microorganisms via a wide variety of parenteral routes such as birth canal, skin, eyes, ears, nostrils, airway, urinary tract, and vagina are also critically important for the establishment of relevant microbial flora, which differ profoundly among these anatomical sites, and 2) it is impossible to administer bacterial suspension orally to neonatal mice immediately after delivery.
using gf and spf mice of the same strain, we have analyzed the gene expression profiles of li and si in previous studies  <cit> .
in the present study, we have extended the target of transcriptome analyses to the spl and liv in addition to the intestines.
spl is a key part of body’s immune system and, therefore, comparison of the gene expression profiles of spl and intestines may help to elucidate the cross-talk between intestinal local immunity and general immunity.
the gene expression profile of liv enzymes, especially those of steroid and xenobiotic metabolism, has been shown to be profoundly affected by intestinal microflora via nutrients and metabolites transported through portal vein  <cit> .
further, translocation of bacteria-derived substances, e.g., lipopolysaccharide to the liv, is also known to have a great impact on the host immune system.
go analysis: overview and gene clustering
firstly, we have summarized the biological impact of enteric microbiota on host organs by mgp analysis on go  categories.
figure  <dig> shows venn diagrams of the numbers of differentially expressed genes and overrepresented go categories of biological process .
results from  <dig> comparisons  revealed alteration of an extremely large number of genes  in the si of spf mice compared with those of 0wexgf and 5wexgf mice.
there were approximately  <dig> to  <dig> commonly regulated probe sets in the li, si, liv, and spl among the  <dig> colonization models.
at least  <dig> probe sets were listed in each group in the li, si, and spl; however, no or only a small number of go categories were overrepresented in most cases with the exception of li in 0wexgf mice, si in spf mice, and spl in 5wexgf mice.
the overrepresented go categories are listed in table  <dig> , table  <dig>  and table  <dig> .
specific alteration of genes for antigen presentation in li of 0exgf mice, and for energy cycle and nervous development in si of spf mice, were identified from these tables.
mgp analysis was performed on go bp categories which contains 4 ~  <dig> genes .
the categories were sorted according to the integrated p value  and the top  <dig> are represented.
the categories related to antigen presentation are noted .
the results of unsupervised hierarchical cluster analysis of gene expression patterns are shown in figure  <dig>  figure  <dig>  figure  <dig> and figure  <dig>  the clusters denoting overrepresented go categories included immune system development/b cell activation/hematopoiesis  and metabolic process/lipid metabolic process .
analysis of signal pathway
next, in order to gain information regarding the molecular events evoked by enteric microbiota, we applied the gene set enrichment analysis program mgp to the signal pathway database of tranthpath.
figure  <dig>  figure  <dig>  and figure  <dig> show the signal pathways overrepresented in li, si and spl.
tlr, rac <dig>  and ifn-α signaling were commonly overrepresented in the li of spf and 0wexgf mice  and epidermal growth factor  and nerve growth factor  signaling were common in the li of 0wexgf and 5wexgf mice.
there was low commonality between 5wexgf and spf.
in 5wexgf mice, specific alteration of gene expression for chemotaxis was indicated.
the lists for si differed from those of li in spf and 0wexgf mice .
however, surprisingly, the list for si seemed essentially the same as that of li in 5wexgf mice .
in all  <dig> groups, the profiles of si were characterized by large scale alteration of egf, platelet-derived growth factor  and ngf/neurotrophin signaling.
overrepresentation of tlr pathway and rac <dig> signaling was observed in 0wexgf si but not in spf si.
in spl, only a few overrepresented signal pathways passed the statistical criteria , but overrepresentation of tlr pathway and rac <dig> signaling in 0wexgf mice and apoptosis pathway signaling in 5wexgf mice was noted.
no pathways reached statistical significance of q <  <dig>  in liv.
investigation of tlr signaling by quantitative rt-pcr
gene set enrichment analysis provides useful information about the statistical significance of changes in the expression of gene sets.
however, the analysis does not provide information on the degree and direction of change of expression for any particular gene.
furthermore, it is unlikely that any existing database would include all of the signaling pathways involved in complex physiological and pathophysiological events.
for examples, "tlr pathway" in the trantpath database seemed to contain relatively little information on tlr <dig> and tlr <dig>  therefore, in order to obtain more information regarding differences in tlr signaling, we investigated the expression of various tlr-related genes by real time rt-pcr analysis in li of the  <dig> groups of mice.
the results are summarized in figure  <dig> in terms of the difference of gene expression compared with spf mice.
the status of enteric microbiota  influenced various tlr-related genes in diverse ways.
the most prominent and specific changes were seen in tlr7/ <dig>  and irf <dig> .
discussion
in our previous study  <cit> , using gf and spf iqi mice and the old version of affimetrix genechip , we compared gene expression profiles of li in gf and spf mice.
a list of genes most affected by the presence of enteric microbiota was generated by the simple welch's t-statistic and enumeration in the order corresponding to the fold-change.
the list clearly indicated that the expression of genes related to the induction and production of type  <dig> ifn, such as irf <dig>  irf <dig>  isgf3g, ifit <dig> and g1p <dig> , was markedly decreased in spf mice compared with gf mice.
in situ hybridization and immunohistochemistry indicated that these ifn-related genes were expressed mainly in lamina proprial cd11b+ cells.
in another study  <cit> , using a synthetic ifn inducer and an herbal medicine that enhances ifn production, we demonstrated that the difference in the steady-state expression level of these ifn-related genes was correlated with the difference in the timing of ifn-α release triggered by the inducer.
in the present study, we used the same mouse strain but a different version of genechip  and employed a completely different strategy and algorithm for bioinformatic analysis.
in place of single gene analysis, we have adopted a gene set enrichment analysis with the primary aim of identifying the most affected gene sets in terms of statistical significance.
the predetermined gene sets of categories, cascades and pathways have been furnished from public  and commercial  databases.
the results indicated that, together with tlr signaling and rac <dig> signaling, the signaling cascade for the production of type  <dig> ifn was most prominently affected in spf and 0wexgf mice.
further, in good accordance with our previous paper  <cit> , the changes in gene expression profiles were dissimilar between si and li of spf mice, while those in gf mice showed substantial commonality between si and li.
the gene expression profile of spf mice was similar to that of 0wexgf, particularly in li, but not to that of 5wexgf mice.
the main differences were observed in tlr, rac <dig>  and ifn-α signaling.
tlrs play an important role in the recognition of microbes by host sentinel cells that contribute to subsequent innate and adaptive immune responses  <cit> .
tlrs recognize molecular patterns specific for microbes to eliminate pathogens and engender commensal colonization of symbiotic bacteria  <cit> .
therefore, the difference in tlr signaling may contribute to the altered immune function of 5wexgf mice.
similarly, 5wexgf mice may have unique characteristics with respect to the ifn-α response, which is critically important not only for protection against infection by virus and certain microorganism, but also for control of autoimmune responses.
only in 5wexgf intestines, both si and li, large scale alteration of chemokine signaling was observed.
in response to microbial stimuli, fetal intestinal epithelial cells release chemokines very rapidly, and this response may be a prerequisite for intestinal tolerance to commensal bacteria  <cit> .
acute chemokine release in response to microbial activation was subsequently maintained for several days in gf mice, but did not occur in mice harboring conventional microbiota due to tolerance acquisition at birth  <cit> .
to clarify how the observed changes in the expression profile of chemokine signaling in 5wexgf mice relates to the failure to acquire immune tolerance, it may be necessary to examine chemokine production in epithelial cells isolated from 5wexgf intestines.
our findings collectively suggest that the encounter with environmental microbiota during the specific time interval within the neonatal period is critically important for the development of normal immune system responsiveness to microorganisms, both commensal and infectious.
previous transcriptomics studies focusing on the role of commensal bacteria  <cit>  differ from the present study with respect to animal species , tissues , ages at the time of sampling , means of manipulation of microbiota , microarray platforms, and statistical methodologies.
however, the results of these studies are similar to ours in several respects such as prominent alteration of the expression profiles for antigen presentation, xenobiotic systems, and ifn signaling.
to shed additional light on the "hygiene hypothesis", a careful comparison of our results with those of schumann et al.
<cit>  would be particularly valuable because in the latter study, enteric bacteria were eradicated by antibiotic treatment of neonatal rats from postnatal day  <dig> to  <dig>  in the study by schumann and colleagues, in addition to significant down-regulation of antigen presentation systems, expression of paneth cell products such as α-defensins, matrylsin, and type iia phospholipase a <dig> decreased and expression of mast cell proteases increased as a result of a drastic reduction in enteric bacteria evoked by antibiotics.
similar changes were observed in the si of gf mice compared with the other  <dig> colonization models in the present study .
furthermore, schumann et al.. have also found that eradication of enteric bacteria had a greater effect, in terms of the number of affected probe sets, in si than in the colon.
clustering analysis of the affected probe sets in the study revealed only a few functional categories; i.e., endocytosis and vesicle-mediated transport in the proximal si, immune response in the distal si, and ion transport processes in the colon.
the present analysis based on mgp identified expression changes for many functional categories, including proton transport and vesicle-mediated transport in the si  and immune responses in the li .
in spite of the extensive differences in species, experimental settings, microarray platforms, and approaches to statistical analysis, a detailed comparative examination of these two transcriptome data sets might provide new insights into the impact of commensal microbiota colonization during the neonatal period.
mgp analysis was performed on go bp categories which contains 4 ~  <dig> genes .
the categories were sorted according to the integrated p value  and top  <dig> are represented.
the categories related to energy cycle are noted .
CONCLUSIONS
in summary, the gene expression profiles of mice with bacterial colonization at different times suggest that the encounter with environmental commensal microbiota during the specific time interval within neonatal period is essential for normal development of the immune system, especially of the li.
microbiota-mediated development of regulatory circuits of tlrs and type i ifn seem to play a particularly important role.
rectification of chemokine expression might also be involved in microbiota-related immunological dysregulation because expression of chemokines was activated specifically in mice colonized at pre-pubertal ages.
the present study provides important insights for clarification and refinement of the so-called "hygiene hypothesis".
mgp analysis was performed on go bp categories which contains 4 ~  <dig> genes .
the categories were sorted according to the integrated p value  and top  <dig> are represented.

 the so-called "hygiene hypothesis" suggests that reduced exposure of children to microbes is associated with increased prevalence of common allergies in developed countries. at birth, the gastrointestinal tract is sterile and the neonatal immune response is characterized by a polarized t helper cytokine profile. during gut colonization by commensal microorganisms, the gut immune system is constantly challenged by a myriad of bacterial and food antigens. gut colonization apparently plays a major role in driving the initial th2-skewed immune response toward a more finely balanced th1/th response, by boosting counterregulatory th immune responses. numerous studies using animal models have suggested the possible involvement of immunoregulatory lymphocytes, e.g., regulatory t cells and/or interleukin- producing b cells and cytokines in intestinal homeostasis, which are driven by the intestinal bacterial burden many epidemiological studies suggest that there is an inverse relationship between infections in early childhood and the subsequent development of allergic diseases. therapy with broad-spectrum antibiotics is frequently performed in pediatric practice and children receiving this therapy within their first year of life are particularly prone to develop allergic diseases later in life. studies in animal models, such as gf animals, have also suggested that microbiota play a critical role in normal development of the immune system. oyama et al. reported that antibiotic use during infancy in mice promotes a shift in the th1/th balance toward th2-dominant immunity. further, they demonstrated that gf mice do not develop oral tolerance, which was restored by microbial reconstitution in neonatal mice but not in older mice. impairment of immune tolerance has been shown to augment disease in various models of allergy and/or autoimmunity, including diabetes onset in nod mice, collagen-induced arthritis, and experimental colitis. taken together, these findings provide new perspectives on the pathogenesis and recurrence of these diseases in the present study, to elucidate the impact of microbes on the immune system during the neonatal period, we performed microarray analysis of li, si, liv, and spl of mice with or without enteric microbiota, and of gf mice reconstituted with microbiota at different ages. in order to extract useful information from the massive amount of gene expression data obtained by microarray assay, we employed a gene set enrichment analysis approach in the present study. this type of analysis uses predetermined aggregations of genes rather than individual genes to assess for coordinated expression in the samples. single gene analysis may miss important effects on signaling because cellular processes often affect sets of genes acting in concert, with moderate effects on the strength of expression. subtle but orchestrated changes of internally-related genes have often been found to be more important than a dramatic increase/decrease of a single gene. therefore, it may be important to evaluate the statistical significance of changes in a gene aggregate or gene set, rather than in an individual gene. since subramanian et al.. initially proposed “gene set enrichment analysis ”, a number of algorithms optimized for this type of analysis have been developed, and in the present study, we employed “metagene profiler ” developed by gupta et al.. unlike most existing methods, the main characteristic of mgp is that it evaluates statistical data for a set of genes independently from data of other gene sets. because of the advantage of the analysis, it is logical to compare the results of tests for the same set of genes observed under different conditions, such as case–control experiments with multiple cases and time-course experiments, because the statistical evidence is evaluated using the same standard. mgp was therefore suitable for the present study, which required multiple cross comparisons of the same samples. this study is part of a research project using gf and spf mice of the iqi strain, which has been established as an inbred strain from icr mice at the central institute for experimental animals. we previously reported some of the results of microarray analyses conducted for the project. these studies indicated that 1) activation of the ifn-α system in li differs significantly between gf and spf mice, and 2) there is a striking commonality in transcriptome profiles between gf li and si, while the profiles of spf li and si share almost nothing in common. the present study, despite the use of different cohorts of animals, different versions of genechip arrays, and completely different strategies and algorithms for bioinformatics analysis from previous studies, has given support to our earlier findings. furthermore, we demonstrated that the lack of neonatal encounter with commensal microorganisms may result in profound alteration of certain signaling pathways including tlrs, rac and type- ifn, which cannot be restored by later exposure to microbiota. these findings may provide important insights into the molecular events underlying the interaction between neonatal immune systems and commensal microbiota.  generation of gnotobiotic mice and targets of microarray analysis in the present study we prepared groups of mice with different status of enteric flora: gf, spf, 0wexgf, and 5wexgf. bacterial reconstitution was executed not by oral administration of fecal suspension but by cohabitation of pregnant gf mice with female spf mice of the same age, because 1) contact with microorganisms via a wide variety of parenteral routes such as birth canal, skin, eyes, ears, nostrils, airway, urinary tract, and vagina are also critically important for the establishment of relevant microbial flora, which differ profoundly among these anatomical sites, and 2) it is impossible to administer bacterial suspension orally to neonatal mice immediately after delivery. using gf and spf mice of the same strain, we have analyzed the gene expression profiles of li and si in previous studies. in the present study, we have extended the target of transcriptome analyses to the spl and liv in addition to the intestines. spl is a key part of body’s immune system and, therefore, comparison of the gene expression profiles of spl and intestines may help to elucidate the cross-talk between intestinal local immunity and general immunity. the gene expression profile of liv enzymes, especially those of steroid and xenobiotic metabolism, has been shown to be profoundly affected by intestinal microflora via nutrients and metabolites transported through portal vein. further, translocation of bacteria-derived substances, e.g., lipopolysaccharide to the liv, is also known to have a great impact on the host immune system. go analysis: overview and gene clustering firstly, we have summarized the biological impact of enteric microbiota on host organs by mgp analysis on go categories. figure  shows venn diagrams of the numbers of differentially expressed genes and overrepresented go categories of biological process. results from comparisons revealed alteration of an extremely large number of genes in the si of spf mice compared with those of 0wexgf and 5wexgf mice. there were approximately to commonly regulated probe sets in the li, si, liv, and spl among the colonization models. at least probe sets were listed in each group in the li, si, and spl; however, no or only a small number of go categories were overrepresented in most cases with the exception of li in 0wexgf mice, si in spf mice, and spl in 5wexgf mice. the overrepresented go categories are listed in table , table and table. specific alteration of genes for antigen presentation in li of 0exgf mice, and for energy cycle and nervous development in si of spf mice, were identified from these tables. mgp analysis was performed on go bp categories which contains 4 ~  genes. the categories were sorted according to the integrated p value and the top are represented. the categories related to antigen presentation are noted. the results of unsupervised hierarchical cluster analysis of gene expression patterns are shown in figure  figure  figure  and figure  the clusters denoting overrepresented go categories included immune system development/b cell activation/hematopoiesis and metabolic process/lipid metabolic process. analysis of signal pathway next, in order to gain information regarding the molecular events evoked by enteric microbiota, we applied the gene set enrichment analysis program mgp to the signal pathway database of tranthpath. figure  figure  and figure  show the signal pathways overrepresented in li, si and spl. tlr, rac and ifn-α signaling were commonly overrepresented in the li of spf and 0wexgf mice and epidermal growth factor and nerve growth factor signaling were common in the li of 0wexgf and 5wexgf mice. there was low commonality between 5wexgf and spf. in 5wexgf mice, specific alteration of gene expression for chemotaxis was indicated. the lists for si differed from those of li in spf and 0wexgf mice. however, surprisingly, the list for si seemed essentially the same as that of li in 5wexgf mice. in all groups, the profiles of si were characterized by large scale alteration of egf, platelet-derived growth factor and ngf/neurotrophin signaling. overrepresentation of tlr pathway and rac signaling was observed in 0wexgf si but not in spf si. in spl, only a few overrepresented signal pathways passed the statistical criteria, but overrepresentation of tlr pathway and rac signaling in 0wexgf mice and apoptosis pathway signaling in 5wexgf mice was noted. no pathways reached statistical significance of q <  in liv. investigation of tlr signaling by quantitative rt-pcr gene set enrichment analysis provides useful information about the statistical significance of changes in the expression of gene sets. however, the analysis does not provide information on the degree and direction of change of expression for any particular gene. furthermore, it is unlikely that any existing database would include all of the signaling pathways involved in complex physiological and pathophysiological events. for examples, "tlr pathway" in the trantpath database seemed to contain relatively little information on tlr and tlr therefore, in order to obtain more information regarding differences in tlr signaling, we investigated the expression of various tlr-related genes by real time rt-pcr analysis in li of the groups of mice. the results are summarized in figure  in terms of the difference of gene expression compared with spf mice. the status of enteric microbiota influenced various tlr-related genes in diverse ways. the most prominent and specific changes were seen in tlr7/ and irf. discussion in our previous study, using gf and spf iqi mice and the old version of affimetrix genechip, we compared gene expression profiles of li in gf and spf mice. a list of genes most affected by the presence of enteric microbiota was generated by the simple welch's t-statistic and enumeration in the order corresponding to the fold-change. the list clearly indicated that the expression of genes related to the induction and production of type ifn, such as irf irf isgf3g, ifit and g1p, was markedly decreased in spf mice compared with gf mice. in situ hybridization and immunohistochemistry indicated that these ifn-related genes were expressed mainly in lamina proprial cd11b+ cells. in another study, using a synthetic ifn inducer and an herbal medicine that enhances ifn production, we demonstrated that the difference in the steady-state expression level of these ifn-related genes was correlated with the difference in the timing of ifn-α release triggered by the inducer. in the present study, we used the same mouse strain but a different version of genechip and employed a completely different strategy and algorithm for bioinformatic analysis. in place of single gene analysis, we have adopted a gene set enrichment analysis with the primary aim of identifying the most affected gene sets in terms of statistical significance. the predetermined gene sets of categories, cascades and pathways have been furnished from public and commercial databases. the results indicated that, together with tlr signaling and rac signaling, the signaling cascade for the production of type ifn was most prominently affected in spf and 0wexgf mice. further, in good accordance with our previous paper, the changes in gene expression profiles were dissimilar between si and li of spf mice, while those in gf mice showed substantial commonality between si and li. the gene expression profile of spf mice was similar to that of 0wexgf, particularly in li, but not to that of 5wexgf mice. the main differences were observed in tlr, rac and ifn-α signaling. tlrs play an important role in the recognition of microbes by host sentinel cells that contribute to subsequent innate and adaptive immune responses. tlrs recognize molecular patterns specific for microbes to eliminate pathogens and engender commensal colonization of symbiotic bacteria. therefore, the difference in tlr signaling may contribute to the altered immune function of 5wexgf mice. similarly, 5wexgf mice may have unique characteristics with respect to the ifn-α response, which is critically important not only for protection against infection by virus and certain microorganism, but also for control of autoimmune responses. only in 5wexgf intestines, both si and li, large scale alteration of chemokine signaling was observed. in response to microbial stimuli, fetal intestinal epithelial cells release chemokines very rapidly, and this response may be a prerequisite for intestinal tolerance to commensal bacteria. acute chemokine release in response to microbial activation was subsequently maintained for several days in gf mice, but did not occur in mice harboring conventional microbiota due to tolerance acquisition at birth. to clarify how the observed changes in the expression profile of chemokine signaling in 5wexgf mice relates to the failure to acquire immune tolerance, it may be necessary to examine chemokine production in epithelial cells isolated from 5wexgf intestines. our findings collectively suggest that the encounter with environmental microbiota during the specific time interval within the neonatal period is critically important for the development of normal immune system responsiveness to microorganisms, both commensal and infectious. previous transcriptomics studies focusing on the role of commensal bacteria differ from the present study with respect to animal species, tissues, ages at the time of sampling, means of manipulation of microbiota, microarray platforms, and statistical methodologies. however, the results of these studies are similar to ours in several respects such as prominent alteration of the expression profiles for antigen presentation, xenobiotic systems, and ifn signaling. to shed additional light on the "hygiene hypothesis", a careful comparison of our results with those of schumann et al. would be particularly valuable because in the latter study, enteric bacteria were eradicated by antibiotic treatment of neonatal rats from postnatal day to in the study by schumann and colleagues, in addition to significant down-regulation of antigen presentation systems, expression of paneth cell products such as α-defensins, matrylsin, and type iia phospholipase a decreased and expression of mast cell proteases increased as a result of a drastic reduction in enteric bacteria evoked by antibiotics. similar changes were observed in the si of gf mice compared with the other colonization models in the present study. furthermore, schumann et al.. have also found that eradication of enteric bacteria had a greater effect, in terms of the number of affected probe sets, in si than in the colon. clustering analysis of the affected probe sets in the study revealed only a few functional categories; i.e., endocytosis and vesicle-mediated transport in the proximal si, immune response in the distal si, and ion transport processes in the colon. the present analysis based on mgp identified expression changes for many functional categories, including proton transport and vesicle-mediated transport in the si and immune responses in the li. in spite of the extensive differences in species, experimental settings, microarray platforms, and approaches to statistical analysis, a detailed comparative examination of these two transcriptome data sets might provide new insights into the impact of commensal microbiota colonization during the neonatal period. mgp analysis was performed on go bp categories which contains 4 ~  genes. the categories were sorted according to the integrated p value and top are represented. the categories related to energy cycle are noted.  in summary, the gene expression profiles of mice with bacterial colonization at different times suggest that the encounter with environmental commensal microbiota during the specific time interval within neonatal period is essential for normal development of the immune system, especially of the li. microbiota-mediated development of regulatory circuits of tlrs and type i ifn seem to play a particularly important role. rectification of chemokine expression might also be involved in microbiota-related immunological dysregulation because expression of chemokines was activated specifically in mice colonized at pre-pubertal ages. the present study provides important insights for clarification and refinement of the so-called "hygiene hypothesis". mgp analysis was performed on go bp categories which contains 4 ~  genes. the categories were sorted according to the integrated p value and top are represented. 
in this study, we examined the transcriptome profiles of the large intestine, small intestine, liver, and spleen of bacterial colonization models—specific pathogen-free mice, ex-germ-free mice with bacterial reconstitution at the time of delivery, and ex-germ-free mice with bacterial reconstitution at 5 weeks of age —and compared them with those of germ-free mice. hundreds of genes were affected in all tissues in each of the colonized models; however, a gene set enrichment analysis method, metagene profiler, demonstrated that the specific changes of gene ontology categories occurred predominantly in 0wexgf li, spf si, and 5wexgf spl, respectively. mgp analysis on signal pathways revealed prominent changes in toll-like receptor - and type interferon -signaling in li of 0wexgf and spf mice, but not 5wexgf mice, while 5wexgf mice showed specific changes in chemokine signaling. the present study provides important new information regarding the molecular mechanisms of the so-called "hygiene hypothesis". hygiene hypothesisgerm-freetoll-like receptortype interferonmetagene profiler epidemiological studies have suggested that the encounter with commensal microorganisms during the neonatal period is essential for normal development of the host immune system. rt-pcr analysis of tlr-related genes showed that the expression of interferon regulatory factor, a crucial rate-limiting transcription factor in the induction of type ifn, prominently decreased in 0wexgf and spf mice but not in 5wexgf and gf mice. basic research involving gnotobiotic mice has demonstrated that colonization at the age of 5 weeks is too late to reconstitute normal immune function.
 epidemiological studies have suggested that the encounter with commensal microorganisms during the neonatal period is essential for normal development of the host immune system. basic research involving gnotobiotic mice has demonstrated that colonization at the age of 5 weeks is too late to reconstitute normal immune function. in this study, we examined the transcriptome profiles of the large intestine, small intestine, liver, and spleen of bacterial colonization models—specific pathogen-free mice, ex-germ-free mice with bacterial reconstitution at the time of delivery, and ex-germ-free mice with bacterial reconstitution at 5 weeks of age —and compared them with those of germ-free mice.  hundreds of genes were affected in all tissues in each of the colonized models; however, a gene set enrichment analysis method, metagene profiler, demonstrated that the specific changes of gene ontology categories occurred predominantly in 0wexgf li, spf si, and 5wexgf spl, respectively. mgp analysis on signal pathways revealed prominent changes in toll-like receptor - and type interferon -signaling in li of 0wexgf and spf mice, but not 5wexgf mice, while 5wexgf mice showed specific changes in chemokine signaling. rt-pcr analysis of tlr-related genes showed that the expression of interferon regulatory factor, a crucial rate-limiting transcription factor in the induction of type ifn, prominently decreased in 0wexgf and spf mice but not in 5wexgf and gf mice.  the present study provides important new information regarding the molecular mechanisms of the so-called "hygiene hypothesis". hygiene hypothesisgerm-freetoll-like receptortype interferonmetagene profiler 