Reproducible Colonization of Germ-Free Mice With the Oligo-Mouse-Microbiota in Different Animal Facilities

The Oligo-Mouse-Microbiota (OMM¹²) is a recently developed synthetic bacterial community for functional microbiome research in mouse models (Brugiroux et al., 2016). To date, the OMM¹² model has been established in several germ-free mouse facilities world-wide and is employed to address a growing variety of research questions related to infection biology, mucosal immunology, microbial ecology and host-microbiome metabolic cross-talk. The OMM¹² consists of 12 sequenced and publically available strains isolated from mice, representing five bacterial phyla that are naturally abundant in the murine gastrointestinal tract (Lagkouvardos et al., 2016). Under germ-free conditions, the OMM¹² colonizes mice stably over multiple generations. Here, we investigated whether stably colonized OMM¹² mouse lines could be reproducibly established in different animal facilities. Germ-free C57Bl/6J mice were inoculated with a frozen mixture of the OMM¹² strains. Within 2 weeks after application, the OMM¹² community reached the same stable composition in all facilities, as determined by fecal microbiome analysis. We show that a second application of the OMM¹² strains after 72 h leads to a more stable community composition than a single application. The availability of such protocols for reliable de novo generation of gnotobiotic rodents will certainly contribute to increasing experimental reproducibility in biomedical research.

Introduction

The mammalian gut is a complex ecosystem, hosting a diverse microbial community that influences normal physiology and disease susceptibility (mainly) through its metabolic activities. The gut microbiome is highly dynamic throughout life and can vary substantially between individuals due to diet, lifestyle and genetic factors (Ursell et al., 2012; Falony et al., 2016). Microbiome community profiling and metagenome-based approaches have recently elucidated how inter-individual microbiome differences correlate with disease states and health (Knight et al., 2017; Hadrich, 2018). However, causal relationships between the microbiota and disease conditions can rarely be deciphered with these methods. Ethical and regulatory issues further limit the possibilities of human intervention studies and robust experimental animal models are urgently needed to research the causality between the gut microbiota and various human diseases.

The laboratory mouse is currently the primary experimental model organism in biomedical research (Taylor et al., 2008). The availability of numerous genetically engineered and mutant mouse strains greatly facilitates functional studies (Eppig et al., 2015). Mice raised in different research institutions, obtained from different vendors or the wild can exhibit profound differences in microbiota composition (Stecher et al., 2010; Wang et al., 2014; Rausch et al., 2016; Sadler et al., 2017; Thiemann et al., 2017), mimicking the inter-personal microbiota variation in human populations. Apart from vendor-specific microbiomes, genotype and environmental conditions such as diet, cage-type, temperature and bedding can profoundly influence microbiota composition and function (Hildebrand et al., 2013). These microbiome differences between genetically identical mouse models have led to the serendipitous discovery of numerous microbiota-dependent disease phenotypes in the past, including immune-cell priming, colitis susceptibility and resistance to infections (Ivanov et al., 2009; Surana and Kasper, 2017; Velazquez et al., 2019). Overall, there is growing evidence that the microbiota is a major confounding factor, which complicates cross-study comparisons and eventually jeopardizes experimental reproducibility (Franklin and Ericsson, 2017). To identify the mechanisms of host genetics-imposed control of the microbiome, adequate experimental design is instrumental (Laukens et al., 2015; McCoy et al., 2017; Mamantopoulos et al., 2018). In this respect, gnotobiology (greek: gnosis: knowledge; bios: life; logos: study) has gained importance as it allows for optimal control of microbiota within and in between animal facilities. Besides, gnotobiology has become an essential method to mechanistically investigate microbiota functioning and to assess causality in disease-associated alterations of gut microbiota composition (Trexler and Reynolds, 1957; Orcutt et al., 1987). Gnotobiotic mice can be generated by microbial reconstitution of germ-free mice with fecal transplants, single organisms and defined mouse- or human-derived microbial consortia (Clavel et al., 2016; McCoy et al., 2017). Several recent studies have focused on the isolation of bacteria from humans and mice and the establishment of culture collections (Hugon et al., 2015; Lagkouvardos et al., 2016; Forster et al., 2019). The most important requirement for the assembly of defined bacterial consortia is the availability of well-characterized and genome-sequenced strains – preferably in public culture collections. Further, state of the art methods to trace and quantify each member of a microbial consortium are needed for quality control and functional studies. Finally, experimental protocols should be optimized to ensure reproducible colonization of germ-free mice with a given microbial consortium.

A majority of studies uses human-derived bacteria to colonize gnotobiotic mice even though phylogenetic differences between human and mouse microbiota (Xiao et al., 2015) may affect microbe-host interaction and long-term stability of colonization. To circumvent this problem, we have developed a model based on twelve bacteria from the murine gut microbiota (Brugiroux et al., 2016; Garzetti et al., 2017). This Oligo-Mouse-Microbiota (OMM¹²) provides several improvements to other defined consortia, including broad phylogenetic diversity, public availability of the strains from the German Type Culture Collection¹ (Brugiroux et al., 2016). Most importantly, the OMM¹² exhibits long-term stability in gnotobiotic mice (Brugiroux et al., 2016). Correspondingly, C57BL/6 mice stably colonized with OMM¹² have recently also been designated stable defined moderately diverse microbiota mice (sDMDMm2) (Li et al., 2015). The Altered Schaedler Flora (ASF), another frequently used model for gnotobiotic research for more than 40 years is also based on murine bacteria (Wymore Brand et al., 2015). Compared to ASF-colonized mice, OMM¹² mice show increased resistance to pathogen colonization (Brugiroux et al., 2016). These characteristics indicate that OMM¹² mice mimic the normal physiology better than ASF mice and are therefore preferable for microbiome intervention studies that investigate microbiome-correlated diseases. Accordingly, several studies have already implemented the OMM¹² as model to systematically probe for a causal role of individual microbes in protection against different pathogens along the lines of Koch’s postulates (Brugiroux et al., 2016; Studer et al., 2016; Herp et al., 2019). To date, the OMM¹² model is used by over 30 research groups world-wide to address fundamental research questions related to microbial ecology, metabolism, mucosal immunology and infection biology (Li et al., 2015; Studer et al., 2016; Uchimura et al., 2018).

