Exercise induction of gut microbiota modifications in obese, non-obese and hypertensive rats
- Bernardo A Petriz†1, 5,
- Alinne P Castro†2, 4,
- Jeeser A Almeida1, 3, 5,
- Clarissa PC Gomes1,
- Gabriel R Fernandes1,
- Ricardo H Kruger2,
- Rinaldo W Pereira1, 3 and
- Octavio L Franco1, 3, 4Email author
© Petriz et al.; licensee BioMed Central Ltd. 2014
Received: 3 December 2013
Accepted: 17 June 2014
Published: 21 June 2014
Obesity is a multifactor disease associated with cardiovascular disorders such as hypertension. Recently, gut microbiota was linked to obesity pathogenesisand shown to influence the host metabolism. Moreover, several factors such as host-genotype and life-style have been shown to modulate gut microbiota composition. Exercise is a well-known agent used for the treatment of numerous pathologies, such as obesity and hypertension; it has recently been demonstrated to shape gut microbiota consortia. Since exercise-altered microbiota could possibly improve the treatment of diseases related to dysfunctional microbiota, this study aimed to examine the effect of controlled exercise training on gut microbial composition in Obese rats (n = 3), non-obese Wistar rats (n = 3) and Spontaneously Hypertensive rats (n = 3). Pyrosequencing of 16S rRNA genes from fecal samples collected before and after exercise training was used for this purpose.
Exercise altered the composition and diversity of gut bacteria at genus level in all rat lineages. Allobaculum (Hypertensive rats), Pseudomonas and Lactobacillus (Obese rats) were shown to be enriched after exercise, while Streptococcus (Wistar rats), Aggregatibacter and Sutturella (Hypertensive rats) were more enhanced before exercise. A significant correlation was seen in the Clostridiaceae and Bacteroidaceae families and Oscillospira and Ruminococcus genera with blood lactate accumulation. Moreover, Wistar and Hypertensive rats were shown to share a similar microbiota composition, as opposed to Obese rats. Finally, Streptococcus alactolyticus, Bifidobacterium animalis, Ruminococcus gnavus, Aggregatibacter pneumotropica and Bifidobacterium pseudolongum were enriched in Obese rats.
These data indicate that non-obese and hypertensive rats harbor a different gut microbiota from obese rats and that exercise training alters gut microbiota from an obese and hypertensive genotype background.
KeywordsGut microbiota Obesity Hypertension Exercise Lactobacillus
Exercise practice is a non-pharmacological treatment for a series of diseases . Along with dietary control, appropriate exercise programs are proposed to treat and attenuate obesity  and its associated cardiovascular disorders such as hypertension . It is known that hypertension and obesity frequently coexist , affecting millions of people worldwide .
Gut microbiota have been recently indicated as having a close relationship with obesity, where the microbiota of an obese subject presents an enhanced ability to harvest energy and accumulate fat . The gut harbors the greatest density of these microorganisms in the body (e.g. ~up to 1.5 kg in the human gut)  with Firmicutes, Bacteriodetes and Actinobacteria constituting the dominant phyla . Moreover, obesity is associated with reduced microbiota diversity at phylum level , seen in rodents and in humans [10, 11].
The gut is a dynamic environment, highly exposed to environmental factors such as diet , antibiotics , pathogens  and lifestyle , which constantly interact with microbial communities. In addition, gut microbiota is also shaped throughout life by host-related factors such as host-genotype . Disturbance within gut microbiota has been reported to influence host susceptibility to pathogens and pathological conditions such as gastrointestinal inflammatory diseases and obesity . Moreover, hypertension induction was also seen to alter gut microbiota .
It has been proposed that dysbiosis and pathologies associated with unbalanced gut microbiota may be prevented or treated with prebiotics, probiotics and fecal microbiota transplantation [19, 20]. In addition, controlled exercise intensities are related to protective effects on the gastrointestinal tract, including a reduced risk of colon inflammation and cancer . It is proposed that exercise may reduce intestinal transit time, diminishing the contact between the colon and cancer-promoting agents .
Recently, exercise has also been shown to induce alterations within microbiota composition [22, 23], which suggests that exercise may be included as a possible therapeutic factor along with diet, prebiotics and other treatments. Since exercise plays a prominent role in metabolic regulation and energy expenditure, it might modulate host-microbiota interaction, affecting the host metabolism. Although these relations are still unknown, exercise may enhance the strategies for obesity control, along with other actions such as microbiota transplant .
Although voluntary exercise was shown to alter microbiota in non-pathological animals [22–24], its effects on gut microbiota still need to be further investigated in pathologic phenotypes and through controlled parameters such as exercise volume and intensity. Therefore, in the present study, we proposed to examine the effect of controlled moderate exercise intensity on gut microbial status in rats with different phenotypes, by using pyrosequencing of 16S RNAs genes from fecal microbiota samples. To our knowledge, this is the first study to use controlled exercise parameters and distinct animal strains with different obesity and hypertension genotypes.
Analyzing 16S rRNA sequences revealed a similar microbiota profile shared between Wistar and Hypertensive rats, with both being divergent from Obese rats. Exercise was shown to enhance bacterial diversity and to alter microbial communities at the species level in all animal lineages. Thus, these data contribute to the emerging knowledge regarding the effect of exercise on gut microbiota, but further studies should be performed to establish the mechanism by which exercise signals in the bacterial community and to determine the impact of these modulations on host homeostasis.
Animals were obtained from the Federal University of São Paulo, Brazil (UNIFESP) and started the experiment at ~18 weeks of age. Three different strains from two different genotypes were used: an obese genotype, homozygous (fa/fa) obese Zucker rats  (Obese rats; n = 5, 389.4 ± 21 g) and a hypertensive genotype, spontaneous hypertensive rats (Hypertensive rats; n = 5, 227.4 ± 29.3 g), a strain obtained by the selective breeding of Wistar-Kyoto rats with high blood pressure . A strain of Wistar rats (WR; n = 5, 223.2 ± 27.3 g) was used as normotensive control for SHR and as a non-obese phenotype  (Additional file 1). All animals were allocated to collective cages according to their lineages, being kept in a 12 h light–dark cycle environment, with food and water ad libitum. All experimental procedures and interventions in the present study, involving animal-welfare, were carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health and were approved by the local ethics committee for standards in animal use, at the Institute of Biological Sciences, University of Brasilia, Brazil (UnBDOC no. 48695/2010); these were also in accordance with international standards. After experimental procedures, all animals were deeply anesthetized with 2% Xylazine (50 mg.kg−1) and 10% Ketamine (80 mg.kg−1) and euthanized by cervical dislocation. Through the entire experiment, all efforts were made to minimize animal suffering.
Before the training period, all animals underwent a familiarization period on a treadmill device (Li 870, Letica Scientific Instruments, Barcelona, Spain) for 2 weeks to reduce external stress. Duration and speed on treadmill were increased progressively (up to 12.5 m.min−1 for Obese rats; 20 m.min−1 for Hypertensive and Wistar rats) as previously described [27, 28]. In addition, blood pressure of Wistar and Hypertensive rats was measured by the tail-cuff method  at the beginning of the experiment to characterize the hypertensive phenotype from the SHR strain (171.4 ± 7.7 mm.Hg−1 for Hypertensive rats and 128 ± 5.9 mmHg−1 for Wistar rats) (Additional file 2). As regards the Obese rats group (obese (fa/fa) Zucker rat), these animals are homozygous for the fa allele and are one of the most common models of genetic obesity where rats normally exhibit hyperlipidemia, hyperinsulinemia, hyperphagia and a significant weight gain by the 3rd to the 5th week of age . In this model, conflicting results were reported on whether these obese Zucker rats are hypertensive compared to the lean control rats, showing that blood pressure is not elevated in this model [30, 31]. However, the obesity phenotype may enhance arterial peripheral resistance and the animal may develop hypertension secondary to obesity and to secondary mechanisms [32, 33]. Considering that obese Zucker rats do not present a regular pattern of blood pressure as reported by different studies, the systolic blood pressure was not measured in this group.
After the adaptation period, all animals were trained for 30 min per day, 5 days per week for 4 weeks. Running intensity was set corresponding to maximal lactate steady state (MLSS), previously identified in obese Zucker rats  and SHR rats . Therefore, for Obese rats, running velocity was set at 12.5 m.min−1 and for Hypertensive and Wistar rats at 20 m.min−1. A new MLSS identification was performed after the fourth week of exercise training to assess each animal’s cardiovascular adaptation.
Blood lactate analysis
Capillary blood samples (10 μL) were collected through a small incision in the distal portion of the tail of the animals at rest and every 5 min during the MLSS test. Capillary blood samples were placed in microtubes (0.6 mL) containing 20 μL of 1% sodium fluoride and stored at −20°C. Analyses were performed by electro-enzymatic method with YSI 1500 Sports (Yellow Springs, OH, USA). The MLSS was considered when there was no increase over 1 mmol.L−1 of blood lactate from 10 to 25 min of exercise tests .
Fecal DNA extraction and barcoded pyrosequencing of the 16S rRNA gene
Fecal content from all animals was individually collected in three replicates before exercise training period and in vivo 24 h after the last session of exercise (Additional file 1). Fecal samples were stored in RNA later® (Life technologies) until samples were frozen at −80°C. Three out of five individual fecal samples from each animal lineage were randomly selected. Fecal microbial DNA was extracted from ~0.25 g using the PowerFecal DNA Isolation Kit (MoBio, Carlsbad, CA, USA) according to the manufacturer’s instructions. The triplicate DNA extractions were not pooled.
The bacterial community partial 16S rRNA gene was amplified with primer pair 787 F-1492R . For pyrosequencing analysis, the primer set was modified, where the forward primer included the Roche 454-A pyrosequencing adapter and a 12-bp barcode (unique to each sample), while the reverse primer included only the Roche 454-B pyrosequencing adapter. The 20 μl reaction mixture for the PCR contained approximately 10 ng of metagenomic fecal DNA, 1X PCR buffer (Invitrogen), 3.0 mM MgCl2, 10 pmol of each primer, 0.25 mM dNTP, and 1.5 U Taq DNA polymerase (Invitrogen). The cycling protocol started with an initial denaturation step of 3 min at 95°C, followed by 25 cycles of denaturation for 30 s at 95°C, annealing for 30 s at 58°C, and extension for 1.40 min at 72°C, followed by a final extension for 7 min at 72°C and cooling to 10°C. Finally, the rRNA amplicons from bacterial communities were purified with the QIAquick PCR Purification Kit (Chatsworth, CA). The concentrations of the rRNA amplicons were measured by Qubit fluorometer (Invitrogen), and subsequently the massively parallel GS FLX Titanium technology was performed in the Roche 454 Life Sciences Corporation, Branford, CT, USA.
Analysis of 16S rRNA sequences
A total of 1,398,681 16S rRNA sequences were obtained by the 454 GS FLX Titanium sequencer. All the 16S amplicons were processed by the quantitative insights into microbial ecology (QIIME) pipeline version 1.6.0-dev . Briefly, all the 16S rRNA amplicons were sorted by their barcodes, and subsequently the reads with a length less than 180 bp, ambiguous sequences, bases with Phred values of <30 and their primers, barcodes and adaptor sequences were removed. The remaining sequences were submitted to Denoiser algorithm  to remove pyrosequencing errors. Operational taxonomic units (OTUs) were clustered at 97% similarity using an ‘open-reference’ OTU picking protocol, where sequences are clustered against the Greengenes database  using Uclust . One of the most-abundant reads from each OTU was aligned using the PyNAST algorithm . The chimeric OTUs were detected with ChimeraSlayer , and taxonomic classifications were assigned with the naïve Bayesian classifier of the Ribosomal Database Project (RDP) classifier  applying 80% of confidence threshold. Shannon indices and observed richness were used to evaluate community richness, and the unweighted Unifrac algorithm was performed to generate principal coordinate plots (PCoA).
Statistical differences between the groups pre-exercise and after exercise were tested using the analysis of similarities (ANOSIM) by permutation of group membership with 999 replicates , and bivariate relationships were measured with Pearson correlations and regression analysis available through QIIME . Statistical tests on the taxonomic differences between samples were calculated by Welch’s t-test combined with Welch’s inverted method for calculating confidence intervals (nominal coverage of 95%), using the Statistical Analysis of Metagenomic Profiles (STAMP) software version 2.0.0 (STAMP) .
The 454 FLX Titanium flowgrams (sff files) have been submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive database, project number: PRJNA246617.
Effect of exercise training on aerobic capacity
Composition of fecal microbiota in rat lineages
After quality filtering, 889,124 out of 1,398,681 sequences were obtained from fecal samples collected pre-exercise and post-exercise, after 4 weeks of moderate exercise training (Additional file 1). An average of 49,951 denoised sequences per animal were obtained (average read length of 524.8) which composed an average of 583.1 distinct observed OTUs. Post-training samples presented a higher Shannon index compared to pre-training samples (6.4 ± 0.5 vs. 6.8 ± 0.2) (Additional file 3). Detailed sequencing information from all individual rats is presented in (Additional file 3).
Bacterial diversity was assessed by rarefaction measure of observed species against the number of sequences per sample, and the observed OTUs were identified with 97% of identity. Here, the rarefaction measure showed that species diversity (Additional file 4: A; Wistar rats; B; Hypertensive rats and C; Obese rats) reached a plateau tendency in all samples as the number of sequences increased, indicating that in the present study more sequences are unlikely to yield many additional species. As demonstrated in Additional file 4, the rarefaction curves revealed that OTU richness in post-exercise fecal samples is more species-rich than those found in the pre-exercise fecal samples. This was further evidenced in hypertensive rats and obese rats samples (Additional file 4B and C).
The relative abundance of the main dominant bacterial phylum from all fecal samples collected before and after exercise training is shown in Additional file 5A. Here, Firmicutes and Bacteroidetes are the most dominant phyla, followed by Proteobacteria. Firmicutes was shown to be enhanced after exercise training (1.1 fold change, p < 0.05) (Additional file 5B), thus being more evidenced in Obese rats (Obese rats; 0.69 ± 0.03 vs. Exercised Obese rats; 0.78 ± 0.04, p < 0.05). On the other hand, Proteobacteria were shown to be 1.8 fold reduced after exercise training (p < 0.05) (Additional file 5C). The Bacteroidetes phylum was shown to be 1.3 fold reduced after exercise only in Wistar rats (Wistar rats; 0.23 ± 0.04 vs. Exercised Wistar rats; 0.17 ± 0.03, p < 0.05).
Composition of bacterial communities before, during and after exercise training
Regarding pre-exercise samples, the proportion of sequences from the Bacteroides acidifaciens species was significantly more abundant in Obese rats than in Wistar and Hypertensive rats (p < 0.05) (Figure 3A). However, an opposite profile is observed in sequences attributed to the Ruminococcus flavefaciens species, where greater abundance is seen in Wistar rats followed by Hypertensive rats, with no abundance seen in Obese rats (p < 0.05) (Figure 3B). After exercise training, the proportion of sequences indicated that all listed species (Streptococcus alactolyticus, Bifidobacterium animalis, Ruminococcus gnavus, Aggregatibacter pneumotropica and Bifidobacterium pseudolongum) were more abundant in Obese rats compared to Wistar and Hypertensive rats lineages (Figure 3C-G respectively). The relative abundance of Streptococcus alactolyticus in Obese rats diverged significantly from Hypertensive and Wistar rats (p < 0.05), while a diminished proportion of sequences was seen in both rat strains (Figure 3C). The Bifidobacterium animalis species was seen to be highly enriched in Obese rats (p < 0.05) and absent in Wistar and Hypertensive rats (Figure 3D). In relation to Ruminococcus gnavus, this species was poor in Wistar and almost absent in Hypertensive rats, but more abundant in Obese rats (Figure 3E). The Aggregatibacter pneumotropica species presented a similar profile to the previous species, being also more abundant in Obese rats compared to Wistar and Hypertensive rats (p < 0.05) (Figure 3F). Lastly, from the Actinobacteria phylum, Bifidobacterium pseudolongum abundance was almost exclusive to Obese rats (p < 0.05), being almost absent in Hypertensive rats and completely absent in Wistar rats group (Figure 3G).
Principal coordinates analysis (PCoA)
Correlation of bacterial abundance and lactate concentration
Several environmental  and host-related factors  are known to affect gut microbiota composition. This dynamic ecosystem, highly susceptible to external agents, has a symbiotic link with host health homeostasis. In this context, imbalanced gut microbiota has been associated with the development of inflammatory gastrointestinal diseases, obesity and altered metabolic status .
Recently, physical activity was shown to modulate gut microbiota in diet induced obesity , and in healthy rodents, altering the microbiota diversity and composition  and increasing n-Butyrate concentration in the cecum . In contrast to some of these previous studies that used PCR-TGGE  and PCR-DGGE of bacterial 16S rRNA genes , here a robust pyrosequencing of the 16S rRNA genes was used, along with controlled exercise training parameters to investigate this relationship in non-pathologic and pathologic rat models.
Rarefaction measurements and Shannon index indicated that exercise training enhances bacterial diversity in non-pathological Wistar rats and in Obese and Hypertensive rats (Additional file 3 and Additional file 4). Here, Firmicutes and Bacteroidetes were found to be the most predominant phyla in all animal lineages (Additional file 5A). This predominance was also seen in mice cecum , and in exercised rats . Considering all rat lineages, exercise was shown to enhance Firmicutes abundance and to diminish Proteobacteria content (Additional file 5B, C). Thus, Firmicutes were more abundant in post-exercise samples compared to pre-exercise in obese rats (Obese rats; 0.69 ± 0.03 vs. Exercised Obese rats; 0.78 ± 0.04, p < 0.05), while Bacteroidetes was shown to be reduced after training only in WR (Wistar rats; 0.23 ± 0.04 vs. Exercised Wistar rats; 0.17 ± 0.03, p < 0.05). Moreover, Bacteroidetes have been reported to be diminished in obese mice , while the ratio of Firmicutes to Bacteroidetes was shown to change in favor of Bacteroidetes in overweight and obese subjects, compared to the lean group . Thus, as previously stated by Harris et al.,  studying gut microbiota and their relation with metabolic disorders revealed that there is no difference between obese and lean individuals at phylum level. However, data here reported showed a significant alteration in bacterial community abundance at phylum as well as at genus levels (Figure 2), which could be associated with the effects of exercise and/or pathological conditions.Furthermore, our study revealed a significant alteration in bacterial community abundance at genus and species level as an effect of exercise and/or pathological stimuli (Figure 2).
In accordance with the PCoA analysis presented in this study (Figure 4), other studies have also reported a distinction between non-obese and obese microbiota from Zuckerfa/fa rats  and ob/ob mice . We also reported that Wistar and Hypertensive rats share a similar microbiota composition (Figure 4). In a similar way, it was reported that rats treated with nitric oxide synthase inhibitor NG-nitro-L-arginine methyl ester (L-NAME) develop hypertension, with a variation in cecal microbiota compared to control normotensive rats . Regarding the effect of exercise, PCoA analysis demonstrated that four weeks of moderate exercise training significantly altered microbiota composition in all rat lineages (Figure 4). According to our results, different exercise volumes, 6 days , and 5  and 12 weeks  of voluntary running exercise, were shown to alter microbiota composition, indicating that the microbial community is affected even by a few days of exercise. Together, these data suggest that besides other well-known factors, exercise may be seen as a potential environmental factor capable of modulating gut microbiota.
Here, exercise was also shown to significantly alter six bacterial genera (Figure 2). Fecal samples were more enriched with Allobaculum (Hypertensive rats), Pseudomonas (Obese rats) and Lactobacillus (Obese rats) after exercise training, while Streptococcus (Wistar rats), Aggregatibacter (Hypertensive rats) and Sutturela (Hypertensive rats) were shown to be more abundant before exercise training was performed (Figure 2). The Lactobacillus genus presented higher abundance after exercise only in Obese rats (13.4% p < 0.05) (Figure 2C). In agreement with our data, the recent study of Queipo-Ortuno et al.  revealed that Lactobacillus was also enhanced with exercise in lean rats, with a longer exercise stimulus (6 weeks). Lactic acid bacteria (LAB), represented in our study by Lactobacillus (enriched after exercise), are associated with the mucosal surface of the small intestine and colon in animals, where they produce lactic acid though homo or heterofermentative metabolism . In this second process, besides lactic acid, CO2, acetic acid and/or ethanol are produced , which may contribute to a more acidic environment . It has been reported that LAB in the gastrointestinal tract leads to positive health benefits with influence on microflora, modulation of mucosal immunity and exclusion of pathogens . The enrichment of Lactobacillus in Obese rats after exercise may have some influence on gastrointestinal acidity trough the production of acidic compounds (e.g. lactic acid, acetic acid); however, this parameter was not measured in the present study. The capillary blood lactic acid was measured in order to establish aerobic capacity and thus to be used as a parameter for adaptation to exercise. Moreover, the lactate produced by Lactobacillus bacteria is converted into butyrate in the gut through bacteria such as B. Coccoides and E. rectale, also found to be enhanced after exercise . Furthermore, butyrate is shown to be related to mucin synthesis and gut epithelium protection . In another study, Matsunomo et al.  showed that exercise altered microbiota and enhanced n-butyrate concentration in rats’ cecum. Therefore, enhanced Lactobacillus found in Obese rats group may possibly have a positive role in the gastrointestinal environment of these animals. It has been reported that obesity-associated gut microbiota is enriched in some species of Lactobacillus (e.g. Lactobacillus reuteri)  while other species (e.g. Lactobacillus gasseri BNR17) are involved in metabolism regulation , presenting anti-obesity effects .
In our study, the Sutterela genus was more abundant before exercise training in Hypertensive rats (Figure 2B). The role of this genus in inflammatory bowel disease has been recently investigated with no relation being found . Nevertheless, more research is needed to understand the relation of Sutturela with exercise and its possible gastrointestinal protective effect.
It is observed that these alterations in genera are not consistent across all host phenotypes. We believe that part of this inconsistency may be related to the biologic differences peculiar to each host genotype used in the study. Obesity has been shown to modulate gut microbiota [10, 11], while no study has shown this yet in a hypertensive phenotype. Furthermore, as a first exploratory study to use different genotypes with different pathologies, it is interesting to note that the gut microbiota of 3 different phenotypes (and 2 genotypes) was possibly altered by an external factor such as exercise.
Seven bacterial species were shown to have a significant differential abundance altered between the three animal lineages (Figure 3). From fecal samples collected pre-exercise, only Bacteroides acidifaciens and Ruminococcus flavefaciens presented a significant differential abundance (Figure 3A and B). While Bacteroides acidifaciens was more enriched in Obese rats compared to Wistar and Hypertensive rats (Figure 3A), Ruminococcus flavefaciens showed an opposite profile, being more enriched in Wistar rats followed by Hypertensive rats and significantly depleted in Obese rats (Figure 3B). B acidifaciens has recently been shown to have an important role in the production of imunoglobulin (IgA) in the large intestine of mice . This production plays an adaptive role in the intestinal mucosal immune system . Since IgA is enhanced in metabolic disorders , the relative abundance of Bacteroides acidifaciens in Obese rats (Figure 3A) may be associated with the role of gut microbiota in the inflammatory signalling peculiar to obesity . Otherwise, Ruminococcus flavefaciens is a cellulolytic bacterium present in the rumen of mammals, and it has been shown to be inhibited by probiotic supplementation (L. acidophilus NCFM) in young rats . Our data indicate that the obesity phenotype from obese rats may suppress this particular species.
In samples collected after training, Streptococcus alactolyticus, Bifidobacterium animalis, Ruminococcus gnavus, Aggregatibacter pnemotropica and Bifidobacterium pseudolongum were all more abundant in Obese rats (Figure 3C-G). Streptococcus alactolyticus and Bifidobacterium animalis species were shown to be present in the gut of obese rats. In contrast to our data, Bifidobacterium is often associated with lean phenotypes ; however, our study showed that this species was completely absent in non-obese Wistar rats and Hypertensive rats (Figure 3D). Regarding the relative abundance of Ruminococcus gnavus in Obese rats, this species is known to have an antibacterial effect and to protect the host from pathogens , also found to be reduced in colon cancer tissue . However, this species was also shown to be enhanced in diverticulitis , which is commonly associated with obesity . Bifidobacterium pseudolongum was another species almost exclusively in Obese rats (Figure 3G). Moreover, the content of this species was shown to be enhanced in obese mice induced by diet and probiotic administration, when compared to a group of mice without probiotic supplementation .
In the present study, the MLSS was used to assess aerobic improvement as a result of four weeks of exercise training at moderate intensity . Thus, after exercise training, a significant reduction in BLC was observed in all rat lineages (Figure 1B), which is associated with an improved aerobic capacity when compared to higher BLC from untrained rats (Figure 1B). The OTUs from Oscillospira and Ruminococcus families were found to be negatively correlated with BLC (Clostridiaceae, R = −0.82, P < 0.01; Bacteroidaceae, R = −0.73 P < 0.01) as were the OTUs from Ruminococcus genus (R = −0.75 P < 0.01). In these three cases, the greatest relative abundance of OTUs in both families is correlated with lower BLC, indicating that exercise training may be favorable to the proliferation of these specific OTUs (Figure 5A, B and D). Otherwise, OTUs from Oscillospira presented a positive correlation with BLC (R = 0.78, P < 0.01) (Figure 5C). The relative abundance within these OTUs is seen to increase when the concentration of lactate goes over ~3.5 mmol.L−1. Since lower BLC during the MLSS test was associated with a more trained status, this result may indicate that exercise training may affect the abundance of the OTUs from this genus.
These findings suggest that exercise training is capable of altering gut microbiota at genus level, with significant alteration in bacterial composition and diversity in obese, non-obese and hypertensive rats. Exercise was shown to enhance the relative abundance of three genera, with Lactobacillus being the most abundant, while another three genera were shown to be more abundant before exercise training (Streptococcus, Aggregatibacter and Sutterella). Non-obese Wistar rats and spontaneously hypertensive rats were shown to share similar microbiota, unlike Obese rats. Rat lineages were also shown to harbor a differential abundance at species level, and six species were shown to be significantly more abundant in obese rats. Two bacterial families (Clostridiaceae and Bacteroidaceae) and two genera (Oscillospira and Ruminococcus) were also shown to significantly correlate with blood lactate accumulation, while exercise was shown to be favorable to the two families and Ruminococcus genus in opposition to Oscillospira. In conclusion, this is the first study to use controlled exercise parameters to assess gut bacterial community modification in three different animal lineages, which may reflect the potential of exercise to alter gut microbial community. However, the effect of exercise on the acidity of lumen or fecal samples was not measured, which limits us in establishing a direct link between exercise and gut alteration by acidic induction. Thus, more studies are necessary to establish these modifications as possible therapeutic implications for obesity or hypertension treatment through the modulation of gut microbiota.
The authors acknowledge grant support from CNPq, CAPES, FAPDF and Catholic University of Brasilia – Brazil.
- WHO: Global Recommendations on Physical Activity for Health. 2010, Geneva, Switzerland: WHO Library Cataloguing-in-Publication DataGoogle Scholar
- Almeida JA BAP, Gomes CP, Araujo RC, Pereira RW, Franco OL: Exercise training at MLSS decreases weight gain and increases aerobic capacity in obese zucker rats. Int J Sports Med. 2014, 35 (3): 199-202.PubMedGoogle Scholar
- Halliwill JR, Buck TM, Lacewell AN, Romero SA: Postexercise hypotension and sustained postexercise vasodilatation: what happens after we exercise?. Exp Physiol. 2013, 98 (1): 7-18.PubMedView ArticleGoogle Scholar
- Mende CW: Obesity and hypertension: a common coexistence. J Clin Hypertens (Greenwich). 2012, 14 (3): 137-138.View ArticleGoogle Scholar
- Haslam DW, James WP: Obesity. Lancet. 2005, 366: 1197-1209.PubMedView ArticleGoogle Scholar
- Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI: An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006, 444 (7122): 1027-1031.PubMedView ArticleGoogle Scholar
- Harris K, Kassis A, Major G, Chou CJ: Is the gut microbiota a new factor contributing to obesity and its metabolic disorders?. J Obes. 2012, 2012: 879151-PubMed CentralPubMedGoogle Scholar
- Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende DR, Li J, Xu J, Li S, Li D, Cao J, Wang B, Liang H, Zheng H, Xie Y, Tap J, Lepage P, Bertalan M, Batto JM, Hansen T, Le Paslier D, Linneberg A, Nielsen HB, Pelletier E, Renault P, et al: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464 (7285): 59-65.PubMed CentralPubMedView ArticleGoogle Scholar
- Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2009, 457 (7228): 480-484.PubMed CentralPubMedView ArticleGoogle Scholar
- Ley RE, Backhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI: Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A. 2005, 102 (31): 11070-11075.PubMed CentralPubMedView ArticleGoogle Scholar
- Ley RE, Turnbaugh PJ, Klein S, Gordon JI: Microbial ecology: human gut microbes associated with obesity. Nature. 2006, 444 (7122): 1022-1023.PubMedView ArticleGoogle Scholar
- Murphy EF, Cotter PD, Healy S, Marques TM, O'Sullivan O, Fouhy F, Clarke SF, O'Toole PW, Quigley EM, Stanton C, Ross PR, O'Doherty RM, Shanahan F: Composition and energy harvesting capacity of the gut microbiota: relationship to diet, obesity and time in mouse models. Gut. 2010, 59 (12): 1635-1642.PubMedView ArticleGoogle Scholar
- Perez-Cobas AE, Gosalbes MJ, Friedrichs A, Knecht H, Artacho A, Eismann K, Otto W, Rojo D, Bargiela R, Von Bergen M, von Bergen M, Neulinger SC, Daumer C, Heinsen FA, Latorre A, Barbas C, Seifert J, Dos Santos VM, Ott SJ, Ferrer M, Moya A: Gut microbiota disturbance during antibiotic therapy: a multi-omic approach. Gut. 2012Google Scholar
- Nelson AM, Walk ST, Taube S, Taniuchi M, Houpt ER, Wobus CE, Young VB: Disruption of the human gut microbiota following Norovirus infection. PLoS One. 2012, 7 (10): e48224-PubMed CentralPubMedView ArticleGoogle Scholar
- Santacruz A, Marcos A, Warnberg J, Marti A, Martin-Matillas M, Campoy C, Moreno LA, Veiga O, Redondo-Figuero C, Garagorri JM, Azcona C, Delgado M, Garcia-Fuentes M, Collado MC, Sanz Y: Interplay between weight loss and gut microbiota composition in overweight adolescents. Obesity (Silver Spring). 2009, 17 (10): 1906-1915.View ArticleGoogle Scholar
- Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, Pettersson S: Host-gut microbiota metabolic interactions. Science. 2012, 336 (6086): 1262-1267.PubMedView ArticleGoogle Scholar
- Fei N, Zhao L: An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice. ISME J. 2013, 7 (4): 880-884.PubMed CentralPubMedView ArticleGoogle Scholar
- Xu J, Ahren IL, Prykhodko O, Olsson C, Ahrne S, Molin G: Intake of blueberry fermented by lactobacillus plantarum affects the Gut microbiota of L-NAME treated rats. Evid Based Complement Alternat Med. 2013, 2013: 809128-PubMed CentralPubMedGoogle Scholar
- Million M, Raoult D: The role of the manipulation of the gut microbiota in obesity. Curr Infect Dis Rep. 2013, 15 (1): 25-30.PubMedView ArticleGoogle Scholar
- Aroniadis OC, Brandt LJ: Fecal microbiota transplantation: past, present and future. Curr Opin Gastroenterol. 2013, 29 (1): 79-84.PubMedView ArticleGoogle Scholar
- De Oliveira EP, Burini RC: The impact of physical exercise on the gastrointestinal tract. Curr Opin Clin Nutr Metab Care. 2009, 12 (5): 533-538.PubMedView ArticleGoogle Scholar
- Matsumoto M, Inoue R, Tsukahara T, Ushida K, Chiji H, Matsubara N, Hara H: Voluntary running exercise alters microbiota composition and increases n-butyrate concentration in the rat cecum. Biosci Biotechnol Biochem. 2008, 72 (2): 572-576.PubMedView ArticleGoogle Scholar
- Queipo-Ortuno MI, Seoane LM, Murri M, Pardo M, Gomez-Zumaquero JM, Cardona F, Casanueva F, Tinahones FJ: Gut microbiota composition in male Rat models under different nutritional status and physical activity and its association with serum leptin and ghrelin levels. PLoS One. 2013, 8 (5): e65465-PubMed CentralPubMedView ArticleGoogle Scholar
- Evans CC, LePard KJ, Kwak JW, Stancukas MC, Laskowski S, Dougherty J, Moulton L, Glawe A, Wang Y, Leone V, Antonopoulos DA, Smith D, Chang EB, Ciancio MJ: Exercise prevents weight gain and alters the gut microbiota in a mouse model of high fat diet-induced obesity. PLoS One. 2014, 9 (3): e92193-PubMed CentralPubMedView ArticleGoogle Scholar
- Zucker LM, Antoniades HN: Insulin and obesity in the Zucker genetically obese rat “fatty”. Endocrinology. 1972, 90 (5): 1320-1330.PubMedView ArticleGoogle Scholar
- Okamoto AK: Development of a strain of spontaneously hypertensive rat. Jap Circ J. 1963, 27: 282-293.PubMedView ArticleGoogle Scholar
- Almeida JA, Petriz BA, Gomes CP, Pereira RW, Franco OL: Assessment of maximal lactate steady state during treadmill exercise in SHR. BMC Res Notes. 2012, 5: 661-PubMed CentralPubMedView ArticleGoogle Scholar
- Almeida JA, Petriz BA, Gomes CP, Rocha LA, Pereira RW, Franco OL: Determination of the maximal lactate steady state in obese Zucker rats. Int J Sports Med. 2013, 34 (3): 214-217.PubMedGoogle Scholar
- Fritz M, Rinaldi G: Blood pressure measurement with the tail-cuff method in Wistar and spontaneously hypertensive rats: influence of adrenergic- and nitric oxide-mediated vasomotion. J Pharmacol Toxicol Methods. 2008, 58 (3): 215-221.PubMedView ArticleGoogle Scholar
- Barringer DL, Bunag RD: Uneven blunting of chronotropic baroreflexes in obese Zucker rats. Am J Physiol. 1989, 256 (2 Pt 2): H417-H421.PubMedGoogle Scholar
- O’Donnell MP, Kasiske BL, Cleary MP, Keane WF: Effects of genetic obesity on renal structure and function in the Zucker rat II. Micropuncture studies. J Lab Clin Med. 1985, 106 (5): 605-610.PubMedGoogle Scholar
- Alonso-Galicia M, Brands MW, Zappe DH, Hall JE: Hypertension in obese Zucker rats. Role of angiotensin II and adrenergic activity. Hypertension. 1996, 28 (6): 1047-1054.PubMedView ArticleGoogle Scholar
- Kurtz TW, Morris RC, Pershadsingh HA: The Zucker fatty rat as a genetic model of obesity and hypertension. Hypertension. 1989, 13 (6 Pt 2): 896-901.PubMedView ArticleGoogle Scholar
- Billat VL, Sirvent P, Py G, Koralsztein JP, Mercier J: The concept of maximal lactate steady state: a bridge between biochemistry, physiology and sport science. Sports Med. 2003, 33 (6): 407-426.PubMedView ArticleGoogle Scholar
- Armougom FDR: Exploring microbial diversity using 16S rRNA high-throughput methods. J Comput Sci Syst Biol. 2009, 2 (1): 74-92.View ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7 (5): 335-336.PubMed CentralPubMedView ArticleGoogle Scholar
- Reeder J, Knight R: Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat Methods. 2010, 7 (9): 668-669.PubMed CentralPubMedView ArticleGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72 (7): 5069-5072.PubMed CentralPubMedView ArticleGoogle Scholar
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26 (19): 2460-2461.PubMedView ArticleGoogle Scholar
- Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, Methe B, DeSantis TZ, Human Microbiome C, Petrosino JF, Knight R, Birren BW: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011, 21 (3): 494-504.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007, 73 (16): 5261-5267.PubMed CentralPubMedView ArticleGoogle Scholar
- Clarke KR: Non-parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993, 18: 117-143.View ArticleGoogle Scholar
- Parks DH, Beiko RG: Identifying biologically relevant differences between metagenomic communities. Bioinformatics. 2010, 26 (6): 715-721.PubMedView ArticleGoogle Scholar
- Natividad JM, Verdu EF: Modulation of intestinal barrier by intestinal microbiota: pathological and therapeutic implications. Pharmacol Res. 2013, 69 (1): 42-51.PubMedView ArticleGoogle Scholar
- Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, Hardt PD: Microbiota and SCFA in lean and overweight healthy subjects. Obesity (Silver Spring). 2010, 18 (1): 190-195.View ArticleGoogle Scholar
- Waldram A, Holmes E, Wang Y, Rantalainen M, Wilson ID, Tuohy KM, McCartney AL, Gibson GR, Nicholson JK: Top-down systems biology modeling of host metabotype-microbiome associations in obese rodents. J Proteome Res. 2009, 8 (5): 2361-2375.PubMedView ArticleGoogle Scholar
- Klaenhammer T, Altermann E, Arigoni F, Bolotin A, Breidt F, Broadbent J, Cano R, Chaillou S, Deutscher J, Gasson M, van de Guchte M, Guzzo J, Hartke A, Hawkins T, Hols P, Hutkins R, Kleerebezem M, Kok J, Kuipers O, Lubbers M, Maguin E, McKay L, Mills D, Nauta A, Overbeek R, Pel H, Pridmore D, Saier M, van Sinderen D, Sorokin A, et al: Discovering lactic acid bacteria by genomics. Antonie Van Leeuwenhoek. 2002, 82 (1–4): 29-58.PubMedGoogle Scholar
- Makarova K, Slesarev A, Wolf Y, Sorokin A, Mirkin B, Koonin E, Pavlov A, Pavlova N, Karamychev V, Polouchine N, Shakhova V, Grigoriev I, Lou Y, Rohksar D, Lucas S, Huang K, Goodstein DM, Hawkins T, Plengvidhya V, Welker D, Hughes J, Goh Y, Benson A, Baldwin K, Lee JH, Diaz-Muniz I, Dosti B, Smeianov V, Wechter W, Barabote R, et al: Comparative genomics of the lactic acid bacteria. Proc Natl Acad Sci U S A. 2006, 103 (42): 15611-15616.PubMed CentralPubMedView ArticleGoogle Scholar
- Peng L, Li ZR, Green RS, Holzman IR, Lin J: Butyrate enhances the intestinal barrier by facilitating tight junction assembly via activation of AMP-activated protein kinase in Caco-2 cell monolayers. J Nutr. 2009, 139 (9): 1619-1625.PubMed CentralPubMedView ArticleGoogle Scholar
- Armougom F, Henry M, Vialettes B, Raccah D, Raoult D: Monitoring bacterial community of human gut microbiota reveals an increase in Lactobacillus in obese patients and Methanogens in anorexic patients. PLoS One. 2009, 4 (9): e7125-PubMed CentralPubMedView ArticleGoogle Scholar
- Xie N, Cui Y, Yin YN, Zhao X, Yang JW, Wang ZG, Fu N, Tang Y, Wang XH, Liu XW, Wang CL, Lu FG: Effects of two Lactobacillus strains on lipid metabolism and intestinal microflora in rats fed a high-cholesterol diet. BMC Complement Altern Med. 2011, 11: 53-PubMed CentralPubMedView ArticleGoogle Scholar
- Kang JH, Yun SI, Park HO: Effects of Lactobacillus gasseri BNR17 on body weight and adipose tissue mass in diet-induced overweight rats. J Microbiol. 2010, 48 (5): 712-714.PubMedView ArticleGoogle Scholar
- Mukhopadhya I, Hansen R, Nicholl CE, Alhaidan YA, Thomson JM, Berry SH, Pattinson C, Stead DA, Russell RK, El-Omar EM, Hold GL: A comprehensive evaluation of colonic mucosal isolates of Sutterella wadsworthensis from inflammatory bowel disease. PLoS One. 2011, 6 (10): e27076-PubMed CentralPubMedView ArticleGoogle Scholar
- Yanagibashi T, Hosono A, Oyama A, Tsuda M, Suzuki A, Hachimura S, Takahashi Y, Momose Y, Itoh K, Hirayama K, Takahashi K, Kaminogawa S: IgA production in the large intestine is modulated by a different mechanism than in the small intestine: Bacteroides acidifaciens promotes IgA production in the large intestine by inducing germinal center formation and increasing the number of IgA+ B cells. Immunobiology. 2013, 218 (4): 645-651.PubMedView ArticleGoogle Scholar
- Peterson DA, McNulty NP, Guruge JL, Gordon JI: IgA response to symbiotic bacteria as a mediator of gut homeostasis. Cell Host Microbe. 2007, 2 (5): 328-339.PubMedView ArticleGoogle Scholar
- Rodriguez-Segade S, Camina MF, Carnero A, Lorenzo MJ, Alban A, Quinteiro C, Lojo S: High serum IgA concentrations in patients with diabetes mellitus: agewise distribution and relation to chronic complications. Clin Chem. 1996, 42 (7): 1064-1067.PubMedGoogle Scholar
- Hotamisligil GS: Inflammation and metabolic disorders. Nature. 2006, 444 (7121): 860-867.PubMedView ArticleGoogle Scholar
- Kaplan CW, Astaire JC, Sanders ME, Reddy BS, Kitts CL: 16S ribosomal DNA terminal restriction fragment pattern analysis of bacterial communities in feces of rats fed Lactobacillus acidophilus NCFM. Appl Environ Microbiol. 2001, 67 (4): 1935-1939.PubMed CentralPubMedView ArticleGoogle Scholar
- Dabard J, Bridonneau C, Phillipe C, Anglade P, Molle D, Nardi M, Ladire M, Girardin H, Marcille F, Gomez A, Fons M: Ruminococcin A, a new lantibiotic produced by a Ruminococcus gnavus strain isolated from human feces. Appl Environ Microbiol. 2001, 67 (9): 4111-4118.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen W, Liu F, Ling Z, Tong X, Xiang C: Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLoS One. 2012, 7 (6): e39743-PubMed CentralPubMedView ArticleGoogle Scholar
- Hansen SG, Skov MN, Justesen US: Two cases of Ruminococcus gnavus bacteremia associated with diverticulitis. J Clin Microbiol. 2013, 51 (4): 1334-1336.PubMed CentralPubMedView ArticleGoogle Scholar
- Strate LL, Liu YL, Aldoori WH, Syngal S, Giovannucci EL: Obesity increases the risks of diverticulitis and diverticular bleeding. Gastroenterology. 2009, 136 (1): 115-122. e111PubMed CentralPubMedView ArticleGoogle Scholar
- Park DY, Ahn YT, Park SH, Huh CS, Yoo SR, Yu R, Sung MK, McGregor RA, Choi MS: Supplementation of Lactobacillus curvatus HY7601 and Lactobacillus plantarum KY1032 in diet-induced obese mice is associated with gut microbial changes and reduction in obesity. PLoS One. 2013, 8 (3): e59470-PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.