Fasting and refeeding triggers specific changes in bile acid profiles and gut microbiota

Article type: Research Article

Keywords: bile acid, biosynthesis, fasting, gut microbiota, reabsorption, refeeding, 胆汁酸, 肠道菌群, 禁食, 复食, 生物合成, 再吸收

Authors: Yi Zhang ⚬, Hongyan Qi, Long Wang, Chunyan Hu ⚬, Aibo Gao, Qihan Wu, Qiaoling Wang, Huibin Lin, Banru Chen, Xingyu Wang, Shuangyuan Wang, Hong Lin, Weiqing Wang ⚬, Yufang Bi ⚬, Jiqiu Wang, Jieli Lu ⚬, Ruixin Liu

Affiliations: Department of Endocrine and Metabolic Diseases, Shanghai Institute of Endocrine and Metabolic Diseases Ruijin Hospital, Shanghai Jiao Tong University School of Medicine Shanghai China; Shanghai National Clinical Research Center for Metabolic Diseases, Key Laboratory for Endocrine and Metabolic Diseases of the National Health Commission of the PR China, Shanghai Key Laboratory for Endocrine Tumor, State Key Laboratory of Medical Genomics Ruijin Hospital, Shanghai Jiao Tong University School of Medicine Shanghai China

License: © 2023 The Authors. Journal of Diabetes published by Ruijin Hospital, Shanghai JiaoTong University School of Medicine and John Wiley & Sons Australia, Ltd. CC BY 4.0 This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Article links: DOI: 10.1111/1753-0407.13356  |  PubMed: 36682739  |  PMC: PMC9934961

Relevance: Moderate: mentioned 3+ times in text

Full text: PDF (9.2 MB)

INTRODUCTION

Nutrient homeostasis, which is maintained by liver during fasting and feeding, is critical for health.ref. jdb13356-bib-0001 Imbalance in nutrient homeostasis is a common trigger for the progression of metabolic diseases including obesity,ref. jdb13356-bib-0002 type 2 diabetes mellitus (T2DM),ref. jdb13356-bib-0002 nonalcoholic fatty liver disease,ref. jdb13356-bib-0003 metabolic syndrome,ref. jdb13356-bib-0004 and cardiovascular disease.ref. jdb13356-bib-0005

Bile acids (BAs) are necessary for nutrient absorption, emulsification, and transport.ref. jdb13356-bib-0006 In hepatocytes, BAs are derived from cholesterol catabolism, a process that involves the activation of two major pathways and at least 17 liver enzymes.ref. jdb13356-bib-0006 The “classical pathway” is initiated by cholesterol 7α‐hydroxylase (CYP7A1), followed by the sterol 12α‐hydroxylase (CYP8B1) and enzymatic action of sterol 27‐hydroxylase (CYP27A1) to produce unconjugated primary BAs (PBAs), including chenodeoxycholic acid (CDCA) and cholic acid (CA). The “alternative pathway” is initiated by CYP27A1 and produces CDCA via oxysterol 7α‐hydroxylase (CYP7B1). In addition to CDCA and CA, mice also produce α muricholic acid (MCA), β MCA, and ursodeoxycholic acid (UDCA).ref. jdb13356-bib-0007 In hepatocytes, most unconjugated PBAs are conjugated to either taurine (95% BAs are conjugated to taurine in mice) or glycine (predominant in human beings) by bile acid: CoA synthetase (BACS) and bile acid‐CoA: amino acid N‐acyltransferase (BAAT) and produce the conjugated PBAs prior to their secretion into intestine.ref. jdb13356-bib-0008 In the gut, microbial enzymes from bacteria metabolize PBAs. Conjugated PBAs are deconjugated via bile salt hydrolases (BSH), 7α‐dehydroxylated, C‐6 epimerized and 6β‐epimerization to form unconjugated secondary BAs (SBAs), including lithocholic acid (LCA), deoxycholic acid (DCA), and ω MCA.ref. jdb13356-bib-0009 Unconjugated SBAs are conjugated to taurine or glycine and then produce conjugated SBAs in liver.ref. jdb13356-bib-0007 After food stimulation, BAs are released into the duodenum, and then both conjugated and unconjugated BAs are reabsorbed in the ileum by the apical sodium dependent BA transporter (ASBT)ref. jdb13356-bib-0010 and organic solute transporter‐α/β (OST‐α/β).ref. jdb13356-bib-0008 BAs are recirculated to the liver through the portal blood, a process known as enterohepatic circulation.ref. jdb13356-bib-0007 Recently, BAs have also emerged as important signaling molecules in the control of energyref. jdb13356-bib-0011 and glucose metabolism.ref. jdb13356-bib-0012 BAs regulate gut microbiome configuration,ref. jdb13356-bib-0007 glycogen synthesis,ref. jdb13356-bib-0013 hepatic gluconeogenesis,ref. jdb13356-bib-0014 intestinal incretin secretion,ref. jdb13356-bib-0015 and energy expenditureref. jdb13356-bib-0016 by binding to the nuclear hormone farnesoid X receptorref. jdb13356-bib-0017 and takeda G protein receptor 5ref. jdb13356-bib-0018 in multiple organs. It has been found that in population‐based studies, impaired BAs metabolism is associated with the progression of metabolic disorder including T2DM,ref. jdb13356-bib-0019 obesity,ref. jdb13356-bib-0020 insulin resistance,ref. jdb13356-bib-0021 and polycystic ovary syndrome.ref. jdb13356-bib-0022 However, how nutrient availability systematically affects BA composition remain unclear. Maintenance of BA homeostasis based on nutrient supply may be a potential therapy to relieve metabolic disorders.