Here, we present a feasibility trial and scrutinize colonization of germ-free C57BL/6 mice with the OMM¹² community. We established a protocol using a standardized OMM¹² inoculum and compared community composition when inoculations were performed repeatedly in the same facility or at different facilities. We show that double inoculation leads to efficient introduction of consortium members in different facilities. Notably, we observed subtle differences in the absolute abundance of some of the strains between different experiments. This also translates to differences in community profiles. Whether this induces functional changes to the community and phenotypic differences in the mice should be a matter of future investigation. Still, care should be taken when comparing results obtained with the model at different facilities. Overall, our study shows that generating mice colonized with the OMM¹² synthetic bacterial consortium shows excellent reproducibility between different animal facilities.

Materials and Methods

Animal Facilities

Five European germ-free rodent facilities participated in this study. All mouse experiments were approved by the local authorities and performed according to the legal requirements. Detailed questionnaires were distributed to record the mouse husbandry conditions in the different breeding facilities (Table 1). Germ-free status was routinely confirmed by aerobic and anaerobic culture as well as Sytox green (Invitrogen) and Gram staining (Harleco) of caecal contents to detect unculturable contaminants. C57BL/6J Agr2^–/– mice were provided by David Erle (Park et al., 2009) and re-derived germ-free from Agr2^± and Agr2^–/– conventional mice as described and colonized with OMM¹² to generate an isobiotic mouse line (Herp et al., 2019). Agr2 encodes a disulfide isomerase that is required for folding and export of the mucin Muc2. For the experiment shown in Figure 1, germ-free Agr2^± mice were used. Agr2^± mice behave like wild-type mice in terms of Muc2 secretion and mucus structure (Bergstrom et al., 2014). For all other experiments, wild-type C57BL/6J mice were used.

TABLE 1

Table 1. Description of housing conditions at Germ-free animal facilities.

FIGURE 1

Figure 1. Dynamics of fecal community composition after inoculation of OMM¹² colonization in germ-free mice. (A) Experimental scheme. Germ-free mice were inoculated with the OMM¹² mixture and kept in a germ-free isolator in facility 4. Total number of mice and collection time points of fecal samples are indicated. (B) Relative abundance of fecal microbiota composition at the indicated time points. Abundance of individual strains is shown as relative abundance and expressed as% of cumulative 16S rRNA gene copy numbers of all OMM¹² strains. One bar corresponds to one mouse. (C) PCoA based on the distance matrix of Bray–Curtis dissimilarity of relative OMM¹² abundance profiles shows the effect of time after inoculation. Points are colored by time (days) after inoculation. Samples taken from two mice are connected to visualize their trajectory during time.

Bacterial Cultivation

All glassware and media used for cultivation were kept under anoxic conditions (3% H₂, rest N₂) in an anaerobic chamber for at least 2 days before the start of the experiment.

Glycerol cryostocks (for preparation see Brugiroux et al., 2016) of individual OMM¹² strains (Table 2) were thawed in a 1% Virkon S (V.P. Products) solution (37°C) and the entire content of the vial was transferred into 100 ml Wheaton glass serum bottles (Sigma) sealed with a butyl rubber stoppers (Geo-Microbial Technologies) containing 10 ml of Anaerobic Akkermansia Medium (AAM; 18.5 g l^–1 brain heart infusion (BHI), 5 g l^–1 yeast extract, 15 g l^–1 trypticase soy broth, 2.5 g l^–1 K₂HPO₄, 1 mg l^–1 haemin, 0.5 g l^–1 glucose, 0.4 g l^–1 Na₂CO₃, 0.5 g l^–1 cysteine hydrochloride, 5 mg l^–1 menadione, 3% complement-inactivated fetal calf serum). These subcultures were gassed (7% H₂, 10% CO₂, 83% N₂) and incubated at 37°C for 24 h.

TABLE 2

Table 2. Growth conditions for OMM¹² strains.

Subsequently, 100 μl of each subculture was transferred into a 100 ml Wheaton glass serum bottle containing 10 ml of AAM. These cultures were gassed (7% H₂, 10% CO₂, 83% N₂) and incubated at 37°C for 1 or 2 days depending on the growth rate (Table 2). Afterward, culture purity of each strain was confirmed by Gram staining and 16S rRNA gene sequencing. The OD₆₀₀ of individual cultures was determined and all cultures were adjusted to the lowest OD₆₀₀ value by dilution. The respective culture volumes of all strains were transferred into a 50 ml Falcon tube under anoxic conditions. For cryopreservation, glycerol supplemented with palladium black crystals (Sigma-Aldrich) was added to these bacterial mixtures [final concentration of 10% (v/v)]. 1 ml aliquots were prepared in 1.5 ml glass vials (Sigma-Aldrich), sealed with butyl-rubber stoppers (Wheaton) and aluminum crimp seals (Sigma-Aldrich). Mixtures were frozen at −80°C within 1 h of preparation. Frozen aliquots were shipped to the different facilities on dry ice and stored at −80°C.

Inoculation of Mice and Fecal Sampling

The frozen OMM¹² mixtures were thawn in a 1% Virkon S (V.P. Produkte) disinfectant solution (37°C) and used for inoculation of germ-free mice in gnotocages or germ-free isolators. In any case, the mixture was used within 30 min after thawing. Mice were inoculated by gavage (50 μl orally, 100 μl rectally). Exposure of the mixture to oxygen was restricted to a short time (up to 5 min). For the double inoculation protocol, inoculation was repeated 72 h after the initial inoculation using the same protocol. To confirm the colonization of the 12 strains, fresh fecal pellets were obtained and frozen at −80°C within 30 min.

gDNA Extraction From Fecal Samples

Fecal samples were shipped on dry ice. DNA extraction and qPCR were performed centralized at the same laboratory to minimize experimental bias, which is known to be introduced by laboratory-specific experimental procedures. Fecal gDNA was either extracted using the QIAamp DNA Stool Mini Kit (Qiagen; time course shown in Figure 1) or a phenol-chloroform based protocol (all other experiments). The QIAamp DNA Stool Mini Kit protocol was performed following the manufacturer’s instructions with the following modifications. An initial bead-beating step using differentially sized beads [Zirkonia beads: 0.5–0.75 mm (BioSpec products) and acid-washed glass beads: <100 μm (Sigma-Aldrich)] was included and 20 mg/ml lysozyme was added to the lysis buffer. gDNA extraction using the phenol-chloroform based protocol was performed as described previously (Herp et al., 2019). The resulting gDNA was purified using the NucleoSpin gDNA clean-up kit (Macherey-Nagel).

Quantitative PCR of Bacterial 16S rRNA Genes

