- Open Access
Treadmill exercise has minimal impact on obesogenic diet-related gut microbiome changes but alters adipose and hypothalamic gene expression in rats
Nutrition & Metabolism volume 17, Article number: 71 (2020)
Exercise has been extensively utilised as an effective therapy for overweight- and obesity-associated changes that are linked to health complications. Several preclinical rodent studies have shown that treadmill exercise alongside an unhealthy diet improves metabolic health and microbiome composition. Furthermore, chronic exercise has been shown to alter hypothalamic and adipose tissue gene expression in diet-induced obesity. However, limited work has investigated whether treadmill exercise commenced following exposure to an obesogenic diet is sufficient to alter microbiome composition and metabolic health.
To address this gap in the literature, we fed rats a high-fat/high-sugar western-style cafeteria diet and assessed the effects of 4 weeks of treadmill exercise on adiposity, diet-induced gut dysbiosis, as well as hypothalamic and retroperitoneal white adipose tissue gene expression. Forty-eight male Sprague-Dawley rats were allocated to either regular chow or cafeteria diet and after 3 weeks half the rats on each diet were exposed to moderate treadmill exercise for 4 weeks while the remainder were exposed to a stationary treadmill.
Microbial species diversity was uniquely reduced in exercising chow-fed rats, while microbiome composition was only changed by cafeteria diet. Despite limited effects of exercise on overall microbiome composition, exercise increased inferred microbial functions involved in metabolism, reduced fat mass, and altered adipose and hypothalamic gene expression. After controlling for diet and exercise, adipose Il6 expression and liver triglyceride concentrations were significantly associated with global microbiome composition.
Moderate treadmill exercise induced subtle microbiome composition changes in chow-fed rats but did not overcome the microbiome changes induced by prolonged exposure to cafeteria diet. Predicted metabolic function of the gut microbiome was increased by exercise. The effects of exercise on the microbiome may be modulated by obesity severity. Future work should investigate whether exercise in combination with microbiome-modifying interventions can synergistically reduce diet- and obesity-associated comorbidities.
Overweight and obesity leads to reductions in physical and mental health, and quality of life , resulting in increased direct and indirect costs to the global economy . Along with gross metabolic changes, obesity is associated with altered fecal microbial species diversity  and composition . Separate studies involving transfer of obese human fecal microbiome samples induced fat gain in naïve mice  and supplementation with Akkermansia muciniphila improved insulin sensitivity and reduced body weight in overweight and obese people , providing some evidence for a potential role of diet- and obesity-associated gut microbiota changes in adiposity and metabolic dysfunction.
Weight loss through lifestyle intervention is an effective strategy for reducing obesity-related comorbidities . One such intervention is moderate exercise, a practical and sustainable approach for people with overweight and obesity : while exercising at this intensity is unlikely to cause weight loss independent of caloric restriction, it confers cardiovascular and metabolic benefits, and assists with weight maintenance [9, 10]. Furthermore, regular exercise is known to improve glucose regulation and insulin sensitivity  as well as reducing cardiovascular disease  and cancer risk , and there is increasing interest in the effects of exercise on the gut microbiota.
The first study to indicate an effect of exercise on fecal microbiome showed that elite, professional athletes exhibited a distinct microbiome composition with increased microbial species diversity . However, since athletes consume a distinct diet from healthy people in the community, further work has been undertaken in rodents to identify the specific effects of exercise on fecal microbiome. A recent systematic review of primarily rodent studies concluded that while there was no consistent effect of exercise on microbial species richness, exercise increases the relative abundance of Firmicutes . The different types of exercise used (forced versus voluntary) have been shown to exert different effects on fecal microbiome composition in mice [16, 17] which may in part explain the inconsistent findings.
Furthermore, there were considerable differences in exercise duration used (6–16 weeks) which may contribute to the range of responses observed. Microbiome compositional changes were observed in mice maintained on a healthy diet following 6 weeks of moderate treadmill exercise  and after 16 weeks in mice fed a high-fat diet . In contrast, eight weeks of low-to-moderate exercise did not confer microbiome compositional changes in mice fed a high-fat diet . Furthermore, most studies examining the effects of treadmill exercise on fecal microbiome in diet-induced obesity used a design where exercise was co-administered with high-fat diet, making it difficult to translate the findings to people, in terms of implementing exercise after a history of unhealthy eating and obesity.
Here we sought to examine whether moderate treadmill exercise in rats could exert benefits to gut microbiome composition following exposure to either a healthy or a high-fat, high-sugar western-style cafeteria diet. We aimed to investigate whether any changes in gut microbiome were associated with altered gene expression in white adipose tissue (WAT) and the hypothalamus, which are known to be affected by both obesogenic diets and exercise.
Materials and methods
Subjects and diet manipulation
This protocol was approved by the Animal Care and Ethics Committee of UNSW Sydney in accordance with the Australian guidelines for the use and care of animals for scientific purposes (National Health and Medical Research Council).
Forty-eight male Sprague-Dawley rats (6–7 weeks, 165-185 g; Animal Resource Centre, Australia) were housed 3/box (18-22 °C; 12 h light/dark) and handled daily for one week while maintained on standard chow (11 kJ/g; Premium Rat and Mouse Maintenance diet; Gordon’s Specialty Stock Feeds, Australia) and water ad libitum.
Following acclimatization, weight-matched groups were randomly allocated to Chow plus water or Cafeteria diet (n = 12 rats, n = 4 cages per group) ad libitum which comprised 10% sucrose solution alongside commercially-produced cakes, cookies and savoury foods  in addition to chow and water. Body weight and 24-h food intake were measured twice weekly and food intake was calculated assuming equal intake per rat in each cage. Body composition was measured during week 7 by EchoMRI-900 (EchoMRI LLC, USA).
Following 3 weeks, half the rats in each dietary condition were allocated to treadmill exercise (generating 4 weight-matched groups: Chow Sedentary (CSed), Chow Exercise (CEx), Cafeteria Sedentary (CafSed) and Cafeteria Exercise (CafEx)) for 4 weeks until the day before sacrifice. Moderate exercise consisted of 10-12 m/min for 45 min, five days a week at zero incline. The first week comprised two training sessions (0–10 m/min) after which time spent at 12 m/min was gradually increased. Sedentary rats were placed in stationary treadmills during the exercise protocol. Exercising rats were closely monitored. Three CafEx rats showed signs of fatigue, distress or vocalisation, and were removed from the treadmill and thereafter exercised at a lower intensity (5-6 m/min for 45 min). These rats were analysed together with moderately exercised CafEx rats.
At 7 weeks of diet, rats were deeply anaesthetized (ketamine/xylazine 15/100 mg/kg intraperitoneally). Body weight, naso-anal length, girth and blood glucose were measured following induction of anesthesia. Rats were decapitated for trunk blood.
The hypothalamus (within coronal block defined by rostro-caudal limits of Circle of Willis) was rapidly dissected and collected. Retroperitoneal and gonadal WAT, and liver were dissected and weighed. One fecal pellet was removed from the distal colon. Tissue and feces were snap frozen in liquid nitrogen and stored at -80 °C.
Protein and triglyceride measurements
Plasma leptin and insulin concentrations were measured using commercial kits (CAT#90040 and CAT#90060, CrystalChem Inc., USA). Plasma and liver triglyceride content were measured spectrophotometrically using triglyceride reagent (Roche Diagnostics Australia Pty Ltd., Australia) at 37 °C alongside a standard curve generated from glycerol standard (G7793-5ML, Sigma-Aldrich Pty Ltd., Australia). Livers were extracted by homogenization in 2:1 chloroform/methanol and incubated overnight; 0.6% NaCl was added and samples were vortexed and centrifuged (1000 g for 10 min). The lower phase was then evaporated at 40 °C under nitrogen gas. Dried extract was re-dissolved in absolute ethanol and measured spectrophotometrically. Retroperitoneal WAT IL-6 content was determined using a DuoSet IL-6 Rat ELISA (DY506, R&D Systems Inc.) following manufacturers’ recommendations.
RNA extraction and gene expression assays
RNA was extracted from hypothalamus and retroperitoneal WAT using TRI Reagent (Sigma-Aldrich Pty Ltd., Australia). Following DNAse I treatment (Catalogue# 42885; Merck, Australia), 1.5 or 2 μg of RNA (WAT and hypothalamus respectively) were reverse transcribed to produce cDNA (High Capacity Reverse Transcriptase Kit; Thermofisher Scientific, USA). Gene expression was assessed using Taqman inventoried gene expression assays (Life Technologies Australia Pty Ltd., Australia; see Supplementary Table 1). Genes of interest were normalized against the geometric mean of the two most stable housekeeping genes (Gapdh and Hprt1 for WAT, Hprt1 and B2m for hypothalamus) identified by the Normfinder package . Analysis of relative gene expression was performed using the ∆∆CT method normalized to an independent calibrator  made from all samples.
Data were analyzed using two-way between-subjects ANOVA, while measures over time were analysed using 3-way mixed ANOVA. Post-hoc pairwise comparisons were performed using a Tukey adjustment where appropriate (THSD) and presented in the associated figures and tables only when p < 0.05. Pearson’s correlations were used to identify associations. All analyses were completed using IBM SPSS Statistics 23 (Australia).
Fecal DNA extraction, microbiome community sequencing and statistical analyses
DNA extraction was performed using the PowerFecal DNA Isolation Kit (MoBio Laboratories, USA). Microbial community composition was assessed by Illumina amplicon sequencing (2 × 250 bp MiSeq chemistry, V4 region, 515F-806R primer pair) using a standard protocol. Sequence data were analyzed using Mothur  using modified commands from MiSeq SOP , including alignment with the SILVA database, singleton removal, chimera checking with UCHIME and classification against the latest RDP training set. Sequence depth was normalized by subsampling to 8346 total clean reads per sample.
Operational taxonomic unit (OTU) correlations were completed using Calypso , with the Benjamini-Hochberg false discovery rate (FDR) procedure  used to control for multiple tests. FDR-corrected DESeq2 was performed using the Phyloseq  R package for the negative binomial Walk test in DESeq2 . OTU abundances were analyzed using SPSS with Kruskal-Wallis tests, followed by non-parametric Bonferroni-Dunn post-hoc testing where appropriate. OTUs of interest were identified using SINA Aligner .
Alpha diversity metrics, distanced-based linear modelling (dbLM), permutational multivariate ANOVA (PERMANOVA), non-metric multidimensional scaling (nMDS) and canonical analysis of principal coordinates were completed using Primer V6 (Primer-E Ltd., Plymouth, United Kingdom ). All Primer analyses utilized a Bray-Curtis similarity matrix constructed at the OTU level.
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) was performed using Galaxy web to predict putative functions (through metagenomic prediction) from the 16S OTU data using Greengenes 13.5 for taxonomic classification . Pathway counts were compared across groups using FDR-corrected Kruskal-Wallis tests followed by non-parametric Bonferroni-Dunn post-hoc testing where appropriate.
Energy intake, body weight and composition, and WAT gene expression
Over the 7-week study, cafeteria-fed rats ate more than twice the energy consumed by chow-fed groups (Fig. 1a; CSed: 18990 kJ/rat; CEx: 19380 kJ/rat; CafSed: 49234 kJ/rat; CafEx: 53485 kJ/rat). When energy intake was stratified into pre- and during-exercise intervention, a significant interaction between time and diet was observed (Fig. 1b; F (1, 12)=9.42, p = 0.01) which appeared to be due to increased cafeteria diet intake while rats were exercising, although this comparison did not reach statistical significance. Significant interactions between time and diet were also identified for fat and carbohydrate intakes (Supplementary Figure 1A; F (1, 12)=5.92, p = 0.032 and Supplementary Figure 1C; F (1, 12)=6.83, p = 0.023 respectively) where both macronutrient intakes increased in the cafeteria-fed rats during exercise. Protein intake, while unaffected by time or exercise, was elevated by cafeteria diet (Supplementary Figure 1B; F (3, 12)=24.06, p < 0.001). Sucrose intake increased over time (Supplementary Figure 1D; F (1, 6)=51.90, p < 0.001) but was not affected by exercise.
All rats gained body weight over time (Fig. 1c) and cafeteria-fed rats gained significantly more than chow-fed controls prior to (F (1, 44)=61.54, p < 0.001) and during (Fig. 1d; F (1, 44)=57.70, p < 0.001) exercise. Exercise significantly reduced weight gain overall (F (1, 44)=8.97, p = 0.004). THSD post-hoc comparisons showed that compared to CSed and CEx, CafSed (p < 0.001 and p < 0.001) and CafEx (p = 0.038 and p < 0.001) gained significantly more weight over the exercise intervention. Of note, CafSed rats gained significantly more weight than CafEx rats over the exercise intervention (p = 0.043).
Body composition data following 4 weeks of treadmill exercise showed that relative fat mass was significantly increased by cafeteria diet (Fig. 1e; F (1, 44)=154.58, p < 0.001; absolute fat mass presented in Table 1) and reduced by exercise (F (1, 44)=6.35, p = 0.015). THSD post-hoc comparisons revealed that relative fat mass in both cafeteria-fed groups were significantly greater than both chow-fed groups. Lean mass however showed only an overall diet effect (Fig. 1f; F (1, 44)=27.34, p < 0.001) with significantly more lean mass in cafeteria-fed rats than chow-fed rats.
Cafeteria diet also increased naso-anal length, girth, plasma insulin and triglycerides, with no effect of exercise (Table 1). Plasma leptin levels were significantly increased by diet (F (1, 43)=185.1, p < 0.001) and reduced by exercise (F (1, 43)=4.64, p = 0.037). Unfasted blood glucose did not differ between groups (Table 1).
In line with plasma leptin, retroperitoneal and epidydimal fat pad weights were significantly greater in cafeteria- than chow-fed rats (F (1, 44)=90.21, p < 0.001 and F (1, 44)=161.41, p < 0.001 for retroperitoneal and epidydimal fat pads respectively), and were significantly reduced by exercise overall (F (1, 44)=5.73, p = 0.021 and F (1, 44)=4.92, p = 0.032 for retroperitoneal and epidydimal fat pads respectively) (See Table 1).
Microbial species diversity and microbiome composition
Microbial species diversity was assessed using Shannon’s diversity index, microbial species richness and microbial species evenness. Shannon’s diversity and evenness were significantly reduced by cafeteria diet overall (F (1, 44)=5.202, p = 0.027; Fig. 2a and f (1, 44)=8.278, p = 0.006; Fig. 2c respectively) and THSD post-hoc comparisons showed that CafEx rats exhibited reduced evenness relative to CEx (p = 0.017) and reduced Shannon’s diversity relative to CSed (p = 0.040). No significant differences were observed for bacterial species richness (Fig. 2b).
Microbiome composition at the OTU level was significantly affected by diet (Pseudo-F (1, 32)=12.17, p = 0.030) and cage (Pseudo-F (3, 32)=1.55, p = 0.006) but not by exercise when assessed using 4-way PERMANOVA (999 permutations) and confirmed with non-metric multidimensional scaling (Fig. 2d). Supplementary Figure 2 shows groups differences in microbiome composition at the Phylum level.
DESeq2 analyses were used to identify OTUs differentially enriched with exercise exposure amongst the top 200 OTUs. While exercise was not associated with differentially expressed OTUs in cafeteria-fed rats, Muribaculum_OTU72 was significantly depleted in CEx rats relative to CSed (Fig. 2e, relative abundance in Fig. 2f; this OTU was originally classified as Akkermansia when aligned with the RDP reference library). Muribaculum_OTU72 was putatively identified as an unknown bacterium from the genus Muribaculum using SINA Aligner (97.6% alignment identity).
Predicted microbiome function
To determine whether the subtle microbiome composition changes observed with exercise affected microbiome function, we inferred microbiome function using PICRUSt. Following an FDR correction, amino acid metabolism (Fig. 2g; H (4)=16.3, p = 0.001), overall energy metabolism (Fig. 2h; H (4)=12.64, p = 0.006), D-glutamine and D-glutamate metabolism (Fig. 2i; H (4)=14.58, p = 0.002) and one carbon pool by folate (Fig. 2j; H (4)=11.44, p = 0.010) exhibited overall group differences. Amino acid metabolism was significantly elevated in CafEx relative to CSed (p < 0.001) and CafSed (p = 0.023) rats while the 3 other processes were significantly elevated in both exercised groups relative to CSed.
WAT and hypothalamic gene expression
Examination of WAT inflammatory signaling and browning genes (Fig. 3a) revealed a significant interaction effect for Ucp1 (F (1, 40)=4.41, p = 0.042) while THSD post-hoc comparisons showed that CafEx rats exhibited elevated Ucp1 expression relative to CafSed (p = 0.024). Lep expression was significantly elevated by cafeteria diet consumption (F (1, 41)=4.98, p = 0.031), while Lepr expression was increased with exercise overall (F (1, 41)=4.75, p = 0.035). No significant differences were observed in pro-inflammatory markers or Adipoq gene expression.
Hypothalamic gene expression was analyzed to determine if 4 weeks of moderate exercise was sufficient to reverse diet-induced changes in expression of genes involved in feeding (Fig. 3b), blood brain barrier integrity and pro-inflammatory signaling (Fig. 3c). A significant interaction effect was observed for Npy gene expression in the hypothalamus (F (1, 41)=7.58, p = 0.009), and THSD showed that CafSed rats exhibited downregulated Npy relative to CSed (p < 0.001), CEx (p < 0.001) and CafEx rats (p < 0.001). No significant differences were observed for Pomc, Agrp, Npy1r, Lepr or Insr. Crh gene expression was significantly increased in cafeteria-fed rats overall (F (1, 41)=14.89, p < 0.001), but was unaffected by exercise (F (1, 41)=1.06, p = 0.309). Post-hoc THSD analysis revealed that CafEx rats exhibited significantly upregulated Crh relative to CSed (p = 0.001) and CEx (p = 0.001).
Both Cln5 and Glut1 were significantly increased in cafeteria-fed rats overall (F (1, 41)=10.96, p = 0.002 and F (1, 41)=5.44, p = 0.025 respectively). While no group differences were apparent for Cln5 using THSD comparisons, Glut1 expression was significantly elevated in CafSed rats relative to CSed controls (p = 0.045). No significant differences were observed in Ocln, Tjp1 or any of the pro-inflammatory genes assessed.
Associations between variables of interest and microbiome
When variables were assessed for their unique contribution to the variance in overall microbiome composition, several adiposity measures, as well as retroperitoneal WAT Il6 gene expression and hypothalamic Crh and Npy expression, were identified as significant predictors of global microbiome composition (Supplementary Table 2). When the contribution of variables of interest to overall microbiome composition was assessed while controlling for diet and treadmill exercise, both liver triglyceride concentration (R2 = 0.031, p = 0.008) and retroperitoneal WAT Il6 gene expression (R2 = 0.026, p = 0.027) were significant predictors of microbiome composition at the OTU level (Table 2, complete model predicts 34.7% of the variance in microbiome composition).
Porphyromonadaceae unclassified_OTU106 was increased with cafeteria diet (Fig. 4a; H (4)=11.42, p = 0.010) and significantly associated with WAT Il6 gene expression (Fig. 4b). Interestingly, this OTU was also negatively correlated with hypothalamic Npy gene expression (Fig. 4c). IL-6 protein in WAT exhibited a significant interaction effect (F (1, 41)=5.310, p = 0.026; Fig. 4d) and was positively associated with Porphyromonadaceae unclassified_OTU106 (Fig. 4e), following a similar, although non-significant, trend to that observed with Il6 gene expression. Like Il6 gene expression, WAT IL-6 content was a significant independent predictor of overall microbiome composition (Supplementary Table 2) but was not a significant predictor in the final model, after controlling for diet and exercise. SINA Aligner putatively identified Porphyromonadaceae unclassified_OTU106 as a strain of Bacteroides eggerthii, a Gram-negative bacterium known to hydrolyze carbohydrates including simple sugars .
We found that 4 weeks of moderate treadmill exercise reduced fat mass and plasma leptin concentrations and altered WAT expression of some adipokine- and metabolism-associated genes. Overall microbiome composition and microbial species diversity was changed by cafeteria diet but not by exercise. However, predicted microbial functions associated with metabolism were increased by exercise. Cafeteria diet-induced changes in hypothalamic Npy and Glut1 gene expression were returned to control levels by exercise.
Exercise induced modest changes in gut microbiome composition that were statistically significant in chow-fed rats only and the relative abundance of OTU72 was significantly reduced in chow-fed rats that exercised (CEx) rats relative to sedentary controls (CSed). While CafEx rats exhibited reduced microbial species diversity, this reduction appeared dependent on cafeteria diet exposure rather than exercise. To date work examining the effects of exercise on microbial species diversity has produced inconsistent results : some rodent studies have shown that exercise is associated with reductions in fecal microbial species richness  while others have reported no such effect .
After 4 weeks of treadmill exercise we did not observe changes in overall microbiome composition. Our data are in line with findings in humans showing that changes in microbiome composition are dependent on obesity status, such that more severe obesity is associated with smaller effects of exercise . Since the cafeteria diet used here  tends to produce a more severe metabolic phenotype than purified high-fat diets [35,36,37,38], which to our knowledge have been used in all studies investigating the interrelationship between exercise and diet on microbiome composition, the diet-induced effects on microbiome composition here may be more resistant to the effects of exercise than previously reported. Additionally, a number of rodent studies report no differences in overall microbiome composition with exercise [17, 19], and there is evidence that this effect may be moderated by age  which may have contributed to the inconsistent findings in the literature.
Predicted microbial functions associated with metabolism, specifically overall energy metabolism, amino acid metabolism, one carbon pool by folate and D-glutamine and D-glutamate metabolism, were increased in exercised rats. This is in line with metagenomic results where fecal microbiome from male elite athletes exhibited increased amino acid biosynthesis and overall energy metabolism relative to sedentary, normal-weight controls . While PICRUSt analysis produces predicted functional data, unlike metagenomic analysis, this is an interesting finding that warrants follow-up to determine if and how exercise shifts the metabolic profile of the gut microbiome and whether any such shift is affected by exercise intensity and duration. Furthermore, confirming these results across a range of diets would be useful to determine whether the shift in microbial function with exercise is modulated by the macro- and micronutrients available.
Here, a moderate exercise intervention reduced fat mass and plasma leptin concentrations, and increased WAT Lepr gene expression in both exercised groups and Ucp1 gene expression in exercised, cafeteria-fed (CafEx) rats uniquely. Increased Ucp1 in WAT depots is a marker of adipocyte beiging , known to be promoted by exercise , and is most likely related to exercise-induced fat loss. While there were no significant effects of exercise on WAT pro-inflammatory gene expression, after controlling for diet and exercise WAT Il6 expression was significantly associated with global microbiome composition.
WAT is one of the major sources of IL-6 in obese humans  and mice [44, 45], which is a key component of the low-grade systemic inflammation observed in overweight and obesity  and is associated with insulin resistance . WAT Il6 expression was strongly associated with the relative abundance of a strain of Bacteroides eggerthii (OTU106), a Gram-negative sugar-scavenging bacterium  which is enriched in obese children relative to normal-weight controls . The associations between WAT Il6 expression, global microbiome composition and OTU106 abundance are therefore likely to be due to the effects of cafeteria diet on both adipose inflammatory processes and the gut microbiome. However, probiotic treatment with a strain of Bifidobacterium in mice fed a high-fat diet reduced WAT macrophage infiltration and plasma IL-6 concentration , indicating that changes to the gut microbiome may contribute to WAT inflammatory signaling and IL-6 production. Further studies determining whether specific bacterial species can modulate WAT IL-6 production are warranted, as interventions that could reduce WAT IL-6 expression in obesity may provide an avenue for preventing insulin resistance and type 2 diabetes.
In contrast to the overall effect of exercise on WAT gene expression, hypothalamic genes disrupted by cafeteria diet were typically normalized with exercise, with no differences observed in chow-fed rats. Npy was downregulated in CafSed rats, as shown previously , but normalized to control levels with exercise, which is in line with other rodent work showing increased hypothalamic NPY mRNA and protein in response to both acute and chronic exercise [51,52,53,54]. Hypothalamic Glut1 gene expression was increased in CafSed rats and reduced to control levels with exercise. This is in contrast to studies showing that acute exercise increased GLUT1 protein expression across the rat brain  and prolonged exercise increases whole brain resting glucose uptake in people . Further work investigating acute and chronic exercise-induced changes in Glut1 expression in the hypothalamus and other brain regions is required.
Cafeteria diet-induced Crh upregulation in the hypothalamus was not normalized by exercise. Crh is transcribed to corticotrophin-releasing hormone (CRH), which activates the hypothalamic-pituitary-adrenal (HPA) axis, key in physiological and behavioral responses to stress. Involuntary exercise is a known stressor , and forced exercise increased activity in CRH neurons relative to both voluntary exercise and sedentary control groups in rats .
Exercise can impact the gut microbiome via several pathways, and there are multiple potential mechanisms by which these gut microbiome changes could subsequently confer health benefits. Exercise may shape microbiome composition and function by increasing host metabolic requirements and shifting metabolite availability. A recent study showed that marathon runners exhibit increased abundance of Veillonella atypica, which rely on lactate as a source of carbon, and inoculation with this species increased running time in mice . Alternatively, exercise is known to activate the HPA axis . The HPA axis appears to exert effects on the gut microbiome, and may be impacted by microbiome modification . Therefore, the effects of exercise on stress behaviors and physiology may be in part due to exercise-induced changes to the gut microbiome, or vice versa.
In summary, short-term (4 weeks) treadmill exercise induced subtle microbiome composition changes in chow-fed rats but did not overcome the microbiome changes induced by cafeteria diet. Despite this, predicted microbiome function was altered by exercise in both chow- and cafeteria-fed groups. Exercise reduced adiposity and altered both WAT and hypothalamic gene expression. WAT Il6 expression was significantly associated with microbiome composition and negatively correlated with the abundance of a strain of Bacteroides eggerthii, bacteria known to scavenge sugars. Future work identifying the circumstances where exercise exerts effects on the microbiome, and the potential causal role of these changes for reducing diet- and obesity-associated comorbidities, is warranted.
Availability of data and materials
The study metadata and sequence data are available in the European Nucleotide Archive under accession number PRJEB36541. The study is currently set on private and will be released upon acceptance. All other data will be made available from the corresponding author on reasonable request.
- Adipoq :
- Agrp :
- B2m :
- Cln5 :
- Crh :
Distance-based linear modelling
False discovery rate
- Gapdh :
Glyceraldehyde 3-phosphate dehydrogenase
- Glut1 :
Glucose transporter 1
- Hprt1 :
Hypoxanthine phosephoribosyltransferase 1
- Il6 :
- Insr :
- Lep :
- Lepr :
Non-metric multidimensional scaling
- Npy :
- Npy1r :
Neuropeptide Y receptor 1
- Ocln :
Operational taxonomic unit
Permutational multivariate ANOVA
Phylogenetic investigation of communities by reconstruction of unobserved states
- Pomc :
- Tjp1 :
Tight junction protein 1
Tukey honest significant difference
White adipose tissue
- Ucp1 :
Uncoupling protein 1
Blüher M. Obesity: global epidemiology and pathogenesis. Nat Rev Endocrinol. 2019;15(5):288–98.
Tremmel M, Gerdtham UG, Nilsson PM, Saha S. Economic Burden of Obesity: A Systematic Literature Review. Int J Environ Res Public Health. 2017;14(4):435.
Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500:541..
Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3..
Ridaura VK, Faith JJ, Rey FE, Cheng J, Duncan AE, Kau AL, et al. Gut microbiota from twins discordant for obesity modulate metabolism in mice. Science (New York, NY). 2013;341(6150):1241214.
Depommier C, Everard A, Druart C, Plovier H, Van Hul M, Vieira-Silva S, et al. Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study. Nat Med. 2019;25(7):1096–103..
Dixon JB. The effect of obesity on health outcomes. Mol Cell Endocrinol. 2010;316(2):104–8.
Bray GA, Frühbeck G, Ryan DH, Wilding JPH. Management of obesity. Lancet. 2016;387(10031):1947–56.
Swift DL, McGee JE, Earnest CP, Carlisle E, Nygard M, Johannsen NM. The effects of exercise and physical activity on weight loss and maintenance. Prog Cardiovasc Dis. 2018;61(2):206–13.
Petridou A, Siopi A, Mougios V. Exercise in the management of obesity. Metabolism. 2019;92:163–9.
Mul JD, Stanford KI, Hirshman MF, Goodyear LJ. Exercise and regulation of carbohydrate metabolism. Prog Mol Biol Transl Sci. 2015;135:17–37.
Fiuza-Luces C, Santos-Lozano A, Joyner M, Carrera-Bastos P, Picazo O, Zugaza JL, et al. Exercise benefits in cardiovascular disease: beyond attenuation of traditional risk factors. Nat Rev Cardiol. 2018;15(12):731–43.
Hojman P. Exercise protects from cancer through regulation of immune function and inflammation. Biochem Soc Trans. 2017;45(4):905–11.
Clarke SF, Murphy EF, O'Sullivan O, Lucey AJ, Humphreys M, Hogan A, et al. Exercise and associated dietary extremes impact on gut microbial diversity. Gut. 2014;63(12):1913–20.
Mitchell CM, Davy BM, Hulver MW, Neilson AP, Bennett BJ, Davy KP. Does exercise Alter gut microbial composition? A systematic review. Med Sci Sports Exerc. 2019;51(1):160–7.
Allen JM, Miller MEB, Pence BD, Whitlock K, Nehra V, Gaskins HR, et al. Voluntary and forced exercise differentially alters the gut microbiome in C57BL/6J mice. J Appl Physiol. 2015;118(8):1059–66.
Lamoureux EV, Grandy SA, Langille MGI. Moderate Exercise Has Limited but Distinguishable Effects on the Mouse Microbiome. mSystems. 2017;2(4):e00006–17.
Kang SS, Jeraldo PR, Kurti A, Miller ME, Cook MD, Whitlock K, et al. Diet and exercise orthogonally alter the gut microbiome and reveal independent associations with anxiety and cognition. Mol Neurodegener. 2014;9:36.
Ribeiro FM, Ribeiro CFA, G ACM, Castro AP, Almeida JA, Franco OL, et al. Limited Effects of Low-to-Moderate Aerobic Exercise on the Gut Microbiota of Mice Subjected to a High-Fat Diet. Nutrients. 2019;11(1):149.
Leigh SJ, Kendig MD, Morris MJ. Palatable western-style cafeteria diet as a reliable method for modeling diet-induced obesity in rodents. JoVE. 2019;153:e60262.
Andersen CL, Jensen JL, Orntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64(15):5245–50.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods (San Diego, Calif). 2001;25(4):402–8.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20.
Zakrzewski M, Proietti C, Ellis JJ, Hasan S, Brion M-J, Berger B, et al. Calypso: a user-friendly web-server for mining and visualizing microbiome-environment interactions. Bioinformatics (Oxford, England). 2017;33(5):782–3.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
McMurdie PJ, Holmes S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PloS one. 2013;8(4):e61217.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Pruesse E, Peplies J, Glöckner FO. SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012;28(14):1823–9.
Clarke KR. Non-parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993;18(1):117–43.
Douglas GM, Maffei VJ, Zaneveld J, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2: An improved and extensible approach for metagenome inference. bioRxiv. 2019:672295.
Wexler HM. Bacteroides: the good, the bad, and the nitty-gritty. Clin Microbiol Rev. 2007;20(4):593–621.
Welly RJ, Liu T-W, Zidon TM, Rowles JL 3rd, Park Y-M, Smith TN, et al. Comparison of diet versus exercise on metabolic function and gut microbiota in obese rats. Med Sci Sports Exerc. 2016;48(9):1688–98.
Allen JM, Mailing LJ, Niemiro GM, Moore R, Cook MD, White BA, et al. Exercise alters gut microbiota composition and function in lean and obese humans. Med Sci Sports Exerc. 2018;50(4):747–57.
Sampey BP, Vanhoose AM, Winfield HM, Freemerman AJ, Muehlbauer MJ, Fueger PT, et al. Cafeteria diet is a robust model of human metabolic syndrome with liver and adipose inflammation: comparison to high-fat diet. Obesity (Silver Spring, Md). 2011;19(6):1109–17.
Higa TS, Spinola AV, Fonseca-Alaniz MH, Evangelista FS. Comparison between cafeteria and high-fat diets in the induction of metabolic dysfunction in mice. Int J Physiol, Pathophysiol Pharmacol. 2014;6(1):47–54.
Buyukdere Y, Gulec A, Akyol A. Cafeteria diet increased adiposity in comparison to high fat diet in young male rats. PeerJ. 2019;7:e6656.
Oliva L, Aranda T, Caviola G, Fernandez-Bernal A, Alemany M, Fernandez-Lopez JA, et al. In rats fed high-energy diets, taste, rather than fat content, is the key factor increasing food intake: a comparison of a cafeteria and a lipid-supplemented standard diet. PeerJ. 2017;5:e3697.
Mika A, Van Treuren W, González A, Herrera JJ, Knight R, Fleshner M. Exercise is More Effective at Altering Gut Microbial Composition and Producing Stable Changes in Lean Mass in Juvenile versus Adult Male F344 Rats. PloS one. 2015;10(5):e0125889.
Barton W, Penney NC, Cronin O, Garcia-Perez I, Molloy MG, Holmes E, et al. The microbiome of professional athletes differs from that of more sedentary subjects in composition and particularly at the functional metabolic level. Gut. 2018;67(4):625–33.
Wu J, Boström P, Sparks Lauren M, Ye L, Choi Jang H, Giang A-H, et al. Beige adipocytes are a distinct type of Thermogenic fat cell in mouse and human. Cell. 2012;150(2):366–76.
Aldiss P, Betts J, Sale C, Pope M, Budge H, Symonds ME. Exercise-induced ‘browning’ of adipose tissues. Metabolism. 2018;81:63–70.
Fried SK, Bunkin DA, Greenberg AS. Omental and subcutaneous adipose tissues of obese subjects release interleukin-6: depot difference and regulation by glucocorticoid. J Clin Endocrinol Metab. 1998;83(3):847–50.
Starr ME, Evers BM, Saito H. Age-associated increase in cytokine production during systemic inflammation: adipose tissue as a major source of IL-6. J Gerontol A Biol Sci Med Sci. 2009;64(7):723–30.
Harkins JM, Moustaid-Moussa N, Chung YJ, Penner KM, Pestka JJ, North CM, et al. Expression of interleukin-6 is greater in preadipocytes than in adipocytes of 3T3-L1 cells and C57BL/6J and Ob/Ob mice. J Nutr. 2004;134(10):2673–7.
Maachi M, Piéroni L, Bruckert E, Jardel C, Fellahi S, Hainque B, et al. Systemic low-grade inflammation is related to both circulating and adipose tissue TNFα, leptin and IL-6 levels in obese women. Int J Obes. 2004;28(8):993–7.
Bastard JP, Maachi M, Van Nhieu JT, Jardel C, Bruckert E, Grimaldi A, et al. Adipose tissue IL-6 content correlates with resistance to insulin activation of glucose uptake both in vivo and in vitro. J Clin Endocrinol Metab. 2002;87(5):2084–9.
Lopez-Contreras BE, Moran-Ramos S, Villarruel-Vazquez R, Macias-Kauffer L, Villamil-Ramirez H, Leon-Mimila P, et al. Composition of gut microbiota in obese and normal-weight Mexican school-age children and its association with metabolic traits. Pediatr Obes. 2018;13(6):381–8.
Moya-Pérez A, Neef A, Sanz Y. Bifidobacterium pseudocatenulatum CECT 7765 Reduces Obesity-Associated Inflammation by Restoring the Lymphocyte-Macrophage Balance and Gut Microbiota Structure in High-Fat Diet-Fed Mice. PloS one. 2015;10(7):e0126976.
Hansen MJ, Jovanovska V, Morris MJ. Adaptive responses in hypothalamic neuropeptide Y in the face of prolonged high-fat feeding in the rat. J Neurochem. 2004;88(4):909–16.
Smith JK. Exercise, Obesity and CNS Control of Metabolic Homeostasis: A Review. Front Physiol. 2018;9:574.
Bi S, Scott KA, Hyun J, Ladenheim EE, Moran TH. Running wheel activity prevents Hyperphagia and obesity in Otsuka long-Evans Tokushima fatty rats: role of hypothalamic signaling. Endocrinology. 2005;146(4):1676–85.
Chen J-X, Zhao X, Yue G-X, Wang Z-F. Influence of acute and chronic treadmill exercise on rat plasma lactate and brain NPY, L-ENK, DYN A1–13. Cell Mol Neurobiol. 2007;27(1):1–10.
Wang J, Chen C, Wang R-Y. Influence of short- and long-term treadmill exercises on levels of ghrelin, obestatin and NPY in plasma and brain extraction of obese rats. Endocrine. 2008;33(1):77–83.
Takimoto M, Hamada T. Acute exercise increases brain region-specific expression of MCT1, MCT2, MCT4, GLUT1, and COX IV proteins. J Appl Physiol(Bethesda, Md : 1985). 2014;116(9):1238–50.
Robinson MM, Lowe VJ, Nair KS. Increased brain glucose uptake after 12 weeks of aerobic high-intensity interval training in young and older adults. J Clin Endocrinol Metab. 2018;103(1):221–7.
Moraska A, Deak T, Spencer RL, Roth D, Fleshner M. Treadmill running produces both positive and negative physiological adaptations in Sprague-Dawley rats. Am J Phys Regul Integr Comp Phys. 2000;279(4):R1321–R9.
Yanagita S, Amemiya S, Suzuki S, Kita I. Effects of spontaneous and forced running on activation of hypothalamic corticotropin-releasing hormone neurons in rats. Life Sci. 2007;80(4):356–63.
Scheiman J, Luber JM, Chavkin TA, MacDonald T, Tung A, Pham L-D, et al. Meta-omics analysis of elite athletes identifies a performance-enhancing microbe that functions via lactate metabolism. Nat Med. 2019;25(7):1104–9.
Duclos M, Tabarin A. Exercise and the Hypothalamo-pituitary-adrenal Axis. Front Horm Res. 2016;47:12–26.
Farzi A, Frohlich EE, Holzer P. Gut microbiota and the neuroendocrine system. Neurotherapeutics. 2018;15(1):5–22.
The authors thank Robyn Lawler at the Biological Resources Centre, UNSW Sydney, for assistance with animal husbandry, Dr. Sabiha Chowdhury, Saroj Khatiwada and Josephine Yu for assistance with dissection. This research includes data processed by the computational cluster Katana supported by Research Technology Services at UNSW Sydney.
The work was supported by an NHMRC project grant (#1126929) to MJM and RFW as well as a Brain Sciences UNSW seed grant to SJL, MJM and RFW. SJL was supported by a Research Training Program scholarship. RME was supported by the project grant PSI2016–77234-R from Ministerio de Economía y Competitividad and by the grant PRX15/00455 from Ministerio de Educación, Cultura y Deporte.
Ethics approval and consent to participate
This protocol was approved by the Animal Care and Ethics Committee of UNSW Sydney (Ethics number 16-67A) in accordance with the Australian guidelines for the use and care of animals for scientific purposes (National Health and Medical Research Council).
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Taq assay probe information. Supplementary Figure 1. Weekly macronutrient intake and sucrose intake over the study. Supplementary Figure 2. Average composition of phyla between experimental groups. Supplementary Table 2. Distance-based linear modelling to determine contributions of metabolic and gene expression measures to the variance observed in microbiota composition at the OTU level.
About this article
Cite this article
Leigh, SJ., Kaakoush, N.O., Escorihuela, R.M. et al. Treadmill exercise has minimal impact on obesogenic diet-related gut microbiome changes but alters adipose and hypothalamic gene expression in rats. Nutr Metab (Lond) 17, 71 (2020). https://doi.org/10.1186/s12986-020-00492-6
- White adipose tissue
- Cafeteria diet