Growing evidence support the crosstalk between host health and gut microbiota. Alterations in gut microbiota composition are associated with various metabolic functions, including energy absorption, glucose homeostasis, and intestinal barrier homeostasis.ref. jdb13356-bib-0023, ref. jdb13356-bib-0024, ref. jdb13356-bib-0025 Several studies have shown that the abundance of Firmicutes, Bacteroidetes, and butyrate‐producing bacteria such as Akkermansia is associated with metabolic disturbances, including T2DM and obesity.ref. jdb13356-bib-0023, ref. jdb13356-bib-0026, ref. jdb13356-bib-0027, ref. jdb13356-bib-0028 Supplementation with Akkermansia improves gut‐barrier function and metabolic endotoxemia, thereby improving the systemic metabolism.ref. jdb13356-bib-0026, ref. jdb13356-bib-0029, ref. jdb13356-bib-0030 Diet pattern (the type, amount, and timing of the meals) is a key factor that restructures the gut microbiota.ref. jdb13356-bib-0031, ref. jdb13356-bib-0032, ref. jdb13356-bib-0033 Sustained caloric restriction (CR) during very‐low‐calorie diet (VLCD),ref. jdb13356-bib-0034 intermittent fasting (IF),ref. jdb13356-bib-0035 and long‐term fastingref. jdb13356-bib-0036 have an impact on the gut microbiome. Schwartzenberg et al. found weight loss caused by VLCD was associated with impaired nutrient absorption and Clostridioides difficile enrichment, and enrichment in Clostridium difficile was associated with decreased BAs levels.ref. jdb13356-bib-0034 Shi et al. demonstrated BAs were potential mediators in the microbe‐host interaction involved in IF intervention.ref. jdb13356-bib-0035 These studies emphasized the importance of diet‐microbiome‐BA interactions in regulating host metabolic homeostasis. However, the consequences of rapid changes in gut microbiota and BAs in response to nutrient supply for host pathophysiology have not been fully explored.

In this study, we contributed a systematic evaluation of changes of BA profiles (in plasma, feces, and liver), gene expression related to BA biosynthesis and intestinal BA reabsorption, as well as gut microbiota in mice fed ad libitum, fasted for 24 h, fasted for 24 h and then fed ad libitum for 24 h. We identified specific and rapid changes in BA components, microbial taxa, BA biosynthesis, and intestinal reabsorption‐related gene expression in response to nutrient availability. Of note, Akkermansia was increased in the fasting state and negatively correlated with plasma unconjugated PBAs, plasma unconjugated SBAs, and glucose levels, whereas it was positively correlated with plasma conjugated SBAs, fecal unconjugated PBAs, and fecal unconjugated SBAs. Hence, our research characterized the BA profiles, gut microbiota, gene expression related to BA biosynthesis, and intestinal reabsorption in fasting and refeeding mice, identifying the rapid changes in response to fasting‐feeding cycle. Our results thus highlighted the important and rapid effects of nutrient supply on BA profiles and gut microbiota.

MATERIALS AND METHODS

Animal care and experimental procedures

Specific‐pathogen‐free‐grade 6‐week‐old male C57BL/6 mice (n = 36) purchased from Shanghai Laboratory Animal Center were housed in groups of four mice per cage in a stable environment with free access to water. Mice were reared without any intervention for 1 week to acclimate to the new environment and then randomly divided into three groups: (a) mice in normal chow diet (NCD) fed ad libitum, (b) mice in NCD starved for 24 h, and (c) mice in NCD starved for 24 h and then fed ad libitum for 24 h (Figure 1A). The NCD was composed of 10% fat, 20% protein, and 70% carbohydrate. Mice were anesthetized with 10% chloral hydrate at the end of the experiment, and blood from the orbital plexus was collected in tubes. Blood samples were centrifuged at 3000 rpm, 4°C for 15 min. Plasma was separated and stored at −80°C for further biochemical testing. Using a blood glucose meter (LifeScan), blood glucose levels in mice were measured. The liver and ileum tissues of each mouse were dissected precisely and stored immediately at −80°C for subsequent analysis. Feces in the colon were collected precisely and stored immediately at −80°C for BA profiling and 16S rRNA gene sequencing. All these procedures were approved by the Animal Care Committee of Shanghai Jiao Tong University School of Medicine.

Bile acids analysis

BAs in plasma, feces, and liver tissue were measured. Detailed analytical methods are provided in the Supplementary Materials. BAs were quantified using ultraperformance liquid chromatography coupled to tandem mass spectrometry (UPLC‐MS/MS) in combination with multiple reactions monitoring methods on the ACQUITY UPLC system (Waters). All calibration standard solutions were mixed at appropriate concentrations and analyzed every 10 samples for quality control.

BAs were divided into four groups according to the degree of conjugation,ref. jdb13356-bib-0008 including (a) unconjugated PBAs: UDCA, α MCA, β MCA, CA, and CDCA; (b) conjugated PBAs: tauroursodeoxycholic acid (TUDCA), tauro‐β‐muricholic acid (T‐β‐MCA), tauro‐α‐muricholic acid (T‐α‐MCA), glycochenodeoxycholic acid (GCDCA), taurochenodeoxycholic acid (TCDCA), glycocholic acid (GCA), and taurocholic acid (TCA); (c) unconjugated SBAs: DCA, LCA, hyocholic acid (HCA), and ω MCA; and (d) conjugated SBAs: taurolithocholic acid (TLCA), taurohyocholic acid (THCA), taurodeoxycholic acid (TDCA), and tauro‐ω‐muricholic acid (T‐ω‐MCA).

Sequencing the 16S rRNA gene and metataxonomic analysis

Using the E.Z.N.A.® Soil DNA Kit (Omega Bio‐tek), total microbial genomic DNA was extracted from feces. DNA quality and concentration were determined by 1.0% agarose gel electrophoresis and NanoDrop ND‐2000 spectrophotometer (Thermo Scientific Inc.). DNA was stored at −80°C for subsequent analysis. The V3‐V4 hypervariable region of the bacterial 16S rRNA gene was amplified by ABI GeneAmp® 9700 polymerase chain reaction (PCR) thermocycler (ABI) with primer pairs 338F (5’‐ACTCCTACGGGAGGCAGCAG‐3′) and 806R (5’‐GGACTACHVGGGTWTCTAAT‐3′).ref. jdb13356-bib-0037 PCR reaction mix included 10 ng template DNA, 0.8 μl each primer (5 μM), 2 μl 2.5 mM dNTP, 0.2 μl BSA, 4 μl 5 × Fast Pfu buffer, 0.4 μl Fast Pfu polymerase, and ddH₂O (20 μl final volume). All samples were amplified in triplicate. PCR amplification conditions were as follows: initial denaturation at 95°C for 3 min, amplification for 27 cycles (each cycle consists of denaturation at 95°C for 30 s, annealing at 55°C for 30 s and extension at 72°C for 45 s), single extension at 72°C for 10 min, and end at 4°C. PCR products were extracted from 2% agarose gel, then purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences) and quantified using Quantus™ Fluorometer (Promega). 16S rRNA gene sequencing data were processed using the Illumina MiSeq PE300/NovaSeq PE250 platform, details are provided in the Supplementary Materials.