Quantitative PCR (qPCR) was performed as described previously (Brugiroux et al., 2016). All samples were analyzed by the same person in a centralized way. Briefly, OMM strain-specific 16S rRNA primers and hydrolysis probes were used for amplification. Standard curves using linearized plasmids containing the 16S rRNA gene sequence of the individual OMM¹² strains were used for absolute quantification of 16S rRNA gene copy numbers of individual strains. Since the fecal weight was not always available, 16S rRNA gene copy numbers were normalized to equal volumes of extracted DNA, assuming that DNA extraction is equally efficient between different samples. We confirmed that there is a linear relationship between stool weight and extracted DNA concentration (Supplementary Figure S1).

Statistical Analysis

For comparison of absolute abundance levels of OMM¹² strains between experiments, Kruskal–Wallis test with Dunn’s multiple comparison test was performed using GraphPad Prism version 5.01 for Windows (GraphPad Software). p-values below 0.05 were considered as statistically significant (^∗p < 0.05, ^∗∗p < 0.01, ^∗∗∗p < 0.001).

The vegist function of the R library vegan version 2.5–4 was employed to obtain Bray–Curtis (BC) dissimilarities between the samples based on relative abundance estimates. Principal coordinate analysis was performed in R using ade4 package and figures were generated using the ggplot2 library. Permutational multivariate analyses of variance (PERMANOVA) were performed in R using the function adonis. We used the function capscale with the variable “facility” as constraint to estimate the effect of the facility on the overall BC variance. Statistical significance of the ordinations as well as confidence intervals for the variance were determined by an ANOVA-like permutation test (functions permutest and anova.cca) with 5,000 permutations. The heatmap (Supplementary Figure S1) was generated using the ComplexHeatmap library on the BC values.

To estimate the within facility stability we calculated the mean of all pairwise BC dissimilarities of community profiles of mice from the same facility. Between-facility estimates were calculated using the mean value of all pairwise distances between the community profile of each mouse of one facility to the community profile of all mice of different facilities. R scripts are available under https://github.com/philippmuench/OligoMM-facilities.

Results

OMM¹² Colonization Dynamics After Oral Inoculation

In order to determine the time necessary for stable community formation, we tested temporal dynamics of the OMM¹² community after oral inoculation to germ-free mice in one facility. Changes in OMM¹² microbial community composition were closely monitored over time in fecal samples for 99 days (Figure 1A). In this experiment, we used two germ-free breeders from a heterozygous C57BL/6J Agr2^± breeding, which were reconstituted with the OMM¹² consortium to generate an isobiotic mouse line (Herp et al., 2019). Fecal community composition, represented as relative abundance profiles quantified by strain-specific qPCR, rapidly changed within the first week after OMM¹² inoculation. Beyond day 7 post inoculation, OMM¹² community stabilized and profiles remained highly similar afterward (Figure 1B). Principal Coordinates Analysis (PCoA) of Bray–Curtis dissimilarities showed a gradual shift with time toward a more similar community (Figure 1C). At early time points, high abundance of Bifidobacterium longum subsp. animalis YL2 and Enterococcus faecalis KB1 was observed, which gradually declined within the first week (Figures 2B,J). Inverse colonization dynamics were seen for several other strains (Figure 2). Notably, Muribaculum intestinale YL27 took 4 days to reach detectable levels in all mice (Figure 2I). Colonization levels for Clostridium innocuum I46 and Akkermansia muciniphila YL44 were constant throughout (Figure 2K). The number of 16S rRNA gene copies of Acutalibacter muris KB18 remained below the detection limit in the majority of fecal samples, which was observed previously in stably colonized OMM¹² mice (Brugiroux et al., 2016). We conclude that the OMM¹² consortium adopts a stable composition between 10 and 20 days post-inoculation. Therefore, we reasoned, that a 3-week colonization phase after inoculation is sufficient to verify successful colonization of OMM¹² in the feasibility trial outlined below.

FIGURE 2

Figure 2. Absolute abundance of individual OMM¹² strains in time course analysis of OMM¹² colonization in germ-free mice. Germ-free mice were inoculated with the OMM¹² mixture and kept in a germ-free isolator. Fecal samples were collected at different time points for microbiota analysis. Absolute abundance of each strain was determined by a strain-specific qPCR assay and is plotted as 16S rRNA gene copy numbers of the individual strains per μl of extracted gDNA: (A) Lactobacillus reuteri I49, (B) Enterococcus faecalis KB1, (C) Blautia coccoides YL58, (D) Clostridium innocuum I46, (E) Flavonifractor plautii YL31, (F) Clostridium clostridioforme YL32, (G) Acutalibacter muris KB18, (H) Bacteroides caecimuris I48, (I) Muribaculum intestinale YL27, (J) Bifidobacterium longum subsp. animalis YL2, (K) Akkermansia muciniphila YL44, (L) Turicimonas muris YL45. Statistical analysis was performed using Kruskal-Wallis test with Dunn’s multiple comparison test (^∗p < 0.05, ^∗∗p < 0.01, ^∗∗∗p < 0.001). Green symbols indicate samples collected <20 days post inoculation. Dotted lines indicate detection limits.

A Single Inoculation Does Not Lead to Reproducible Introduction of Muribaculum intestinale in Different Germ-Free Mouse Facilities

Next, we investigated, whether germ-free mice can be reproducibly associated with OMM¹² in different germ-free mouse facilities across four participating institutions in Germany and Switzerland. Germ-free C57Bl/6J mice (n = 2–5) at different facilities were orally inoculated once with the same batch of inoculum, to avoid variations introduced by differences in the inoculum. We obtained feces from the animals and different time points (day 10–72) post-inoculation (Figure 3A). Facility-specific characteristics are outlined in Table 1.

FIGURE 3

Figure 3. Colonization dynamics of OMM¹² mice in four different germ-free facilities after single-dose inoculation. (A) Experimental scheme. Germ-free C57BL/6J mice were inoculated with the OMM¹² mixtures and kept in germ-free isolators or gnotocages at four different animal facilities (1-1, 2, 3, 4-1); the number of mice and collection time points of fecal samples are indicated. (B) Fecal microbiota composition at the different time points, displayed as relative abundance and expressed as the fraction of cumulated 16S rRNA gene copy numbers. One bar corresponds to one mouse. (C) PCoA based on Bray–Curtis dissimilarity (relative abundances) between samples obtained from mice in different facilities. Points are colored by facility.

