Unveiling the effect of host genetics on gut microbiome using high (FVB) vs. low (C57BL/6) PXR/CAR expressers
Figure 2 shows the alpha diversity (chao1 index), which quantifies the richness of taxa (Fig. 2A), as well as the beta diversity, which quantifies the differences in microbiome among different groups (Fig. 2B), in male and female C57BL/6 and FVB/NJ mice at adolescent and adult ages. FVB/NJ mice had higher microbial richness compared to their C57BL/6 counterparts in both developmental ages and sexes (Fig. 2A). Principle coordinate analysis (PCoA) of the weighted UniFrac measurements showed distinct separations between the microbiomes of FVB/NJ and C57BL/6 mouse strains at both ages and sexes (Fig. 2B).
As shown in Fig. 3A, Lactococcus, which may be associated with anti-inflammatory pathways, was specifically colonized in C57BL/6 mice but not in FVB/NJ mice [56]. Conversely, the pro-inflammatory Acinetobacter [57, 58], and E. dolichum, which is enriched in the Western diet [59], were explicitly colonized in FVB/NJ mice but not in C57BL/6 mice. The most abundant taxa and all differentially regulated taxa in the feces of FVB and C57BL/6 mice are shown in Supplemental Fig. 2 and Supplemental Fig. 3. Interestingly, the abundance of Paraprevotellaceae Prevotella sp. was higher in FVB/NJ mice than in C57BL/6 mice at both ages and sexes, except for adolescent females. To note, Prevotella has been shown to exhibit increased inflammatory properties as demonstrated by enhanced release of inflammatory mediators from immune cells [60], and the abundance of Paraprevotellaceae family has been shown to be higher in IBD [61]. The Turicibacter genus was also higher in both FVB/NJ adolescent males and females as compared to their C57BL/6 counterparts, and this pro-inflammatory taxon has been shown to be positively linked to TNFα levels [62].
As shown in Supplemental Fig. 4, to predict the metagenomic functional content of the gut microbiome in FVB/NJ and C57BL/6 mice, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) was performed based on the QIIME output. Interestingly, between the two mouse strains, males exhibited higher numbers of predicted functional changes than females at both ages. Pathways enriched were classified into 5 main groups, namely cellular processes, environmental information processing, genetic information processing, metabolism, and organismal systems. Metabolism was the largest category for all groups, with 78 pathways enriched in adolescent males, 104 pathways enriched in adult males, 4 pathways enriched in adolescent females, and 6 enriched in adult females. Overall, the predicted pathways were decreased in C57BL/6 mice and increased in FVB/NJ mice.
The pro-inflammatory microbial signature in FVB mice led us to hypothesize that the basal differences in the gut microbiome of the two mouse strains make FVB/NJ mice more prone to pro-inflammatory signaling under basal conditions. To test this hypothesis, cytokines were quantified from C57BL/6 and FVB/NJ mice to assess how strain differences affect inflammation (Fig. 3B). Overall, pro-inflammatory cytokines were higher in feces of FVB/NJ mice. In particular, MCP-1 had the most significant difference in concentration in all four (4) groups, followed by IL-12, and IL-6. TNFα was higher only in adult female FVB/NJ mice. To note, the anti-inflammatory cytokine IL-10, which is known to increase during inflammation to suppress IL-12 [63], was also higher in FVB/NJ than C57BL/6 mice in males and adult females, indicating a compensatory response.
In summary, profound basal differences in gut microbiome composition were observed between FVB/NJ and C57BL/6 mouse strains. Specifically, we showed that under basal conditions, FVB mice, which had higher basal PXR and CAR expression, have higher richness of the gut microbiome, with a pro-inflammatory microbial signature, corresponding to multiple elevated pro-inflammatory cytokines in feces. This basal configuration may prevent invasions by pathogens through enhanced immune surveillance [64].
Understanding the necessity of host nuclear receptors PXR and CAR on the constitutive regulation of the gut microbiome
16S rRNA gene sequencing was conducted on feces collected over a 24-hour period of adolescent and adult aged WT, PXR-null, CAR-null, and PXR-CAR-double null male and female mice (all in the C57BL/6 background, n = 5 per group), to determine the necessity of host nuclear receptors on the composition and function of the gut microbiome. As shown in Fig. 4A, the alpha diversity of PXR and CAR single or double knockout mice tended to have greater richness than WT controls in all groups, and this trend was especially prominent for PXR-CAR double null mice. The only exception was in adult males; CAR-null mice had the highest richness as compared to the other genotypes. As shown in Fig. 4B, regarding the beta diversity, all four genotypes of mice (WT, PXR-null, CAR-null, and PXR-CAR-double null) exhibited distinct separations among their microbial communities at both ages and sexes. The microbial separations indicate that PXR and CAR are essential and unique modulators of the gut microbiome.
Individual variations contributing to these differences are displayed in the heatmap in Supplemental Fig. 5. In total, 63 taxa were significantly different between WT and PXR-null, CAR-null, and PXR-CAR-null mice (Supplemental Fig. 5). Figure 5 shows the microbial compositional changes at the species level of the top 15 most abundant bacteria as quantified by the percentage of OTUs for the four genotypes of mice (note: the 15th category (Other) includes all other taxa summed together). In all exposure groups, the predominant phyla were Bacteroidetes, Firmicutes, and Proteobacteria. The family S24-7 was lower in all PXR-null and PXR-CAR-double null groups compared to WT feces, as well as in adult CAR-null males; this taxon has been shown to be decreased during colitis [65]. Bacteria in the family Helicobacteraceae, which are positively correlated with IBD [66], increased in relative abundance in PXR-null samples in all groups, as well as both male CAR-null age groups, and CAR-null adult female samples. The pro-inflammatory helicobacter was also higher in the absence of PXR and CAR, whereas the anti-inflammatory A. muciniphila was lower in the absence of PXR and CAR (supplemental Fig. 5). Allobaculum sp., which is a microbial biomarker that is inversely associated with obesity [67], was higher in all PXR-CAR-null groups, except for adult females. Most notably, as shown in Fig. 6, Lactobacillus, which is known to carry bile salt hydrolase (BSH) activity for BA-deconjugation [68, 69], was markedly up-regulated in PXR/CAR-double null mice, and this bacteria also tended to be higher in single receptor gene null mice. Anaerostipes sp. decreased in relative abundance in all nuclear receptor-deficient mice of both ages and sexes. This microbe has been shown to have anti-inflammatory effect through producing butyrate [70]. In contrast, the pro-inflammatory Sutterella [71], which is enriched in IBD, was up-regulated in PXR- and CAR-null mice and further up-regulated in the double knockout mice. To note, Sutterella is also known to be bile-resistant [72].
The functional predictions of the gut microbiome in the absence of PXR and CAR, as shown in Supplemental Fig. 6, were done using PICRUSt. Females exhibited markedly more functional content predictions than their male counterparts at both developmental ages sampled. In both female ages, absence of PXR caused an increase in the greatest number of pathways, with the opposite being true for CAR-null samples. The PXR-CAR-null group appeared to respond to combinatorial effects of the absence of PXR or CAR, as no change in the abundance of these pathways was seen for these females. The opposite is true for the samples from male mice in that PXR-CAR-null in both ages, and CAR-null in adult males had a marked increase in pathways compared to PXR-null and WT samples. Therefore, sex and host nuclear receptor status affect the function of the gut microbiome in mice.
Because there was a marked increase in the BA-deconjugating Lactobacillus in the PXR/CAR-double null mice (and to a lesser extent tended to be higher in the single null mice), we hypothesized that this would result in a reduction in conjugated BAs in feces of the PXR/CAR-null mice. To test this hypothesis, we conducted LC-MS based targeted metabolomics of all major BAs in mice (Fig. 7 and Supplemental Figs. 7–10). As expected, as shown in Fig. 7, the most abundant conjugated BAs in mice, namely T-CA, T-αMCA, T-βMCA, and T-ωMCA, tended to be lower in feces of the PXR/CAR single and double null mice, and the trend was most predominant in the PXR/CAR double null mice. Other minor T-conjugated secondary BAs, such as T-HDCA T-UDCA, were also down-regulated in the knockout mice in an age and sex-dependent manner (Supplemental Figs. 7–10). We also observed an increase in the major unconjugated secondary BA LCA in adolescent PXR/CAR double null mice of both sexes (Supplemental Figs. 7 and 9). However, several other unconjugated secondary BAs were lower in adult PXR/CAR null mice. The inconsistency of the regulatory patterns of the unconjugated secondary BAs patterns in feces may be due to the hydrophobic nature of these BA species, leading to increased body retention and lack of excretion into the fecal compartment. Because unconjugated BAs are more pro-inflammatory than conjugated BAs [9], enhanced BA de-conjugation may prime the host for inflammation related diseases.
Comparison between mouse and human PXR in regulating the gut microbiome
To compare the role of mouse and human PXR genes on the composition and function of the gut microbiome, 16S rRNA gene sequencing was conducted on feces collected over a 24-hour period of adolescent and adult aged wild type (WT) and hPXR-TG male and female mice (n = 5 per group). The overall microbial richness was similar between WT and hPXR-TG mice, except for adolescent males, as evidenced by a higher microbial richness in hPXR-TG mice (Fig. 8A). As shown in Fig. 8B, in all 4 comparisons, hPXR-TG and WT mice exhibited distinct separations between their microbial communities as measured by beta diversity (weighted uniFrac).
In total, 27 taxa with an abundance above 0.00005 were significantly different between WT and hPXR-TG mice (Supplemental Fig. 11). Figure 9 shows the top 15 species level compositional changes of the gut microbiome as quantified by % OTUs. In all groups, the predominant phyla were Bacteroidetes, Firmicutes, Tenericutes, and Verrucomicrobia, plus Proteobacteria in both adult groups and adolescent males. The most predominant difference was a relative increase in Prevotella in hPXR-TG mice.
As shown in Supplemental Fig. 12, PICRUSt was used to predict the effect of host nuclear receptor genotype on the metagenomic functional content of the gut microbiome in mice. Adult male hPXR-TG mice were the only group with altered functional pathways. These 21 pathways are: lysine degradation, styrene degradation, aminobenzoate degradation, sulfur metabolism, limonene and pinene degradation, atrazine degradation, biosynthesis of siderophore group nonribosomal peptides, biosynthesis of unsaturated fatty acids, arachidonic acid metabolism, retinol metabolism, chlorocyclohexane and chlorobenzene degradation, bacterial invasion of epithelial cells, proximal tubule bicarbonate reclamation, metabolism of cofactors and vitamins, ion channels, flavonoid biosynthesis, cytochrome P450-mediated xenobiotic metabolism, caprolactam degradation, tryptophan metabolism, and cell motility and secretion. Notably, as a group, all of these pathways decreased in the adult male hPXR-TG mice, with some individual variation. Therefore, species specificity of PXR (i.e. mPXR vs. hPXR) affects the predicted functional differences of the gut microbiome, in a sex- and age-specific manner.
As shown in Fig. 10, BA profiles were different between WT and hPXR-TG mice. Specifically, the largest increase in relative concentration occurred with the secondary BA DCA in hPXR-TG adult male mice, and this trend was seen in the other groups as well. The primary BA CA, which is the precursor of DCA, also tended to be higher in hPXR-TG mice at all 4 comparisons. Conversely, the major secondary BA T-ωMCA tended to be lower in hPXR-TG mice in all 4 comparisons, whereas its unconjugated form ωMCA also tended to be lower in female hPXR-TG mice. Regarding other minor BAs in feces, HDCA was lower in both hPXR-TG adult mice groups, with this trend observed in the adolescent mice, and UDCA was higher in adolescent male and adult female hPXR-TG mice, with this trend continuing in the other groups as well.