RNA isolation and real‐time polymerase chain reaction

Total RNA was extracted from liver and ileum tissues using Eastep Super Total RNA Extraction Kit (Promega) and 1 μg of RNA was transcribed following the manufacturer’s instruction by using the Reverse Transcription System Kit (Promega), which produced high‐quality complementary DNA. A real‐time PCR was performed on LightCyclerH 480 (Roche) using SYBR Green II Master (Takara). All mRNA quantification data were normalized to the housekeeping 36B4 gene. All samples were run in duplicate in a single 384‐well reaction plate, and data of the reaction were collected and analyzed according to the ∆∆CT method. Primers used in this study are listed in the Supplementary Materials.

Statistical analysis

Results are expressed as mean ± SEM. BA data were analyzed using one‐way analysis of variance (ANOVA) test followed by Bonferroni’s post hoc test. Significant differences between the two groups were assessed by unpaired two‐tailed Student’s t test. All of these data were tested for the normality before using Student’s t test. For nonnormal distributions, nonparametric statistical analysis such as Mann–Whitney U‐test was used. GraphPad prism version 5.01 (GraphPad software) and SPSS Statistics version 23.0 (IBM SPSS) were used for statistical analyses. Using Spearman’s correlation analysis, correlations between BA levels, glucose level, and abundance of gut microbiota were performed, and the post hoc correction was performed using the false discovery rate method. Differences between experimental groups were considered significant at p < .05. The Majorbio Cloud platform (https://cloud.majorbio.com) was used to carry out bioinformatic analysis of gut microbiota. See the Supplementary Materials for detailed analytical methods.

RESULTS

Fasting and refeeding modify bile acid composition

We first measured BA levels in plasma, feces, and liver tissues from the libitum fed mice, fasting mice, and refeeding mice (Figure 1A) using UPLC‐MS/MS. The baseline fecal BA levels (Figure S1) were comparable among the three groups. We found that the mice fasted for 24 h showed a significant decrease in plasma unconjugated PBAs (eg, α MCA, β MCA, UDCA, CA, and CDCA) and unconjugated SBAs (eg, DCA and HCA, despite ꙍ MCA and LCA showing a decreased trend) and a significant increase in conjugated SBAs (eg, TDCA), when compared with the libitum fed mice. It is noteworthy that these changes were completely recovered in mice refed for another 24 h. Nevertheless, plasma conjugated PBAs did not change in fasting‐refeeding cycle (Figure 1B,C).

Next, we found that fecal unconjugated PBAs (eg, α MCA, β MCA, UDCA, and CA), unconjugated SBAs (eg, ꙍ MCA, DCA, LCA, and HCA), and conjugated SBAs (T‐ꙍ‐MCA, TDCA, and TLCA) were significantly elevated in the fasting mice compared with the libitum fed mice. Although the level of conjugated PBAs subgroup did not change in fasting mice, T‐β‐MCA, TUDCA, and GCA were significantly increased. Similarly, these changes were overall recovered in mice after refeeding, back to the levels of those in the libitum fed mice (Figure 1D,E). Of note, unconjugated PBAs and unconjugated SBAs showed an opposite changing pattern whereas conjugated SBAs shared similar changing pattern in feces and plasma in response to the fasting‐refeeding cycle.

In liver, compared with the libitum fed mice, the levels of unconjugated PBAs subgroup and conjugated PBAs subgroup did not change in the fasting mice, whereas changes in individual BAs were observed. In the fasting mice, level of UDCA in unconjugated PBAs were significantly increased, whereas TUDCA, T‐α‐MCA, and T‐β‐MCA in conjugated PBAs were significantly decreased. The alterations of these BAs were reversed in the refeeding mice (Figure 1F,G). Unconjugated SBAs (eg, ꙍ MCA and DCA) were significantly increased, and conjugated SBAs (eg, T‐ꙍ‐MCA and TLCA) were significantly decreased in the fasting mice. These alterations were also recovered in mice after refeeding for 24 h (Figure 1F,G). Taken together, these results demonstrated a rapid and specific change in individual BA species and BA subgroups in plasma, feces, and liver in response to fasting‐refeeding status.

To determine which fecal BAs associated the most with plasma BAs under fasting‐refeeding condition, we further examined the correlation between BAs in plasma and feces (Figure S2). Levels of β MCA, CA, UDCA, α MCA, DCA, ꙍ MCA, HCA, and LCA in plasma were negatively correlated with their levels in feces, whereas levels of T‐ꙍ‐MCA, TDCA, and TLCA in plasma were positively correlated with their levels in feces (Figure S2).

The mRNA levels of genes related to bile acid biosynthesis and intestinal reabsorption under fasting‐refeeding cycle

Liver biosynthesis and intestinal reabsorption play a critical role in the regulation of enterohepatic circulation of BAs, BA pool size, and BA composition.ref. jdb13356-bib-0006 We next analyzed the mRNA expression levels of BA biosynthesis‐related genes in liver and BA reabsorption‐related genes in ileum. Figure 2A summarized the BA biosynthesis pathways in liver. In hepatocytes, “classical pathway” involves CYP7A1, 3β‐hydroxy‐∆⁵‐C₂₇‐steroid dehydrogenase (HSD3B7), CYP8B1, aldo‐keto reductase family 1 member D1 (AKR1D1) and CYP27A1, and “alternative pathway” involves CYP27A1 and CYP7B1. A species‐specific sterol 6β‐hydroxylase (CYP2C70) is required for the conversion of MCAs from CDCA in mouse liver. Unconjugated PBAs are subsequently conjugated to glycine or taurine by BAAT and BACS (Figure 2A).

In liver, the expression of CYP7A1, CYP7B1, CYP2C70, and AKR1D1 were significantly decreased in the fasting mice compared with the libitum fed mice, whereas the expression of CYP8B1 and BAAT were significantly increased, and all of these changes were completely recovered in the refeeding mice. The expression of CYP27A1 and BACS remained unchanged in the fasting‐refeeding cycle (Figure 2B). In ileum tissues, mice fasted for 24 h showed a significant decrease in the expression of OST‐α, OST‐β, and ASBT (intestinal reabsorption‐related genes) when compared with the libitum fed mice. The alterations in the expression of these genes were reversed in the refeeding mice (Figure 2C). These results together demonstrated that BAs biosynthesis in liver and intestinal reabsorption in ileum were suppressed under fasting status, whereas refeeding ameliorated this effect.

Gut microbiota features in fasting and refeeding

Using 16S rRNA gene sequencing, the characteristics of gut microbiota in libitum fed mice, fasting mice, and refeeding mice were analyzed. We first examined the baseline microbial features among the three groups and found that bacterial diversity (Shannon index, Figure S1) and the community structure (Figure S1) of gut microbiota were comparable among three groups at baseline. We further examined the gut microbiota of three groups with different nutrient status. Significantly, principal coordinate analysis (PCoA) indicated a distinct gut microbiota structure among the three groups (p = .001) (Figure 3A). Further, fasting mice had a higher bacterial diversity than libitum fed mice (Shannon index, p = .029, one‐way ANOVA test followed by Bonferroni’s post hoc test), whreeas in refeeding mice, bacterial diversity was again decreased (Shannon index, p for refeeding mice and libitum fed mice = 1.011 × 10⁻⁷, p for refeeding mice & fasting mice = 2.98 × 10⁻⁹, one‐way ANOVA test followed by Bonferroni’s post hoc test) (Figure 3B). Venn diagram showed that 326 operational taxonomic units (OTUs) were common to the three groups. 12 OTUs were unique for the libitum fed mice, 19 OTUs were unique for the fasting mice, whereas 45 OTUs were unique for the refeeding mice (Figure 3C). A list of the unique enriched bacterial genera among three groups is shown in the Supplementary Material (Figure S3). Overall, these results indicated that gut microbiota changed rapidly in response to food availability.

We used the linear discriminant analysis effect size (LEfSe) tool to identify specific communities within different groups (Figure 3D).ref. jdb13356-bib-0038 In the libitum fed mice, six groups of bacteria were detected to be significantly enriched, namely Actinobacteria (the class, including its order Bifidobacteriales, its family Bifidobacteriaceae, and its genus Bifidobacterium), RF39 (the order, including its family and genus), Erysipelotrichales (the order, including its family Erysipelotrichaceae), Clostridia_UCG‐014 (the order, including its family and genus), Christensenellales (the order, including its family Christensenellaceae and its genus Christensenellaceae_R‐7_group), and Eubacteriales (the order, including its family and genus) (Figure 3D). In the fasting mice, eight groups of bacteria were significantly enriched, namely Bacteroidetes (the phylum, including its class Bacteroidia, its order Bacteroidales, its family Muribaculaceae and Bacteroidaceae, its genus Bacteroides and Muribaculum), Verrucomicrobiota (the phylum, including its class Verrucomicrobiae, its order Verrucomicrobiales, its family Akkermansiaceae, and its genus Akkermansia), Vampirivibrionia (the class, including its order, family, and genus), Proteobacteria (the phylum, including its class Gammaproteobacteria, its order Burkholderiales, its family Sutterellaceae, and its genus Parasutterella), Peptococcales (the order, including its family and genus), Clostridia (from order to genus), Defluviitaleaceae (the family, including its genus), and Oscillospirales (the order, including its family Oscillospiraceae, UCG‐010 and Eubacterium_coprostanoligenes_group, its genus Oscillibacter, NK4A214_group, Colidextribacter, Candidatus_Soleaferrea and Anaerotruncus) (Figure 3D). In the refeeding mice, Firmicutes (the phylum, including its class Bacilli, its order Lactobacillales, its family Lactobacillaceae, and its genus Lactobacillus) were significantly enriched (Figure 3D).

Figures 3E depicted the microbial composition of each group at genus level. Of note, libitum fed mice and refeeding mice were dominated by the Lactobacillus. Fasting for 24 h increased the relative abundance of norank_f_Muribaculaceae and Akkermansia and decreased the relative abundance of Lactobacillus and Bifidobacterium, compared with the libitum fed mice. These alterations were reversed in mice after refeeding for 24 h (Figure 3E). The top 15 abundant genera that differed among three groups were shown in Figure 3F. Of note, the relative abundances of seven dominant genera differed significantly between libitum fed mice and fasting mice. Lactobacillus (29.49% versus 3.44%, p < .001) and Bifidobacterium (5.04% versus 0.09%, p < .001) were significantly decreased, whereas norank_f_Muribaculaceae (13.00% versus 35.83%, p < .001), Akkermansia (6.25%versus 23.43%, p < 0.001), Muribaculum (0.59% versus 3.24%, p < .001), Parabacteroides (0.48% versus 2.14%, p < .05), and norank_f_Eubacterium_coprostanoligenes_group (0.11% versus 2.69%, p < .001) were significantly increased in fasting mice, compared with libitum fed mice. These alterations were also reversed in mice after refeeding for 24 h (Figure 3F). Taken together, these results illustrated that rapid change in nutrient supply profoundly affected the composition of gut microbiota in mice.

Associations of BAs and plasma glucose with gut microbiota in fasting and refeeding

These results demonstrated that the fasting‐refeeding cycle changed gut microbial community and BA profiles. The effect of nutrient supply on microbial communities may potentially affect BAs and glucose characteristics. Therefore, we further investigated whether gut microbial community, BAs, and glucose are correlated during the fasting‐refeeding cycle (Figure 4). Spearman’s correlation analysis and redundancy analysis/canonical correlation analysis (RDA/CCA) revealed that the levels of plasma BAs (Figure 4A,B), fecal BAs (Figure 4C,D), and liver BAs (Figure 4E,F) were correlated with the top 20 abundant gut microbes at genus level. Lactobacillus and Bifidobacterium were positively correlated with plasma unconjugated PBAs (β MCA, α MCA, UDCA, CDCA, and CA) and unconjugated SBA levels (LCA, DCA, HCA, and ꙍ MCA), whereas norank_f_Muribaculaceae and Akkermansia were negatively correlated with unconjugated PBAs and unconjugated SBA levels. Akkermansia showed different correlation patterns with conjugated SBA levels. It was positively correlated with the levels of T‐ꙍ‐MCA, TDCA, and TLCA and negatively correlated with THCA levels (Figure 4A). More importantly, plasma glucose level was positively correlated with Lactobacillus and Bifidobacterium, and it was negatively correlated with norank_f_Muribaculaceae and Akkermansia (Figure 4A). RDA/CCA revealed that the microbial community structure was closely related to BAs and glucose. After removal of the redundant variables 11 plasma BAs and blood glucose were chosen for RDA/CCA. Plasma TDCA (p = .001), GCA (p = .043), and glucose (p = .001) were significantly associated with the gut microbiota structure (Figure 4B). In feces, unconjugated PBAs levels were negatively correlated with Lactobacillus and positively correlated with norank_f_Muribaculaceae, Akkermansia and Bifidobacterium. Unconjugated SBAs were negatively correlated with Lactobacillus and positively correlated with norank_f_Muribaculaceae, Akkermansia (Figure 4C). After removal of the redundant variables, 11 fecal BAs were chosen for RDA/CCA. Fecal GCDCA (p = .001), DCA (p = .002), and LCA (p = .001) were significantly associated with the gut microbiota structure (Figure 4D). In liver, α MCA and UDCA were negatively correlated with Lactobacillus, and they were positively correlated with norank_f_Muribaculaceae. The level of UDCA was positively correlated with Akkermansia. Akkermansia was negatively correlated with T‐β‐MCA, T‐α‐MCA, and TUDCA and positively correlated with UDCA, DCA and ꙍ MCA (Figure 4E). After removal of the redundant variables, 11 liver BAs were chosen for RDA/CCA. Liver TUDCA (p = .018), TCDCA (p = .018), TDCA (p = .002), and TLCA (p = .001) were significantly associated with the gut microbiota structure (Figure 4F). These results demonstrated the important and rapid effects of nutrient supply on the correlations between BA profile, glucose, and gut microbiota.

DISCUSSION

Previous studies support the crosstalk among gut microbiota, BAs, and nutrient supply in host health homeostasis, and changes of gut microbiota and BA profiles are involved in the progression of metabolic diseases such as T2DM and obesity.ref. jdb13356-bib-0006, ref. jdb13356-bib-0007 Thus, elucidating the specific changes of gut microbiota and BA profile in response to nutrient supply will provide new targets for preventing and treating metabolic disorders. In the current study, we characterized BA profile, gut microbiota and BA biosynthesis, and intestinal BA reabsorption‐related gene expression in response to rapid change of nutrient supply.

To the best of our knowledge, for the first time, we contributed a systematic assessment of BA profiles (in plasma, feces, and liver), gut microbiota and BA biosynthesis and intestinal BA reabsorption‐related genes expression changes triggered by fasting‐refeeding cycle in mice. We found that in the fasting mice, plasma unconjugated PBAs and plasma unconjugated SBAs decreased whereas plasma conjugated SBAs, fecal unconjugated PBAs, fecal unconjugated SBAs, fecal conjugated SBAs, and liver unconjugated SBAs were increased, compared with the libitum fed mice. All these changes in BA composition were completely recovered after refeeding. BA biosynthesis in liver and intestinal reabsorption were suppressed under fasting status, whreeas refeeding ameliorated the decrease in biosynthesis and reabsorption of BAs due to fasting. Correspondingly, comparing with the gut microbiota in the libitum fed mice, Akkermansia, Parabacteroides, Muribaculum, norank_f_Muribaculaceae, and norank_f_Eubacterium_coprostanoligenes_group were largely increased in the fasting mice whereas Lactobacillus and Bifidobacterium were decreased. Also, these microbial changes were recovered under refeeding. More importantly, Akkermansia was negatively correlated with plasma unconjugated PBAs, plasma unconjugated SBAs, and glucose levels, whereas it was positively correlated with plasma conjugated SBAs, fecal unconjugated PBAs, and fecal unconjugated SBAs. Glucose contributed the most to the changes of gut bacterial community structure, followed by plasma GCA and TDCA. These results highlighted the important and rapid effects of nutrient supply on BA profile, gut microbiota, BA biosynthesis, and intestinal reabsorption.

In our results, compared with the ad libitum fed mice, the levels of plasma unconjugated PBAs and unconjugated SBAs were decreased in the fasting mice, whereas fecal unconjugated PBAs, unconjugated SBAs, and conjugated SBAs were increased in the fasting mice. These BAs showed an opposite changing pattern in feces and plasma in response to the fasting‐refeeding cycle. This might be because of the alterations in hepatic biosynthesis and intestinal reabsorption in the fasting mice, which play a critical role in the regulation of enterohepatic circulation of BAs, BA pool size, and BA composition.ref. jdb13356-bib-0001 Consistently, we found that overnight fasting repressed, whereas refeeding induced the expression of CYP7A1 (the rate‐limiting enzyme in the classic pathway of BA synthesis, determining the rate of BA synthesis), CYP7B1 (involved in alternative pathway), AKR1D1 (involved in both classic pathway and alternative pathway), and CYP2C70 (CDCA to α MCA conversion) in liver. The reduction in BA biosynthesis gene expression in liver may contribute to the reduced plasma unconjugated PBAs and unconjugated SBAs after fasting. On the other hand, we detected a decreased expression of intestinal reabsorption‐related genes (OST‐α, OST‐β, and ASBT) in the ileum of fasting mice, which represents an attenuated enterohepatic circulation after fasting and may result in more BA accumulation in intestine and feces. The suppression of hepatic BA biosynthesis and intestinal BA reabsorption may explain the opposite changing pattern of BAs in plasma and feces. In summary, the alterations of BA profiles in response to rapid nutrient supply probably reflects a combination of decreased liver biosynthesis, attenuated enterohepatic circulation, reduced gallbladder emptying of PBAs owing to decreased nutrient consumption,ref. jdb13356-bib-0039 and the changes in the abundance of BA‐metabolizing taxa (changes in gut microbiota related to 7α‐dehydroxylase and BSH³⁵ such as Lactobacillus and Bifidobacterium as uncovered in our study).

CR intervention induces a long‐term reduction in nutrition supply. A clinical trial including 80 postmenopausal women who were overweight or obese showed 8‐week VLCD intervention increased HCA level and decreased GCDCA, GDCA, DCA TCA, and TDCA levels in feces, compared with a 4‐week weight maintenance.ref. jdb13356-bib-0034 Mice with 25% CR (25% less than the average daily calorie intake) for 14 days induced changes of BA profiles characterized by decreased TCA and TDCA levels in feces and increased CDCA, TUDCA, and TCA levels in plasma, compared with the ad libitum fed mice.ref. jdb13356-bib-0040 The aforementioned studies emphasized a relatively longer‐term nutrient restriction induced significant changes in BA profiles. Setchell et al. demonstrated that serum unconjugated BAs (including LCA, UDCA, DCA, CDCA, and CA) increased after breakfast but returned to fasting level in the absence of lunch in two healthy subjects.ref. jdb13356-bib-0041 These results were consistent with our findings in fasting and refeeding mice. For the first time, we observed a considerable number of BA species’ changes in plasma, feces, and liver in response to rapid fasting‐refeeding cycle. Compared with preceding observations in sustained CR intervention studies,ref. jdb13356-bib-0034, ref. jdb13356-bib-0040 a more substantial amounts of changes were found in overnight fasting status. The differentially changed BA profiles between overnight fasting and long‐term CR could be explained by the following: (a) fasting is a complete food depletion whereas CR is part of food depletion; (b) after long‐term intervention, compensated regulations may be induced (including regulation of key enzymes); and (c) different gut microbiome and phenotypic changes may also affect BA profiles.

Consistent with the effects of sustained CRref. jdb13356-bib-0042, ref. jdb13356-bib-0043 and IF,ref. jdb13356-bib-0044, ref. jdb13356-bib-0045, ref. jdb13356-bib-0046 our study showed overnight fasting significantly increased α‐diversity of gut microbiota. Our results demonstrated that 24‐h fasting stimulated shifts in the abundance of several bacteria in feces, and these alterations were reversed upon refeeding. Of note, Lactobacillus and Bifidobacterium decreased whereas Akkermansia, Parabacteroides, Muribaculum, Eubacterium, and Muribaculaceae increased in mice fasted for 24 h, and these alterations were also found to be reversed in mice refed for 24 h. Lactobacillus and Bifidobacterium were involved in BA metabolism including BA deconjugation via BSH and 7α‐dehydroxylation via 7α‐dehydroxylase.ref. jdb13356-bib-0008 The decreased levels of plasma unconjugated SBAs and increased levels of plasma conjugated SBAs matched well with decreased abundances of Lactobacillus and Bifidobacterium in fasting mice presented in our study. Thus, the bacterial contribution to the observed BA phenotype in response to rapid fasting‐refeeding cycle is likely connected with the regulation of the conversion of conjugated BAs. More importantly, our study showed overnight fasting also induced a marked, but reversible change in the relative abundances of Muribaculum, Parabacteroides, Eubacterium, and Muribaculaceae, which have not been reported before. Bacteria in the family Muribaculaceae and genus Muribaculum produce propionate as the final product of fermentationref. jdb13356-bib-0047 and respond positively to acarbose in mice.ref. jdb13356-bib-0048 Parabacteroides was found to alleviates obesity by producing succinate and SBAs in mice.ref. jdb13356-bib-0049 Previous studies have shown that Eubacterium, a representative family of butyrate‐producing species, can alleviate intestinal inflammation, T2DM and obesity symptoms by producing short‐chain fatty acids.ref. jdb13356-bib-0050, ref. jdb13356-bib-0051 However, detailed connections among fasting‐induced changes in host BAs, gut microbiota, and metabolism homeostasis remain to be explored in the future.

Of note, we found Lactobacillus and Bifidobacterium decreased whereas Akkermansia increased in mice fasted for 24 h, and these alterations were reversed in refeeding mice, reflecting the rapid changes of these taxa to nutrient supply. Previous studies showed the abundance of Akkermansia was increased whereas the abundance of Bifidobacterium was decreased in wild type C57BL/6 mice fed NCD after IF and CR intervention,ref. jdb13356-bib-0042, ref. jdb13356-bib-0044, ref. jdb13356-bib-0046 which were consistent with our results. In terms of the effect of diets on the abundance of Lactobacillus, previous studies have shown inconsistent results, with some studies suggesting that CR and IF interventions led to an increase in the abundance of Lactobacillus ref. jdb13356-bib-0042, ref. jdb13356-bib-0043, ref. jdb13356-bib-0044, ref. jdb13356-bib-0045 and others showing that IF led to a decrease in the Lactobacillus.ref. jdb13356-bib-0045 Liu et al. found that mice underwent IF (alternative‐day feed deprivation, 15 cycles total) had a significant increase in the abundance of Lactobacillus and Akkermansia compared with the control mice.ref. jdb13356-bib-0044 Zhang et al. investigated the effects of three CR regimens (20% CR, 40% CR, and 60% CR) for 4 weeks on the gut microbiome composition in mice, and they found that CR dose‐dependently reduced the relative abundances of Bifidobacterium and remarkably increased the relative abundances of Lactobacillus and Akkermansia compared with the mice fed ad libitum.ref. jdb13356-bib-0042 Pan et al. reported a unique Lactobacillus‐predominated microbial community was attained in mice after 2‐week 30% CR.ref. jdb13356-bib-0043 Zhang et al. showed 30% CR led to a significant increase in the Lactobacillus, whereas 5:2 IF (2‐day fasting followed by a 5‐day ad libitum period) regimen led to a significant reduction in the Lactobacillus.ref. jdb13356-bib-0045 Li et al. evaluated the effects of three IF regimens for 1 month, that is, 12‐, 16‐, and 20‐h fasting per day, respectively, on the gut microbiome, and they reported that 16‐h fasting regimen led to an increased level of Akkermansia, and these effects disappeared after the cessation of fasting.ref. jdb13356-bib-0046

We noticed that changes of microbiota taxa especially Lactobacillus, Bifidobacterium, and Akkermansia were closely related to the changes of BAs and glucose during the fasting‐refeeding cycle. Lactobacillus and Bifidobacterium were proved to be related to BA metabolism in previous studies.ref. jdb13356-bib-0008 Akkermansia has been proven to improve obesity, T2DM, hepatic steatosis, intestinal inflammation, and different cancers in mice.ref. jdb13356-bib-0052 In our study, Akkermansia was negatively correlated with plasma unconjugated PBAs, plasma unconjugated SBAs, and glucose levels, whereas it was positively correlated with plasma conjugated SBAs, feces unconjugated PBAs, and feces unconjugated SBAs. We highlighted the potential impact of Akkermansia on regulating BA metabolism during the fasting‐refeeding cycle. Several small sample size clinical studies have also demonstrated the change of BAs is closely related to the change of gut microbiota under CR diet intervention. Schwartzenberg et al. found VLCD‐induced Clostridium difficile enrichment was associated with decreased BAs levels.ref. jdb13356-bib-0034 A prospective cohort study including 10 obese postmenopausal women showed VLCD‐induced microbiota taxa change, especially Bifidobacterium, has great impact on the SBAs metabolism pathway.ref. jdb13356-bib-0053 They also found positive correlations between Akkermansia and allolithocholic acid (an isomer of LCA).ref. jdb13356-bib-0053

Our correlation analysis demonstrated that gut microbiota showed closer associations with plasma BAs than fecal BAs, suggesting changes in plasma BA profile were most closely related to the changes in gut microbiota during rapid fasting‐refeeding cycle. As BAs circulate, the majority (~ 95%) are actively absorbed across the distal small intestine wall.ref. jdb13356-bib-0006 A small amount (~ 5%) escape ileal bile salt transport and enter the large intestine. In the large intestine primarily, but also in the small intestine, BAs are subjected to chemical modification by bacteria.ref. jdb13356-bib-0006 Sayin et al. comprehensively mapped the BA profiles of C57BL/6 mice, and they found the serum BA profile most closely resembled the BAs in the distal small intestine rather than that in the colon and feces, which probably reflects the fact that most BAs are reabsorbed from the ileum.ref. jdb13356-bib-0054 Thus, the BA profiles in feces were affected by many factors such as intestinal reabsorption rate and gut microbiota, which may explain the relatively weak association between fecal BA and gut microbiota. Future studies are needed to elucidate this phenomenon.

Furthermore, we found that the BA changes in plasma, feces, and liver during the fasting‐refeeding cycle may also contribute differently to microbial community. Plasam glucose contributed the most to the changes of gut microbiota, followed by plasma GCA and TDCA. LCA in feces and TDCA in liver contributed the most to gut microbiota remodeling. Thus, our study demonstrated that during rapid fasting‐refeeding cycle, there existed an interaction between BAs and gut microbiota, which could conceivably have a significant impact in the dietary intervention‐resulted metabolic outcome. A comprehensive assessment of BA profile‐gut microbiota crosstalk during the fasting‐refeeding cycle could provide insights into possible pathways that define the pathogenesis of metabolic diseases.

In summary, our work clearly demonstrated that overnight fasting led to significant, but reversible changes in BA profile in plasma, feces, and liver. Overnight fasting also induced reversible alterations in the gut microbiome structure. BA biosynthesis in liver and intestinal BA reabsorption were suppressed under fasting status, whereas refeeding ameliorated these changes. More broadly, this study emphasized the important and rapid impact of nutrient supply on BA profile and gut microbiota, allowing the mechanistic dissection of complex relationships between the host, microbiome, and BAs metabolism.

DISCLOSURE

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

References

1. Putting the balance back in diet.. Cell., 2015. [DOI | PubMed]
2. Impact of nutrient overload on metabolic homeostasis.. Nutr Rev., 2018. [DOI | PubMed]
3. Impact of nutritional changes on nonalcoholic fatty liver disease.. Nutrients., 2019. [DOI | PubMed]
4. Modeling diet‐induced metabolic syndrome in rodents.. Mol Nutr Food Res., 2020. [DOI | PubMed]
5. Diet quality, change in diet quality and risk of incident CVD and diabetes.. Public Health Nutr., 2020. [DOI | PubMed]
6. Bile acids in glucose metabolism in health and disease.. J Exp Med., 2018. [DOI | PubMed]
7. Intestinal crosstalk between bile acids and microbiota and its impact on host metabolism.. Cell Metab., 2016. [DOI | PubMed]
8. Bile acid‐microbiota crosstalk in gastrointestinal inflammation and carcinogenesis.. Nat Rev Gastroenterol Hepatol., 2018. [DOI | PubMed]
9. A top‐down systems biology view of microbiome‐mammalian metabolic interactions in a mouse model.. Mol Syst Biol., 2007. [PubMed]
10. The solute carrier family SLC10: more than a family of bile acid transporters regarding function and phylogenetic relationships.. Naunyn Schmiedebergs Arch Pharmacol., 2006. [PubMed]
11. Role of bile acids and bile acid receptors in metabolic regulation.. Physiol Rev., 2009. [DOI | PubMed]
12. Prevention and treatment of obesity, insulin resistance, and diabetes by bile acid‐binding resin.. Diabetes., 2007. [PubMed]
13. The farnesoid X receptor modulates hepatic carbohydrate metabolism during the fasting‐refeeding transition.. J Biol Chem., 2005. [PubMed]
14. Bile acids regulate gluconeogenic gene expression via small heterodimer partner‐mediated repression of hepatocyte nuclear factor 4 and Foxo1.. J Biol Chem., 2004. [DOI | PubMed]
15. Role of bile acids and GLP‐1 in mediating the metabolic improvements of bariatric surgery.. Gastroenterology., 2019. [DOI | PubMed]
16. Lowering bile acid pool size with a synthetic farnesoid X receptor (FXR) agonist induces obesity and diabetes through reduced energy expenditure.. J Biol Chem., 2011. [DOI | PubMed]
17. Farnesoid X receptor deficiency improves glucose homeostasis in mouse models of obesity.. Diabetes., 2011. [DOI | PubMed]
18. TGR5‐mediated bile acid sensing controls glucose homeostasis.. Cell Metab., 2009. [DOI | PubMed]
19. Association of Serum Bile Acids Profile and Pathway Dysregulation with the risk of developing diabetes among normoglycemic Chinese adults: findings from the 4C study.. Diabetes Care., 2021. [DOI | PubMed]
20. A dysregulated bile acid‐gut microbiota axis contributes to obesity susceptibility.. EBioMedicine., 2020. [DOI | PubMed]
21. Effect of insulin resistance on monounsaturated fatty acid levels: a multi‐cohort non‐targeted metabolomics and mendelian randomization study.. PLoS Genet., 2016. [DOI | PubMed]
22. Gut microbiota‐bile acid‐interleukin‐22 axis orchestrates polycystic ovary syndrome.. Nat Med., 2019. [DOI | PubMed]
23. Akkermansia muciniphila improves metabolic profiles by reducing inflammation in chow diet‐fed mice.. J Mol Endocrinol., 2017. [PubMed]
24. The influence of the gut microbiome on obesity, metabolic syndrome and gastrointestinal disease.. Clin Transl Gastroenterol., 2015. [DOI | PubMed]
25. Gut microbiota orchestrates energy homeostasis during cold.. Cell., 2015. [DOI | PubMed]
26. Cross‐talk between Akkermansia muciniphila and intestinal epithelium controls diet‐induced obesity.. Proc Natl Acad Sci USA., 2013. [DOI | PubMed]
27. Akkermansia muciniphila secretes a glucagon‐like peptide‐1‐inducing protein that improves glucose homeostasis and ameliorates metabolic disease in mice.. Nat Microbiol., 2021. [DOI | PubMed]
28. Pediatric obesity is associated with an altered gut microbiota and discordant shifts in firmicutes populations.. Environ Microbiol., 2017. [DOI | PubMed]
29. A purified membrane protein from Akkermansia muciniphila or the pasteurized bacterium improves metabolism in obese and diabetic mice.. Nat Med., 2017. [DOI | PubMed]
30. Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof‐of‐concept exploratory study.. Nat Med., 2019. [DOI | PubMed]
31. Impact of microbial transformation of food on health ‐ from fermented foods to fermentation in the gastro‐intestinal tract.. Curr Opin Biotechnol., 2011. [DOI | PubMed]
32. The microbiome as a circadian coordinator of metabolism.. Endocrinology., 2020. [DOI]
33. Time‐limited diets and the gut microbiota in cardiometabolic disease.. J Diabetes., 2022. [DOI | PubMed]
34. Caloric restriction disrupts the microbiota and colonization resistance.. Nature., 2021. [DOI | PubMed]
35. Restructuring the gut microbiota by intermittent fasting lowers blood pressure.. Circ Res., 2021. [DOI | PubMed]
36. Unique and shared responses of the gut microbiota to prolonged fasting: a comparative study across five classes of vertebrate hosts.. FEMS Microbiol Ecol., 2014. [DOI | PubMed]
37. Denitrifying sulfide removal process on high‐salinity wastewaters in the presence of Halomonas sp.. Appl Microbiol Biotechnol., 2016. [DOI | PubMed]
38. Numerical analysis of grassland bacterial community structure under different land management regimens by using 16 S ribosomal DNA sequence data and denaturing gradient gel electrophoresis banding patterns.. Appl Environ Microbiol., 2001. [PubMed]
39. Gallbladder motility and gallstone formation in obese patients following very low calorie diets. Use it (fat) to lose it (well).. Int J Obes Relat Metab Disord., 1998. [PubMed]
40. Caloric restriction increases levels of taurine in the intestine and stimulates taurine uptake by conjugation to glutathione.. J Nutr Biochem., 2021. [DOI | PubMed]
41. Diurnal changes in serum unconjugated bile acids in normal man.. Gut., 1982. [PubMed]
42. Caloric restriction, friend or foe: effects on metabolic status in association with the intestinal microbiome and metabolome.. J Agric Food Chem., 2022. [DOI | PubMed]
43. Predominant gut lactobacillus murinus strain mediates anti‐inflammaging effects in calorie‐restricted mice.. Microbiome., 2018. [DOI | PubMed]
44. Intermittent fasting reshapes the gut microbiota and metabolome and reduces weight gain more effectively than melatonin in mice.. Front Nutr., 2021. [DOI | PubMed]
45. The effect of calorie intake, fasting, and dietary composition on metabolic health and gut microbiota in mice.. BMC Biol., 2021. [DOI | PubMed]
46. The effects of daily fasting hours on shaping gut microbiota in mice.. BMC Microbiol., 2020. [DOI | PubMed]
47. Sequence and cultivation study of Muribaculaceae reveals novel species, host preference, and functional potential of this yet undescribed family.. Microbiome., 2019. [DOI | PubMed]
48. Genomes assembled from metagenomes suggest genetic drivers of differential response to acarbose treatment in mice.. mSphere., 2021. [DOI | PubMed]
49. Parabacteroides distasonis alleviates obesity and metabolic dysfunctions via production of succinate and secondary bile acids.. Cell Rep., 2019. [DOI | PubMed]
50. The gut microbiota, bacterial metabolites and colorectal cancer.. Nat Rev Microbiol., 2014. [DOI | PubMed]
51. Gut microbes from the phylogenetically diverse genus and their various contributions to gut health.. Gut Microbes., 2020. [DOI | PubMed]
52. Akkermansia muciniphila: paradigm for next‐generation beneficial microorganisms.. Nat Rev Gastroenterol Hepatol., 2022. [DOI | PubMed]
53. Fecal microbiota and bile acid interactions with systemic and adipose tissue metabolism in diet‐induced weight loss of obese postmenopausal women.. J Transl Med., 2018. [DOI | PubMed]
54. Gut microbiota regulates bile acid metabolism by reducing the levels of tauro‐beta‐muricholic acid, a naturally occurring FXR antagonist.. Cell Metab., 2013. [DOI | PubMed]