Relative OMM¹² abundance profiles of mice from different facilities showed, at large, high similarity (Figure 3B). PCoA based on Bray–Curtis dissimilarity (relative abundances) between samples showed that community composition was overall similar with the exception of facility 2, which clustered separately (Figure 3C). We used canonical analysis of principal coordinates (CAP; Anderson and Willis, 2003) to estimate the influence of the facility on the beta diversity. CAP analysis constrained by the facility revealed that the facility explains 35% of the overall variance of Bray–Curtis dissimilarity between samples from different facilities (p < 0.001). Based on PERMANOVA analysis of Bray–Curtis dissimilarities, facility “2” was clearly distinguishable from the other three facilities (Figure 3C and Supplementary Table S1).

Absolute abundance, as determined by qPCR revealed that 7 of the 12 species were detected in all mice at the different facilities at comparable levels (Figure 4). The absolute abundance of Lactobacillus reuteri I49, E. faecalis KB1, B. longum subsp. animalis YL2 and A. muris KB18 varied substantially between the facilities and was below the detection limit in most samples (Figures 4A,B,G,J). Additionally, Muribaculum intestinale YL27 was only detectable in facility “2” (Figure 4I). This may explain the notable different community profile of mice from this facility (Figure 3C and Supplementary Table S1). We reasoned that single inoculation of mice might not be sufficient to ensure reliable colonization of M. intestinale YL27.

FIGURE 4

Figure 4. Colonization dynamics of OMM¹² mice in four different germ-free facilities after single-dose inoculation: absolute abundance of individual OMM¹² strains. Absolute abundance of each strain was determined using a strain-specific qPCR assay for the experiment described in Figure 3, Data are plotted as 16S rRNA gene copy numbers of the individual strains per μl of extracted gDNA: (A) Lactobacillus reuteri I49, (B) Enterococcus faecalis KB1, (C) Blautia coccoides YL58, (D) Clostridium innocuum I46, (E) Flavonifractor plautii YL31, (F) Clostridium clostridioforme YL32, (G) Acutalibacter muris KB18, (H) Bacteroides caecimuris I48, (I) Muribaculum intestinale YL27, (J) Bifidobacterium longum subsp. animalis YL2, (K) Akkermansia muciniphila YL44, (L) Turicimonas muris YL45. Statistical analysis was performed using Kruskal-Wallis test with Dunn’s multiple comparison test (^∗p < 0.05, ^∗∗p < 0.01, ^∗∗∗p < 0.001). Green points indicate samples collected <20 days post inoculation.

Double-Dose Application Increases Reproducibility of OMM¹² Colonization in Different Germ-Free Mouse Facilities

Since single application of the OMM¹² mixture did not lead to reproducible colonization of OMM¹² strains across different germ-free mouse facilities, we aimed to amend the protocol. Previous work indicated that consecutive inoculations might increase the chance of successful introduction of oxygen-sensitive members of a consortium (Becker et al., 2011) (Taconic protocols). Therefore, we modified the initial inoculation protocol and applied the mixture twice with 72 h in-between inoculations. This time, three facilities participated and two independent trials were performed at facility “4.” We found that relative OMM¹² abundance profiles of mice from different facilities and trials were rather uniform (Figure 5B). PCoA based on Bray–Curtis dissimilarity (relative abundances) between samples showed that community composition was overall similar between facilities (Figure 5C). CAP analysis constrained by the facility revealed that in this trial, the facility explains 20% of the overall variance in Bray–Curtis dissimilarity of the data (p < 0.001). In this trial, no obvious clustering was apparent between samples from different facilities (Figure 5C), yet based on PERMANOVA analysis of Bray–Curtis dissimilarities, some differences between facilities were observed (Supplementary Table S2). Muribaculum intestinale YL27 was reliably detected in all mice colonized in the different facilities (Figure 6I), which is a substantial improvement compared to the single-dose experiment. However, levels of Lactobacillus reuteri I49, E. faecalis KB1, B. longum subsp. animalis YL2 and A. muris KB18 still varied between the tested facilities and were below detection limit in some samples (Figure 6).

FIGURE 5

Figure 5. Colonization dynamics of OMM¹² mice in four different germ-free facilities after double-dose inoculation reveals high reproducibility. (A) Experimental scheme. Germ-free C57BL/6 mice were inoculated twice with the OMM¹² mixtures and kept in germ-free isolators or gnotocages at different facilities (1-2, 5, 4-2, 4-3). In case of facility “1” and “4,” inoculations were done on several independent occasions. The number of mice and collection time points of fecal samples are indicated. (B) Fecal microbiota composition at the different time points. Microbiota composition is shown as relative abundance and expressed as the fraction of cumulated 16S rRNA gene copy numbers. One bar corresponds to one mouse. (C) PCoA based on the distance matrix of Bray–Curtis dissimilarity of relative OMM¹² abundance profiles shows samples obtained from mice in different facilities. Points are colored by facility.

FIGURE 6

Figure 6. Absolute abundance of OMM¹² strains in different germ-free facilities after double-dose inoculation. Absolute abundance of each strain was determined using a strain-specific qPCR assay for the experiment described in Figure 5, Data are plotted as 16S rRNA gene copy numbers of the individual strains per μl of extracted gDNA: (A) Lactobacillus reuteri I49, (B) Enterococcus faecalis KB1, (C) Blautia coccoides YL58, (D) Clostridium innocuum I46, (E) Flavonifractor plautii YL31, (F) Clostridium clostridioforme YL32, (G) Acutalibacter muris KB18, (H) Bacteroides caecimuris I48, (I) Muribaculum intestinale YL27, (J) Bifidobacterium longum subsp. animalis YL2, (K) Akkermansia muciniphila YL44, (L) Turicimonas muris YL45. Statistical analysis was performed using Kruskal–Wallis test with Dunn’s multiple comparison test (^∗p < 0.05, ^∗∗p < 0.01, ^∗∗∗p < 0.001). Green color indicate samples collected <20 days post inoculation.

Finally, we assessed variations in overall community profiles within and between the facilities for the two different trials. To this end, we compared Bray–Curtis dissimilarities between samples of the same facility to dissimilarities between samples of different facilities. An overview of the pair-wise compositional Bray–Curtis dissimilarity of relative community profiles of mice between and across facilities and studies (single-dose and double-dose) is shown in Supplementary Figure S2. For the first trial (single-dose inoculation), the mean of pairwise Bray–Curtis dissimilarity values of the community profile of mice housed in the same facility is lower than the mean pairwise Bray–Curtis dissimilarity between samples from mice located in different facilities. The mean Bray–Curtis dissimilarity for within-and between facility are 0.12 ± 0.04 (mean ± SD) and 0.18 ± 0.04, respectively; p-value = 0.07; paired two-sided t-test (Figure 7A). For the double-dose trial, the mean Bray–Curtis dissimilarity was lower (0.16 ± 0.06 and 0.19 ± 0.03 for within and between facilities, respectively; p-value = 0.25; paired two-sided t-test, Figure 7B).

FIGURE 7

Figure 7. OMM¹² community profiles within and between different facilities for the two trials. The mean of pairwise BC dissimilarity values of the community profiles of mice housed in the same facility (within facility analysis) was plotted against the mean pairwise BC dissimilarity values of the community profile of mice located in different facilities for single-dose inoculation trial for (A) Single-dose inoculation and (B) double-dose inoculation trial. For single-dose inoculation, the mean BC dissimilarity is 0.12 ± 0.04 (mean ± SD) and 0.18 ± 0.04, within and between the facilities, respectively (p-value = 0.07, paired two-sided t-test). For double-dose inoculation, the mean BC dissimilarity is 0.16 ± 0.06 and 0.19 ± 0.03 for within and between the facilities, respectively (p-value = 0.25; paired two-sided t-test).

Discussion

Variations in the gut microbiota within and between animal facilities can be a major factor accounting for the lack of reproducibility of animal models of human biology and disease (Franklin and Ericsson, 2017). Several guidelines were established to optimally control for microbiota differences when using different genetic strains of experimental mice within an animal unit (McCoy et al., 2017; Wullaert et al., 2018). Gnotobiotic mouse models based synthetic microbial consortia are becoming increasingly popular. In particular, these models offer the opportunity to generate “isobiotic” mice, which may significantly enhance experimental reproducibility across different institutions (Macpherson and McCoy, 2015).

This study reports on the first comparative inter-facility trial conducted to evaluate and optimize the effectiveness of a protocol for colonization of germ-free mice with a synthetic bacterial community. The majority of current protocols for introducing synthetic bacterial communities to germ-free mice use mixtures generated from fresh pure bacterial cultures for inoculations (Faith et al., 2011; Desai et al., 2016; Gomes-Neto et al., 2017). These protocols require a sophisticated cultivation setup, including devices for anaerobic bacterial cultivation. Our protocol overcomes these limitations by generating frozen aliquots of a mixture of the strains, which can be distributed, thawed and directly applied. We note that in our trial, frozen aliquots of bacterial mixtures remain viable for at least 18 months at −80°C. Another common method for generating gnotobiotic mice is by co-housing with a colonized donor animal (Geuking et al., 2011), a setup which requires no expertise for bacterial cultivation. When using this approach, it should be considered that serial passage of a bacterial community within the same facility or in between facilities, e.g., through breeding, might promote genomic diversification of the individual community members by consecutive rounds of within-host selection (Robinson et al., 2018). The evolved bacterial community would differentiate genetically and functionally from the parental strains over time. To date, the degree and temporal course of within-host evolution of microbial communities is not known and experimental data are only available for evolution of individual bacterial populations in the gut (Leatham et al., 2005; De Paepe et al., 2011; Barroso-Batista et al., 2014). Based on the analysis of individual bacterial populations, it is expected that mutants emerge rapidly and are selected based on improved competitiveness within days. In order to retain genetic identity of a minimal bacterial consortium, it is advisable to regenerate gnotobiotic mouse lines every 18–24 months using original cultures. This period is the result of cost-benefit considerations, keeping the degree of genomic diversification within acceptable boundaries but at the same time, minimizing costs and experimental efforts.

In this trial, we used aliquots of the same batch of frozen OMM¹² mixed cultures for the inoculations. This allowed us to compare efficiency of inoculation between different facilities using the same starting material. Thereby the inoculum could be eliminated as confounding factor in microbial community establishment in the gut. The gut microbiome of mice can be significantly influenced by husbandry-related factors, such as type of laboratory animal diet (Hildebrandt et al., 2009; Ooi et al., 2014) water (Sofi et al., 2014), housing effects, genetic background (Deloris Alexander et al., 2006; Hildebrand et al., 2013) and a wide range of other environmental and stress-related factors (Bangsgaard Bendtsen et al., 2012). Many of these variables are likely to differ across germ-free animal facilities and account for the small but measurable differences in OMM¹² community composition observed in our study. Further, due to facility-specific differences in procedures and experimental protocols, it was not possible to match sex, age and number of inoculated mice and obtain fecal samples at matched time points post inoculation. This may account for part of the facility-dependent differences and should be optimized in future trials.

Our results suggest that a double-dose application with a 72 h interval improves the engraftment of Muribaculum intestinale YL27. Representatives of the Muribaculaceae family (former S24-7) are highly diverse and dominant members of the murine gut microbiota (Lagkouvardos et al., 2019). M. intestinale YL27, the first cultured representative, is strictly anaerobic and genome-based prediction indicates the potential to degrade complex carbohydrates (Lagkouvardos et al., 2016; Lee et al., 2019). In our study, M. intestinale YL27 showed slow colonization dynamics after oral inoculation compared to all other OMM¹² strains disclose that expansion time can take up to 20 days. The oxidation/reduction (O/R) potential is increased in germ-free mice (+200 mV) compared to conventional mice but becomes reduced in response to colonization with a complex microbiota (−200 mV) (Celesk et al., 1976). As M. intestinale YL27 is oxygen-sensitive, we reason that high O/R potential in the germ-free mouse gut may inhibit its expansion early after inoculation. Although a recent study showed that the luminal contents of germ-free mice can chemically consume oxygen (e.g., via lipid oxidation reactions), the gut lumen of germ-free and antibiotic-treated mice may also exhibit increased luminal oxygen concentration compared to mice colonized with a complex microbiota (Friedman et al., 2018). Oxygen-tolerant members of the microbiota are among the first colonizers in a germ-free environment, after antibiotic treatment or in the course of intestinal colonization of the neonatal gut. They are thought to consume oxygen and anaerobic electron acceptors and to reduce the O/R potential sufficiently for oxygen-sensitive strains to colonize (Reese et al., 2018). Presumably, colonization dynamics of the OMM¹² consortium are subject to similar principles: Some members of the consortium (Enterococcus faecalis, Bifidobacterium longum) that are also among the early colonizers of the human neonatal gut predominate at early colonization stages. Obligate anaerobic bacteria (Bacteroides acidifaciens, Clostridiales) follow with a delay of 3 days. We assume that a second dose of OMM¹² administered when the O/R potential has already been lowered by the initial colonizers increases the chance of successful engraftment of obligate anaerobes such as Muribaculum intestinale.

In summary, our study demonstrates that germ-free mice in different facilities can be reproducibly associated with the OMM¹² synthetic bacterial community, employing a protocol using previously frozen aliquots of a mixture of the strains. Furthermore, application of the consortium at two consecutive time points increased the chance of community engraftment. We envision that guidelines and validated protocols for generation of gnotobiotic models based on synthetic microbial communities will contribute to optimizing experimental reproducibility in this intense area of research.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by the Regierung von Oberbayern, Kantonales Veterinäramt Zürich, and the Lower Saxony State Office for Consumer Protection and Food Safety (LAVES).

Author Contributions

BS and CE conceived and designed the experiments. CE, DR, PM, MBe, MBa, ES, MS, DS, and AL performed the experiments. CE and PM analyzed the data. JF, AB, and BS contributed the materials and analysis tools. BS coordinated the project and wrote the original draft. All authors reviewed and edited the draft of the manuscript. Correspondence and requests for materials should be addressed to BS.

Funding

This research was supported by the German Research Foundation (DFG) Priority Program SPP1656 (Grant Numbers STE 1971/4-2 and BL 953/52) and the CRC1371 to BS, the German Center of Infection Research (DZIF), and the Center for Gastrointestinal Microbiome Research (CEGIMIR).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Eric Franzosa for advice on the data analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02999/full#supplementary-material

Footnotes

1. ^ https://www.dsmz.de/

References

Anderson, M. J., and Willis, T. J. (2003). Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology 84, 511–525. doi: 10.1890/0012-9658(2003)084[0511:CAOPCA]2.0.CO;2

CrossRef Full Text | Google Scholar

Bangsgaard Bendtsen, K. M., Krych, L., Sorensen, D. B., Pang, W., Nielsen, D. S., Josefsen, K., et al. (2012). Gut microbiota composition is correlated to grid floor induced stress and behavior in the BALB/c mouse. PLoS One 7:e46231. doi: 10.1371/journal.pone.0046231

PubMed Abstract | CrossRef Full Text | Google Scholar

Barroso-Batista, J., Sousa, A., Lourenco, M., Bergman, M. L., Sobral, D., Demengeot, J., et al. (2014). The first steps of adaptation of Escherichia coli to the gut are dominated by soft sweeps. PLoS Genet. 10:e1004182. doi: 10.1371/journal.pgen.1004182

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, N., Kunath, J., Loh, G., and Blaut, M. (2011). Human intestinal microbiota: characterization of a simplified and stable gnotobiotic rat model. Gut Microbes 2, 25–33. doi: 10.4161/gmic.2.1.14651

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergstrom, J. H., Berg, K. A., Rodriguez-Pineiro, A. M., Stecher, B., Johansson, M. E., and Hansson, G. C. (2014). AGR2, an endoplasmic reticulum protein, is secreted into the gastrointestinal mucus. PLoS One 9:e104186. doi: 10.1371/journal.pone.0104186

PubMed Abstract | CrossRef Full Text | Google Scholar

Brugiroux, S., Beutler, M., Pfann, C., Garzetti, D., Ruscheweyh, H. J., Ring, D., et al. (2016). Genome-guided design of a defined mouse microbiota that confers colonization resistance against Salmonella enterica serovar Typhimurium. Nat. Microbiol. 2:16215. doi: 10.1038/nmicrobiol.2016.215

PubMed Abstract | CrossRef Full Text | Google Scholar

Celesk, R. A., Asano, T., and Wagner, M. (1976). The size pH, and redox potential of the cecum in mice associated with various microbial floras. Proc Soc. Exp. Biol. Med. 151, 260–263. doi: 10.3181/00379727-151-39187

PubMed Abstract | CrossRef Full Text | Google Scholar

Clavel, T., Lagkouvardos, I., Blaut, M., and Stecher, B. (2016). The mouse gut microbiome revisited: from complex diversity to model ecosystems. Int. J. Med. Microbiol. 306, 316–327. doi: 10.1016/j.ijmm.2016.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

De Paepe, M., Gaboriau-Routhiau, V., Rainteau, D., Rakotobe, S., Taddei, F., and Cerf-Bensussan, N. (2011). Trade-off between bile resistance and nutritional competence drives Escherichia coli diversification in the mouse gut. PLoS Genet. 7:e1002107. doi: 10.1371/journal.pgen.1002107

PubMed Abstract | CrossRef Full Text | Google Scholar

Deloris Alexander, A., Orcutt, R. P., Henry, J. C., Baker, J. Jr., Bissahoyo, A. C., and Threadgill, D. W. (2006). Quantitative PCR assays for mouse enteric flora reveal strain-dependent differences in composition that are influenced by the microenvironment. Mamm. Genome 17, 1093–1104. doi: 10.1007/s00335-006-0063-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Desai, M. S., Seekatz, A. M., Koropatkin, N. M., Kamada, N., Hickey, C. A., Wolter, M., et al. (2016). A dietary fiber-deprived gut microbiota degrades the colonic mucus barrier and enhances pathogen susceptibility. Cell 167:1339-1353.e21. doi: 10.1016/j.cell.2016.10.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Eppig, J. T., Blake, J. A., Bult, C. J., Kadin, J. A., Richardson, J. E., and Mouse Genome Database, (2015). The Mouse Genome Database (MGD): facilitating mouse as a model for human biology and disease. Nucleic Acids Res. 43, D726–D736.

Google Scholar

Faith, J. J., McNulty, N. P., Rey, F. E., and Gordon, J. I. (2011). Predicting a human gut microbiota’s response to diet in gnotobiotic mice. Science 333, 101–104. doi: 10.1126/science.1206025

PubMed Abstract | CrossRef Full Text | Google Scholar

Falony, G., Joossens, M., Vieira-Silva, S., Wang, J., Darzi, Y., Faust, K., et al. (2016). Population-level analysis of gut microbiome variation. Science 352, 560–564. doi: 10.1126/science.aad3503

PubMed Abstract | CrossRef Full Text | Google Scholar

Forster, S. C., Kumar, N., Anonye, B. O., Almeida, A., Viciani, E., Stares, M. D. M., et al. (2019). A human gut bacterial genome and culture collection for improved metagenomic analyses. Nat. Biotechnol. 37, 186–192. doi: 10.1038/s41587-018-0009-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Franklin, C. L., and Ericsson, A. C. (2017). Microbiota and reproducibility of rodent models. Lab Anim. 46, 114–122. doi: 10.1038/laban.1222

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedman, E. S., Bittinger, K., Esipova, T. V., Hou, L., Chau, L., Jiang, J., et al. (2018). Microbes vs. chemistry in the origin of the anaerobic gut lumen. Proc. Natl. Acad. Sci. U.S.A. 115, 4170–4175. doi: 10.1073/pnas.1718635115

PubMed Abstract | CrossRef Full Text | Google Scholar

Garzetti, D., Brugiroux, S., Bunk, B., Pukall, R., McCoy, K. D., Macpherson, A. J., et al. (2017). High-quality whole-genome sequences of the oligo-mouse-microbiota bacterial community. Genome Announc. 5:e758-17. doi: 10.1128/genomeA.00758-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Geuking, M. B., Cahenzli, J., Lawson, M. A., Ng, D. C., Slack, E., Hapfelmeier, S., et al. (2011). Intestinal bacterial colonization induces mutualistic regulatory T cell responses. Immunity 34, 794–806. doi: 10.1016/j.immuni.2011.03.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes-Neto, J. C., Kittana, H., Mantz, S., Segura Munoz, R. R., Schmaltz, R. J., Bindels, L. B., et al. (2017). A gut pathobiont synergizes with the microbiota to instigate inflammatory disease marked by immunoreactivity against other symbionts but not itself. Sci. Rep. 7:17707. doi: 10.1038/s41598-017-18014-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hadrich, D. (2018). Microbiome research is becoming the key to better understanding health and nutrition. Front. Genet. 9:212. doi: 10.3389/fgene.2018.00212

PubMed Abstract | CrossRef Full Text | Google Scholar

Herp, S., Brugiroux, S., Garzetti, D., Ring, D., Jochum, L. M., Beutler, M., et al. (2019). Mucispirillum schaedleri antagonizes salmonella virulence to protect mice against colitis. Cell Host Microbe 25:681-694.e8. doi: 10.1016/j.chom.2019.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hildebrand, F., Nguyen, T. L., Brinkman, B., Yunta, R. G., Cauwe, B., Vandenabeele, P., et al. (2013). Inflammation-associated enterotypes, host genotype, cage and inter-individual effects drive gut microbiota variation in common laboratory mice. Genome Biol. 14:R4. doi: 10.1186/gb-2013-14-1-r4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hildebrandt, M. A., Hoffmann, C., Sherrill-Mix, S. A., Keilbaugh, S. A., Hamady, M., Chen, Y. Y., et al. (2009). High-fat diet determines the composition of the murine gut microbiome independently of obesity. Gastroenterology 137, e1–e2. doi: 10.1053/j.gastro.2009.08.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Hugon, P., Dufour, J. C., Colson, P., Fournier, P. E., Sallah, K., and Raoult, D. (2015). A comprehensive repertoire of prokaryotic species identified in human beings. Lancet Infect. Dis. 15, 1211–1219. doi: 10.1016/S1473-3099(15)00293-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Ivanov, I. I., Atarashi, K., Manel, N., Brodie, E. L., Shima, T., Karaoz, U., et al. (2009). Induction of intestinal Th17 cells by segmented filamentous bacteria. Cell 139, 485–498. doi: 10.1016/j.cell.2009.09.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Knight, R., Callewaert, C., Marotz, C., Hyde, E. R., Debelius, J. W., McDonald, D., et al. (2017). The microbiome and human biology. Annu. Rev. Genomics Hum. Genet. 18, 65–86.

Google Scholar

Lagkouvardos, I., Lesker, T. R., Hitch, T. C. A., Galvez, E. J. C., Smit, N., Neuhaus, K., et al. (2019). Sequence and cultivation study of Muribaculaceae reveals novel species, host preference, and functional potential of this yet undescribed family. Microbiome 7:28. doi: 10.1186/s40168-019-0637-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Lagkouvardos, I., Pukall, R., Abt, B., Foesel, B. U., Meier-Kolthoff, J. P., Kumar, N., et al. (2016). The mouse intestinal bacterial collection (miBC) provides host-specific insight into cultured diversity and functional potential of the gut microbiota. Nat. Microbiol. 1:16131.

Google Scholar

Laukens, D., Brinkman, B. M., Raes, J., De Vos, M., and Vandenabeele, P. (2015). Heterogeneity of the gut microbiome in mice: guidelines for optimizing experimental design. FEMS Microbiol. Rev. 40, 117–132. doi: 10.1093/femsre/fuv036

PubMed Abstract | CrossRef Full Text | Google Scholar

Leatham, M. P., Stevenson Gauger, E. J., Krogfelt, K. A., Lins, J. J., Haddock, T. L., et al. (2005). Mouse intestine selects nonmotile flhDC mutants of Escherichia coli MG1655 with increased colonizing ability and better utilization of carbon sources. Infect. Immun. 73, 8039–8049. doi: 10.1128/iai.73.12.8039-8049.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, K. S., Palatinszky, M., Pereira, F. C., Nguyen, J., Fernandez, V. I., Mueller, A. J., et al. (2019). An automated Raman-based platform for the sorting of live cells by functional properties. Nat. Microbiol. 4, 1035–1048. doi: 10.1038/s41564-019-0394-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Limenitakis, J. P., Fuhrer, T., Geuking, M. B., Lawson, M. A., Wyss, M., et al. (2015). The outer mucus layer hosts a distinct intestinal microbial niche. Nat. Commun. 6:8292. doi: 10.1038/ncomms9292

PubMed Abstract | CrossRef Full Text | Google Scholar

Macpherson, A. J., and McCoy, K. D. (2015). Standardised animal models of host microbial mutualism. Mucosal Immunol. 8, 476–486. doi: 10.1038/mi.2014.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Mamantopoulos, M., Ronchi, F., McCoy, K. D., and Wullaert, A. (2018). Inflammasomes make the case for littermate-controlled experimental design in studying host-microbiota interactions. Gut Microbes 9, 374–381. doi: 10.1080/19490976.2017.1421888

PubMed Abstract | CrossRef Full Text | Google Scholar

McCoy, K. D., Geuking, M. B., and Ronchi, F. (2017). Gut microbiome standardization in control and experimental mice. Curr. Protoc. Immunol. 117, 2311–23113. doi: 10.1002/cpim.25

PubMed Abstract | CrossRef Full Text | Google Scholar

Ooi, J. H., Waddell, A., Lin, Y. D., Albert, I., Rust, L. T., Holden, V., et al. (2014). Dominant effects of the diet on the microbiome and the local and systemic immune response in mice. PLoS One 9:e86366. doi: 10.1371/journal.pone.0086366

PubMed Abstract | CrossRef Full Text | Google Scholar

Orcutt, R. P., Gianni, F. J., and Judge, R. J. (1987). Development of an “altered Schaedler flora” for NCI gnotobiotic rodents. Microecol. Ther. 17:59.

Google Scholar

Park, S. W., Zhen, G., Verhaeghe, C., Nakagami, Y., Nguyenvu, L. T., Barczak, A. J., et al. (2009). The protein disulfide isomerase AGR2 is essential for production of intestinal mucus. Proc. Natl. Acad. Sci. U.S.A. 106, 6950–6955. doi: 10.1073/pnas.0808722106

PubMed Abstract | CrossRef Full Text | Google Scholar

Rausch, P., Basic, M., Batra, A., Bischoff, S. C., Blaut, M., Clavel, T., et al. (2016). Analysis of factors contributing to variation in the C57BL/6J fecal microbiota across German animal facilities. Int. J. Med. Microbiol. 306, 343–355. doi: 10.1016/j.ijmm.2016.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Reese, A. T., Cho, E. H., Klitzman, B., Nichols, S. P., Wisniewski, N. A., Villa, M. M., et al. (2018). Antibiotic-induced changes in the microbiota disrupt redox dynamics in the gut. eLife 7:e35987. doi: 10.7554/eLife.35987

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, C. D., Klein, H. S., Murphy, K. D., Parthasarathy, R., Guillemin, K., and Bohannan, B. J. M. (2018). Experimental bacterial adaptation to the zebrafish gut reveals a primary role for immigration. PLoS Biol. 16:e2006893. doi: 10.1371/journal.pbio.2006893

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadler, R., Singh, V., Benakis, C., Garzetti, D., Brea, D., Stecher, B., et al. (2017). Microbiota differences between commercial breeders impacts the post-stroke immune response. Brain Behav. Immun. 66, 23–30. doi: 10.1016/j.bbi.2017.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sofi, M. H., Gudi, R., Karumuthil-Melethil, S., Perez, N., Johnson, B. M., and Vasu, C. (2014). pH of drinking water influences the composition of gut microbiome and type 1 diabetes incidence. Diabetes Metab. Res. Rev. 63, 632–644. doi: 10.2337/db13-0981

PubMed Abstract | CrossRef Full Text | Google Scholar

Stecher, B., Chaffron, S., Kappeli, R., Hapfelmeier, S., Freedrich, S., Weber, T. C., et al. (2010). Like will to like: abundances of closely related species can predict susceptibility to intestinal colonization by pathogenic and commensal bacteria. PLoS Pathog. 6:e1000711. doi: 10.1371/journal.ppat.1000711

PubMed Abstract | CrossRef Full Text | Google Scholar

Studer, N., Desharnais, L., Beutler, M., Brugiroux, S., Terrazos, M. A., Menin, L., et al. (2016). Functional intestinal bile acid 7alpha-Dehydroxylation by clostridium scindens associated with protection from clostridium difficile infection in a gnotobiotic mouse model. Front. Cell Infect. Microbiol. 6:191. doi: 10.3389/fcimb.2016.00191

PubMed Abstract | CrossRef Full Text | Google Scholar

Surana, N. K., and Kasper, D. L. (2017). Moving beyond microbiome-wide associations to causal microbe identification. Nature 552, 244–247. doi: 10.1038/nature25019

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, K., Gordon, N., Langley, G., and Higgins, W. (2008). Estimates for worldwide laboratory animal use in 2005. Altern. Lab. Anim. 36, 327–342. doi: 10.1177/026119290803600310

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiemann, S., Smit, N., Roy, U., Lesker, T. R., Galvez, E. J. C., Helmecke, J., et al. (2017). Enhancement of IFNgamma production by distinct commensals ameliorates Salmonella-Induced Disease. Cell Host Microbe 21:682-694.e5. doi: 10.1016/j.chom.2017.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Trexler, P. C., and Reynolds, L. I. (1957). Flexible film apparatus for the rearing and use of germfree animals. Appl. Microbiol. 5, 406–412.

Google Scholar

Uchimura, Y., Fuhrer, T., Li, H., Lawson, M. A., Zimmermann, M., Yilmaz, B., et al. (2018). Antibodies set boundaries limiting microbial metabolite penetration and the resultant mammalian host response. Immunity 49:545-559.e5. doi: 10.1016/j.immuni.2018.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ursell, L. K., Clemente, J. C., Rideout, J. R., Gevers, D., Caporaso, J. G., and Knight, R. (2012). The interpersonal and intrapersonal diversity of human-associated microbiota in key body sites. J. Allergy Clin. Immunol. 129, 1204–1208. doi: 10.1016/j.jaci.2012.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Velazquez, E. M., Nguyen, H., Heasley, K. T., Saechao, C. H., Gil, L. M., Rogers, A. W. L., et al. (2019). Endogenous Enterobacteriaceae underlie variation in susceptibility to Salmonella infection. Nat. Microbiol. 4, 1057–1064. doi: 10.1038/s41564-019-0407-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Linnenbrink, M., Kunzel, S., Fernandes, R., Nadeau, M. J., Rosenstiel, P., et al. (2014). Dietary history contributes to enterotype-like clustering and functional metagenomic content in the intestinal microbiome of wild mice. Proc. Natl. Acad. Sci. U.S.A. 111, E2703–E2710. doi: 10.1073/pnas.1402342111

PubMed Abstract | CrossRef Full Text | Google Scholar

Wullaert, A., Lamkanfi, M., and McCoy, K. D. (2018). Defining the impact of host genotypes on microbiota composition requires meticulous control of experimental variables. Immunity 48, 605–607. doi: 10.1016/j.immuni.2018.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wymore Brand, M., Wannemuehler, M. J., Phillips, G. J., Proctor, A., Overstreet, A. M., Jergens, A. E., et al. (2015). The altered schaedler flora: continued applications of a defined murine microbial community. ILAR J. 56, 169–178. doi: 10.1093/ilar/ilv012

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, L., Feng, Q., Liang, S., Sonne, S. B., Xia, Z., Qiu, X., et al. (2015). A catalog of the mouse gut metagenome. Nat. Biotechnol. 33, 1103–1108. doi: 10.1038/nbt.3353

PubMed Abstract | CrossRef Full Text | Google Scholar